database for fungus ?

Question

[Picasa]

In order to find some fungi in my dataset,

Can you recommend a complete fungi database ?

Thanks

Answer 1

[5heikki]

ftp://ftp.ncbi.nlm.nih.gov/genomes/genbank/fungi/

Answer 2

[Emily]

http://fungi.ensembl.org/index.html

Answer 3

[Picasa]

Thanks but do I need to download each species ? or there is a global database. (Im looking for fasta actually)

[GenoMax]

In future use the "Add comment" button when you are providing supplementary information or comments.
You would need to get the fasta genome files from the link @5heikki provided. Similar link for Ensembl is here. It may be possible to recover these from refseq_genome blast database (if you have that handy) and the blastdbcmd tool.

[5heikki]

Something like below to download the latest fungal assemblies from GenBank into a dir:

wget -q -O- ftp://ftp.ncbi.nlm.nih.gov/genomes/genbank/fungi/assembly_summary.txt | cut -f 20 -d $'\t' | awk '{FS="/"}{print $0"/"$6"_genomic.fna.gz"}' | xargs -n 1 wget

[Picasa]

Thanks for your help, can you explain your command ??

[5heikki]

wget -q -O- ftp://ftp.ncbi.nlm.nih.gov/genomes/genbank/fungi/assembly_summary.txt

Download ftp://ftp.ncbi.nlm.nih.gov/genomes/genbank/fungi/assembly_summary.txt and output to STDOUT.

cut -f 20 -d $'\t'

Parse the 20th column (columns separated by tabs) from STDIN and output to STDOUT (20th column of the assembly_summary.txt file has the base ftp urls for the genome assemblies).

awk '{FS="/"}{print $0"/"$6"_genomic.fna.gz"}'

Print the base ftp url and append: slash, the 6th field of the base ftp url, and _genomic.fna.gz. Output to STDOUT (fields separated by slashes, e.g. ftp://ftp.ncbi.nlm.nih.gov/genomes/all/GCA_001600535.1_JCM_30696_assembly_v001 becomes ftp://ftp.ncbi.nlm.nih.gov/genomes/all/GCA_001600535.1_JCM_30696_assembly_v001/GCA_001600535.1_JCM_30696_assembly_v001_genomic.fna.gz). As far as I know, this is a valid url construction method to all latest GenBank assemblies.

xargs -n 1 wget

Provide each constructed ftp url to wget one by one from STDIN:

Commands are glued together by pipes, i.e. "|"