Hi
See
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.
Hi
See
CONFOLD/confold.pl
Line 1459 in 6b29b50
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.