今天给大家介绍一下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实例分析相关的内容可以搜索亿速云之前的文章或者浏览下面的文章进行学习哈!相信小编会给大家增添更多知识,希望大家能够支持一下亿速云!
免责声明:本站发布的内容(图片、视频和文字)以原创、转载和分享为主,文章观点不代表本网站立场,如果涉及侵权请联系站长邮箱:is@yisu.com进行举报,并提供相关证据,一经查实,将立刻删除涉嫌侵权内容。