Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • westeros
    Junior Member
    • Sep 2013
    • 4

    Bowtie2 mapping and paired-end reads

    Hi all,

    I have started using bowtie2 to map 100bp Illumina reads to a reference genome. Our ultimate goal is to have gene expression estimates.

    I have two questions:

    1) Bowtie2 use a score to define if a read is aligned or not. What is the best setting to call a read uniquely aligned ? Most papers using small reads allow a few mismatches and the threshold is pretty easy to grasp. But how do we handle it using scores ?

    2) What is the best way to use paired-end mappings and read counts ? Is it correct to count the two mates as two reads mapped to a gene ?

    Thanks a lot.
  • arcolombo698
    Senior Member
    • Nov 2013
    • 142

    #2
    Score Threshold

    Hello. I am interested in the similar topic regarding score threshold, but I am not sure how to set it.

    I would guess the following from the manual

    For an alignment to be considered "valid" (i.e. "good enough") by Bowtie 2, it must have an alignment score no less than the minimum score threshold. The threshold is configurable and is expressed as a function of the read length. In end-to-end alignment mode, the default minimum score threhsold is -0.6 + -0.6 * L, where L is the read length. In local alignment mdoe, the default minimum score threshold is 20 + 8.0 * ln(L), where L is the read length. This can be configured with the --score-min option. For details on how to set options like --score-min that correpond to functions, see the section on setting function options.


    so if you know the length is 100 bp, you can calculate the min score using the above info.


    however is the 100bp the FRAGMENT LENGTH of the template used for the sequencing machine? or the actual read itself? how do you calculate the length of read? I am not sure if this implies the fragment length used to initially prepare the sequence reads. This is important because after the sequencer outputs the reads, normally we have a distribution of read lengths.

    Comment

    • mastal
      Senior Member
      • Mar 2009
      • 666

      #3
      It means the read length.

      See the section of the manual under the heading Setting Function Options,




      it tells you how to specify the parameters.

      Comment

      • dpryan
        Devon Ryan
        • Jul 2011
        • 3478

        #4
        The length is the length of the read and has nothing to do with the fragment length. If you have a 50bp read and L,-0.6,-0.6, then the minimum score is -31.6.

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Pathogen Surveillance with Advanced Genomic Tools
          by seqadmin




          The COVID-19 pandemic highlighted the need for proactive pathogen surveillance systems. As ongoing threats like avian influenza and newly emerging infections continue to pose risks, researchers are working to improve how quickly and accurately pathogens can be identified and tracked. In a recent SEQanswers webinar, two experts discussed how next-generation sequencing (NGS) and machine learning are shaping efforts to monitor viral variation and trace the origins of infectious...
          03-24-2025, 11:48 AM
        • seqadmin
          New Genomics Tools and Methods Shared at AGBT 2025
          by seqadmin


          This year’s Advances in Genome Biology and Technology (AGBT) General Meeting commemorated the 25th anniversary of the event at its original venue on Marco Island, Florida. While this year’s event didn’t include high-profile musical performances, the industry announcements and cutting-edge research still drew the attention of leading scientists.

          The Headliner
          The biggest announcement was Roche stepping back into the sequencing platform market. In the years since...
          03-03-2025, 01:39 PM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, 03-20-2025, 05:03 AM
        0 responses
        49 views
        0 reactions
        Last Post seqadmin  
        Started by seqadmin, 03-19-2025, 07:27 AM
        0 responses
        57 views
        0 reactions
        Last Post seqadmin  
        Started by seqadmin, 03-18-2025, 12:50 PM
        0 responses
        50 views
        0 reactions
        Last Post seqadmin  
        Started by seqadmin, 03-03-2025, 01:15 PM
        0 responses
        201 views
        0 reactions
        Last Post seqadmin  
        Working...