bac-genomics-scripts icon indicating copy to clipboard operation
bac-genomics-scripts copied to clipboard

cdd2cog and RPSBLAST 2.2.31+ report parse

Open mberacochea opened this issue 8 years ago • 0 comments

Hi,

Thanks for the tool, really useful!.

There seems to be a change on the output of RPS-Blast, an example line:

fig|992418.4.peg.8018 gnl|CDD|223110 26.42 159 100 6 20 165 11 165 4e-13 65.6 45

The subject id format changed, and the current regex on cdd2cog doesn't work. Changing the line 302 to my $pssm_id = $1 if $line[1] =~ /^gnl\|CDD\|(\d+)/; # get PSSM-Id from the subject hit fix the issue.

Regards

mberacochea avatar Feb 04 '18 23:02 mberacochea