Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Alignment with BWA

    Hello!

    Im using this BWA to align some reads

    Code:
    bwa index -p my_reference my_reference.fasta
    bwa mem -t 4 my_reference input.fastq > output.sam 2
    After what i found that in the position im interesteed in there are some bugs, or something like this. As shown on the picture i've attached in some reads motif TACCCGAGG in mutant reads algned with 3M1D6M with 1 mismatch, but in others as 9M with 1 mismatch. How can i make BWA to align all such reads as 9M?
    Attached Files

  • #2
    Perhaps I'm confused, but it looks like there are 3 classes of reads:

    TACTCAAAGG (normal)
    TACCCGAAGG (mutant 1)
    TACCCGAGG (mutant 2)

    The second class of mutant has a deletion in it, so of course it has to be aligned with a deletion somewhere.

    Comment


    • #3
      Originally posted by Brian Bushnell View Post
      Perhaps I'm confused, but it looks like there are 3 classes of reads:

      TACTCAAAGG (normal)
      TACCCGAAGG (mutant 1)
      TACCCGAGG (mutant 2)

      The second class of mutant has a deletion in it, so of course it has to be aligned with a deletion somewhere.
      So is there any way to align (mutant 2) as 6M1D3M with 2 mismatches rather than 3M1D6M with 2 mismatches, or I have to edit Sam file manually after each alignment? Because real mutations are T4C and A6G, but frequant sequance mistake 6delA in mutant reads disturb alignment and pileup consequantly and makes score of real mutations lower in SNP calling.

      Comment


      • #4
        You will have to manually edit the cigar strings. It's possible that if you align to the reverse-complement of the reference, the deletion will be put in the alternate location, though.

        Comment

        Latest Articles

        Collapse

        • seqadmin
          The Impact of AI in Genomic Medicine
          by seqadmin



          Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
          Yesterday, 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
        43 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 02-21-2024, 08:52 AM
        0 responses
        59 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 02-20-2024, 08:57 AM
        0 responses
        47 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 02-14-2024, 09:19 AM
        0 responses
        65 views
        0 likes
        Last Post seqadmin  
        Working...
        X