Skip to content

count_ss_match uses wrong file #1

@johnlamb

Description

@johnlamb

Hi
See

sub count_ss_match{

In the count_ss_match function, $file_sec is sent in as a parameter(in both instances when the function is used it is the secondary structure file that is sent in). However, this file is never referenced and %residue_ss1 is instead created from $file_fasta.
This causes the later comparison to compare amino acid vs secondary structure characters which will only work in cases of H, E and C and will hence give a wrong secondary structure match count.
I am happy to do a pull request with a correction if you can confirm it is a bug.

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions