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
          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
        • SEQadmin2
          New Genomics Technologies Take Aim at Long-Standing Limits
          by SEQadmin2


          Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

          We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
          ...
          09-28-2026, 10:25 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by SEQadmin2, Yesterday, 11:13 AM
        0 responses
        16 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, Yesterday, 10:48 AM
        0 responses
        17 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 10-15-2026, 10:52 AM
        0 responses
        35 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 10-14-2026, 05:39 AM
        0 responses
        42 views
        0 reactions
        Last Post SEQadmin2  
        Working...