This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Where can I download "whole bacteria database"?

I need a database for whole bacteria genome to mapping. And It should be unique.

I have try the ncbi ftp://ftp.ncbi.nlm.nih.gov/genomes/genbank/bacteria/assembly_summary.txt but this class is not unique.

rna-seq
• 5,329 views • link
written by felix.kuo.1211 • felix.kuo.1211

How would you define unique? All the proteins are exactly similar? 95% identity of 95% of the proteome? 80%? Do you see where I'm going here?

• 1 views • link
written by Asaf Asaf

So you want to remove 100% identity. Do you need the transcripts or whole genomes? I think that there is no easy way to select a subset of the assemblies but you can download all of them and in a smart way screen duplicates (maybe using mummer or similar software)

• 1 views • link
written by Asaf Asaf

I need whole bacteria genomes. Its there any way to download it?

• 1 views • link
written by felix.kuo.1211 • felix.kuo.1211

1 answer

You might find this post useful.

• 1 views • link
written by Sej Modha Sej Modha

Ya! I did use this. But its not unique. It has many repeat.

• 1 views • link
written by felix.kuo.1211 • felix.kuo.1211

There is no non-redundant database of bacterial genomes. You would need to make one yourself. While there may be a few repeats many of the genomes are likely different strains of a particular species and so may appear redundant.

• 1 views • link
written by GenoMax GenoMax

Are you only looking for "refseq" genomes (ftp://ftp.ncbi.nlm.nih.gov/genomes/refseq/assembly_summary_refseq.txt actual data will be in ftp://ftp.ncbi.nlm.nih.gov/genomes/refseq/bacteria/ ) or all bacterial genomes?

In any case this would involve getting a list of available genomes following @5heikki's receipe (C: Download All The Bacterial Genomes From Ncbi ). Then parsing that list to make a non-redundant set. Downloading the fasta genomes for that set and making blast indexes.

• 1 views • link
written by GenoMax GenoMax

I did use the @5heikki's. And how to do "Then parsing that list to make a non-redundant set."? I want to use bowtie2.

• 1 views • link
written by felix.kuo.1211 • felix.kuo.1211

Log in to answer this question.

More posts like this