07 不同 vcf 注释结果转 maf

ANNOVAR

ANNOVAR 注释后拿到的结果,还需要做一个处理,就是在每一个 ${id}.hg38_multianno.txt文件中增加两列,为 Tumor_Sample_BarcodeMatched_Norm_Sample_Barcode,而最后的一列 Otherinfo 我们暂时丢弃掉(如果你需要就不要丢弃),添加上这两列信息之后再把所有的 vcf 文件合并成一个 annovar_merge.vcf。处理的方法如下:

cat config | while read id
do grep -v '^Chr' ./7.annotation/annovar/${id}.hg38_multianno.txt | cut -f 1-20 

07 不同 vcf 注释结果转 maf相关推荐

  1. c++注释快捷键_JAVA编程中你一定要掌握的“快捷键”

    导读 相信很多编程小白刚开始的时候,都对大神充满了,膜拜.因为他们可以随随便便敲出很多行代码,重点是他们没有动鼠标!这时候就有人问了,我怎么不动鼠标啊? 今天,我就要和大家分享一些很实用的Java快捷 ...

  2. 突变位点注释软件snpEff 使用详解

    通常在获得SNP位点以后,我们肯定想知道这些突变位点在基因组上出现的位置.突变的类型以及危害程度.在此,小编给大家推荐一种变异位点注释软件--snpEff. 1. 软件下载与安装 该软件的开发基于ja ...

  3. 13 对发生拷贝数变异的基因进行 KEGG 注释

    13 对发生拷贝数变异的基因进行 KEGG 注释 数据准备 首先还是同样的读入数据,进行一定的处理.我们同样用 VEP 注释后的 maf 文件,然后取出需要用到的几列 rm(list=ls()) op ...

  4. R包——maftools可视化神器

    参考:http://bioconductor.org/packages/release/bioc/vignettes/maftools/inst/doc/maftools.html 目录 介绍 准备 ...

  5. 生信初学者必备的基础知识

    如果你是一个生信初学者,又或者你是一个学临床的,为了发文章开始学生信,学了点数据挖掘,GEO,TCGA什么的,但是对很多专有名词不理解,对很多流程或者步骤云里雾里,可以看看这个教程:生物信息学最佳实践 ...

  6. R语言bioconductor包—maftools的使用

    maftools包分析和可视化大规模测序研究中的突变注释格式(MAF)文件.该软件包提供了各种功能来执行癌症基因组学中最常用的分析,并以最小的工作量创建功能丰富的可定制可视化. 1. 数据读入,生成M ...

  7. 版本控制 Git RPM打包

    版本控制 版本库 典型的客户/服务器系统 版本库是版本控制中心的核心 任意数量客户端 客户端通过写数据库分享代码 版本控制系统 集中式版本控制系统 开发者之间共用一个仓库(repository) 所有 ...

  8. Knowledge based | 如何获取某类肿瘤中所有已被报道的相关突变基因?

    当你的肿瘤高通量测序样本越多,或没有病人自身配对的癌旁与血液等非体细胞突变受累组织的配对时,那么WES或WGS所获得的体细胞突变数量将会非常多,所涉及的基因也很多,以至于突变Maf图无法阅读 (基因太 ...

  9. 使用 QTLtools 进行 PCA 分析

    3 使用 QTLtools 进行 PCA 分析 QTLtools 工具可以进行基因型的PCA分析,也可以进行表型的PCA分析,以下教程分别针对基因型和表型的PCA进行介绍: 3.1 对基因型进行PCA ...

最新文章

  1. 退出python命令行-退出python命令
  2. java web 的标准目录结构(zz)
  3. ubuntu11.04服务器安装
  4. MVC设计模式学习总结
  5. 基于FPGA的IIR滤波器设计
  6. dubbo-admin管理平台搭建
  7. 三国历史的教训读后感_历史上伟大的创新者的教训,你永远都不会忘记
  8. python进阶教程之异常处理
  9. nginx简介(轻量级开源高并发web服务器:大陆使用者百度、京东、新浪、网易、腾讯、淘宝等)(并发量5w)(一般网站apache够用了,而且稳定)...
  10. CodeVS 1576 最长严格上升子序列 (DP)
  11. nodejs/pomelo 使用 mongodb 连接 mongo时 出现
  12. Eclipse 修改项目名称
  13. 在linux中安装字体
  14. Codeforces 1102F Elongated Matrix 状压dp
  15. [转][darkbaby]任天堂传——失落的泰坦王朝(中)
  16. 小鸡啄米之React组件内部的API
  17. 编写.bat执行jar包
  18. Ubuntu 扩大/home磁盘分区
  19. alooa是华为什么型号_alooa是华为什么型号 华为alooa型号是华为畅享7Plus
  20. linux串口termios

热门文章

  1. IPV6 长度为46位
  2. 移动跨平台ReactNative动画组件Animated【14】
  3. 三阶魔方与四阶魔方总结
  4. 拜尔滤色拜尔滤色镜_如何在Windows 10上启用滤色器以更清晰地阅读屏幕
  5. (三)DDD上下文映射图——老师,我俩可是纯洁的男女关系!
  6. Apple Sign in with Apple(苹果授权登录PHP)
  7. 集合中的compare方法返回值说明
  8. 分析Ajax爬取今日头条街拍图片
  9. linux igb网卡,网卡驱动程序igb ixgbe
  10. Hadoop单节点设置