Remove contaminants by aligning to two genomes and selecting reads from one genome using samtools view
I have sequencing reads from a non-human primate sequencing project, but also have some human contamination. To remove the contaminants I would like to align the reads to both human and non-human reference genomes then select only the reads which primary alignment is to the non-human primate genome. I know that particular regions can be sub-set using samtools view. However the non-human genome has over 3 million scaffolds is there a way to select all reads by genome rather than by scaffold?
• 3,371 views
•
link
0 answers
No answers yet.
Log in to answer this question.
Use BBSplit.sh from BBMap to bin the reads: How to remove contamination from NGS data
http://seqanswers.com/forums/showthread.php?t=41288
Ok thanks for your suggestion.