Vegetable databases usage
Good morning, I'm using all nt database to analyse the vegetables. However, in some cases we have animal mixed with vegetables and want only to detect the vegetable species, normally, we would only expect to detect the vegetables only. But this is normal since we are using the all nt. So now, could anyone point me out, the way to detect only vegetables? Thanks for the help.
• 769 views
•
link
0 answers
No answers yet.
Log in to answer this question.
Cross posted at StackExchange
https://bioinformatics.stackexchange.com/questions/11384/vegetable-databases-usage
It is generally unadvisable to cross post the same question in multiple communities, as it leads to duplication of effort.
You should add a lot more information here, because your question is wildly unclear.
What are you working on? Do you really mean "meat", or perhaps "animal"?!
yes, i meant that, sorry!
Do you have a list of target species and know if the genomes are available? If you have a list of vegetables with genomes available, you can index the genomes and use fastqscreen for initial contamination screen. You can also use any target genome indices and run against fastqscreen to know the contamination for first few sequences in the fastq or randomly chosen reads in the fasta (AFAIK).
C: Restricting database by taxa in local blast searches
As asked this question does not make sense. What do you mean by
we have animal mixed with vegetables and want only to detect the vegetable species? There are many basic metabolism genes that would be common to animals/plants. I assumevegetable = plants?How are you restricting your
ntsearches? My answer in this thread: A: blastx local terminal would be the way to restrict your blast+ searches.Thank you genomax, this really helped me!
What data do you have? Reads?
You could run it through Kraken or Centrifuge and have the reads classified. You can then just retain the ones which have
Plantaehits. This isn't going to be fast though, even with a small-ish dataset.Hello anasofiamoreira94!
We believe that this post does not fit the main topic of this site.
Lets have information concentrated in that SE thread that received an answer already.
For this reason we have closed your question. This allows us to keep the site focused on the topics that the community can help with.
If you disagree please tell us why in a reply below, we'll be happy to talk about it.
Cheers!
Hello ATpoint, thank you and sorry for the trouble.