Entering edit mode
3.8 years ago
poecile.pal
▴
50
Good afternoon,
I am trying to get sequence from genbank, using sequinr library.
library("seqinr")
choosebank("genbank")
A.thaliana <- query("A.thaliana", "sp=Arabidopsis thaliana AND K=rbcL")
This results in an error:
Error in readLines(socket, n = nelem, ok = FALSE): too few lines read in readLines
Traceback:
1. query("A.thaliana", "sp=Arabidopsis thaliana AND K=rbcL")
2. readLines(socket, n = nelem, ok = FALSE)
Please tell me what this may be related to and how to get the sequence? Thank you in advance!
Best regards, Poecile
Thank you very much! It works.