SMNCopyNumberCaller icon indicating copy to clipboard operation
SMNCopyNumberCaller copied to clipboard

Low coverage and confidence

Open freconfarmer opened this issue 2 years ago • 2 comments

Amazing little program!

A few questions

  1. The difference from raw to mapped was 81% from my Nebula results. There was a warning before finishing that stated low coverage, confidence low. Would you trust this result? Smn 1=1 Smn2 =1 Would it still be safe to say is a carrier with warning and low confidence?
  2. The confidence was set at 1 for all. What is a perfect score?
  3. I used a Bam file that was converted from the Cram using WGS Extract. Would the Cram file with reference produce a better result?

Thanks for your time!

Conrad

freconfarmer avatar Dec 05 '22 00:12 freconfarmer

Hi Conrad,

Are you using a full WGS bam? What is your sequencing coverage? The warning may suggest something suspicious. If you like, you could send me your .json output file, either attaching it here or emailing me at [email protected].

Thanks, Xiao

xiao-chen-xc avatar Dec 05 '22 04:12 xiao-chen-xc

Hi Xiao,

I am using a 30X WGS bam from Nebula. It was a conversion from the cram by WGS Extract. I'm not sure of which reference genome it used, but it appears as "hs38(Chr), rCRS, 195 SNs" on the stats. Sequencing coverage was 81% with an average read of 29X. I will send you an email of the .json file. For some reason it isn't being recognized to copy.

Thanks

Conrad

freconfarmer avatar Dec 05 '22 13:12 freconfarmer