search

smithlabcode / falco

A C++ drop-in replacement of FastQC to assess the quality of sequence read data
https://falco.readthedocs.io
GNU General Public License v3.0
96 stars 10 forks source link

Adapter content plot shows wrong read lengths #37

Closed shigita closed 2 years ago

shigita commented 2 years ago

In the version 1.0.0, the plot shows read length of 489 bp, even though the actual length is 100 bp.

image

The command, fastq file, and adapter list to reproduce this issue are as follows:

falco -a qc_adaptors_list.txt example.fq.gz

qc_adaptors_list.txt example.fq.gz

andrewdavidsmith commented 2 years ago

@shigita This should be fixed in commit 78eac8d44b67c04005542855d70aafe1879a608c and I'm able to get the image with the range to 100 bp using your files and command (thanks for those!!). Be aware: adaptor enrichment can only be calculated up to 91 positions, if the smallest adapter is 10bp and the read length is 100, but the most recent commit plots all the way to 100, instead of 91, which is reasonable since the plot is cumulative. I'm closing the issue, but happy to reopen if needed.