Commit 12b6c560 authored by Eric Coissac's avatar Eric Coissac

change tRNA and CDS annotations

parent ba00e56f
......@@ -371,7 +371,7 @@ blastx \
n = split($0,parts,",")
j = 1
for (i = 1; i <= n; i++) {
if (length(line) + length(parts[i]) > 79) {
if (length(line) + length(parts[i]) > 78) {
print line ","
line = "FT "
j = i
......
......@@ -160,11 +160,11 @@ function Unk(s) {
gname = (Ngene == 1 ? GeneName : GeneName "_" ++Igene)
locus = ""
Feature("gene", GeneLocation())
QQualifier("gene", gname)
QQualifier("locus_tag", locus)
if (FrameShift)
QQualifier("pseudogene","unknown")
# Feature("gene", GeneLocation())
#QQualifier("gene", gname)
#QQualifier("locus_tag", locus)
#if (FrameShift)
# QQualifier("pseudogene","unknown")
Feature("CDS", CdsLocation())
SQualifier("codon_start", 1)
......@@ -180,7 +180,7 @@ function Unk(s) {
}
QQualifier("product", Product)
QQualifier("inference", "similar to DNA sequence:" Simil)
QQualifier("inference", "org.annot -- detect pass:" PassType ":" PassInfo)
# QQualifier("inference", "org.annot -- detect pass:" PassType ":" PassInfo)
if (FrameShift==0) {
if (match(Translat,/\*/)>0) {
QQualifier("pseudogene","unknown")
......
......@@ -168,7 +168,7 @@ function emblTRNA(geneid,trna,loc,anti,intron,notes,seq) {
print "FT /gene=\""trna"\"";
print "FT /anticodon=\""anti"\"";
print "FT /product=\""product"("anti")\"";
print "FT /inference=\"Aragorn-1.2.38\"";
# print "FT /inference=\"Aragorn-1.2.38\"";
if (notes!="-")
print "FT /note=\""notes"\"";
......
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