Identify dovetail read-pairs from a bam/sam file

Hi, Is there a way to identify read pairs which are aligned as "dove-tail" by bowtie2? I've used the flag --dovetail while mapping with bowtie2.

alignment

• 2,761 views • link

updated 10.1 years ago by Pierre Lindenbaum

• written 10.1 years ago by neja •

1 answer

Using samjs: https://github.com/lindenb/jvarkit/wiki/SamJS

java -jar dist/samjs.jar -e 'record.getReadPairedFlag() && !record.getReadUnmappedFlag() && !record.getMateUnmappedFlag() && record.getReferenceName().equals(record.getMateReferenceName()) && record.getAlignmentStart()<= record.getMateAlignmentStart() && record.getMateAlignmentStart()<=record.getAlignmentEnd()'  input.bam

you might consider the strand of the read/mate too.

• 1 views • link

written 10.1 years ago by Pierre Lindenbaum

More posts like this

Cut & run - bowtie2 dovetail option
written 5.8 years ago by blur

Hi, I am working on a CUT & RUN library for the first time. The methods state: "Dovetail alignment policy. Bowtie210 aligns each mate of …
Plotting range/type of discordance in alignment data
written 6.8 years ago by rbronste

I was wondering if there was a way, in a Bowtie2 output, to plot the range and type of discordance? The really obvious kind when …
Bowtie2 reports no secondary alignments
written 7.7 years ago by goodez

This has been bugging me for awhile now. I know there are multiple mapping reads in my alignments, but bowtie2 does not assign the 256 …
extracting orphan reads from bam file
written 8.2 years ago by Assa Yeroslaviz

Hi, I have a bam file from WGS sequencing, wehre I try to identify the integration site(s) of a transgene insertion in the mouse genome. …
very-sensitive-local bowtie2 or bwa_mem
written 8.9 years ago by Illusion •

I'm now analysising targeted sequencing pair-end reads (approximately 100M per sample) to find edit pattern of crispr. But how to do mapping and alignment? bowtie2 …
bowtie2 summary looks very weird. Neither concordant nor discordant.
written 8.9 years ago by zhangzhenyusz •

Hi everybody, My bowtie2 summary looks so weird that most read pairs are aligned neither concordantly nor discordantly: output: 4385453 reads; of these: 4385453 (100.00%) …
unmapped reads by Tophat can be aligned by Bowtie2
written 9.0 years ago by zhangzhenyusz •

I used Tophat to map my pair-end reads to hg19. I sorted my unmapped.bam file, converted it to sam file and extracted pairs into two …
How to identify bowtie2 unique mapping using sam flag value?
written 10.5 years ago by jinxinhao1988

<p>I used BWA before,after mapping to genome or cDNA,I can filter unique mapping reads by sam flag value,but BWA mapping efficiency is not as good …
Please help me with bowtie2 output and samtools report
written 10.6 years ago by seta

Hi friends, I have made a de novo transcriptome assembly, followed by mapping back read to the assembly using bowtie2 with this command: bowtie2 -x …
Unaligned Reads From Bowtie2
written 12.8 years ago by deepthithomaskannan

<p>Hi all , I am using bowtie2 for mapping illumina paired end mapping. Is there any option for getting the reads which are aligned neither …