gggenes icon indicating copy to clipboard operation
gggenes copied to clipboard

Perfect alignment of gene island

Open mweberr opened this issue 4 years ago • 1 comments

Hallo , first of all this is a absolutely wonderful visualisation tool for gene maps. I want to draw gene maps with differing composition of genes due to mutations and frame shifts. This means that in each line different genes are missing. Still, I want the maps to be aligned according to overlapping genes.

Screenshot from 2019-07-05 14-56-16

Here all lines are aligned to genE, like in the example, but not line 1. Which should aligned to protD. How it would be possible to extend this alignment feature ?

mweberr avatar Jul 05 '19 13:07 mweberr

I'm afraid there's currently no good way to do this with gggenes and make_alignment_dummies(). Ultimately I think the solution will be to build in some kind of auto-alignment feature that attempts to find the optimal alignment across all the molecules. I've played around with implementing this feature but haven't yet found an algorithm that works in a reasonable amount of time.

wilkox avatar Jul 12 '19 06:07 wilkox