CW_mgatk.read.RdThis function allows you to read raw data from XX/final folder, the output from mitoV
CW_mgatk.read(path, Processed = F)The XX/final folder, the output from mitoV
Boolean variable (Default F), if true directly readRDS("VariantsGTSummary.RDS") or, generate and saveout "VariantsGTSummary.RDS"
this returns depth which is a list of 4 df (Total/VerySensitive/Sensitive/Specific), each is a genotype summary
WD<-"/lab/solexa_weissman/cweng/Projects/MitoTracing_Velocity/SecondaryAnalysis/Donor01_CD34_1_Multiomekit/MTenrichCombine/Enrich/CW_mgatk/final"
DN1CD34_1.VariantsGTSummary<-CW_mgatk.read(WD,Processed =T)