Question: paired reads aligned at the same genomic position
0
gravatar for irfanwustl
4 weeks ago by
irfanwustl10
irfanwustl10 wrote:

We know sometimes paired reads can align at the same position. In such a case, is there any simple command to remove one of the paired reads from the bam file?

bam alignment • 99 views
ADD COMMENTlink modified 4 weeks ago by MatthewP780 • written 4 weeks ago by irfanwustl10
1
gravatar for MatthewP
4 weeks ago by
MatthewP780
China
MatthewP780 wrote:

Paired reads align at same position probably because some insert fragments is too short that even shorter than your reads, you can use fastp to get fragment length distribution report, so you can see is there some fragments too short.
I don't think you need to remove one of them, you can see bowtie2 manual "dovetail" part, aligner can handle such case.

ADD COMMENTlink written 4 weeks ago by MatthewP780

Actually, I want to count how many reads aligned to a specific position (like depth) from a python code. That's why I was trying to remove one read.

ADD REPLYlink modified 29 days ago • written 29 days ago by irfanwustl10
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: 1585 users visited in the last hour