Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Can BWA align reads to references that are even shorter?

    We have obtained a lot of RNA-seq reads from HiSeq with 100bp in length. We are trying to align them to a splice junction database which only contain sequences that are 62bp in length. Is this possible? Are there any parameter values that need to be tweaked? Would soft clipping automatically mark both hanging 5' and 3' read fragments?

    Hope some BWA guru can answer the above questions. Thanks!

  • #2
    I'm trying to pop this post up for myself. Nobody had any experience with the above challenge?

    Comment


    • #3
      Sorry, never tried it. I used SOAP for bacterial alignments, and I know it seg faults when reference sequence is smaller than reads.

      I recommend just trying it. If it doesn't work, maybe pad the 62 base junctions with N's.

      Comment


      • #4
        I know I had some reads that were over-hanging from the end of the contigs (when mapping RNA seq onto CDS sequences) and it worked fine with BWA.

        Comment


        • #5
          On second thought, I think the way bwa works is it concatenates all the reference sequences into one long file. So if a read falls off the edge of one, and lands on another, that won't stop mapping.

          However, what this means is that if your read starts a few bases upstream of reference junction B, and Junction A was concatenated to junction B such that the start of B is right after the end of A, samtools is going to report that your read is in junction A, becuase that's where it started. I think there will be a flag, or the MQ will be zero, something like that, so it would be catchable, but no easily, and it won't be easy to figure out. It's probably easier to pad your reference with n's so that this won't happen.

          Comment

          Latest Articles

          Collapse

          • seqadmin
            Genetic Variation in Immunogenetics and Antibody Diversity
            by seqadmin



            The field of immunogenetics explores how genetic variations influence immune responses and susceptibility to disease. In a recent SEQanswers webinar, Oscar Rodriguez, Ph.D., Postdoctoral Researcher at the University of Louisville, and Ruben Martínez Barricarte, Ph.D., Assistant Professor of Medicine at Vanderbilt University, shared recent advancements in immunogenetics. This article discusses their research on genetic variation in antibody loci, antibody production processes,...
            11-06-2024, 07:24 PM
          • seqadmin
            Choosing Between NGS and qPCR
            by seqadmin



            Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...
            10-18-2024, 07:11 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, 11-01-2024, 06:09 AM
          0 responses
          29 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 10-30-2024, 05:31 AM
          0 responses
          21 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 10-24-2024, 06:58 AM
          0 responses
          26 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 10-23-2024, 08:43 AM
          0 responses
          57 views
          0 likes
          Last Post seqadmin  
          Working...
          X