Question: One of two mates listed as "properly paired" in Bowtie2
0
gravatar for aredash
10 weeks ago by
aredash0
aredash0 wrote:

I have been running bowtie2 on paired end sequences and noticed something funny when I tried to filter out reads that weren't properly paired. I seemed to be getting orphan reads left behind after filtering with samtools view -f 2. After digging around, I noticed that after the initial alignment, I was getting reads that were flagged 163 and 113.

NS500418:295:HWG3TBGXX:4:11405:20193:10248 163 chr1 1624577 255

NS500418:295:HWG3TBGXX:4:11405:20193:10248 113 chr1 1677965 255

Running samtools fixmate changes the flags to 163 and 81 , but that still leaves the bit 2 unaffected.

As I understand it, it seems like this means that one of the two is "read mapped in proper pair" as indicated by flag 2, but the other one isn't. Am I failing to understand something about what flag2 means?

Here is the bowtie command I ran:

bowtie2 -k1 -N1 -p16 -x hg19 -1 R1.fastq -2 R2.fastq -S file.sam

Any thoughts or help would be appreciated! Thanks!

alignment • 177 views
ADD COMMENTlink modified 10 weeks ago • written 10 weeks ago by aredash0
1
gravatar for Devon Ryan
10 weeks ago by
Devon Ryan73k
Freiburg, Germany
Devon Ryan73k wrote:

Looks like you found a bug, please report it on github.

BTW, you used the "RNA-seq" tag, so I assume you're using this for that purpose. If so, please be aware that this is not an appropriate tool for mammalian RNAseq data. Please use STAR or another splice-aware aligner instead.

ADD COMMENTlink written 10 weeks ago by Devon Ryan73k

Thanks Devon. I also removed the rna-seq tag.

ADD REPLYlink written 10 weeks ago by aredash0
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: 1146 users visited in the last hour