温馨提示×

温馨提示×

您好,登录后才能下订单哦!

密码登录×
登录注册×
其他方式登录
点击 登录注册 即表示同意《亿速云用户服务条款》

vcf文件annovar的注释结果绘制瀑布图maftools实例分析

发布时间:2022-01-20 13:49:12 来源:亿速云 阅读:656 作者:柒染 栏目:开发技术

今天给大家介绍一下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实例分析相关的内容可以搜索亿速云之前的文章或者浏览下面的文章进行学习哈!相信小编会给大家增添更多知识,希望大家能够支持一下亿速云!

向AI问一下细节

免责声明:本站发布的内容(图片、视频和文字)以原创、转载和分享为主,文章观点不代表本网站立场,如果涉及侵权请联系站长邮箱:is@yisu.com进行举报,并提供相关证据,一经查实,将立刻删除涉嫌侵权内容。

AI