Просмотр исходного кода

更新 'tasks/arriba.wdl'

24
master
biochenglinliu 3 лет назад
Родитель
Сommit
5dd594fa3c
1 измененных файлов: 2 добавлений и 2 удалений
  1. +2
    -2
      tasks/arriba.wdl

+ 2
- 2
tasks/arriba.wdl Просмотреть файл



mkdir ./output/ mkdir ./output/
STAR \ STAR \
--runThreadN 32 \
--runThreadN 24 \
--genomeDir ${STAR_INDEX_DIR} \ --genomeDir ${STAR_INDEX_DIR} \
--genomeLoad NoSharedMemory \ --genomeLoad NoSharedMemory \
--readFilesIn ${fastq1} ${fastq2} \ --readFilesIn ${fastq1} ${fastq2} \
-t /arriba_v2.1.0/database/known_fusions_hg38_GRCh38_v2.1.0.tsv.gz \ -t /arriba_v2.1.0/database/known_fusions_hg38_GRCh38_v2.1.0.tsv.gz \
-p /arriba_v2.1.0/database/protein_domains_hg38_GRCh38_v2.1.0.gff3 -p /arriba_v2.1.0/database/protein_domains_hg38_GRCh38_v2.1.0.gff3
samtools index -@ 32 ./output/${sample_id}.Aligned.sortedByCoord.out.bam ./output/${sample_id}.Aligned.sortedByCoord.out.bam.bai
samtools index -@ 24 ./output/${sample_id}.Aligned.sortedByCoord.out.bam ./output/${sample_id}.Aligned.sortedByCoord.out.bam.bai
Rscript /arriba_v2.1.0/draw_fusions.R \ Rscript /arriba_v2.1.0/draw_fusions.R \
--fusions=./output/${sample_id}_fusions.tsv \ --fusions=./output/${sample_id}_fusions.tsv \

Загрузка…
Отмена
Сохранить