Question: Remove contaminants by aligning to two genomes and selecting reads from one genome using samtools view
0
gravatar for aylward.megan
3.2 years ago by
aylward.megan0 wrote:

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?

ADD COMMENTlink written 3.2 years ago by aylward.megan0
1

Use BBSplit.sh from BBMap to bin the reads: How to remove contamination from NGS data
http://seqanswers.com/forums/showthread.php?t=41288

ADD REPLYlink modified 3.2 years ago • written 3.2 years ago by genomax71k

Ok thanks for your suggestion.

ADD REPLYlink modified 3.2 years ago • written 3.2 years ago by aylward.megan0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1761 users visited in the last hour