SizeSeqFasta {GeneR} | R Documentation |
Give the size of a sequence from a fasta/embl/gbk file.
sizeSeqFasta(file,name=NA) sizeSeqEmbl(file,name=NA, index="ix") sizeSeqGbk(file,name=NA)
file |
File name of the bank |
name |
Name of the sequence in the bank (default: first sequence) |
index |
character. "ix": Name is accn of sequence. "id" : Le name is id of sequence. |
Size of sequence, NULL if error.
P.Durosay(C). L.Cottret(R)
## WARNING: Because NCBI changed the web service the seqNcbi() function ## needs to connect to, the function is broken. Please contact ## the maintainer of the GeneR package about this. Thanks! if (interactive()) { seqNcbi("NC_004718",file="NC_004718.gbk",type="G") sizeSeqGbk("NC_004718.gbk") # Or: sizeSeqGbk("NC_004718.gbk","NC_004718") # Idem with fasta format seqNcbi("NC_004718",file="NC_004718.fa",type="F") sizeSeqFasta("NC_004718.fa") # Or: sizeSeqFasta("NC_004718.fa","gi|30271926|ref|NC_004718.3|") } ## Not run: # With Embl download.file("http://bioinfo.hku.hk/sars/AY291451.seq", destfile="AY291451.seq") sizeSeqEmbl("AY291451.seq") ## End(Not run)