With the current implementation, each GenomeRegion object has an attribute coords, which stores the start and end coordinates of the region. coords is set to the global coordinate (relative to the start of the genome) as defined in the reference coordinates file.
Since align and retrieve rely on local and global coordinates to find overlapping regions and generate output, I think it would be best to store each GenomeRegion'sglobal_coords and local_coords.
This change would further simplify the logic in align and retrieve and would resolve issue #20.
With the current implementation, each
GenomeRegion
object has an attributecoords
, which stores the start and end coordinates of the region.coords
is set to the global coordinate (relative to the start of the genome) as defined in the reference coordinates file.Since
align
andretrieve
rely on local and global coordinates to find overlapping regions and generate output, I think it would be best to store eachGenomeRegion's
global_coords
andlocal_coords
.This change would further simplify the logic in
align
andretrieve
and would resolve issue #20.