Possible bug in perlscript formatNcbiData.pl

[e099...@gmail.com]

Hi,
First of all I wanted to say thanks for the script formatNcbiData.pl published to split microbial data from NCBI db, I also wanted to report a little bug that I think I've found in the perl script used when clustering NCBI db data for 5s/16s/23s. The counter of 5S is not incremented in the if statement when checking for 5S:

[formatNcbiData.pl (line 104)]
---------------------------------------
} elsif($r->[0] eq '5S') {
unless(exists $derepHash{'5S'}->{$md5}) {
print OUT5 '>'.$db.'5S_'.$count5."\n".$tmpseq."\n";
$derepHash{'5S'}->{$md5}->{c} = $count5;
$derepHash{'5S'}->{$md5}->{l} = length($tmpseq);
#-->> $count16++; <<--#
}
$derepHash{'5S'}->{$md5}->{i}->{$id} = [@taxs];
---------------------------------------

Greetings,
Xavi

[Robert Schmieder]

Thanks Xavi for the bug report! I fixed the 5S counter and updated the version available on SourceForge.