Question: Bam2Cfg Parameter Description
0
gravatar for ryan.hartmaier
4.5 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 bam2cfg.pl - 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 • 980 views
ADD COMMENTlink modified 7 months ago by lionair0 • written 4.5 years ago by ryan.hartmaier0
2
gravatar for Chris Miller
4.5 years ago by
Chris Miller18k
Washington University in St. Louis, MO
Chris Miller18k 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 4.5 years ago by Chris Miller18k

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

ADD REPLYlink written 4.5 years ago by ryan.hartmaier0
0
gravatar for lionair
7 months ago by
lionair0
Malaysia
lionair0 wrote:

Hi,

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 bam2cfg.pl -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 bam2cfg.pl. However in the explanation of the parameters, it stated -q is minimum mapping quality, not alternative mapping quality.

http://breakdancer.sourceforge.net/breakdancermax.html

-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.

Thanks.

ADD COMMENTlink written 7 months ago by lionair0
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: 1676 users visited in the last hour