Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • coverageBed: command not found

    Dear all,

    I am mapping one end mRNA sequencing on TOPHAT with human reference genome hg19.

    It looks like script worked partially and It has generated .BAM file of 7 million MB. But geneexpression.txt are empty.

    Here is the error file:


    [Fri Oct 28 18:13:55 2011] Beginning TopHat run (v1.3.0)
    -----------------------------------------------
    [Fri Oct 28 18:13:55 2011] Preparing output location /home/ah644/tophat_output/intestinal/mRNA/C1055/
    [Fri Oct 28 18:13:55 2011] Checking for Bowtie index files
    [Fri Oct 28 18:13:55 2011] Checking for reference FASTA file
    [Fri Oct 28 18:13:55 2011] Checking for Bowtie
    Bowtie version: 0.12.7.0
    [Fri Oct 28 18:13:55 2011] Checking for Samtools
    Samtools Version: 0.1.17
    [Fri Oct 28 18:13:55 2011] Generating SAM header for /home/ah644/hg19andhs37/hg19
    [Fri Oct 28 18:13:56 2011] Preparing reads
    format: fastq
    quality scale: phred33 (default)
    [Fri Oct 28 18:13:56 2011] Reading known junctions from GTF file
    Left reads: min. length=76, count=151037396
    [Fri Oct 28 19:33:26 2011] Mapping left_kept_reads against hg19 with Bowtie
    [Fri Oct 28 21:46:41 2011] Processing bowtie hits
    [Fri Oct 28 23:43:18 2011] Mapping left_kept_reads_seg1 against hg19 with Bowtie (1/3)
    [Sat Oct 29 00:38:57 2011] Mapping left_kept_reads_seg2 against hg19 with Bowtie (2/3)
    [Sat Oct 29 01:34:48 2011] Mapping left_kept_reads_seg3 against hg19 with Bowtie (3/3)
    [Sat Oct 29 02:21:31 2011] Searching for junctions via segment mapping
    [Sat Oct 29 05:17:18 2011] Retrieving sequences for splices
    [Sat Oct 29 05:22:36 2011] Indexing splices
    [Sat Oct 29 05:23:14 2011] Mapping left_kept_reads_seg1 against segment_juncs with Bowtie (1/3)
    [Sat Oct 29 05:48:31 2011] Mapping left_kept_reads_seg2 against segment_juncs with Bowtie (2/3)
    [Sat Oct 29 06:11:28 2011] Mapping left_kept_reads_seg3 against segment_juncs with Bowtie (3/3)
    [Sat Oct 29 06:28:37 2011] Joining segment hits
    [Sat Oct 29 08:53:02 2011] Reporting output tracks
    -----------------------------------------------
    Run complete [16:42:32 elapsed]
    /opt/torque/mom_priv/jobs/80964.bulldogn-rocks.wss.yale.edu.SC: line 15: coverageBed: command not found


    Pipeline:

    #!/bin/bash
    # usage: qsub -v SAMPLE="name" intestinaltophat1lane.pbs
    #PBS -l walltime=100:00:00
    #PBS -m ae
    #PBS -M [email protected]
    ### Set the queue name, given to you when you get a reservation.
    #PBS -l nodes=1: ppn=8
    ROOT=/home/ah644/projects
    cd $ROOT/$SAMPLE
    ROOTB=/home/ah644/tophat_output/
    TOPHATOUTPUT=$ROOTB$SAMPLE
    mkdir $TOPHATOUTPUT
    SEQFILES=`ls *_sequence.txt`
    tophat -p 8 -r 50 --microexon-search -G /home/ah644/hg19andhs37/Homo_sapiens.GRCh37.62.gtf -o $TOPHATOUTPUT /home/ah644/hg19andhs37/hg19 $SEQFILES
    coverageBed -abam $TOPHATOUTPUT/accepted_hits.bam -b /home/ah644/hg19andhs37/hg19ensembl.BED > $TOPHATOUTPUT/geneexpression.txt
    cut -f 1,2,3,4,13,14,15,16 $TOPHATOUTPUT/geneexpression.txt > $TOPHATOUTPUT/temp
    cut -f 5 /home/ah644/hg19andhs37/hg19ensemblgeneinfo.txt > $TOPHATOUTPUT/temp2
    paste $TOPHATOUTPUT/temp $TOPHATOUTPUT/temp2 > $TOPHATOUTPUT/reducedgeneexpression.txt
    rm $TOPHATOUTPUT/temp*



    please help me in correcting the script.
    many thanks,
    Adam

  • #2
    You probably don't have coverageBed in your $PATH variable (Google it to find out more) You can add it in the script (Google "path variable")

    Comment


    • #3
      I have BEDtool software installed and this is how it shows.

      -bash-3.2$ cd programs/
      -bash-3.2$ cd BEDTools-Version-2.13.4/
      -bash-3.2$ ls
      bin data genomes LICENSE Makefile obj README.rst RELEASE_HISTORY src



      I also have BEDtools in my path.

      # .bash_profile

      # Get the aliases and functions
      if [ -f ~/.bashrc ]; then
      . ~/.bashrc
      fi

      # User specific environment and startup programs

      PATH=$PATH:$HOME/bin:/home/ah644:/home/ah644/programs/cufflinks-1.0.3.Linux_x86_64:/home/ah644/programs/samtools-0.1.17:/home/ah644/programs/samtools-0.1.17/misc:/home/ah644/programs/tophat-1.3.0.Linux_x86_64:/home/ah644/programs/BEDTools-Version-2.13.4:/home/ah644/programs/bowtie-0.12.7:/home/ah644/programs/maq-0.7.1:/home/ah644/programs/maq-0.7.1/scripts:/home/ah644/programs/plink-1.07-x86_64/

      export PATH



      but I am still getting the same error.
      Last edited by ibn.adam; 10-29-2011, 02:10 PM.

      Comment


      • #4
        The coverageBed executable is in the bin subdirectory of BEDTools, so you probably need to add /bin to the BEDTools entry in your path

        Comment


        • #5
          you mean I should change /home/ah644/programs/BEDTools-Version-2.13.4 to /home/ah644/programs/BEDTools-Version-2.13.4/bin

          Comment


          • #6
            Yes, that's what I meant.

            Comment


            • #7
              its working.
              thanks a lot and warm regards.
              Adam

              Comment

              Latest Articles

              Collapse

              • seqadmin
                Strategies for Sequencing Challenging Samples
                by seqadmin


                Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                03-22-2024, 06:39 AM
              • seqadmin
                Techniques and Challenges in Conservation Genomics
                by seqadmin



                The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                Avian Conservation
                Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                03-08-2024, 10:41 AM

              ad_right_rmr

              Collapse

              News

              Collapse

              Topics Statistics Last Post
              Started by seqadmin, Yesterday, 06:37 PM
              0 responses
              8 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, Yesterday, 06:07 PM
              0 responses
              8 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-22-2024, 10:03 AM
              0 responses
              49 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-21-2024, 07:32 AM
              0 responses
              66 views
              0 likes
              Last Post seqadmin  
              Working...
              X