Merge pull request #223 from golobor/parse_refactor

Parse refactor
open2c · Apr 22, 2024 · 9d99660 · 9d99660
2 parents 6000a2d + f6cc28d
commit 9d99660
Show file tree

Hide file tree

Showing 3 changed files with 185 additions and 162 deletions.
diff --git a/pairtools/cli/parse.py b/pairtools/cli/parse.py
@@ -123,7 +123,7 @@
 @click.option(
     "--walks-policy",
     type=click.Choice(["mask", "5any", "5unique", "3any", "3unique", "all"]),
-    default="mask",
+    default="5unique",
     help="the policy for reporting unrescuable walks (reads containing more"
     " than one alignment on one or both sides, that can not be explained by a"
     " single ligation between two mappable DNA fragments)."

diff --git a/pairtools/lib/pairsam_format.py b/pairtools/lib/pairsam_format.py
@@ -100,6 +100,10 @@
     "dist_to_3",
     "seq",
     "mismatches",  # Format: "{ref_letter}:{mut_letter}:{phred}:{ref_position}:{read_position}"
+    "read_side",
+    "algn_idx",
+    "same_side_algn_count"
+
 ]
 
 DTYPES_EXTRA_COLUMNS = {
@@ -115,4 +119,7 @@
     "dist_to_3": int,
     "seq": str,
     "mismatches": str,
+    "read_side": int,
+    "algn_idx": int,
+    "same_side_algn_count": int,
 }