fedarko / strainFlye

Pipeline for analyzing (rare) mutations in metagenome-assembled genomes
BSD 3-Clause "New" or "Revised" License
8 stars 1 forks source link

Including (some) rare mutations from the decoy contig in the "fdr fix" BCF? #54

Open fedarko opened 1 year ago

fedarko commented 1 year ago

Might be nice, especially if the FDR estimates used for fixing were based on a decoy context. (For example, if the FDR estimates used were based on rare CP2 mutations in the decoy contig, then we are making the assumption that all rare CP2 mutations in the decoy contig are false -- but we have not said anything about rare non-CP2 mutations in the decoy contig.)

There are various ways we could try to "rescue" these mutations, although most of these would probably start to cause problems with the validity of the FDR estimation logic in the first place:

Proobably not worth worrying about for now.