Conversion of a GenBank format file into a glimmer-like one. Eukaryotic version.
This function reads a file in GenBank format and converts the features corresponding to CDS (Coding Sequences) into a format similar to glimmer program output. This function is specifically made for eukaryotic sequences, i.e. with introns.
gbk2g2.euk(gbkfile = system.file("sequences/ame1.gbk", package ="seqinr"), g2.coord = "g2.coord")
gbkfile |
The name of the GenBank file |
g2.coord |
The name of the output file |
This function returns the coordinates of the exons annotated in the GenBank format file.
A data frame with three columns will be written to the g2.coord
file. The first column corresponds to the name of the gene, given in the
GenBank file through the /gene
feature. The second and third
column contain the start and the stop position of the exon.
J.R. Lobry, A. Necşulea
citation("seqinr")
## Not run: gbk2g2.euk()
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.