首页
学习
活动
专区
工具
TVP
发布
精选内容/技术社群/优惠产品,尽在小程序
立即前往
  • 您找到你想要的搜索结果了吗?
    是的
    没有找到

    满分室间质评之GATK Somatic SNV+Indel+CNV+SV(下)性能优化

    #此处是原先Manta分析SV的步骤一,生成runWorkflow.py,因为这一不步速度很快,所以串行执行 rm -f ${result}/${sn}/runWorkflow.py python ${tools.manta} \ --normalBam ${result}/${sn}NC_marked.bam \ --tumorBam ${result}/${sn}_marked.bam \ --referenceFasta ${refs.hum} \ --exome \ --callRegions /opt/ref/projects/Illumina_pt2.bed.zip \ --runDir ${result}/${sn} # 对bam文件碱基质量校正的第二步,Normal & Tumor并行处理 ${tools.gatk} ApplyBQSR \ --bqsr-recal-file ${result}/${sn}_recal.table \ -L ${refs.interval} \ -R ${refs.hum} \ -I ${result}/${sn}_marked.bam \ -O ${result}/${sn}_bqsr.bam & ​ ​ ${tools.gatk} ApplyBQSR \ --bqsr-recal-file ${result}/${sn}NC_recal.table \ -L ${refs.interval} \ -R ${refs.hum} \ -I ${result}/${sn}NC_marked.bam \ -O ${result}/${sn}NC_bqsr.bam & ​ #原先QC步骤,获取insert size,Normal & Tumor并行 ${tools.gatk} CollectInsertSizeMetrics \ -I ${result}/${sn}_marked.bam \ -O ${result}/${sn}_insertsize_metrics.txt \ -H ${result}/${sn}_insertsize_histogram.pdf & ​ ​ ${tools.gatk} CollectInsertSizeMetrics \ -I ${result}/${sn}NC_marked.bam \ -O ${result}/${sn}NC_insertsize_metrics.txt \ -H ${result}/${sn}NC_insertsize_histogram.pdf & ​ # 运行manta SV分析 python ${result}/${sn}/runWorkflow.py -m local -j ${envis.threads} & ​ # 运行cnvkit CNV分析 ${tools.cnvkit} batch \ ${result}/${sn}_marked.bam \ --normal ${result}/${sn}NC_marked.bam \ --method hybrid \ --targets ${refs.bed} \ --annotate /opt/ref/refFlat.txt \ --output-reference ${result}/${sn}_reference.cnn \ --output-dir ${result}/ \ --diagram \ -p 0 & ​ #samtools统计测序深度 ${tools.samtools} depth -b ${refs.bed} ${result}/${sn}_marked.bam > ${result}/${sn}_marked.depth & ${tools.samtools} depth -b ${refs.bed} ${result}/${sn}NC_marked.bam > ${result}/${sn}NC_marked.depth & #samtools统计比对信息 ${tools.samtools} flagstat --threads ${envis.threads} ${result}/${sn}_marked.bam > ${result}/$

    01

    运维平台的建设思考-元数据管理(四)(r8笔记第16天)

    对于服务器的一些信息,如果数据量大了之后总是感觉力不从心,需要了解,但是感觉得到的这些信息不够清晰明了。 比如我们得到一台服务器,需要知道最基本的硬件配置,内存情况,磁盘空间情况,哪些磁盘空间问题需要关注,哪些磁盘空间问题可以忽略,swap的使用情况 如何,服务器的操作系统版本,内核版本,上面运行有几个实例,是否启用了ASM,甚至服务器运行了多少天呢,这些信息看起来非常琐碎,也可以通过脚本得 到,但是一直以来感觉都是比较笼统模糊。 今天使用shell脚本进行了简单的改进。 我们来看看基本的效果情况。有了这些

    015
    领券