enforce the annotation of full-length rrna genes only

Former-commit-id: 0669e643dea2391c8370246cd3182b6328c068dc
Former-commit-id: 7820b800405b2134f2f0c2e27c4ef1ec73e1dc0d
This commit is contained in:
2015-10-12 07:26:18 -03:00
parent 01c5880ace
commit 921682f75a

View File

@ -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