Skip to content
Snippets Groups Projects
Verified Commit c6865206 authored by Laurent Modolo's avatar Laurent Modolo
Browse files

hisat2: remove samtools report

parent 170a745c
No related branches found
No related tags found
No related merge requests found
......@@ -77,10 +77,8 @@ process mapping_fastq {
-1 ${reads[0]} \
-2 ${reads[1]} 2> \
${file_prefix}_ht2_mapping_report.txt \
| samtools view -@ ${task.cpus} -bS - 2>> \
${file_prefix}_ht2_samtools_view_report.txt \
| samtools sort -@ ${task.cpus} -o ${file_prefix}.bam - 2>> \
${file_prefix}_ht2_samtools_sort_report.txt
| samtools view -@ ${task.cpus} -bS - \
| samtools sort -@ ${task.cpus} -o ${file_prefix}.bam
if grep -q "Error" ${file_prefix}_ht2_mapping_report.txt; then
exit 1
......@@ -93,11 +91,8 @@ process mapping_fastq {
-x ${index_id} \
-U ${reads} 2> \
${file_prefix}_ht2_mapping_report.txt \
| samtools view -@ ${task.cpus} -bS - 2>> \
${file_prefix}_ht2_samtools_view_report.txt \
| samtools sort -@ ${task.cpus} -o ${file_prefix}.bam - 2>> \
${file_prefix}_ht2_samtools_sort_report.txt
| samtools view -@ ${task.cpus} -bS - \
| samtools sort -@ ${task.cpus} -o ${file_prefix}.bam
if grep -q "Error" ${file_prefix}_ht2_mapping_report.txt; then
exit 1
fi
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment