I'm interested in obtaining coding sequences of my favourite gene in all individuals from the 1000 genomes project (and similar endeavours). I use GATK to get the right subset of variants, vcf-consensus to map these variants onto the reference genome and finally samtools to extract the individual exons. This works fine if the variants are SNPs but if there are any indels, this changes the coordinates of exons and I end up getting the wrong region. Is there any generic way of remapping genomic coordinates to account for the changes created by indels?
Question: Remapping genomic coordinates to account for indels
2.7 years ago by
jergosh • 50
jergosh • 50 wrote:
ADD COMMENT • link •
Please log in to add an answer.
Powered by Biostar version 2.3.0
Traffic: 1351 users visited in the last hour