Kevin2015的个人博客分享 http://blog.sciencenet.cn/u/Kevin2015

博文

全外显子组生信分析流程-14-annovar_to_maftools

已有 3799 次阅读 2019-3-27 17:15 |个人分类:全外显子项目|系统分类:科研笔记

#!/usr/bin/bash

purpuse: combine annovar results and convert it to MAFtools format

#combination
infolder="/home/zhanghl/data_3/workshop/yelianhua/wes/annovar_result"
outfolder="/home/zhanghl/data_3/workshop/yelianhua/wes/annovar_result"

for sample in `cat /home/zhanghl/data_3/workshop/yelianhua/wes/sample_list`;do 
cut -f '1-10' $infolder/${sample}.hg19_multianno.txt | sed '1d'|sed "s/$/\t${sample}/" >>  $outfolder/all_annovar
done

grep -P "\texonic\t" $outfolder/all_annovar > $outfolder/all_annovar2
/home/zhanghl/local/perl_install/bin/perl  /home/zhanghl/supporting_softwares/to-maftools/to-maftools.pl $outfolder/all_annovar2         #将文件转换为maf格式 



#convertion
require(maftools)
var.annovar.maf <- annovarToMaf(annovar = "all_annovar3", Center = 'CSI-NUS', refBuild = 'hg19',tsbCol = 'Tumor_Sample_Barcode', table = 'refGene')
write.table(var.annovar.maf,file="var_annovar_maf",quote= F,sep="\t",row.names=F)
var_maf = read.maf(maf ="var_annovar_maf")




https://blog.sciencenet.cn/blog-2609994-1169970.html

上一篇:全外显子组生信分析流程-13-table_annovar
下一篇:安装Ubuntu 18.04.2 LTS流程
收藏 IP: 159.226.149.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-2-29 10:14

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部