We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent 1d544aa commit 9eab57eCopy full SHA for 9eab57e
1 file changed
simulating/generate_molecules.py
@@ -97,7 +97,7 @@ def generate_molecules_from_bed(genome_file_path,
97
chr,pos = get_chr_pos_from_abs_postion(position=random.randint(1,genome_length), lengths=chrom_lengths)
98
gene = 'no_gene'
99
start = pos - math.floor(random.normalvariate(mu=mol_len/2, sigma=15))
100
- end = pos + mol_len
+ end = start + mol_len
101
if start <= 0 or end > len(genome[chr]):
102
continue
103
molecules.append((chr, gene, start, end))
0 commit comments