HRGV / phyloFlash

phyloFlash - A pipeline to rapidly reconstruct the SSU rRNAs and explore phylogenetic composition of an illumina (meta)genomic dataset.
GNU General Public License v3.0
77 stars 25 forks source link

Improved single end support #102

Closed rpetit3 closed 4 years ago

rpetit3 commented 4 years ago

Running phyloFlash with single end reads led to the following error:

[12:28:23] Reading remappings to summarize taxonomy of unmapped reads
[12:28:23] Reading results of remapping to SPAdes results
[12:28:23] total fwd reads remapping:8900
Illegal division by zero at /local/home/rpetit3/repos/phyloFlash/phyloFlash.pl line 1961.

This was due to not counting remapped reads for single end reads due to the flags being different then the paired end reads.

Also, if CHAO becomes the string 'n.d.' the following warning is produced.

Argument "n.d." isn't numeric in sprintf

Added logic to only use sprintf if $chao1 is not equal to 'n.d.'