Question: Chimeric read insert length in SAM $9
0
gravatar for Yijun Tian
8 months ago by
Yijun Tian20
Medical College of Wisconsin
Yijun Tian20 wrote:

Hello everyone,

I am recently doing a script to dissect chimeric alignment from BAM file, here is my code:

#!/bin/bash
samtools=`which samtools`
bn=$(basename $1 .bam)
echo processing $bn
$samtools view -@ 20 -h -b -f97 -F22 $1 > $bn.flag97.bam
$samtools view -@ 20 -h -b -f145 -F42 $1 > $bn.flag145.bam
$samtools view -@ 20 -h -b -f161 -F22 $1 > $bn.flag161.bam
$samtools view -@ 20 -h -b -f81 -F42 $1 > $bn.flag81.bam
echo merging tmp BAM files
$samtools merge $bn.uns.pp.bam $bn.flag97.bam $bn.flag145.bam $bn.flag161.bam $bn.flag81.bam
echo sorting final BAM file
$samtools view -H $bn.uns.pp.bam -o SamHeader
$samtools view $bn.uns.pp.bam |awk '$7=="=" && $9^2>150^2  {print $0}' > SamChimeric
cat SamHeader SamChimeric |$samtools view -bS > $bn.Chimeric.bam
$samtools sort $bn.Chimeric.bam > $bn.Chimeric.sorted.bam

I have a question about the insert length. Since I was looking for some alignment located in circular DNA, so I assume the insert lengths for such read pair are incalculable. But seeing from the BAM file, the insert length are all given for each alignment:

$ samtools view s4_E.PE.sorted.Chimeric.sorted.bam | cut -f 9 | wc -l
454883
$ samtools view s4_E.PE.sorted.Chimeric.sorted.bam | cut -f 9 | grep -v '\*' | wc -l
454883

How does aligner give the insert size? (I'm using bowtie2)

alignment • 198 views
ADD COMMENTlink written 8 months ago by Yijun Tian20
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: 1445 users visited in the last hour