How to generated unaligned .sam file using bowtie2
0
0
Entering edit mode
2.3 years ago
wangxiong • 0

I hope to retain unaligned read for further circRNA identification, so I used the following code to generate unaligned reads using bowtie2.

bowtie2 -x ~/database/genome/UCSC/mm10/index/bowtie2/mm10\ -U ~/project/m6A/Mapping/bowtie2/sample1.fq.gz -p 8\ -S ~/project/m6A/Mapping/bowtie2/sample1.sam\ 2>~/project/m6A/Mapping/bowtie2/sample1.log\ && samtools view -hbuS ~/project/m6A/Mapping/bowtie2/sample1.sam | samtools view -hf 4 - | samtools view -Sb - | samtools sort -@ 6 -o ~/project/m6A/Mapping/bowtie2/sample1.unmapped.bam -

I don't know whether the sample1.unmapped.sam only contained the unaligned reads.

Thanks.

bowtie2 unaligned reads • 556 views
ADD COMMENT
0
Entering edit mode

Please see: http://bowtie-bio.sourceforge.net/bowtie2/manual.shtml

...and in there the "output options" section. Explains what you need.

ADD REPLY
0
Entering edit mode

I used "--un ~/project/m6A/Mapping/bowtie2/Sample1.unmapped.sam"

I don't know wether this code generated the the unaligned reads

ADD REPLY

Login before adding your answer.

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