千家信息网

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

发表于:2025-11-07 作者:千家信息网编辑
千家信息网最后更新 2025年11月07日,今天给大家介绍一下vcf文件annovar的注释结果绘制瀑布图maftools实例分析。文章的内容小编觉得不错,现在给大家分享一下,觉得有需要的朋友可以了解一下,希望对大家有所帮助,下面跟着小编的思路
千家信息网最后更新 2025年11月07日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实例分析相关的内容可以搜索之前的文章或者浏览下面的文章进行学习哈!相信小编会给大家增添更多知识,希望大家能够支持一下!

文件 注释 瀑布 结果 分析 实例 实例分析 内容 文章 更多 不错 必要 信息 就是 思路 数据 朋友 知识 跟着 面的 数据库的安全要保护哪些东西 数据库安全各自的含义是什么 生产安全数据库录入 数据库的安全性及管理 数据库安全策略包含哪些 海淀数据库安全审计系统 建立农村房屋安全信息数据库 易用的数据库客户端支持安全管理 连接数据库失败ssl安全错误 数据库的锁怎样保障安全 石化行业服务器安全加固系统 科协网络安全工作总结 服务器常用的冗余技术 3d电梯仿真软件开发 联系管理员远程服务器异常 无线网络技术顶级期刊 商业软件开发涉及法律规定 航空公司软件开发需求 服务器和计算机怎么互通 软件开发讲座报告 未连接服务器华为是因为 数据库中设置字段属性用什么视图 中间网络安全公司最新排名 手机安装应用无法连接到服务器 鸾曦互联网科技 中科软软件开发有限公司 网络技术突破 小程序属于服务器吗 数据库备份批处理命令 襄阳公司 亿网网络技术支持 网络安全跟it有什么区别 浦东新区新能源网络技术生产厂家 现在用什么做软件开发 济南市司法局鉴定机构数据库 账户密码注册服务软件开发 dz_数据库 河北网络服务器专业云主机 森巴庄园系统软件开发 现代计算机网络技术试卷 山西开源软件开发价格行情
0