Bug Bulletin: The GenomeLocPArser error in SplitNCigarReads has been fixed; if you encounter it, use the latest nightly build.
Service note: the GATK website will be undergoing maintenance over the next few hours. Some pages may be unavailable for that time. All docs remain available in the forum.

Local reference for UG SNP-calling

BlueBlue Posts: 24Member
edited March 2013 in Ask the GATK team

The Unified Genotyper calls SNPs relative to the specified, publicly-available reference assembly.

How can I call SNPs (in many samples, with UG) relative to an in-house individual, which I have sequenced at high-coverage?

My current solution is to perform a de novo assembly on the in-house reference individual using e.g. Velvet, and then simply use the fasta as the reference for UG.

Can the publicly-available reference assembly still be useful here for speeding up the mapping and filling-in missing parts ?

My organism is Drosophila melanogaster.

Post edited by Blue on

Best Answer

Sign In or Register to comment.