Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • BWA Combined Reference vs sequential alignment

    Hi,

    I am currently developing a pipeline for novel viral detection and I am using BWA for filtering the mapped sequnces.

    I recognized during my tests that I have a different number of unmapped reads after I align my sample against the whole genome once, in contrast to if I align each chromosome seperately and pipe the output.

    Is this causing BWA heuristically concept, that I get a different number of mapped sequences by using the same references? Is this normal?

    What is better, to align chromosome by chromosome or to align the whole genome at once?

    Would be pelased if anyone could help.

  • #2
    If you align against each chromosome individually, you will have some reads that ought to align to a different chromosome being forced to align to the chromosome you gave it. So the same read will align to multiple chromosomes, and would be counted multiple times. Does that explain what you are seeing?

    It's definately better to give the aligner the whole reference sequence, so the aligner can find the best match, the real match for each read, rather than the best match of the sequences you gave it. Sure, it takes more memory, and more time, but that's what it takes.

    Comment

    Latest Articles

    Collapse

    • seqadmin
      The Impact of AI in Genomic Medicine
      by seqadmin



      Article Coming Soon......
      Today, 02:07 PM
    • 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

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 02-23-2024, 04:11 PM
    0 responses
    31 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 02-21-2024, 08:52 AM
    0 responses
    46 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 02-20-2024, 08:57 AM
    0 responses
    36 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 02-14-2024, 09:19 AM
    0 responses
    63 views
    0 likes
    Last Post seqadmin  
    Working...
    X