Announcement

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

  • ehlin
    replied
    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

    Leave a comment:


  • swbarnes2
    replied
    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.

    Leave a comment:


  • ehlin
    started a topic bwa aln for long gaps

    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?

Latest Articles

Collapse

  • seqadmin
    Multiomics Techniques Advancing Disease Research
    by seqadmin


    New and advanced multiomics tools and technologies have opened new avenues of research and markedly enhanced various disciplines such as disease research and precision medicine1. The practice of merging diverse data from various ‘omes increasingly provides a more holistic understanding of biological systems. As Maddison Masaeli, Co-Founder and CEO at Deepcell, aptly noted, “You can't explain biology in its complex form with one modality.”

    A major leap in the field has
    ...
    02-08-2024, 06:33 AM
  • seqadmin
    The 3D Genome: New Technologies and Emerging Insights
    by seqadmin


    The study of three-dimensional (3D) genomics explores the spatial structure of genomes and their role in processes like gene expression and DNA replication. By employing innovative technologies, researchers can study these arrangements to discover their role in various biological processes. Scientists continue to find new ways in which the organization of DNA is involved in processes like development1 and disease2.

    Basic Organization and Structure
    Understanding...
    01-22-2024, 03:25 PM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, Today, 08:57 AM
0 responses
9 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-14-2024, 09:19 AM
0 responses
42 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-12-2024, 03:37 PM
0 responses
402 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-09-2024, 03:36 PM
0 responses
646 views
0 likes
Last Post seqadmin  
Working...
X