浏览代码

revise samtools sort -m

tags/v0.1.0
stead99 4 年前
父节点
当前提交
7c829c5944
共有 3 个文件被更改,包括 6 次插入6 次删除
  1. +2
    -2
      defaults
  2. +1
    -1
      tasks/fastp.wdl
  3. +3
    -3
      tasks/samtools.wdl

+ 2
- 2
defaults 查看文件

"stringtie_docker": "registry.cn-shanghai.aliyuncs.com/pgx-docker-registry/stringtie:v1.3.4", "stringtie_docker": "registry.cn-shanghai.aliyuncs.com/pgx-docker-registry/stringtie:v1.3.4",
"stringtie_cluster": "OnDemand bcs.a2.large img-ubuntu-vpc", "stringtie_cluster": "OnDemand bcs.a2.large img-ubuntu-vpc",
"ballgown_docker": "registry.cn-shanghai.aliyuncs.com/pgx-docker-registry/pgx-ballgown:0.0.1", "ballgown_docker": "registry.cn-shanghai.aliyuncs.com/pgx-docker-registry/pgx-ballgown:0.0.1",
"ballgown_cluster": "OnDemand bcs.a2.large img-ubuntu-vpc",
"ballgown_cluster": "OnDemand bcs.ps.g.large img-ubuntu-vpc",
"count_docker": "registry.cn-shanghai.aliyuncs.com/pgx-docker-registry/count:v1.0", "count_docker": "registry.cn-shanghai.aliyuncs.com/pgx-docker-registry/count:v1.0",
"count_cluster": "OnDemand bcs.a2.large img-ubuntu-vpc",
"count_cluster": "OnDemand bcs.ps.g.large img-ubuntu-vpc",
"insert_size":"8000", "insert_size":"8000",
"count_length": "150" "count_length": "150"
} }

+ 1
- 1
tasks/fastp.wdl 查看文件

## Trim ## Trim
if [ "${trim_adapter}" != 'true' ]; then if [ "${trim_adapter}" != 'true' ]; then
cp ${read1} ${sample_id}_R1.fq.gz cp ${read1} ${sample_id}_R1.fq.gz
cp ${read2} ${sample_id}_R2.fq.gz
cp ${read2} ${sample_id}_R2.fq.gz
else else
fastp --thread $nt -l ${length_required} -q ${qualified_quality_phred} -u ${length_required1} --adapter_sequence ${adapter_sequence} --adapter_sequence_r2 ${adapter_sequence_r2} --detect_adapter_for_pe --trim_front1 ${trim_front1} --trim_tail1 ${trim_tail1} --max_len1 ${max_len1} --trim_front2 ${trim_front2} --trim_tail2 ${trim_tail2} --max_len2 ${max_len2} -i ${read1} -I ${read2} -o ${sample_id}_R1.fq.gz -O ${sample_id}_R2.fq.gz -j ${sample_id}.json -h ${sample_id}.html fastp --thread $nt -l ${length_required} -q ${qualified_quality_phred} -u ${length_required1} --adapter_sequence ${adapter_sequence} --adapter_sequence_r2 ${adapter_sequence_r2} --detect_adapter_for_pe --trim_front1 ${trim_front1} --trim_tail1 ${trim_tail1} --max_len1 ${max_len1} --trim_front2 ${trim_front2} --trim_tail2 ${trim_tail2} --max_len2 ${max_len2} -i ${read1} -I ${read2} -o ${sample_id}_R1.fq.gz -O ${sample_id}_R2.fq.gz -j ${sample_id}.json -h ${sample_id}.html
fi fi

+ 3
- 3
tasks/samtools.wdl 查看文件

nt=$(nproc) nt=$(nproc)
/opt/conda/bin/samtools view -@ $nt -bS ${sam} > ${bam} /opt/conda/bin/samtools view -@ $nt -bS ${sam} > ${bam}
/opt/conda/bin/samtools sort -@ $nt -m 1000000000 ${bam} -o ${sorted_bam} /opt/conda/bin/samtools sort -@ $nt -m 1000000000 ${bam} -o ${sorted_bam}
/opt/conda/bin/samtools index -@ $nt ${sorted_bam}
/opt/conda/bin/samtools stats -@ $nt ${sorted_bam} > ${samstats}
/opt/conda/bin/samtools stats -@ $nt -i ${insert_size} ${sorted_bam} |grep ^IS|cut -f 2- > ${ins_size}
/opt/conda/bin/samtools index ${sorted_bam}
/opt/conda/bin/samtools stats ${sorted_bam} > ${samstats}
/opt/conda/bin/samtools stats -i ${insert_size} ${sorted_bam} |grep ^IS|cut -f 2- > ${ins_size}
>>> >>>


runtime { runtime {

正在加载...
取消
保存