Question: Bam2Cfg Parameter Description
gravatar for ryan.hartmaier
7.7 years ago by
United States
ryan.hartmaier0 wrote:

I've been looking through all the Breakdancer docs for more detailed description of the options in - specifically for what -s represents. The 'minimal mean insert size' seems a little ambiguous. If it was 'minimum insert size' or 'mean insert size' I'd know exactly what it means, but 'minimal insert size' doesn't make any sense yet it seems that it likely needs to be modified from the default value for large insert libraries. Not sure how critical this is but the documentation for the bam2cfg is very limited. Thanks for any help/clarification.

breakdancer • 1.8k views
ADD COMMENTlink modified 3.8 years ago by lionair10 • written 7.7 years ago by ryan.hartmaier0
gravatar for Chris Miller
7.7 years ago by
Chris Miller21k
Washington University in St. Louis, MO
Chris Miller21k wrote:

I'm looking through the code here and it seems as though for each library in your bam file, bam2cfg says "what's the mean insert size of this library?". Then it checks to see if that mean insert size is below the minimum threshold, (set with -s). If so, it skips that library. So really, it is properly named, although I understand your confusion.

ADD COMMENTlink written 7.7 years ago by Chris Miller21k

Thanks Chris. I'm too focused on my single library bam file - that make sense.

ADD REPLYlink written 7.7 years ago by ryan.hartmaier0
gravatar for lionair
3.8 years ago by
lionair10 wrote:


I also have a question on Bam2Cfg Parameter Description. What is the difference between minimum mapping quality and alternative mapping quality? If lets say I'd like to use minimum mapping quality 30, do I need to include -m ( Using mapping quality instead of alternative mapping quality) or -q (Minumum mapping quality [35]) or both?

perl -m -q 30 bam_file ?

It's very confusing because if say I use minimum mapping quality, for BreakDancerMax on the other hand there is -q parameter ( minimum alternative mapping quality [35]). And there is no parameter for minimum mapping quality like However in the explanation of the parameters, it stated -q is minimum mapping quality, not alternative mapping quality.

-q INT minimum alternative mapping quality [35]

"The -q specifies the MAQ mapping quality threshold and can be used to skip reads that are not confidently mapped."

Hopefully someone can help me shed some light on this.


ADD COMMENTlink written 3.8 years ago by lionair10
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1695 users visited in the last hour