why bowtie2 doesn’t report longer than 500bp inserts?
1
1
Entering edit mode
6 weeks ago
lidan0001 ▴ 50

I used The 9th column of .bam file to get the insert size of atac-seq data.

samtools view bam/${FILE}"_mm10.bam" | awk '$9>0' | cut -f 9 | sort | uniq -c | sort -b -k2,2n > insert_size_plotting/${FILE}".txt"

But I found all of the inserts are less than 500bp, From the insert size plot, I can see there should be longer than 500bp inserts, but why bowtie2 doesn’t report >500bp inserts? enter image description here Thanks

bowtie2 BAM • 288 views
ADD COMMENT
3
Entering edit mode
6 weeks ago
acvill ▴ 110

By default, bowtie2 expects a maximum fragment length of 500 bp for paired-end alignments. To adjust this, change the -X or --maxins parameter.

http://bowtie-bio.sourceforge.net/bowtie2/manual.shtml

ADD COMMENT
0
Entering edit mode

I am wondering can I recalculate the insert size from .bam files? Thanks a lot

ADD REPLY

Login before adding your answer.

Traffic: 932 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6