Commit 0a9d14ac authored by Sander Bollen's avatar Sander Bollen
Browse files

gatk starts counting from 1

parent 391ca51a
......@@ -35,7 +35,7 @@ def split_genome(ref, approx_n_chunks=100):
chunk_size = tot_size//approx_n_chunks
chunks = []
for chrom_name, chrom_value in fa.records.items():
pos = 0
pos = 1
while pos <= len(chrom_value):
end = pos+chunk_size
if end <= len(chrom_value):
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment