1,102
edits
Changes
From Genome Analysis Wiki
Jump to navigationJump to searchSequencing Workshop Analysis of Indels (view source)
Revision as of 16:47, 16 June 2014
, 16:47, 16 June 2014→Normalization
It is usually useful to examine the call sets against known data sets for the passed variants.
It is usually useful to examine the call sets against known data sets for the passed variants.
vt profile_indels -g /net/fantasia/home/atks/ref/vt/grch37/indel.reference.txt -r /net/fantasia/home/atks/ref/vt/grch37/hs37d5.fa run/final/all.genotypes.bcf -i 22:36000000-37000000 -f "PASS"
vt profile_indels -g indel.reference.txt -r hs37d5.fa all.genotypes.bcf -i 22:36000000-37000000 -f "PASS"
data set
data set
We perform the same analysis for the failed variants again, the relatively low overlap with known data sets imply a reasonable tradeoff in sensitivity and specificity.
We perform the same analysis for the failed variants again, the relatively low overlap with known data sets imply a reasonable tradeoff in sensitivity and specificity.
vt profile_indels -g /net/fantasia/home/atks/ref/vt/grch37/indel.reference.txt -r /net/fantasia/home/atks/ref/vt/grch37/hs37d5.fa run/final/all.genotypes.bcf -i 22:36000000-37000000 -f "~PASS"
vt profile_indels -g indel.reference.txt -r hs37d5.fa all.genotypes.bcf -i 22:36000000-37000000 -f "~PASS"
data set
data set
To normalize and remove duplicate variants:
To normalize and remove duplicate variants:
vt normalize mills.genotypes.bcf -r ~/ref/vt/grch37/hs37d5.fa | vt mergedups - -o mills.normalized.genotypes.bcf
vt normalize mills.genotypes.bcf -r hs37d5.fa | vt mergedups - -o mills.normalized.genotypes.bcf
and you will observe that 3994 variants had to be left aligned and 1092 variants were removed.
and you will observe that 3994 variants had to be left aligned and 1092 variants were removed.
