[Seqinr-forum] searching for sequences from Aspergillus nidulans in 'genbank'

Coghlan, Avril A.Coghlan at ucc.ie
Sat Nov 14 18:58:44 CET 2009


Dear SeqinR forum,

Today I used SeqinR to retrieve sequences from the fungus Aspergillus
nidulans from the ACNUC 'genbank' database, using the commands:
> choosebank("genbank")
> query("anidulans","SP=aspergillus nidulans")
> anidulans$nelem
[1] 18948
This means that there were 18948 sequences from Aspergillus nidulans
found. 

As far as I understand it, the ACNUC 'genbank' database corresponds to
the NCBI Nucleotide database, is that right?

I also did a search directly of the NCBI Nucleotide database on the NCBI
website for Aspergillus nidulans sequences, by going to
http://www.ncbi.nlm.nih.gov/nucleotide/ and searching for "Aspergillus
nidulans"[ORGN]. That search found 29119 nucleotide sequences (12271 of
which are ESTs). 

I am wondering why there is a difference between the search that I did
of the ACNUC 'genbank' database, and on the NCBI Nucleotide Database
website?
I don't think it can be due to the ACNUC database missing some sequences
recently submitted to NCBI, as the ACNUC website says that the ACNUC
'genbank' database was very recently updated, on Nov 13, 2009 (from
http://pbil.univ-lyon1.fr/cgi-bin/get_relnum?db=GenBank&ident=1929541324
). 

I will be very grateful for your advice, as I would like to use the
SeqinR library for a bioinformatics practical for students, and want to
make sure I understand how it works. 

Regards,
Avril

Avril Coghlan
Department of Microbiology
University College Cork
Ireland





More information about the Seqinr-forum mailing list