Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • lj907
    Junior Member
    • Jul 2012
    • 9

    #1

    sh: samtools: command not found (Tophat)

    Dear all,

    I have been using tophat for my RNA-seq data. The command line I'm using is the standard: ./tophat -p 8 -G /path_to_genes.gtf/ -o <output_folder> /path_to_genome_index/ XXX.fastq

    Tophat runs happily, finds bowtie and samtools, checks for bowtie index, Fasta reference, generates SAM header... all goes well and I have 4 outputs: accepted_hits.bam deletions.bam insertions.bam and junctions.bam... However, when the sam version of accepted_hits.bam is going to be written, samtools fails and I get a message saying sh:samtools: command not found...
    Having gone through with the bam files till the very end (cufflinks then cuffdiff and R), I seem to get the analysis running but I don't understand why samtools fails during the mapping.... Tophat clearly recognises it at first.. and ./samtools <options> work nicely for other applications.
    Just to add I have used samtools-0.1.17.tar.bz2 (I'm on a iMac) untar ed it and cp ed it into my working path (I haven't got permission for /usr/bin), I have also cp ed libbam.a into my working directory...
    Can anyone help me
    Cheers
    lj907
  • flxlex
    Moderator
    • Nov 2008
    • 412

    #2
    Moved to a more appropriate forum

    Comment

    • maubp
      Peter (Biopython etc)
      • Jul 2009
      • 1544

      #3
      Perhaps you can add your working directory to the $PATH environment variable, if that is where your copy of samtools is?

      Personally, I use $HOME/bin for things I install under my own account, and put that on the $PATH. Most Linux tools can be installed in this way, if they use a configure script try adding --prefix=$HOME to this.

      Comment

      • lj907
        Junior Member
        • Jul 2012
        • 9

        #4
        Cheers Peter, I'm copying executables into $HOME/bin

        Comment

        Latest Articles

        Collapse

        • SEQadmin2
          16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
          by SEQadmin2





          NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification,...
          10-30-2026, 06:46 AM
        • SEQadmin2
          Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
          by SEQadmin2



          Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

          Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
          ...
          10-14-2026, 09:55 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by SEQadmin2, Today, 07:59 AM
        0 responses
        13 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 11-05-2026, 09:25 AM
        0 responses
        32 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 11-03-2026, 11:54 AM
        0 responses
        45 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 10-28-2026, 09:57 AM
        0 responses
        169 views
        0 reactions
        Last Post SEQadmin2  
        Working...