This is a test version of Biostars. For the public version, visit https://www.biostars.org.
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?

samtools contamination multiple genome alignment

0 answers

No answers yet.

Log in to answer this question.