Question: How can I remove Indels from a VCF file that are actually microsatellite variants?
0
gravatar for anp375
23 months ago by
anp375160
anp375160 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 • 948 views
ADD COMMENTlink modified 23 months ago by Pierre Lindenbaum122k • written 23 months ago by anp375160
1
gravatar for Pierre Lindenbaum
23 months ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum122k 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 23 months ago by Pierre Lindenbaum122k

Thank you. Hopefully I get the BAMs soon.

ADD REPLYlink written 22 months ago by anp375160
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: 2311 users visited in the last hour