Rsamtools : separating mapped and unmapped reads from BAM or SAM files (for a newbie)
1
1
Entering edit mode
nginet ▴ 10
@nginet-15664
Last seen 6.6 years ago

Dear all,

As a start I must confess that I am not really trained in bioinformatics. I managed to use bowtie2 and some samtools embedded in RStudio to align Illumina reads to a reference genome. I wish now to write into 2 separate files (to be converted in fastq) the mapped and unmapped reads. I need to perform a de novo assembly with the unmapped reads. Can anyone suggest a script that I can copy and adapt to my needs ?

Many thanks

Nico

R sequencing alignment rsamtools mapping • 3.4k views
ADD COMMENT
0
Entering edit mode
@james-w-macdonald-5106
Last seen 8 hours ago
United States

https://www.biostars.org/p/56246/
 

ADD COMMENT

Login before adding your answer.

Traffic: 550 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