SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Novoalign: No pairs found with simNGS data oiiio Bioinformatics 4 11-20-2012 10:55 AM
visualizing genomic data bogdan Bioinformatics 2 11-10-2011 04:42 PM
Why of number of F3 and F5 reads differ a lot of my solid data chenyao Bioinformatics 20 08-19-2011 07:32 AM
Seeking statistics on genomic data Fixee Bioinformatics 7 07-03-2011 10:44 PM
calling snps with samtools on novoalign data rcorbett Bioinformatics 1 02-05-2010 06:44 PM

Reply
 
Thread Tools
Old 10-26-2011, 12:45 AM   #1
lorddoskias
Junior Member
 
Location: Mars

Join Date: Jul 2011
Posts: 3
Default Novoalign with a lot of genomic data

For our usecase we intend to load Novoalign with the reference genomes for ALL Bacteria, Human, Viral, Fungi, Mouse, Honey Bee, Tazmanian Devil, Gorilla, Microbat, Megabat, Chicken and Swine.

After playing for a couple of hours with novoalign I was left with the impression that every specie should have a separate index file - that is 1 index for a particular virus which currently results in 2419 index files so after adding all the bacterial/archaeal and reamining species I will probably have around 10k or even more index files which basically would mean that if I have an arbitrary sequence which I want to align (without knowledge of what might be in it) I will have to invoke novoalign N times where N is the number of index files created?

If this is the case I was wondering if it is possible to multiplex all my source sequences into a single index file BUT retain the ability to identify to what specie/organism my sequence has mapped or this can be done purely by having separate species/organisms in separate index files?
lorddoskias is offline   Reply With Quote
Old 03-14-2012, 08:56 AM   #2
kopardev
Member
 
Location: VA, USA

Join Date: Oct 2011
Posts: 18
Default

Of course, you can do this. Just merge the fna files into a big fna file and create a novoalign index for it. You can always extract the gi numbers from alignment and get back the species/organism info with some scripting. Beware that novoindex does not allow the index to be bigger than 20GB in size. I once tried to do a novoindex for the entire nt, and I had to manually split it up into 3 fna files and so I had 3 indexes and had to invoke novoalign 3 times. 3 times for the entire nt .... no big deal!
kopardev is offline   Reply With Quote
Reply

Tags
alignment, classiffication, novoalign

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 01:23 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO