Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • miaom
    Member
    • Oct 2012
    • 12

    How to set Tophat insert size

    correctly setting 2 parameters in Tophat:
    --mate-inner-dist
    --mate-std-dev
    should be very important, but i have no idea how to set them without running bowtie (or bowtie2) -> samtools -> picard tools first, which seems a little bit redundant for Tophat.

    Does any one have suggestions?
    Thanks in advance!!
  • gringer
    David Eccles (gringer)
    • May 2011
    • 845

    #2
    Those are precisely the options that you apply to the tophat command line to change the expected insert size distribution for Tophat. Tophat takes as arguments the bowtie index (and optionally transcriptome GTF file) together with raw reads in FASTQ format, calling bowtie/bowtie2 internally. You shouldn't need to directly run the other programs in order to use Tophat:

    Code:
    Usage: tophat [options]* <genome_index_base> \
      <reads1_1[,...,readsN_1]> [reads1_2,...readsN_2]
    Perhaps you could rephrase your question. What have you tried (command line is useful), and why do you think it doesn't give you the correct results?
    Last edited by gringer; 12-01-2013, 03:28 PM.

    Comment

    • miaom
      Member
      • Oct 2012
      • 12

      #3
      I think the parameters (--mate-inner-dist, --mate-std-dev) are for Tophat only. Bowtie does not need to specify them.
      Running bowtie can generate SAM files, which contain information about fragment size. Without pre-running bowtie (or maybe tophat), how can I know the values for --mate-inner-dist and --mate-std-dev to feed in Tophat?

      Comment

      • gringer
        David Eccles (gringer)
        • May 2011
        • 845

        #4
        Those are experimental values. You should have a good idea of your paired-end distance before doing any mapping. If you can get a good enough idea of that based on genomic (or transcriptomic) mapping distance by running bowtie, maybe you don't need to be running tophat at all....

        Comment

        Latest Articles

        Collapse

        • SEQadmin2
          Cancer Drug Resistance: The Lingering Barrier to Rising Survival
          by SEQadmin2



          Cancer survival rates have significantly increased in the last few decades in the United States, reaching a combined 70% 5-year survival rate by 2021. Behind this number, there are years of research to find new therapies, drug targets, and early detection methods. But there is one core challenge that keeps slowing down these advances, and it’s about drug resistance.

          There is no single reason why many patients don’t respond to treatment as expected. Cancer is...
          Today, 05:17 AM
        • GATTACAT
          Reply to Nine Things a Sample Prep Scientist Thinks About Before Sequencing
          by GATTACAT
          Love this - good data definitely starts from good input, and poor input can only give relatively poor data. I particularly like the mention of Nanodrop/absorbance based methods for quantification. It's such a toss up if you'll get an accurate reading or what amounts to a randomly generated number, and a lot of library/sequencing related issues can be traced back to poor quant.
          07-01-2026, 11:43 AM
        • SEQadmin2
          Nine Things a Sample Prep Scientist Thinks About Before Sequencing
          by SEQadmin2


          I’m not a sequencing expert. I’m a purification scientist who uses NGS to evaluate workflows my group develops. With this perspective, we think about the sample first and the NGS workflow second. The sequencer is an exceptionally honest reporter, but it can only report on what you give it, so whether you get clean, interpretable data from an NGS workflow is largely determined before you begin.

          Here are nine questions we think about, in roughly the order they matter, before...
          06-18-2026, 07:11 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by SEQadmin2, Today, 10:08 AM
        0 responses
        5 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, Yesterday, 11:05 AM
        0 responses
        7 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 07-02-2026, 11:08 AM
        0 responses
        30 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 06-30-2026, 05:37 AM
        0 responses
        28 views
        0 reactions
        Last Post SEQadmin2  
        Working...