
<html><head><meta http-equiv="Content-Type" content="text/html; charset=windows-1252"></head><body><div><div style="font-family: Calibri,sans-serif; font-size: 11pt;">I'm trying to get some specific annotations using seqinr. I recently started using R and also, too seqinr. So, the question is: I want to get some specific fields of the annotations section using first prepgetannot and then getAnnots. I tell you what i'm doing: <br><br>First (until here no problem)<br>s <- choosebank("genbank")<br>cavia <- query("cavia", "TID=10140 AND O=mitochondrion", socket = s$socket)<br><br>At this moment I can use getAnnots to get some annotations of cavia like getAnnots(cavia$req[1:2]) but my intention is for example to filter this and get only the definition or the organism. <br><br>I suppose that i could do this with the prepgetannot function but, really, honestly I don't know how to do that. My first idea was using prepgetannot and then getAnnots but this don't seems to work out. This is what I do. <br><br>prepgetannots(what = "OR", setfor = c("scan", "getannots"), socket = s$socket, verbose = FALSE)<br><br>And then use the getAnnots() function like before but I get the same result, that is, the complete annotation with all the fields. Obviously I'm doing something wrong but I cant realize what's my mistake. Thanks, and sorry for my english, this is my first time posting a question here. If you want I can post the output of the annotations but is simply the entire annotations of the accesions that I request<br><br>Thanks!!<br><br>Federico<br></div></div></body></html>