Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • keycard
    Junior Member
    • Dec 2010
    • 1

    #1

    A question about tophat -g/--max-multihits option

    If the --max-multihits option of tophat is set to 1, am I going to get the "best" alignment for the read? For example, will tophat pick the alignment for a read based on the best mapping score or the minimun number of mismatches? I am wondering how tophat works if there are multiple hits and the --max-multihits option is set to 1. Thanks.
  • brentp
    Member
    • Apr 2010
    • 72

    #2
    i think you need to use -g , not --max-multihits due to a bug in the arg parsing in tophat.
    but if you use -g 1, then you will get the best alignment and send -m 1 to bowtie.
    this: http://bowtie-bio.sourceforge.net/ma...e-options-best
    describes how the "best" is determined. it looks at mismatches, then quality.

    Comment

    • Boel
      Member
      • Oct 2009
      • 62

      #3
      I don't believe that is correct that setting --max-multihits to 1 gives you the BEST alignment for the read. According to the tophat manual --max-multihits 1 will trow away all reads that have more than 1 valid alignment. Hence you will get back sam files that only contain the reads that have exactly one alignment if this option is set to 1.

      Instructs TopHat to allow up to this many alignments to the reference for a given read, and suppresses all alignments for reads with more than this many alignments. The default is 40.

      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
      54 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 10-22-2026, 10:48 AM
      0 responses
      40 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 10-15-2026, 10:52 AM
      0 responses
      46 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 10-14-2026, 05:39 AM
      0 responses
      46 views
      0 reactions
      Last Post SEQadmin2  
      Working...