SEED analyses in MEGAN: How to connect Seed database with blast data

0

Entering edit mode

7.5 years ago

Bumble ▴ 10

I would like to use SEED to identify functions in a metegenomic dataset. I have downloaded the latest SEED mapping file here, translated sequence data to proteins, performed a blastp against RefSeq protein files downloaded from NCBI and imported the data into MEGAN. Taxonomy assignment works, but SEED does not show any hits. I assume that my sequences cannot be connected to the Seed database.

The blast command I use is

blastp -db Prot_DB.faa -query SampleProt.fasta -evalue 1e-5 -max_target_seqs 1 -outfmt 5 -out ProtBlastResults.txt -num_threads 12

What am I doing wrong? And is there a detailed tutorial on how to prepare data for SEED analyses in MEGAN? Any help would be greatly appreciated!

Metagenome MEGAN SEED blast • 2.3k views

ADD COMMENT • link updated 7.4 years ago by Biostar 20 • written 7.5 years ago by Bumble ▴ 10