Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • broken pipe during tophat analysis

    Hi Does anyone know what problem with my analysis?

    ../../Program/tophat/tophat -p 12 ../../Database/bowtie-index/human_g1k_v37 ../../Data/MCC/MCC3154.fastq

    [Mon Nov 28 11:39:27 2011] Beginning TopHat run (v1.3.1)
    -----------------------------------------------
    [Mon Nov 28 11:39:27 2011] Preparing output location ./tophat_out/
    [Mon Nov 28 11:39:27 2011] Checking for Bowtie index files
    [Mon Nov 28 11:39:27 2011] Checking for reference FASTA file
    [Mon Nov 28 11:39:27 2011] Checking for Bowtie
    Bowtie version: 0.12.7.0
    [Mon Nov 28 11:39:27 2011] Checking for Samtools
    Samtools Version: 0.1.18
    [Mon Nov 28 11:39:27 2011] Generating SAM header for ../../Database/bowtie-index/human_g1k_v37
    [Mon Nov 28 11:39:29 2011] Preparing reads
    format: fastq
    quality scale: phred33 (default)
    Left reads: min. length=31, count=216599
    [Mon Nov 28 11:39:48 2011] Mapping left_kept_reads against human_g1k_v37 with Bowtie
    [Mon Nov 28 11:40:12 2011] Processing bowtie hits
    [Mon Nov 28 11:42:18 2011] Mapping left_kept_reads_seg1 against human_g1k_v37 with Bowtie (1/17)
    ....
    [Mon Nov 28 11:44:45 2011] Mapping left_kept_reads_seg17 against human_g1k_v37 with Bowtie (17/17)
    [Mon Nov 28 11:44:47 2011] Searching for junctions via segment mapping
    [Mon Nov 28 11:47:04 2011] Retrieving sequences for splices
    [Mon Nov 28 11:51:45 2011] Indexing splices
    [Mon Nov 28 11:51:49 2011] Mapping left_kept_reads_seg1 against segment_juncs with Bowtie (1/17)
    ....
    [Mon Nov 28 11:53:01 2011] Mapping left_kept_reads_seg17 against segment_juncs with Bowtie (17/17)
    [Mon Nov 28 11:53:05 2011] Joining segment hits
    sort: write failed: standard output: Broken pipe
    sort: write error
    [Mon Nov 28 11:55:11 2011] Reporting output tracks
    -----------------------------------------------

  • #2
    Copied and pasted from my previous answer to debug failed TopHat runs: http://seqanswers.com/forums/showthread.php?t=14140

    When debugging failed TopHat runs, you can also try to run the individual commands from the command line yourself.

    In the logs/ sub-directory of your output directory there should be a file called run.log which shows the commands that TopHat runs. There are also other log files in there - look to see if you can find errors (especially in the newest one or two files). Otherwise try running each command by itself from the directory you started TopHat in and see if and where empty files are made. The tmp stuff needed to do this is usually kept if a run fails.

    Chris

    Comment


    • #3
      it drops here

      /Volumes/Voyager/Package/Program/tophat/samtools view ./tophat_out/tmp/left_kept_reads.unspl.bam | sort -k1,1n -m -T./tophat_out/tmp/ - ./tophat_out/tmp/left_kept_reads_3192.sam.merge_sort.fifo | /Volumes/Voyager/Package/Program/tophat/samtools view -S -b - > ./tophat_out/tmp/left_kept_reads.candidates_and_unspl.bam


      after

      /Volumes/Voyager/Package/Program/tophat/samtools view -S -b - > ./tophat_out/tmp/left_kept_reads.candidates_and_unspl.bam


      What is wrong with samtools view?

      [samopen] SAM header is present: 84 sequences.
      [sam_read1] reference '6012' is recognized as '*'.
      Parse warning at line 223: mapped sequence without CIGAR
      Parse error at line 223: sequence and quality are inconsistent
      Abort trap

      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