search

DEIB-GECO / GMQL

GMQL - GenoMetric Query Language
http://www.bioinformatics.deib.polimi.it/geco/
Apache License 2.0
18 stars 11 forks source link

PROJECT - coordinate region_update #28

Closed marcomass closed 7 years ago

marcomass commented 7 years ago

When PROJECT is used to modify region coordinates, despite the shema is not changed, the modified coordinates are set in two additional attributes, not described in the schema.

For testing, use the following query and compare the two materialize datasets: TEAD4_rep_broad = SELECT(dataType == 'ChipSeq' AND view == 'Peaks' AND setType == 'exp' AND antibody_target == 'TEAD4' AND (cell == 'ECC-1')) HG19_ENCODE_BROAD; MATERIALIZE TEAD4_rep_broad INTO TEAD4_rep_broad; TEAD4_rep_broad_less = PROJECT(region_update: left AS left + 100, right AS right - 100) TEAD4_rep_broad; MATERIALIZE TEAD4_rep_broad_less INTO TEAD4_rep_broad_less;

The same apply also when using start/stop; in addition in this case (only) the stop attribute is added in the schema, which should not be. Check the following: TEAD4_rep_broad_right = PROJECT(region_update: start AS start + 100, stop AS stop - 100) TEAD4_rep_broad; MATERIALIZE TEAD4_rep_broad_right INTO TEAD4_rep_broad_right;

marcomass commented 7 years ago

@akaitoua @pp86 Hi, could you please give top priority to solve this new issue? It is blocking the student who for his thesis needs this back to work as before the last update. Thanks!