I just got hold of the iSAAC aligner from Illumina, but it seems to use a very different type of command line than most of the other aligners I have used, particularly when specifying a FASTQ. Does anyone have a command line that has been working for them?
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
-
Thanks I'm downloading it now, although I wonder if this download will finish before my indexing command...
Even though it looks like the downloadable ref is based on the same UCSC hg19 I am using, I can't help but feel a little loss of experimental control when it comes to comparisons against the other software I ran
Comment
-
After downloading the reference I'm now hitting this error:
2013-06-10 17:52:33 [0x11d36a0] __gnu_parallel:sort for 64000000 elements
terminate called recursively
terminate called recursively
./isaac-run.sh: line 5: 40176 Aborted (core dumped) isaac-align --reference-genome IsaacIndex/sorted-reference.xml -b isaac/ --base-calls-format fastq -m 12 --verbosity 3
Has anyone put a successful run through with Isaac yet?
Comment
-
Has anyone managed to get "isaac-align" to work stand-alone on the command line and if so on what OS?
"Isaac-align"er is one of the components of HiSeq Analysis Software (HAS) and I have managed to get it to work when run as a part of a HAS run (but only under SGE, no luck with LSF).
Trying to get it to work stand-alone has been an exercise (Oiiio: I could not get it to work with the options you noted, on individual files etc). I thought I had managed to to it to work using the following command line (on a rocks v. 5.4.2 cluster)
Code:isaac-align -b /Flow_cell_folder/Data/Intensities/BaseCalls -r /path_to/ISAAC_genomes/hg19/Homo_sapiens/UCSC/hg19/Sequence/IsaacIndex/sorted-reference.xml -o /output_dir -m 48
Code:Error: boost::exception: Throw location unknown (consider using BOOST_THROW_EXCEPTION) Dynamic exception type: boost::exception_detail::clone_impl<boost::exception_detail::current_exception_std_exception_wrapper<std::bad_alloc> > std::exception::what: std::bad_alloc [boost::tag_original_exception_type*] = St9bad_alloc
Last edited by GenoMax; 07-12-2013, 05:56 AM.
Comment
-
Hi GenoMax,
if you managed to run it under HAS, why don't you use the command line generated by HAS? That command line can be found either on the standard output, or at the beginning of the log file "Logging/Isaac.stderr.txt" that you will find in the HAS output directory.
Regarding the specific behaviour that you observed when trying to run directly from the command line, there is clearly an issue with the memory settings, but it is difficult to diagnose without more information. Could you provide the output of:
Code:ulimit -a tail -n 9999 /proc/sys/vm/* cat /proc/meminfo
Thanks,
Come
Comment
-
Originally posted by craczy View PostHi GenoMax,
if you managed to run it under HAS, why don't you use the command line generated by HAS? That command line can be found either on the standard output, or at the beginning of the log file "Logging/Isaac.stderr.txt" that you will find in the HAS output directory.
Come
If that is not possible (or is not in the scope of isaac package) then it would be best to have that spelled out clearly.
I am not sure if you are @Illumina/UK (part of the isaac development team?). If you are, then we could take this conversation temporarily offline to omit thread-clogging info that may be specific to problem I am noticing.
Comment
-
I will have to look through my results from back last year. I have since used HAS several time though only for human data.
Before you start be aware that ISAAC/HAS was designed to take over an entire server (using *all* hardware resources on it). It is not going to run on a cluster.
Comment
-
We had the command on this GCAT page work:
Comment
-
Is anyone running current git version? It still seems to be a problem running isaac-aligner on its own.
Currently it is always stuck here with memory issues (well, not necessarily the amount, I've tried to run on 512GB machines, nor limits).
isaac part of HAS is quite old.
Comment
-
@Sven: Are you getting an error or does it seem to just hang?
I have not successfully run isaac-aligner as a standalone app in the past (going to upgrade to the latest version from git now, thanks for the heads-up) since there was no support available.
Comment
Latest Articles
Collapse
-
by seqadmin
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...-
Channel: Articles
04-22-2024, 07:01 AM -
-
by seqadmin
Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...-
Channel: Articles
04-04-2024, 04:25 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Yesterday, 08:47 AM
|
0 responses
12 views
0 likes
|
Last Post
by seqadmin
Yesterday, 08:47 AM
|
||
Started by seqadmin, 04-11-2024, 12:08 PM
|
0 responses
60 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
60 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
||
Started by seqadmin, 04-10-2024, 09:21 AM
|
0 responses
54 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 09:21 AM
|
Comment