I want to make alternate references using phased variants. (Really I just want the sequences of two genes and the space between them, <10 kb). My individuals are diploid, which is why I needed to phase the variants.
However, the alternate reference programs I've found, like GATK AlternateReference, vcf2fasta, and bcftools, don't seem to take diploidy into account. They just stick in all variants and throw out one variant randomly if there are two at the same locus. And even if there was a program that did, these are not long, contiguous haplotypes so I can see how it would be difficult :\
Does anyone know of any solution to this? Anything better than imposing the variants manually? I guess I could do quick assemblies of only the reads mapped to my region but I'd be giving up the higher quality of my GATK-called variants.
However, the alternate reference programs I've found, like GATK AlternateReference, vcf2fasta, and bcftools, don't seem to take diploidy into account. They just stick in all variants and throw out one variant randomly if there are two at the same locus. And even if there was a program that did, these are not long, contiguous haplotypes so I can see how it would be difficult :\
Does anyone know of any solution to this? Anything better than imposing the variants manually? I guess I could do quick assemblies of only the reads mapped to my region but I'd be giving up the higher quality of my GATK-called variants.