Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • xhuister
    Member
    • Apr 2010
    • 41

    #1

    TopHat Error: you must set the mean inner distance between mates with -r?

    I tried to use TOPHAT on Arabidopsis mRNA-seqs, they are not pair end reads. I've run the example data you provided, so the installation is all right.
    For my own data, I've built bowtie index "tair9_bwt" for Arabidopsis genome sequences. I have 6 fastq files of mRNA-seqs. I want to use TOPHAT to find new transcript boundaries in intergenic regions.

    I used the following command to test two of the six fq files, it could generate the 3 output files successfully.
    tophat ~/bowtie-0.12.5/tair9_bwt/tair9_bwt tail.fastq,test1.fastq

    But when I tried to use options like -i -I or -G or other options, it failed:
    -----------------------
    $ tophat ~/bowtie-0.12.5/tair9_bwt/tair9_bwt tail.fastq,test1.fastq -i 30 -I 15000
    Error: you must set the mean inner distance between mates with -r
    Traceback (most recent call last):
    File "/home/wux3/mrnaseq/tophat/bin/tophat", line 1635, in ?
    sys.exit(main())
    File "/home/wux3/mrnaseq/tophat/bin/tophat", line 1535, in main
    exit(1)
    TypeError: 'str' object is not callable
    ----------------------

    The description of -r option is:
    -r/--mate-inner-dist <int> This is the expected (mean) inner distance between mate pairs. For, example, for paired end runs with fragments selected at 300bp, where each end is 50bp, you should set -r to be 200. There is no default, and this parameter is required for paired end runs.

    But the mRNA-seqs I have are not paired end, anybody knows how to set this option on my data? And, why the "-r" option is not neccessary when no options used in the commandline?

    Thank you in advance.
  • xhuister
    Member
    • Apr 2010
    • 41

    #2
    Now, it is OK. Just used options before the fastq... -.-

    tophat -G TAIR9_GFF3_tophat -i 30 -I 15000 ~/bowtie-0.12.5/tair9_bwt/tair9_bwt tail.fastq,head.fastq
    Last edited by xhuister; 04-28-2010, 01:08 PM.

    Comment

    • maximilianh
      Member
      • Oct 2009
      • 15

      #3
      thanks for posting this, just saved me a lot of frustrations...

      Comment

      Latest Articles

      Collapse

      • SEQadmin2
        16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
        by SEQadmin2




        NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification, detection...
        Today, 06:46 AM
      • 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

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by SEQadmin2, 10-28-2026, 09:57 AM
      0 responses
      23 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 10-28-2026, 09:51 AM
      0 responses
      15 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 10-22-2026, 11:13 AM
      0 responses
      68 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 10-22-2026, 10:48 AM
      0 responses
      48 views
      0 reactions
      Last Post SEQadmin2  
      Working...