Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Realigning unmapped reads to a second reference

    I am looking for suggestions for the following issue/task:

    I am using a reference genome for which one of my scaffolds is an incomplete (e.g. is a much smaller subset of one of the genomic regions that I am after based on the complete genome of a closely related species). All other scaffolds are relatively complete for my focal species.

    What I would like to do is first align my Illumina paired-end reads to my focal species' reference (all scaffolds) and then take the reads that align to my short scaffold PLUS the discarded reads and realign them to the more complete version of this particular region from the closely-related model organism in order to pick up any additional reads that might cover that might align to that region.

    I think the approach outlined below will work, but I'm looking for feedback and answers to the questions in my comments before I proceed with my big files!

    # After initial trimming, alignment with bwa and removal of duplicates with Piccard, get reads aligning to target scaffold based on bed file
    # ??? does removing duplicates with Piccard keep unpaired reads ???

    samtools view -L scaffoldName myAlign.bam -b > myScaffold.bam

    # get discarded reads (reads that are paired but unmapped or with a mate that is unmapped)
    # ??? Are these the right flags for this??? I'm most uncertain about the 0x0001 flag ???

    samtools view -F -0x0001 -f 0x0004 -f 0x0008 -b myAlign.bam -b > discarded.bam

    # merge target aligned and discarded

    samtools merge targetDiscarded.bam

    # use Picard's SamToFastq to convert merged bam to fastq
    # ??? based on the above the unpaired fastq should be empty right ???

    java -Xmx2g SamToFastq.jar INPUT=targetDiscarded.bam FASTQ=r1.fastq SECOND_END_FASTQ=r2.fastq UNPAIRED_FASTQ=unpaired.fastq

    # align to second species ref with bwa mem

    bwa mem -M modelSpeciesRef r1.fastq r2.fastq unpaired.fastq > finalAlignment.sam

    ??? I know there is also a way to do this using the aln and sampe function in samtools...how does this compare?

    Thanks in advance for any advice/help

  • #2
    Why not just concatenate your reference with the model organism's scaffold and map to the combined fasta? I guess I'm not really sure what your goal is.

    Comment


    • #3
      Because we do have part of the region from our target species, I was hoping to make use of this information...In theory if the target and model are close enough it shouldn't matter but I do think it would be good to try it both ways...

      Comment


      • #4
        You can use BBSplit to map to two references simultaneously. In your case, you could set "ref=organism.fa,scaffold.fa" and then the reads would just go to whichever they match best. If you use the flag "ambig2=all" then any read that maps equally well to both references would be assigned to both of them, while if you set "ambig2=best" then it would preferentially map to the first reference, which is probably what you want.

        Comment


        • #5
          Unfortunately I'm constrained to use BWA as my project is part of a larger project with the alignment step of the pipeline already in place...

          I'm really asking if there are any flaws in my thinking as to what the script I have indicated is doing...(e.g. the different flags and file conversions)...

          Thanks

          Comment

          Latest Articles

          Collapse

          • seqadmin
            Best Practices for Single-Cell Sequencing Analysis
            by seqadmin



            While isolating and preparing single cells for sequencing was historically the bottleneck, recent technological advancements have shifted the challenge to data analysis. This highlights the rapidly evolving nature of single-cell sequencing. The inherent complexity of single-cell analysis has intensified with the surge in data volume and the incorporation of diverse and more complex datasets. This article explores the challenges in analysis, examines common pitfalls, offers...
            06-06-2024, 07:15 AM
          • seqadmin
            Latest Developments in Precision Medicine
            by seqadmin



            Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

            Somatic Genomics
            “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
            05-24-2024, 01:16 PM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, Today, 07:24 AM
          0 responses
          9 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, Yesterday, 08:58 AM
          0 responses
          11 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 06-12-2024, 02:20 PM
          0 responses
          16 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 06-07-2024, 06:58 AM
          0 responses
          184 views
          0 likes
          Last Post seqadmin  
          Working...
          X