|
|
@@ -17,7 +17,7 @@ task arriba{ |
|
|
|
|
|
|
|
mkdir ./output/ |
|
|
|
STAR \ |
|
|
|
--runThreadN 16 \ |
|
|
|
--runThreadN 24 \ |
|
|
|
--genomeDir ${STAR_INDEX_DIR} \ |
|
|
|
--genomeLoad NoSharedMemory \ |
|
|
|
--readFilesIn ${fastq1} ${fastq2} \ |
|
|
@@ -50,7 +50,7 @@ task arriba{ |
|
|
|
-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 |
|
|
|
|
|
|
|
samtools index -@ 16 ./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 \ |
|
|
|
--fusions=./output/${sample_id}_fusions.tsv \ |