Not sure if this is the right forum for this.
Basically, I'm having some issues trying to find formatting guidelines for genbank. I'm trying to build a parser for genbank files, and looking at the lambda virus genome here there is a loc string complement(<23231).
What exactly does this mean? Does it refer to the complement between coordinate 1 and 23231? Also, is the reverse possible? i.e. complement(>23231).
•
written
by
jamietmorton •
1 answer
you can find good and ample info on the spec of the genbank (embl/ddbj) formats here
On your specific question: the '<'' (or '>') signs in the coordinate structure points to 'true coordinate is beyond/ before this position . Eg. if a gene is incompletely annotated (or because the genomic sequence is missing) you can denote it as <100 , meaning that you know the gene is truncated and that the correct start is located somewhere upstream of position 100 on that sequence
Log in to answer this question.
-
RNA-seq using rasflow - issues with annotation file
written by cdeantoneo31 •
-
Samtools ampliconclip hard clip not giving same results as soft clip
written by Kevin •
-
How to compare Sars-Cov-2 RefSeq to Influenza A RefSeq
written by bwillkie •
-
Does Trimmomatic merge paired-end reads into "single-end" files?
written by smo •
-
How do I create a gff(3) file from a genebank file?
written by torkel.loman •
-
RSEM with prokaryotes
written by csmatyi1 •
-
Generating consensus sequences without reference genome
written by thenameisizaak •
-
NCBI: how to obtain genomic context for arbitrary annotation release?
written by tsukanoffkirill •
-
RNA-Seq library types
written by compucurious •
-
How Do I View Gene Annotations For Non-Standard (E.G., Virus) Genome In IGV?
written by Obi Griffith
NCBI Genbak uses flat file format. Details can be found in the following link. Also, check the biopython GenBank module.
https://www.ncbi.nlm.nih.gov/Sitemap/samplerecord.html