This is a test version of Biostars. For the public version, visit https://www.biostars.org.
SEED analyses in MEGAN: How to connect Seed database with blast data

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

Hi, did you ever figure this issue out?

0 answers

No answers yet.

Log in to answer this question.