This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Align sequences without loading subject into RAM

Are there any tools for global alignment where the entire db is not loaded into ram? Blast is the only one I could find but it is too slow. I am looking to align RNA-seq reads to NCBI nt (~70gb) and I do not have that kind of ram.

Thanks!

alignment blast
• 3,047 views • link
updated by Ram Ram • written by tucanj tucanj

Where are you getting a 70Gb transcriptome from? o_O Unless you want to map to multiple species...?

• 0 views • link
written by John John

Want to do metatranscriptomic analysis by mapping to NCBI nt. I don't want to use DIAMOND as part of MEGAN because I have ncRNA as well.

• 0 views • link
written by tucanj tucanj

1 answer

Try Clustal Omega

I'm not sure about the way it loads dbs into RAMs, but it provides some pretty good features for global alignment tasks.

It might work for you to run it locally for sequences larger than 2000 or 2MB size by dividing the db into several seperate jobs against your set of RNA reads.

Good luck :)

• 0 views • link
updated by Ram Ram • written by Bara'a Bara'a

My understanding was that Clustal Omega was for aligning sequences against each other and not pairwise to a reference database. Is that correct?

• 0 views • link
written by tucanj tucanj

@tucanj ... you might be right, I'm not sure if there's any workaround in Clustal Omega command line regarding this.

I will check that for you :)

• 0 views • link
updated by Ram Ram • written by Bara'a Bara'a

Hello there @tucanj :)

This is the reply I got from the Clustal team.

Clustal Team Reply Regarding Aligning RNA-reads to NCBI nt

Apparently, what you were asking for is possible using Clustal Omega BUT NOT for RNA-reads unfortunately.

Hope that helps you one day in some other cases.

Good Luck :)

• 0 views • link
updated by Ram Ram • written by Bara'a Bara'a

Log in to answer this question.

More posts like this