From Genome Analysis Wiki
Jump to: navigation, search
Did I find interesting variants?
Did you see a variant at the position?
$GC/bin/tabix $OUT/vcfs/chr22/chr22.filtered.vcf.gz 22:36661906 | head -1
22 36661906 . A G 18 PASS DP=409;MQ=59;NS=62;AN=124;AC=2;AF=0.013827;AB=0.4065;AZ=-0.5287;
GT:DP:GQ:PL 0/0:4:28:0,12,65
Let's check the sequence data to confirm that the variant really exists
$GC/bin/samtools tview $IN/bams/HG01242.recal.bam $REF/human.g1k.v37.chr22.fa
* Type 'g' to go to a specific position
* Type 22:36661906 to move to the position
* Press arrows to move between positions
* Press 'b' if you want to color by base quality
* Press '?' for more help

Navigation menu