Question: One of two mates listed as "properly paired" in Bowtie2
0
gravatar for aredash
2.0 years 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 • 627 views
ADD COMMENTlink modified 2.0 years ago • written 2.0 years ago by aredash0
1
gravatar for Devon Ryan
2.0 years ago by
Devon Ryan91k
Freiburg, Germany
Devon Ryan91k 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 2.0 years ago by Devon Ryan91k

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

ADD REPLYlink written 2.0 years 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: 2205 users visited in the last hour