Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • tophat2 problem: No such file: left_kept_reads.m2g_um.candidates_and_unspl.bam

    I'm running into this problem while trying to get tophat2 to work. I'm running this on the test data that came with the software.

    I did get an accepted_hits.bam but it only contains the header.
    [2012-04-13 12:12:14] Beginning TopHat run (v2.0.0)
    -----------------------------------------------
    [2012-04-13 12:12:14] Checking for Bowtie
    Bowtie version: 2.0.0.5
    [2012-04-13 12:12:14] Checking for Samtools
    Samtools version: 0.1.18.0
    [2012-04-13 12:12:14] Checking for Bowtie index files
    [2012-04-13 12:12:14] Checking for Bowtie index files
    [2012-04-13 12:12:14] Checking for reference FASTA file
    [2012-04-13 12:12:14] Generating SAM header for /scratch/akurdogl/tophatForB37/Homo_sapiens.GRCh37.62
    format: fastq
    quality scale: phred33 (default)
    [2012-04-13 12:12:33] Reading known junctions from GTF file
    [2012-04-13 12:12:51] Preparing reads
    left reads: min. length=75, count=100
    right reads: min. length=75, count=100
    [2012-04-13 12:12:51] Using pre-built transcriptome index..
    [2012-04-13 12:13:30] Mapping left_kept_reads against transcriptome ensembl.63.genes with Bowtie2
    [2012-04-13 12:14:03] Mapping right_kept_reads against transcriptome ensembl.63.genes with Bowtie2
    [2012-04-13 12:14:36] Converting left_kept_reads.m2g to genomic coordinates (map2gtf)
    [2012-04-13 12:14:51] Converting right_kept_reads.m2g to genomic coordinates (map2gtf)
    [2012-04-13 12:15:05] Resuming TopHat pipeline with unmapped reads
    [2012-04-13 12:15:05] Mapping left_kept_reads.m2g_um against Homo_sapiens.GRCh37.62 with Bowtie2
    [2012-04-13 12:15:10] Mapping left_kept_reads.m2g_um_seg1 against Homo_sapiens.GRCh37.62 with Bowtie2 (1/3)
    [2012-04-13 12:15:13] Mapping left_kept_reads.m2g_um_seg2 against Homo_sapiens.GRCh37.62 with Bowtie2 (2/3)
    [2012-04-13 12:15:17] Mapping left_kept_reads.m2g_um_seg3 against Homo_sapiens.GRCh37.62 with Bowtie2 (3/3)
    [2012-04-13 12:15:21] Mapping right_kept_reads.m2g_um against Homo_sapiens.GRCh37.62 with Bowtie2
    [2012-04-13 12:15:25] Mapping right_kept_reads.m2g_um_seg1 against Homo_sapiens.GRCh37.62 with Bowtie2 (1/3)
    [2012-04-13 12:15:28] Mapping right_kept_reads.m2g_um_seg2 against Homo_sapiens.GRCh37.62 with Bowtie2 (2/3)
    [2012-04-13 12:15:32] Mapping right_kept_reads.m2g_um_seg3 against Homo_sapiens.GRCh37.62 with Bowtie2 (3/3)
    [2012-04-13 12:15:36] Searching for junctions via segment mapping
    [2012-04-13 12:17:53] Retrieving sequences for splices
    [2012-04-13 12:20:18] Indexing splices
    Warning: Encountered reference sequence with only gaps
    Warning: Encountered reference sequence with only gaps
    Warning: Encountered reference sequence with only gaps
    Warning: Encountered reference sequence with only gaps
    [2012-04-13 12:22:06] Mapping left_kept_reads.m2g_um_seg1 against segment_juncs with Bowtie2 (1/3)
    [2012-04-13 12:23:18] Mapping left_kept_reads.m2g_um_seg2 against segment_juncs with Bowtie2 (2/3)
    [2012-04-13 12:24:29] Mapping left_kept_reads.m2g_um_seg3 against segment_juncs with Bowtie2 (3/3)
    [2012-04-13 12:25:39] Joining segment hits
    [2012-04-13 12:27:55] Mapping right_kept_reads.m2g_um_seg1 against segment_juncs with Bowtie2 (1/3)
    [2012-04-13 12:29:08] Mapping right_kept_reads.m2g_um_seg2 against segment_juncs with Bowtie2 (2/3)
    [2012-04-13 12:30:19] Mapping right_kept_reads.m2g_um_seg3 against segment_juncs with Bowtie2 (3/3)
    [2012-04-13 12:31:30] Joining segment hits
    [2012-04-13 12:33:46] Reporting output tracks
    Traceback (most recent call last):
    File "/home/akurdogl/bin/tophat2/bin/tophat", line 3778, in ?
    sys.exit(main())
    File "/home/akurdogl/bin/tophat2/bin/tophat", line 3754, in main
    os.remove(m)
    OSError: [Errno 2] No such file or directory: '/scratch/akurdogl/tophatTest/test_data/ref37test//tmp/left_kept_reads.m2g_um.candidates_and_unspl.bam'

  • #2
    Check site tophat site for updated version of tophat2, same version number but different. Is the test set really only 100 reads? Might be insufficient?

    Comment


    • #3
      I have similar errors,

      [2012-05-24 16:47:41] Joining segment hits
      [FAILED]
      Error running 'long_spanning_reads':Loading insertions...Error: malformed insertion coordinate record

      Originally posted by Jon_Keats View Post
      Check site tophat site for updated version of tophat2, same version number but different. Is the test set really only 100 reads? Might be insufficient?

      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
      10 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, Yesterday, 06:07 PM
      0 responses
      9 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
      67 views
      0 likes
      Last Post seqadmin  
      Working...
      X