BWA option about the insert size
Hi,
I am using BWA to do alignment and estimate the insert size of my data and the excepted insert size length is 2000bp,so I use the "sampe" command to do alignment,but I found "-a" option is available,and the description of this option is:
Maximum insert size for a read pair to be considered being mapped properly. Since 0.4.5, this option is only used when there are not enough good alignment to infer the distribution of insert sizes. [500]
I use this option and set it to "-a 2000", but I don't know clearly what impact will this option have on my result?Should I use this option?
Thanks.
• 3,837 views
•
link
written
by
xieduo9394 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Bwa sampe - BGI
written by lorena9132 •
-
bwa fail to infer insert size
written by 0range •
-
SSPACE for scaffolding
written by deepti1rao •
-
RNASEQ preprocessing and insert size
written by Omics data mining
-
Best-guessing insert size and stddev in MISO using paired-end data versus running MISO single ended
written by floris.barthel •
-
Flexibility of inner distance parameter for Tophat
written by acorella •
-
long insert size estimation
written by Damian Kao
-
Very large inferred maximum insert size!
written by Nikleotide
-
Estimating Insert Size From Paired End Data.
written by GouthamAtla
-
Bwa Sampe Very Large Insert Size
written by Bioscientist
If this is standard illumina sequencing it is unlikely that your inserts are 2kb in size. If you have a reference available you can easily figure out the insert size using
bbmap.shfrom BBMap suite.bbmap.sh ref=ref.fasta in1=r1.fq in2=r2.fq ihist=ihist.txtWithout reference, and much faster, is BBMerge (what a surprise that it's part of BBMap!), with a little command line addition to catch output in a text file and allow it to be output to screen also:
This will only work when there are reads (I am not sure how many need to) that overlap.