Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • bwa aln for long gaps

    Hello everyone, and thank you in advance for your help.

    I am trying to align Illumina Ion Torrent reads to a reference genome, and these reads are in the TCR region. This means that there may be long gaps in between reads due to VDJ recombination (such as 1 read that is supposed to be from 7:1-100 and then 7:800-900, etc.) I wanted bwa aln to account for this by adjusting gap extension penalty to 0, and allowing 600,000 gap extensions to cover the whole region. It doesn't seem to be working, does anyone have any insights?

  • #2
    I don't think that bwa is going to map reads like that.

    It would have to be splitting the read up and doing two alignments, and it doesn't do that. The longest indel I've seen it was was about 10 bases.

    bwa would tell you that you would have no coverage in the gap, and it might align some reads only by soft clipping the overhang, so you could figure it out from that, given enough time, but bwa is never going to give you an alignment with a 800 bp gap.

    Comment


    • #3
      Originally posted by swbarnes2 View Post
      I don't think that bwa is going to map reads like that.

      It would have to be splitting the read up and doing two alignments, and it doesn't do that. The longest indel I've seen it was was about 10 bases.

      bwa would tell you that you would have no coverage in the gap, and it might align some reads only by soft clipping the overhang, so you could figure it out from that, given enough time, but bwa is never going to give you an alignment with a 800 bp gap.
      I see, that is unfortunate... but better that I know, since I tried a few more times with SAM files with 0 bytes, I'll keep that in mind. Thank you so much for your help!

      -Edwin

      Comment

      Latest Articles

      Collapse

      • seqadmin
        Exploring the Dynamics of the Tumor Microenvironment
        by seqadmin




        The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
        07-08-2024, 03:19 PM
      • seqadmin
        Exploring Human Diversity Through Large-Scale Omics
        by seqadmin


        In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
        06-25-2024, 06:43 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by seqadmin, Yesterday, 07:20 AM
      0 responses
      23 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 07-16-2024, 05:49 AM
      0 responses
      38 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 07-15-2024, 06:53 AM
      0 responses
      44 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 07-10-2024, 07:30 AM
      0 responses
      41 views
      0 likes
      Last Post seqadmin  
      Working...
      X