[Bioperl-l] Should coords be adjusted after removing alignment columns?
Chris Fields
cjfields at uiuc.edu
Tue Aug 14 17:19:59 UTC 2007
On Aug 14, 2007, at 12:13 PM, Sendu Bala wrote:
...
>
> Yes. For extra clarity, my thinking is that the correct answer is:
>
> a/6-19 tcgatcatcatc
> b/30-42 tcgatcatcatc
> c/51-63 tcgatcatcatc
> ...
> What you say in the rest of your discussion is valid and deserves
> some thought/discussion, but for now just getting the start and end
> correct, ignoring any issues with internal residues, seems like a
> no-brainer.
>
> For my own purposes that is all I need; having removed gaps I only
> need the start and end so I can take that region from each sequence
> and do a new alignment (for example).
It might be worth addressing the split location issue in the bug
report before it gets lost in the ether. Or maybe start a new one as
an enhancement request.
> BTW. Either my patch isn't quite perfect or there's another related
> bug I'm still tracking down. I'll commit when I've solved that,
> unless someone points out any mistakes in my thinking.
Sounds fine by me.
chris
More information about the Bioperl-l
mailing list