worklog——pindel,GATK,delly使用

pindel:

需要文件:sort过的bam文件;.fa的参考文件;config.txt文件(主要是用于标准化bam文件);bam.bai文件;.fa.fai文件将这些文件均放在一个目录A下

步骤:

①进入文件A,在目录A中输入命令:

./pindel -f xxx.fa -i xxxx_config.txt  -c [chromosome]  -o [outputfilename]

②则在A目录下会出现输出文件

GATK:

工作原理:

(一)GATK三大步骤:原始数据的处理——变异检测——初步分析

(二)前期准备:①软件:须安装picard,bwa,htslib

②需要文件:fastq文件

(三)原始数据处理:

①用bwa mem对fastq文件进行过滤比对——得到sam文件

java -jar /home/lsw/xijiaoda/GATK/GenomeAnalysisTK.jar -T CountReads -R exampleFASTA.fasta -I exampleBAM.bam

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容

  • Spring Cloud为开发人员提供了快速构建分布式系统中一些常见模式的工具(例如配置管理,服务发现,断路器,智...
    卡卡罗2017阅读 134,969评论 19 139
  • 《ijs》速成开发手册3.0 官方用户交流:iApp开发交流(1) 239547050iApp开发交流(2) 10...
    叶染柒丶阅读 5,330评论 0 7
  • 第十章 使用序列数据 生物信息学的核心问题之一是处理大量的(通常定义糟糕或模糊)文件格式。久而久之,一些特定的简单...
    yangliunk1987阅读 5,090评论 3 53
  • linux资料总章2.1 1.0写的不好抱歉 但是2.0已经改了很多 但是错误还是无法避免 以后资料会慢慢更新 大...
    数据革命阅读 12,237评论 2 33
  • 报怨形体不美;从零开始;没有把握住自己;围城,生机勃勃的学生过几年变成平庸的人,不断的设计起点,不断人突破围城,...
    D051飞鹰阅读 181评论 0 0