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

博文

全外显子组生信分析流程-11-Mutect2

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

#!/usr/bin/bash

#purpose: identify SNV and InDels from paired tumor and normal samples
#set env
work_dir=/home/zhanghl/data_3/workshop/practice1_ovary
reference=/home/zhanghl/supporting_files/Homo_sapiens_glk_v37/human_b37_bundle/human_g1k_v37.fasta  
bwa_dir=/home/zhanghl/supporting_softwares/bwa-0.7.10
picard_dir=/home/zhanghl/supporting_softwares/picard
gatk=/home/zhanghl/supporting_softwares/GenomeAnalysisTK/gatk-4.0.8.1/gatk
cosmic=/home/zhanghl/supporting_files/Homo_sapiens_glk_v37/human_b37_bundle/b37_cosmic_v54_120711.vcf
snp=/home/zhanghl/supporting_files/Homo_sapiens_glk_v37/human_b37_bundle/dbsnp_138.b37.vcf
indel=/home/zhanghl/supporting_files/Homo_sapiens_glk_v37/human_b37_bundle/1000G_phase1.indels.b37.vcf

mutect2_function(){
    local normal_sample_name=$1
    local tumor_sample_name=$2
    local input=$3
    local output=$4
    
     $gatk Mutect2 \
     -R $reference \
     -I  ${input}/${tumor_sample_name}.sorted.dedup.recal.bam \
     -I  ${input}/${normal_sample_name}.sorted.dedup.recal.bam \
     -tumor ${tumor_sample_name} \
     -normal ${normal_sample_name} \
     -O ${output}/${tumor_sample_name}.mutect2.vcf
     
     $gatk FilterMutectCalls -V  ${output}/${tumor_sample_name}.mutect2.vcf -O  ${output}/${tumor_sample_name}.mutect2.filter.vcf
     }

start=$(date +%s)
time mutect2_function \
SRR5141045 \
SRR5141046 \
/input_dir \
/output_dir
end=$(date +%s)
dur=$((end-start))
echo $dur/60




https://wap.sciencenet.cn/blog-2609994-1169203.html

上一篇:全外显子组生信分析流程-10-depth_and_coverage_bedtools
下一篇:全外显子组生信分析流程-12-Control-Freec
收藏 IP: 159.226.149.*| 热度|

0

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

数据加载中...

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

GMT+8, 2024-3-1 13:36

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部