Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • A question on Illumina paired-end reads alignment. Merging from different samples

    Hi folks, I have a question on how to better merge paired-end reads coming from different samples. Basically, we do align genome resequencing fastq with bwa mem, and compress the resulting .sam file into a .bam.

    If we have to merge two different samples, there are basically two ways. The first is to merge the fastqs and align the resulting file, the second is to make use of samtools merge to merge the .bam files.

    My concern is whether the two procedures are equally valid, or there is some relevant difference in the outcome.

    I think that everything revolves around the functioning of the Burrows Wheeler Transform Alignment. I have broadly understood the application of the BWT and of the indexing, but I still wonder if the number of reads affects the results of the alignment, or each read is aligned independently.

    Can anyone give me more insights on this?

  • #2
    Originally posted by linudz View Post
    Hi folks, I have a question on how to better merge paired-end reads coming from different samples. Basically, we do align genome resequencing fastq with bwa mem, and compress the resulting .sam file into a .bam.

    If we have to merge two different samples, there are basically two ways. The first is to merge the fastqs and align the resulting file, the second is to make use of samtools merge to merge the .bam files.

    My concern is whether the two procedures are equally valid, or there is some relevant difference in the outcome.
    Either way should be fine after you take the "note" below into account.

    Note: When you are referencing "merging samples" are you referring to technical replicates of the same sample? Merging would be appropriate only in that case. If "samples" are true biological replicates then you would want to keep them separate for downstream analysis.

    but I still wonder if the number of reads affects the results of the alignment, or each read is aligned independently.

    Can anyone give me more insights on this?
    Each read pair is independently aligned to the reference so there is no effect of the amount of data on actual alignments.

    Comment


    • #3
      Thank you very much for your answer. This is going to be done on replicates of the same sample, so it should apply. Thanks again

      Comment


      • #4
        Originally posted by linudz View Post
        Thank you very much for your answer. This is going to be done on replicates of the same sample, so it should apply. Thanks again
        Just want to re-emphasize. Are these technical replicates or biological? Merging data is only appropriate for technical replicates.

        Comment


        • #5
          Because technical replicates could have different insert sizes, and the average (or estimated) insert size can have an effect on alignment, you may get different results even for technical replicates between merging prior to mapping versus merging after mapping. In many cases this difference is trivial and in general mapping the reads independently is not guaranteed to give you an optimal answer with respect to insert size anyway, since there is an order dependency, but it's something to be aware of.

          Comment

          Latest Articles

          Collapse

          • seqadmin
            Techniques and Challenges in Conservation Genomics
            by seqadmin



            The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

            Avian Conservation
            Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
            03-08-2024, 10:41 AM
          • 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...
            02-26-2024, 02:07 PM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, 03-14-2024, 06:13 AM
          0 responses
          32 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 03-08-2024, 08:03 AM
          0 responses
          71 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 03-07-2024, 08:13 AM
          0 responses
          80 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 03-06-2024, 09:51 AM
          0 responses
          68 views
          0 likes
          Last Post seqadmin  
          Working...
          X