bwa mem比对、去除PCR重复和统计最佳实践
获取软件的安装位置可以通过which bwa / which samtools 查看软件所在路径。
一、编写配置文件
config.ini 配置文件与程序位于同一目录,内容如下:
[Reference]
hg19 = /public/analysis/pipeline/reference/hg19/ucsc.hg19.fa
[Bin]
chrome_length = /public/analysis/pipeline/input/bin/chr_length
[Software]
bwa = ~/anaconda3/envs/bin/bwa
samtools = ~/anaconda3/envs/bin/samtools
chr_length文件内容如下:
chr1 249250621
chr2 243199373
chr3 19802
本文详细介绍了使用BWA进行比对,去除PCR重复,并进行统计的最佳实践。通过配置文件、主程序代码及运行指令,生成了.bam、.rmdup.bam、.coverage和.log文件。
订阅专栏 解锁全文
452

被折叠的 条评论
为什么被折叠?



