Commit 0669e643 by Eric Coissac

enforce the annotation of full-length rrna genes only

parent 865e2626
......@@ -44,11 +44,12 @@ pushTmpDir ORG.rrna
{loc=$5".."$6} \
($1 ~ /_RC$/) { \
loc="complement("loc")"} \
{match($1,"_..*S");\
(full==1) {match($1,"_..*S");\
rrna=substr($1,RSTART+1,RLENGTH-1);\
print "FT rRNA " loc; \
print "FT /gene=\"rrn"rrna"\""
print "FT /product=\""rrna" ribosomal RNA\""\
print "FT /product=\""rrna" ribosomal RNA\"";\
full=0
}'
popTmpDir
......
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