+ #TODO: add sequence info
+ #cistematic.core.retrieveSequence(genome, chrom, start, stop, sense="F", db="")
+ #ooooo. let's do this with PyGr
+
+ #TODO: add quality - max for each base, but default is to not include just like current
+ # rdsToFastq() sets everything to max
+