Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • TopHat "junction" error

    Hi I have some wey strange error.
    I run ~200 samples with tophat2 and 7 do not run with following error:
    ___________________________________________________________________
    tophat2 -p 8 --read-mismatches 4 -o thaln_HV ~/HV_consensus.fa unmapped.fq

    [2012-10-09 12:01:43] Beginning TopHat run (v2.0.0)
    -----------------------------------------------
    [2012-10-09 12:01:43] Checking for Bowtie
    Bowtie version: 2.0.0.6
    [2012-10-09 12:01:43] Checking for Samtools
    Samtools version: 0.1.18.0
    [2012-10-09 12:01:43] Checking for Bowtie index files
    [2012-10-09 12:01:43] Checking for reference FASTA file
    [2012-10-09 12:01:43] Generating SAM header for HV_consensus.fa
    format: fastq
    quality scale: phred33 (default)
    [2012-10-09 12:01:43] Preparing reads
    left reads: min. length=18, count=515457
    Warning: short reads (<20bp) will make TopHat quite slow and take large amount of memory because they are likely to be mapped to too many places
    [2012-10-09 12:01:58] Mapping left_kept_reads against HV_consensus.fa with Bowtie2
    [2012-10-09 12:02:40] Mapping left_kept_reads_seg1 against HV_consensus.fa with Bowtie2 (1/3)
    [2012-10-09 12:02:47] Mapping left_kept_reads_seg2 against HV_consensus.fa with Bowtie2 (2/3)
    [2012-10-09 12:02:54] Mapping left_kept_reads_seg3 against HV_consensus.fa with Bowtie2 (3/3)
    [2012-10-09 12:02:59] Searching for junctions via segment mapping
    Warning: junction database is empty!
    [2012-10-09 12:03:07] Retrieving sequences for splices
    [2012-10-09 12:03:07] Indexing splices
    Warning: Empty fasta file: 'thaln_HV/tmp/segment_juncs.fa'
    Warning: All fasta inputs were empty
    Error: Encountered internal Bowtie 2 exception (#1)
    Command: /bin/bowtie2-2.0.0-beta6/bowtie2-build thaln_HV/tmp/segment_juncs.fa thaln_HV/tmp/segment_juncs
    [FAILED]
    Error: Splice sequence indexing failed with err =1
    _____________________________________________

    strange thing is if i use --no-novel-junc option it runs ... even though i actually do not specify GTF file .....
    my reference is mixed fasta file for metagenomic study and it works fine for other samples. So i assume something is wrong with my *fq file, but they look fine, and reads phred quality is >= 30

    anyone can help with the mystery?

Latest Articles

Collapse

  • seqadmin
    Recent Innovations in Spatial Biology
    by seqadmin


    Spatial biology is an exciting field that encompasses a wide range of techniques and technologies aimed at mapping the organization and interactions of various biomolecules in their native environments. As this area of research progresses, new tools and methodologies are being introduced, accompanied by efforts to establish benchmarking standards and drive technological innovation.

    3D Genomics
    While spatial biology often involves studying proteins and RNAs in their...
    01-01-2025, 07:30 PM
  • seqadmin
    Advancing Precision Medicine for Rare Diseases in Children
    by seqadmin




    Many organizations study rare diseases, but few have a mission as impactful as Rady Children’s Institute for Genomic Medicine (RCIGM). “We are all about changing outcomes for children,” explained Dr. Stephen Kingsmore, President and CEO of the group. The institute’s initial goal was to provide rapid diagnoses for critically ill children and shorten their diagnostic odyssey, a term used to describe the long and arduous process it takes patients to obtain an accurate...
    12-16-2024, 07:57 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 01-09-2025, 04:04 PM
0 responses
11 views
0 likes
Last Post seqadmin  
Started by seqadmin, 01-09-2025, 09:42 AM
0 responses
18 views
0 likes
Last Post seqadmin  
Started by seqadmin, 01-08-2025, 03:17 PM
0 responses
28 views
0 likes
Last Post seqadmin  
Started by seqadmin, 01-03-2025, 11:18 AM
1 response
45 views
1 like
Last Post Tonia
by Tonia
 
Working...
X