Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • fabrice
    Member
    • Oct 2009
    • 86

    #1

    Tophat 2.0.9 only output chr1-7 in junctions.bed

    The log file segment_juncs.log

    segment_juncs v2.0.9 (4112M)
    ---------------------------
    [samopen] SAM header is present: 84 sequences.
    Loading reference sequences...
    Loading chr1...done
    Loading chr2...done
    Loading chr3...done
    Loading chr4...done
    Loading chr5...done
    Loading chr6...done
    Loading chr7...done
    >> Performing segment-search:
    found 0 potential split-segment junctions
    found 0 potential small deletions
    found 0 potential small insertions
    >> Performing coverage-search:
    Map covers 110216535 bases
    Map covers 108217907 bases in sufficiently long segments
    Map contains 1983953 good islands
    98509451 are left looking bases
    98505476 are right looking bases
    Collecting potential splice sites in islands
    reporting synthetic splice junctions...
    Examining donor-acceptor pairings in chr1
    Examining donor-acceptor pairings in chr2
    Examining donor-acceptor pairings in chr3
    Examining donor-acceptor pairings in chr4
    Examining donor-acceptor pairings in chr5
    Examining donor-acceptor pairings in chr6
    Examining donor-acceptor pairings in chr7
    found 4066401 potential junctions
    Reported 4066401 total potential splices
    Reporting 0 potential deletions...
    Reporting 0 potential insertions...
    Reporting potential fusions...
  • fabrice
    Member
    • Oct 2009
    • 86

    #2
    cat tophat.log

    [2013-11-01 17:39:45] Beginning TopHat run (v2.0.9)
    -----------------------------------------------
    [2013-11-01 17:39:45] Checking for Bowtie
    Bowtie version: 2.1.0.0
    [2013-11-01 17:39:45] Checking for Samtools
    Samtools version: 0.1.19.0
    [2013-11-01 17:39:45] Checking for Bowtie index files (genome)..
    [2013-11-01 17:39:45] Checking for reference FASTA file
    [2013-11-01 17:39:45] Generating SAM header for /home/fabrice/genome/bowtie/human/hg19
    format: fastq
    quality scale: phred33 (default)
    [2013-11-01 17:40:15] Preparing reads
    left reads: min. length=12, max. length=41, 32543791 kept reads (935808 discarded)
    Warning: short reads (<20bp) will make TopHat quite slow and take large amount of memory because they are likely to be mapped in too many places
    Warning: you have only one segment per read.
    If the read length is greater than or equal to 45bp,
    we strongly recommend that you decrease --segment-length to about half the read length because TopHat will work better with multiple segments
    [2013-11-01 17:56:46] Mapping left_kept_reads to genome hg19 with Bowtie2
    [2013-11-01 18:09:16] Searching for junctions via segment mapping
    [2013-11-01 18:42:05] Retrieving sequences for splices
    [2013-11-01 18:44:37] Indexing splices
    [2013-11-02 21:57:29] Mapping left_kept_reads_unmapped to genome segment_juncs with Bowtie2 (1/1)
    [2013-11-04 06:42:31] Joining segment hits
    [2013-11-04 06:50:16] Reporting output tracks
    -----------------------------------------------
    [2013-11-04 07:46:44] A summary of the alignment counts can be found in /home/fabrice/projects/analysis/tophat/align_summary.txt
    [2013-11-04 07:46:44] Run complete: 2 days 14:06:58 elapsed

    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, 10-22-2026, 11:13 AM
    0 responses
    60 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 10-22-2026, 10:48 AM
    0 responses
    42 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 10-15-2026, 10:52 AM
    0 responses
    75 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 10-14-2026, 05:39 AM
    0 responses
    74 views
    0 reactions
    Last Post SEQadmin2  
    Working...