Question: insert size via bamtools stats inner or outer ?
1
gravatar for Siva
2.6 years ago by
Siva10
France
Siva10 wrote:

Hi guys, 

I've a doubt about the meaning of "insert size" from bamtools stats.
I've done this : "bamtools stats -in $1 -insert" input is a bam produced with read1.fastq and read2.fastq (bwa), in each file reads' length is 100 b

and got an insert size of 241 so i supposed it's outer size = R1+R2+insert 

 

**********************************************
Stats for BAM file(s):
**********************************************
Total reads:       820130832 ...

Average insert size (absolute value): 309.402
Median insert size (absolute value): 241

 

Advance thanks, 

Siva

insert bamtools stat • 1.1k views
ADD COMMENTlink modified 2.6 years ago by Zev.Kronenberg11k • written 2.6 years ago by Siva10

Relevant resource: http://seqanswers.com/forums/showthread.php?t=15511

ADD REPLYlink written 2.6 years ago by Ram12k
0
gravatar for Zev.Kronenberg
2.6 years ago by
United States
Zev.Kronenberg11k wrote:

test it:

 

// if the alignment is downstream of mate then you can get the end position of the read.

if(al.Position > al.MatePosition){

int outter = abs(al.MatePosition - al.GetEndPosition());

int inner = al.InsertSize;

if(outter != inner){

cerr << "insert length is inner" << end;

}

 

ADD COMMENTlink modified 2.6 years ago • written 2.6 years ago by Zev.Kronenberg11k
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: 1305 users visited in the last hour