Determine my species from Isolates assemblies
0
0
Entering edit mode
2.3 years ago
davidmaimoun ▴ 50

Hello,

If I have a genome assemblies from isolates, and I want to know what it is. More accurately, I need to know the isolates are Neisseria or H.Influenza, or not. I would like to query them against Refseq database to see if there are matches. Blast doesn't works good if I query a whole genome right? What is the better way to do that?

Thank you!

blast • 1.0k views
ADD COMMENT
1
Entering edit mode

As a first approach try TYGS

ADD REPLY
0
Entering edit mode

Thank you Andres

I tried it but it is very slow

ADD REPLY
0
Entering edit mode

Hi David,

how many genomes do you have?

ADD REPLY
0
Entering edit mode

I have 8, We need to use tool that give an answer rapidly to know what to do with the patients The TYGS gave me a pretty good answer though, thank you for the tip. It take time but identified the right specie

ADD REPLY
1
Entering edit mode

May I ask why you are running assemblies first, if speed is paramount? You could e.g. run StrainPhlAn3 directly on the raw reads...you would need to run it with the --add_viruses parameter to include Influenza, though.

If all that interests you, is if Neisseria and/or Influenza or neither is contained in an isolate, you can certainly speed up your analysis by building a custom reference set with just those and a few controls and then use a kmer-based alignment-free classification.

ADD REPLY
1
Entering edit mode

I thought it would world only on assemblies.

Thank you very much for your help Matthias, this tool is great!

ADD REPLY

Login before adding your answer.

Traffic: 1679 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6