Remove contaminants by aligning to two genomes and selecting reads from one genome using samtools view
0
1
Entering edit mode
8.4 years ago

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 • 2.9k views
ADD COMMENT
2
Entering edit mode
ADD REPLY
0
Entering edit mode

Ok thanks for your suggestion.

ADD REPLY

Login before adding your answer.

Traffic: 1721 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6