Conversion of a GenBank format file into a glimmer-like one
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.
gbk2g2(gbkfile = "https://pbil.univ-lyon1.fr/datasets/seqinr/data/ct.gbk", g2.coord = "g2.coord")
gbkfile |
The name of the GenBank file |
g2.coord |
The name of the output file in glimmer-like format |
Partial CDS (either 5' or 3') and join in features are discarded.
The input file is returned invisibly.
J.R. Lobry
citation("seqinr")
oriloc
which uses glimmer-like files,
gbk2g2.euk
for eukaryotic sequences with introns.
## Not run: # need internet connection suppressWarnings(gbk2g2(g2.coord = "gbk2g2.test")) res <- read.table("gbk2g2.test") head(res) stopifnot(nrow(res) == 892) ## End(Not run)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.