3,045
edits
Changes
From Genome Analysis Wiki
Jump to navigationJump to searchSeqShop: Calling Your Own Genome, December 2014 (view source)
Revision as of 16:34, 12 December 2014
, 16:34, 12 December 2014→Exonic Variants NOT found by 1000G
zcat $OUT/wgs.filtered.rsid.anno.exon.vcf.gz | grep "EXTFILTER=NA,NA" | grep -v -w "0/0" |grep -v rs | perl -lane 'print "$_" if ( /ANNO=Stop_Gain/)' |grep -w PASS
zcat $OUT/wgs.filtered.rsid.anno.exon.vcf.gz | grep "EXTFILTER=NA,NA" | grep -v -w "0/0" |grep -v rs | perl -lane 'print "$_" if ( /ANNO=Stop_Gain/)' |grep -w PASS
Want to see this from the BAM file? Use samtools tview:
$GC/bin/samtools tview $SAMPLE/output/bams/$SAMPLE.recal.bam $GC/gotcloud.ref/human.g1k.v37.fa
Use 'g' & enter the Chr:Pos
If you want to know predicted functional significance of a particular variant, you can search by
If you want to know predicted functional significance of a particular variant, you can search by
