Question: How can I remove Indels from a VCF file that are actually microsatellite variants?
0
gravatar for anp375
17 months ago by
anp375150
anp375150 wrote:

I don't have access to the BAM files that were used to generate the VCF. I only have the VCF file.

How do I go through indels in the VCF file and remove sequence repeat variants?

Thanks.

microsatellite indel • 808 views
ADD COMMENTlink modified 17 months ago by Pierre Lindenbaum118k • written 17 months ago by anp375150
1
gravatar for Pierre Lindenbaum
17 months ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum118k wrote:

Two options: if the positions of the microsat is known, use a BED file of those microsat. e.g:

http://hgdownload.cse.ucsc.edu/goldenPath/hg38/database/stsMap.txt.gz

or using GATK https://software.broadinstitute.org/gatk/documentation/tooldocs/current/org_broadinstitute_gatk_tools_walkers_annotator_VariantAnnotator.php and the module HomopolymerRun ( https://software.broadinstitute.org/gatk/documentation/tooldocs/current/org_broadinstitute_gatk_tools_walkers_annotator_HomopolymerRun.php )

ADD COMMENTlink written 17 months ago by Pierre Lindenbaum118k

Thank you. Hopefully I get the BAMs soon.

ADD REPLYlink written 17 months ago by anp375150
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1504 users visited in the last hour