Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • mrfox
    Senior Member
    • Aug 2010
    • 103

    #1

    A few questions on Tophat

    Hi everyone,
    This is my first post. :-) I am a fresh user of Tophat and here I would like to ask you guys a few questions which confuse me a lot.
    I run Tophat 1.0.13 on my Illumina RNA-Seq data (75nt paired end, fly samples). I actually run Tophat without the option -g1 (called run1 below) and with -g1(called run2) respectively. I suppose that the uniquely mapped reads from run1 is the same set from run2. Thus I extract the mapped reads with the field "NH:i:1" from run1 and compared them with all the mapped reads from run2. However, I found that there are a lot more such reads in run1 that does not appear run2. So my question is:
    1) If I do not use "-g1", then how can I identify uniquely mapped reads? I guess my understanding that "NH:i:1" means "uniquely mapped" might not be correct.
    2)Then I run Tophat on 4 mouse samples, without or with "-g1". This time, I found that for two of the samples, in "accept_hits.sam", I did not find the field "NH", while I did find this field in another two.

    Hope I have addressed the questions clearly. Thank you all.:-)
  • mrfox
    Senior Member
    • Aug 2010
    • 103

    #2
    And in 2), my question is: how can I know when and which these optional fields will be contained in accpet_hits.sam? Thank you.

    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
    48 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 10-22-2026, 10:48 AM
    0 responses
    36 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...