vcf文件annovar的注釋結果繪制瀑布圖maftools實例分析

今天給大家介紹一下vcf文件annovar的注釋結果繪制瀑布圖maftools實例分析。文章的內容小編覺得不錯，現在給大家分享一下，覺得有需要的朋友可以了解一下，希望對大家有所幫助，下面跟著小編的思路一起來閱讀吧。

vcf文件annovar的注釋結果繪制瀑布圖maftools分析

1.vcf文件annovar注釋:

table_annovar.pl  154.raw.somatic.vcf.gz   humandb/hg38/ -buildver hg38 -out  154 -remove -protocol \
 refGene,cosmic70,nci60,esp6500siv2_all,clinvar_20210501,1000g2015aug_all,1000g2015aug_eas,1000g2015aug_sas,avsnp150,gwasCatalog,ljb26_all,cytoBand,dgvMerged,phastConsElements100way,genomicSuperDups -operation g,f,f,f,f,f,f,f,f,f,f,r,r,r,r -nastring . -vcfinput

2. 提取必要的信息

 for i in *.hg38_multianno.txt
  do
      sample=`echo $i|awk -F '.' '{print $2}'`
      cut -f '1-10' $i|sed '1d'|sed "s/$/\t${sample}/">>all_sample.txt
  done
  
  sed -i '1s/^/Chr\tStart\tEnd\tRef\tAlt\tFunc.refGene\tGene.refGene\tGeneDetail.refGene\tExonicFunc.refGene\tAAChange.refGene\tTumor_Sample_Barcode\n/' all_sample.txt

3. 讀入數據，利用maftools繪圖

library(maftools)
  var_maf= annovarToMaf(annovar = "all_sample.txt", 
                              Center = 'NA', 
                              refBuild = 'hg38', 
                              tsbCol = 'Tumor_Sample_Barcode', 
                              table = 'refGene',MAFobj =T,
                              sep = "\t")

  plotmafSummary(maf = var_maf, rmOutlier = TRUE, addStat = 'median')

oncoplot(maf = var_maf, top = 30, fontSize = 12 ,showTumorSampleBarcodes = F )

以上就是vcf文件annovar的注釋結果繪制瀑布圖maftools實例分析的全部內容了，更多與vcf文件annovar的注釋結果繪制瀑布圖maftools實例分析相關的內容可以搜索億速云之前的文章或者瀏覽下面的文章進行學習哈！相信小編會給大家增添更多知識,希望大家能夠支持一下億速云！