Mapping locations (develop branch) doesn't seem to work
simonvh opened this issue · 2 comments
simonvh commented
>>> from genomepy.annotation import Annotation
>>> a = Annotation("hg38")
>>> a.map_locations(a.bed, to="Ensembl").head()
index start end name score strand thickStart thickEnd itemRrgb blockCount blockSizes blockStarts ensembl-name
0 chr1 17368 17436 ENST00000619216.1 0 - 17436 17436 0 1 68, 0, 1
1 chr1 29553 31097 ENST00000473358.1 0 + 31097 31097 0 3 486,104,122, 0,1010,1422, 1
2 chr1 30266 31109 ENST00000469289.1 0 + 31109 31109 0 2 401,134, 0,709, 1
3 chr1 30365 30503 ENST00000607096.1 0 + 30503 30503 0 1 138, 0, 1
4 chr1 34553 36081 ENST00000417324.1 0 - 36081 36081 0 3 621,205,361, 0,723,1167, 1
These are still chr
names.
simonvh commented
Oh! I now noticed that it is in the last column!
Shouldn't we replace the first column with the new name?
siebrenf commented
yea, looks like its indeed a bug.