Hi all,
I've downloaded the latest version of bowtie2 (bowtie2-2.3.0-linux-x86_64.zip) and try to run it on my personal account on the cluster. Based on manual, I unzipped the file, changed directory to the bowtie2 folder and type "bowtie2", which the error says: "command not found". I tried with "./bowtie2", which the below error appeared.
/gfs-dist/seta/software/bowtie2-2.3.0/bowtie2-align-s: error while loading shared libraries: libtbb.so.2: cannot open shared object file: No such file or directory
(ERR): Description of arguments failed!
Exiting now ...
Could you please help me what is the problem and how to solve it?
Thank you
Thank you very much. As the link said "TBB to give superior thread scaling", could you please let me know if the speed of Bowtie without TBB library is low? However, the output of command "locate tbb" is
/usr/include/freetype2/freetype/ftbbox.h
, so it seems to be installed, yes?Thanks
I guess you could download the Bowtie2 source code and then change the makefile so that it points to your location of the TBB library. But at the moment I'am not sure what to change.
And by the way even without the TBB library Bowtie2 is a very fast aligner. But it always depends on what you are going to do. For small mappings that will not contain millions of reads it will make no difference since you will not recognize a gain in speed.
Hi seta,
Just curious, is Bowtie2 faster than BBMap for this data?
Hi Brian,
Unfortunately, I didn't compare the speeds of Bowtie2 and BBMap. Here, I need Bowtie2 for running tophat.
Don't use
TopHat
unless there is a pressing need. EvenTopHat
developers ask users to useHISAT2
, successor of TopHat.Thank you for your advice. No pressing. I just map reads on the genome (barley) for finding differentially expressed genes and some alternative splicing events. However, the HISAT2 manual says "HISAT2 is a fast and sensitive alignment program for mapping next-generation sequencing reads (whole-genome, transcriptome, and exome sequencing data) against the general human population (as well as against a single reference genome)", it means the program isn't appropriate for other genomes!? I'm working on barley (plant)
Could you please kindly tell me your idea about STAR especially when compared to HISAT2?
It's fine for any genome.