bac-genomics-scripts
bac-genomics-scripts copied to clipboard
cdd2cog and RPSBLAST 2.2.31+ report parse
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