ensure we fail when conversion of hic pairs fails, don't track mis-sorted reads to save memory

Author: skoren

src/Snakefiles/8-hicPipeline.sm

@@ -386,8 +386,8 @@ rule transformBWA:
 cd 8-hicPipeline
 
 cat > ./transform_bwa.sh <<EOF
-#!/bin/sh
-set -e
+#!/usr/bin/env bash
+set -e -o pipefail
 {params.SAMTOOLS} view -q 1 ../{input.bwa_mapping} | {PYTHON} {VERKKO}/scripts/parse_sam_pairs.py  > ../{output.byread_mapping}
 EOF
 

src/scripts/parse_sam_pairs.py

@@ -33,25 +33,15 @@ def print_results(names):
 
 name = ""
 names = [ ]
-seen = {}
-out_of_order = 0
 
 for line in input_stream:
    line=line.split()
    if name == "":
       name = line[0]
    if name != line[0]:
-      seen[name] = 1
       print_results(names)
       name = line[0]
       names = [ ]
-   if name in seen:
-      print("Warning: read %s already seen but encountered it again, please confirm your bam file is sorted by read."%(name), file=sys.stderr)
-      out_of_order += 1
    names.append("%s\t%s"%(line[0], line[2]))
 
-if out_of_order > 1000:
-   print("Error: encountered too many unsorted reads (%d), exiting. Please confirm the input bam is sorted by read."%(out_of_order), file=sys.stderr)
-   sys.exit(1)
-
 print_results(names)