Question: Does Genome MuSiC use soft-clipped bases in its coverage calculation?
gravatar for kmegq
24 days ago by
kmegq0 wrote:


I am wondering whether MuSiC uses soft-clipped bases in its coverage and background mutation rate calculations.

I have generated somatic variant calls from WES tumor/normal pairs using MuTect2 with the --dontUseSoftClippedBases option, and I want to be sure that the coverage used for the bmr calculation matches the coverage that was actually used for calling variants. I'm not sure if I need to trim these bases out of the bam files before passing them to MuSiC, or if they would be ignored anyway.

Thank you for your help!

Best regards, Kate

ADD COMMENTlink modified 21 days ago by Chris Miller18k • written 24 days ago by kmegq0
gravatar for Chris Miller
21 days ago by
Chris Miller18k
Washington University in St. Louis, MO
Chris Miller18k wrote:

I believe that the answer is no, as it's using the samtools libraries, and those do not return soft-clipped bases by default, IIRC. That said, it'd be easy to verify. Extract a dozen reads from a region of your bam file, choose a position that has both soft-clipped and non-soft-clipped bases, then run the tool just on that bam and that site. Also, the calcRoiCovg code is here, if it helps:

ADD COMMENTlink written 21 days ago by Chris Miller18k
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: 1588 users visited in the last hour