MaestSi / MetONTIIME

A Meta-barcoding pipeline for analysing ONT data in QIIME2 framework
GNU General Public License v3.0
78 stars 17 forks source link

setting minimum sequence lenght when running MetONTIIME directly from fastq.gz files #38

Closed dloukovi closed 2 years ago

dloukovi commented 2 years ago

Dear all, I have ready fastq.gz files from ONT guppy basecaller. I have run them in MetONTIIME and everything seems ok. However, I want to filter the sequences based on read length (e.g. keep sequences above 1000bp). Can I add a parameter in the command in order to do this??

Thanks a lot!!

MaestSi commented 2 years ago

Hi, you can prefilter them using something like: zless reads.fastq.gz | NanoFilt -l 1000 | gzip > reads.filtered.fastq.gz Let me know if it works! Simone

dloukovi commented 2 years ago

Thanks Simone!!!!