Skip to content

Commit 0c20014

Browse files
authored
Merge pull request #52 from sspatel02/master
Update _extract_allc.py - fix BED formatting from allC file
2 parents 314ccf5 + 471b704 commit 0c20014

1 file changed

Lines changed: 15 additions & 3 deletions

File tree

ALLCools/_extract_allc.py

Lines changed: 15 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -84,7 +84,13 @@ def _extract_bed5_format(allc_line_list):
8484
# only chrom, pos, pos, mc, cov
8585
# mc and cov is binarized
8686
allc_line_list[4], allc_line_list[5] = binary_count(int(allc_line_list[4]), int(allc_line_list[5]))
87-
allc_line_list = [allc_line_list[i] for i in [0, 1, 1, 4, 5]]
87+
allc_line_list = [
88+
allc_line_list[0], # chr
89+
int(allc_line_list[1]) - 1, # convert position to 0-based for start
90+
int(allc_line_list[1]), # keep position 1-based for end
91+
allc_line_list[4], # mC
92+
allc_line_list[5], # cov
93+
]
8894
return "\t".join(map(str, allc_line_list)) + "\n"
8995

9096
else:
@@ -95,8 +101,14 @@ def _extract_allc_format(allc_line_list):
95101

96102
def _extract_bed5_format(allc_line_list):
97103
# only chrom, pos, pos, mc, cov
98-
allc_line_list = [allc_line_list[i] for i in [0, 1, 1, 4, 5]]
99-
return "\t".join(allc_line_list) + "\n"
104+
allc_line_list = [
105+
allc_line_list[0], # chr
106+
int(allc_line_list[1]) - 1, # convert position to 0-based for start
107+
int(allc_line_list[1]), # keep position 1-based for end
108+
allc_line_list[4], # mC
109+
allc_line_list[5], # cov
110+
]
111+
return "\t".join(map(str, allc_line_list)) + "\n"
100112

101113
out_format = str(out_format).lower()
102114
if out_format == "allc":

0 commit comments

Comments
 (0)