change tRNA and CDS annotations
Former-commit-id: 12b6c5605f57940e215643b80c93ffbb48d5406e Former-commit-id: 18663d59e90e6d35b029d9087b66723487b8db1d
This commit is contained in:
@ -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"\"";
|
||||
|
||||
|
Reference in New Issue
Block a user