From Genome Analysis Wiki
Jump to navigationJump to search
177 bytes removed
, 09:22, 4 December 2013
Line 83: |
Line 83: |
| | | |
| <div class=" mw-collapsible mw-collapsed"> | | <div class=" mw-collapsible mw-collapsed"> |
− | Merges duplicate variants by position with the option of considering alleles. (This just discards the duplicate variant that appears later in the VCF file)
| + | Discovers variants from reads in a BAM file. |
| <div class="mw-collapsible-content"> | | <div class="mw-collapsible-content"> |
| </div> | | </div> |
Line 89: |
Line 89: |
| | | |
| <div class=" mw-collapsible mw-collapsed"> | | <div class=" mw-collapsible mw-collapsed"> |
− | #merge duplicate variants and save output in mills.merged.vcf | + | #discover variants from NA12878.bam and write to stdout |
− | vt mergedups mills.vcf -o mills.merged.vcf
| + | vt discover -b NA12878.bam -s NA12878 -r hs37d5.fa -i 20 -v snps,indels,mnps |
− | vt discover -b NA12878.bam -s NA12878 -r ~/ref/genome/hs37d5.fa -i 20 -v snps,indels,mnps | |
| <div class="mw-collapsible-content"> | | <div class="mw-collapsible-content"> |
| usage : vt discover [options] | | usage : vt discover [options] |