Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Partially mapped reads

    Hi there,
    I am mapping illumina reads on a reference by using the bwa tool. I was wondering whether there is an option or an alternative aligner that can allow me to highlight partially mapped reads at the begin and end of the reference sequence. Say for example I have 150 bp reads and the first 50 bp of my reference genome align with the last 50 bases of one of my read. In this case it is reasonable to think the the remaining 100 base on the read lie on the 5' position of the reference genome (I guess!)
    your help is really appreciated!

  • #2
    Take some of those reads, blast them and see what comes up. If the first 100 bp matches a known target, then you know your answer. Also check for contamination by adapter sequences.

    Comment


    • #3
      Thanks for the advice. Actually I have 31 gbites reads (a resequencing project). for this reason I would like to avoid using blast because of speed issues. Any other ideas with BW aligners?
      cheers

      Comment


      • #4
        I'm not saying BLAST everything. Take a couple of reads that have this weird behavior and blast them. The goal is simply to see whether the unmapped portion matches anything else in the genome or if its something completely different.

        Comment


        • #5
          I probably did not underline my problem properly, Sorry about that. I have a piece of genome (my reference sequence) and I would like to undercover the sequence on the 5' and 3' portion of that sequence. If I manage to blast a 50 bp read on the 5' position than the remaining 100bp of the read will be flanking my reference genome. Let's say I would like to do an "in silico" genome walking. In order to do so I have to bast all the reads which will take a lot of time. Actually considering only a small portion of the reference genome may help me to solve the time issue probably.
          thanks a lot for your help anyway!

          Comment


          • #6
            The new bra-mem algorithm for longer reads should be able to report local alignments (i.e, not requiring the whole read to match the reference). If not
            you can try bowtie2 with the --local option and look for reads with partial alignments in the beginning / end of your reference.

            Comment


            • #7
              That's great advice! I will try straight away!
              thanks a lot!

              Comment

              Latest Articles

              Collapse

              • seqadmin
                Current Approaches to Protein Sequencing
                by seqadmin


                Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
                04-04-2024, 04:25 PM
              • seqadmin
                Strategies for Sequencing Challenging Samples
                by seqadmin


                Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                03-22-2024, 06:39 AM

              ad_right_rmr

              Collapse

              News

              Collapse

              Topics Statistics Last Post
              Started by seqadmin, 04-11-2024, 12:08 PM
              0 responses
              13 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 04-10-2024, 10:19 PM
              0 responses
              17 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 04-10-2024, 09:21 AM
              0 responses
              14 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 04-04-2024, 09:00 AM
              0 responses
              43 views
              0 likes
              Last Post seqadmin  
              Working...
              X