Cleaning FASTQ files containing a contamination sequence.
Using the technology [censored] we have found a sequence (n = ~80bp) which look like a contaminant in our 2x150 fastq files.
I would like to trim the whole sequence of the reads containing some significative parts of this contaminant.
What would be the best tool to achieve this ? cutadapt ? whith which params ?
Thanks,
P.
• 2,890 views
•
link
2 answers
bbduk.sh from BBMap with literal=your_contaminant_seq option. Significant part is how many bases? Many other options for bbduk to play with.
• 0 views
•
link
I recently discovered AdapterRemoval https://github.com/MikkelSchubert/adapterremoval, give it a try.
• 0 views
•
link
Log in to answer this question.