Bug Bulletin: we have identified a bug that affects indexing when producing gzipped VCFs. This will be fixed in the upcoming 3.2 release; in the meantime you need to reindex gzipped VCFs using Tabix.

Local reference for UG SNP-calling

BlueBlue Posts: 17Member
edited March 2013 in Ask the 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

  • Geraldine_VdAuweraGeraldine_VdAuwera Posts: 5,202 admin
    Answer ✓

    It sounds like you want to assemble your in house reference genome using the public reference genome as a scaffold rather than doing it entirely de novo.

    Genome assembly isn't our department so I can't advise you on this, but I'd recommend asking somewhere like SeqAnswers what is the best aligner to do this.

Sign In or Register to comment.