Hi,
In order to see if any gaps are closed, I am trying to map cDNA query sequences to an assembly. I have been using exonerate with the affine:overlap model - I have found this method to be very slow.
Are there any alternatives? Ideally I'd like to have a standard looking blast style output showing query/target sequnce alignment, percent identity etc., rather than using an assembly program such as AMOS minimus.
Thanks,
Jon.
•
written
by
jomaco
1 answer
exonerate --model est2genome --bestn 5 --showtargetgff 1 --showcigar 1 --hspfilter 1 --subopt 0 --dnawordlen 16
Using the above settings reduces the running time for a single query by a factor of 10 compared to default for my test case. The most efficient reduction is by increasing the wordlength from 12 to 16, sacrificing some sensitivity. Maybe this works also for your case?
Log in to answer this question.
-
Calculating E-value from Exonerate and choosing minimum percentage identity.
written by 8d8fc19e •
-
Exonerate fails to find match
written by Marina •
-
Polishing PACBIO corrected assembly using ONT and ILLUMINA data (.fastq)
written by K •
-
How to get percent identity from bam file ?
written by David
-
how do I get a wonderful genome draft? : gap closed step
written by vmicrobio
-
How to remove transcripts that have poor alignment scores in exonerate analysis
written by Ginsea Chen
-
Is there a way to find the percent similarity just like percent identity in BLAST?
written by Neha shri •
-
Issues with using Exonerate 2.2.0
written by ykel •
-
Error message exonerate, negative HSP score.
written by Myggan •
-
Global Alignment Of Short Reads
written by Pasta •
I am trying to align ESTs using method est2genome, but it appears to me that exonerate just is very slow. With standard settings it took over 80 minutes to align a single EST against a 600Mb genome, but I have 100.000 to align.