We've moved!
This site is now read-only. You can find our new documentation site and support forum for posting questions here.
Be sure to read our welcome blog!

RNA-editing sites

I am walking myself through the suggested variant call analysis flowthrough and have noticed that in various steps (e.g. UnifiedGenotyper and VariantRecalibrator), one should provide vcf files containing known SNPs.
In the case of UnifiedGenotyper, I assume that omitting the --dbSNP option and adding -contamination dbSNPs,vcf would allow me to keep only unknown (possible RNA-editing) variants.
I am also at loss on how to deal with the above problem in the case of VariantCallRecalibrator and other possible steps downstream UnifiedGenotyper necessary to efficiently obtain high quality variants.

Best Answer


Sign In or Register to comment.