This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Do I have to remove overlapping reads from paire_end data before Metahplan?

I have paired_end files from shoutgun metagenomics analysis (251 bp). Before starting with Metaphan, I run fastqc and fastq_screen to check how my files are. I used KneadData to delete the human genome, and now it is ok. (I also notice that all my files do not pass the “Per Base Sequence Content.” Is this a problem? (All the other control is OK.)

Should I also have to delete overlapping reads between R1 and R2? How can I do it?

Thanks Michela

overlapping paired_ends metaphaln
• 995 views • link
updated by GenoMax GenoMax • written by michelafrancesconi9 • michelafrancesconi9

For the “Per Base Sequence Content”, is really hard to pinpoint what the actual problem is without looking at the graph itself. Is it localised within certain region of the read (leading or trailing)? Also, I don't think it is harmful to keep overlapping R1/R2, it just tells you that your fragments are rather short.

• 0 views • link
written by biofalconch biofalconch

Thank you for the answer! The problem is always located in the first 10/15 bp

enter image description here

• 0 views • link
written by michelafrancesconi9 • michelafrancesconi9

0 answers

No answers yet.

Log in to answer this question.

More posts like this