Difference between revisions of "Vt"
From Genome Analysis Wiki
Jump to navigationJump to searchLine 34: | Line 34: | ||
vt normalize <input.vcf> -r <sequence.fa> -o <output.vcf> | vt normalize <input.vcf> -r <sequence.fa> -o <output.vcf> | ||
<div class="mw-collapsible-content"> | <div class="mw-collapsible-content"> | ||
− | # | + | #normalize variants, send to standard out and remove duplicates. |
− | + | vt normalize mills.vcf -r seq.fa -o - | vt merge_duplicate_variants - -o mills.normalized_merged.vcf | |
− | vt normalize mills.vcf -r seq.fa -o | ||
− | |||
</div> | </div> | ||
</div> | </div> |
Revision as of 17:21, 21 October 2013
Introduction
vt is a variant tool set that discovers short variants from Next Generation Sequencing data. The features are being rolled out to github as major rewriting is being undertaken.
Location
The source files are housed in github.
git clone https://github.com/atks/vt.git
Common options
-i multiple intervals in <seq>:<start>-<end> format delimited by commas.
-I multiple intervals in <seq>:<start>-<end> format listed in a text file line by line.
-o defines the out file which and has the STDOUT set as the default. You may modify the STDOUT to output the binary version of the format.
Programs
Normalization
Normalize variants in a VCF file.
Normalized variants may have their positions changed; in such cases, the normalized variants are reordered and output in an ordered fashion. The local reordering takes place over a window of 10000 base pairs.
vt normalize <input.vcf> -r <sequence.fa> -o <output.vcf>
#normalize variants, send to standard out and remove duplicates. vt normalize mills.vcf -r seq.fa -o - | vt merge_duplicate_variants - -o mills.normalized_merged.vcf
Merge duplicate variants
Merges duplicate variants by position with the option of considering alleles. (This just discards the duplicate variant that appears later in the VCF file)
Options: -i, --input-vcf <string> : Input VCF file -o, --output-vcf <string> : Output VCF file [-] -p, --merge-by-position : Merge by position [false]
Example: e.g. vt merge_duplicate_variants -i 8904indels.dups.genotypes.vcf -o out.vcf e.g. vt merge_duplicate_variants -p -i 8904indels.dups.genotypes.vcf -o out.vcf
Maintained by
This page is maintained by Adrian