Question: Identify dovetail read-pairs from a bam/sam file
0
gravatar for neja
3.7 years ago by
neja30
Japan
neja30 wrote:

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 • 1.3k views
ADD COMMENTlink modified 3.7 years ago by Pierre Lindenbaum127k • written 3.7 years ago by neja30
1
gravatar for Pierre Lindenbaum
3.7 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum127k wrote:

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.

ADD COMMENTlink written 3.7 years ago by Pierre Lindenbaum127k
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: 1178 users visited in the last hour