diff --git a/bin/mapped_2hic_dnase.py b/bin/mapped_2hic_dnase.py index ff593666f03c5b82928ced06b941be296d6169ea..c417907a4483bd33bc6a6e0ceeb752d5c3d32121 100755 --- a/bin/mapped_2hic_dnase.py +++ b/bin/mapped_2hic_dnase.py @@ -362,7 +362,6 @@ if __name__ == "__main__": if cur_handler is not None: if not r1.is_unmapped and not r2.is_unmapped: - ##reorient reads to ease duplicates removal or1, or2 = get_ordered_reads(r1, r2) or1_chrom = samfile.get_reference_name(or1.reference_id) diff --git a/bin/mapped_2hic_fragments.py b/bin/mapped_2hic_fragments.py index cc0e40b472dfb764ba1bf2f550c00c79bce7bd3f..a1298716ce83c4e99931534357c03634fd5d8301 100755 --- a/bin/mapped_2hic_fragments.py +++ b/bin/mapped_2hic_fragments.py @@ -619,7 +619,6 @@ if __name__ == "__main__": r2_chrom = None if r1_resfrag is not None or r2_resfrag is not None: - interactionType = get_interaction_type(r1, r1_chrom, r1_resfrag, r2, r2_chrom, r2_resfrag, verbose) dist = get_PE_fragment_size(r1, r2, r1_resfrag, r2_resfrag, interactionType) cdist = get_cis_dist(r1, r2) diff --git a/bin/mergeSAM.py b/bin/mergeSAM.py index 82ab8c34d979240d12de57b1f2510a54a79d2ca9..d670ec2c6a94b89035ab417b4146806099c2c5aa 100755 --- a/bin/mergeSAM.py +++ b/bin/mergeSAM.py @@ -71,7 +71,6 @@ def get_read_name(read): def sam_flag(read1, read2, hr1, hr2): - f1 = read1.flag f2 = read2.flag