R包——maftools可视化神器

最新推荐文章于 2023-05-18 17:19:06 发布

VIP文章 Doris_xixi

最新推荐文章于 2023-05-18 17:19:06 发布

阅读量1.5w

点赞数 6

分类专栏： R语言

本文链接：https://blog.csdn.net/Doris_xixi/article/details/81204250

版权

参考：http://bioconductor.org/packages/release/bioc/vignettes/maftools/inst/doc/maftools.html

#读入annovar文件转换为maf——annovarToMaf

绘制maf文件的摘要——plotmafSummary

介绍

With advances in Cancer Genomics, Mutation Annotation Format (MAF) is being widely accepted and used to store somatic variants detected. The Cancer Genome Atlas Project has sequenced over 30 different cancers with sample size of each cancer type being over 200. Resulting data consisting of somatic variants are stored in the form of Mutation Annotation Format. This package attempts to summarize, analyze, annotate and visualize MAF files in an efficient manner from either TCGA sources or any in-house studies as long as the data is in MAF format.

随着癌症基因组学的进步，突变注释格式（MAF）被广泛接受并用于存储检测到的体细胞变体。癌症基因组图谱项目对30多种不同的癌症进行了测序，每种癌症类型的样本量超过200种。由体细胞变体组成的结果数据以MAF格式形式存储。只要数据采用MAF格式，该软件包就会尝试从TCGA源或任何内部研究中有效地汇总，分析，注释和可视化MAF文件。

准备

使用前要先将文件转换为maf格式,对于VCF格式文件，可以使用vcf2maf进行格式转换

maf文件包含的内容：

Mandatory fields: Hugo_Symbol, Chromosome, Start_Position, End_Position, Reference_Allele, Tumor_Seq_Allele2, Variant_Classification, Variant_Type and Tumor_Sample_Barcode.
Recommended optional fields: non MAF specific fields containing VAF (Variant Allele Frequecy) and amino acid change information.

格式转换

#将突变结果进行注释，得到txt文件
for i in *.somatic.anno;do perl ~/software/Desktop/annovar/table_annovar.pl $sra_file /home/yang.zou/database/humandb_new/ -buildver hg19 -out variants --otherinfo -remove -protocol ensGene -operation g -nastring NA -outfile;done

#然后将所有.hg19_multianno.txt文件添加一列填入文件名前缀并将所有txt文件拼接成一个文件，提取出含有外显子的信息
for i in *.hg19_multianno.txt;do sed '1d' $

最低0.47元/天解锁文章

Doris_xixi

关注

6
点赞
踩
27

收藏

觉得还不错? 一键收藏
3
评论
R包——maftools可视化神器

参考：http://bioconductor.org/packages/release/bioc/vignettes/maftools/inst/doc/maftools.html目录介绍准备格式转换总体分析框架maftools安装正式处理#读入annovar文件转换为maf——annovarToMafDescriptionUsageArgume...
复制链接

扫一扫