I have been able to use exonerate to map a large set of transctiptome-derived cDNA sequences (~10,000) to an assembled genome (~300 Mb). But I am only able to use the default parameters for exonerate. For instance, if I were to specity anything using parameters --bestn and --percent, I do not get any output, and the program seems to be stuck in its paths. If I include --showquerygff or --showtargetgff, the program runs with no problem.
I am using the exonerate version 2.2.0-x86_64 on a linux cluster, after securing a node with 16 processors.
Any help is appreciated.
(My first post to BioStars)
2 answers
Can you try this version:
http://www.ebi.ac.uk/~guy/exonerate/exonerate-2.4.0.tar.gz
I have experienced the same issue in the past. My assumption was that hits need to be stored in memory when you choose to do this type of filtering, and it could take a long time for output to be written to a file. Though, it's curious you get output when combining other options, so it may be a bug. My workaround has been to just set the alignment constraints and formatting you want, then process the output with a script. This may have been fixed in version 2.4, but if you are comfortable with writing a script, then it may be faster just to take that approach (and the results shouldn't be any different).
Yes, For --bestn option, I think you may be right: hits may be stored in memory first before being written out. However, the result for setting either of percent and bestn was the same for version 2.2: no output for a 24 hr job!
Our admin installed the version 2.4. Issue with --bestn still persists (which is OK), but --percent is working fine.
Thanks for your responses!
That is what I was referring to when I said I've seen the same issue. You will have to wait a long time, more than a day. Running the program with no filtering, and then post-processing the results is much faster, though maybe not as convenient for everyone. It does allow you to inspect all your results though, like hit number, not just the top hits.
Log in to answer this question.
-
How to convert exonerate C4 alignment output to gff?
written by Junior •
-
Transferring Gene Annotations from Assembly with Large Duplications
written by Charles Warden
-
Best strategy starting from allelic chromosomes instead of reference plus VCF
written by alfonsosaera •
-
FastTree error while constructing tree
written by Moses
-
Alignment to fastg assembly graph
written by grp2009 •
-
Minia: Pool allocation failed for 68766984 bytes (kmers alloc), mainbuffer is null?.
written by Beuss
-
Discovar De Novo assembly - "Attempt to allocate memory failed"
written by maxwhjohn1988
-
SPAdes insert size
written by treitlis •
-
Error message exonerate, negative HSP score.
written by Myggan •
