torognes / vsearch

Versatile open-source tool for microbiome analysis
Other
671 stars 125 forks source link

fastq_stripleft when the resulting length is null? #533

Closed frederic-mahe closed 1 year ago

frederic-mahe commented 1 year ago

What happens when the --fastq_stripleft value is equal to or greater than the length of a sequence? This is a question that was sent to me directly. I submit it as a GitHub issue for future referencing.

The answer is that if a sequence is equal to or shorter than the --fastq_stripleft value, vsearch will discard it. The number of discarded sequences is reported by vsearch at the end of the process. Tests were added to cover that particular aspect (https://github.com/frederic-mahe/vsearch-tests/commit/d4b860615e070cedf41939f3b3fa6145d2bceeb9). Documentation will be updated, effectively closing that issue.