biojulia / bigbed.jl Goto Github PK

I'm trying to write data to a BigBed file on each human chromosome, I defined my writer as such :

writer = BigBed.Writer(file, [(chr, length(genome[chr])) for chr in chrs])

With chrsbeing ["1", "2", ...]. Then I'm looping on chrs and do some write operations, but I'm getting a :

ArgumentError: disordered intervals

Because in the writer the chromosomes are getting reordered as : ["1", "10", "11"].

I managed to get around it by getting the chromosomes in the right order with:

ochrs = collect(values(writer.chromnames))[sortperm(collect(keys(writer.chromnames)))]

Would it be possible to have the writer keep the ordering it's given ? (maybe using an OrderedDict) or is there a good reason why it gets reordered ?

Alternatively a chromlist method to get the chromosome in the right order would help.