Commit dac4fb73 authored by Eric Coissac's avatar Eric Coissac

Close #15. Actually the bug in intron location was related to a

misinterpretation of the aragorn output format. Now tRNA, and intron
location are coherent with most of the locations extracted from genbank
file with one or 2 base pairs of difference.
parent 476e4070
......@@ -112,7 +112,8 @@ function emblTRNA(geneid,trna,loc,anti,intron,seq) {
intron=substr(intron,2,l-2);
split(intron,intronpos,",");
ib=intronpos[1];
ie=intronpos[2];
il=intronpos[2];
ie=ib+il-1;
match(loc,"[0-9][0-9]*");
gb=substr(loc,RSTART,RLENGTH);
if (complement==0) {
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment