Kurt-Hetrick / CIDR_SOMATIC_WES_PIPELINES

currently for paired tumor - normal WES samples processed with TWIST UMIs
MIT License
1 stars 0 forks source link

mutect filter output statistics #26

Open Kurt-Hetrick opened 1 year ago

Kurt-Hetrick commented 1 year ago

this is an output summary generated from FilterMutectCalls to consider how to implement into QC summaries. note that this might be different if using multiple normals and/or tumors for a person and joint calling them.

Ln prior of deletion of length 10=-17.294720484381504

Ln prior of deletion of length 9=-17.926711851530523

Ln prior of deletion of length 8=-16.854922493940755

Ln prior of deletion of length 7=-20.72326583694641

Ln prior of deletion of length 6=-16.15525049889726

Ln prior of deletion of length 5=-16.47073888358528

Ln prior of deletion of length 4=-16.873737536879357

Ln prior of deletion of length 3=-15.687966683068682

Ln prior of deletion of length 2=-15.609472346803154

Ln prior of deletion of length 1=-13.166812085363665

Ln prior of SNV=-10.75115164828984

Ln prior of insertion of length 1=-14.513247842301208

Ln prior of insertion of length 2=-16.15980876481199

Ln prior of insertion of length 3=-17.321601219579268

Ln prior of insertion of length 4=-18.93694436843823

Ln prior of insertion of length 5=-16.703017272192877

Ln prior of insertion of length 6=-16.5923855941863

Ln prior of insertion of length 7=-16.59974944689006

Ln prior of insertion of length 8=-20.72326583694641

Ln prior of insertion of length 9=-16.471378376015544

Ln prior of insertion of length 10=-20.72326583694641

Background beta-binomial cluster=weight = 0.2155, alpha = 2.93, beta = 20.50

High-AF beta-binomial cluster=weight = 0.0002, alpha = 10.00, beta = 1.01

Binomial cluster=weight = 0.7842, mean = 0.029

threshold=0.583

fdr=0.133

sensitivity=0.804

filter FP FDR FN FNR weak_evidence 78.72 0.02 177.72 0.03 strand_bias 137.06 0.03 83.13 0.02 contamination 55.17 0.01 8.12 0.0 normal_artifact 58.72 0.01 21.88 0.0 orientation 406.83 0.08 652.94 0.13 slippage 0.94 0.0 0.48 0.0 haplotype 44.4 0.01 52.34 0.01 germline 0.25 0.0 0.3 0.0