Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Using Pilon after MaSuRCA Assembly

    So, lots of good threads on using Pilon...but my question is specifically about the use of Pilon after using MaSuRCA for assembly. MaSuRCA specifically states that the unmodified, raw reads be used without any sort of manipulation:

    "IMPORTANT! Do not use third party tools to pre-process the Illumina data before providing it to MaSuRCA, unless you are absolutely sure you know exactly what the preprocessing tool does. Do not do any trimming, cleaning or error correction. This will likely deteriorate the assembly."

    Again, no problem for the assembly portion. BUT....when I move to the Pilon polishing step(s), do I want to try and use the raw reads as they are for polishing?!? Seems that, depending on how Pilon interprets that data, it might introduce some of the low quality aspects back into the assembly? Or does the mapping of the raw to the assembly to create the .bam resolve any such possibilities?

  • #2
    Personally, I would do quality and adapter trimming (and maybe error-correction) of the raw Illumina reads as MaSuRCA pretty much does this from what I understand and use the quality and adapter trimmed reads for polishing with Pilon.

    Comment


    • #3
      Kind of leaning that way for sure. I dug around a bit in the MaSuRCA documentation and have identified and located the "corrected" PE data created during assembly:

      pe.cor.fa - error corrected PE reads. The ordering of the reads is arbitrary, but the pairs are guaranteed to appear together. No quality scores

      So, the file looks more like merged PE data and is now a .fasta, though I don't know if Pilon was using the quality scores or not...

      I figure it's worth a try anyway....maybe compare it with the raw PE data and see what looks better in a Quast analysis...
      Last edited by jpummil; 10-02-2019, 06:27 AM.

      Comment


      • #4
        Haha! Well, THAT didn't work! BWA did OK with the indexing, etc...but Pilon quickly bailed:

        Suppressed: java.lang.IllegalStateException: Inappropriate call if not paired read

        I'll go back and manually quality trim and error correct the raw files.

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Genetic Variation in Immunogenetics and Antibody Diversity
          by seqadmin



          The field of immunogenetics explores how genetic variations influence immune responses and susceptibility to disease. In a recent SEQanswers webinar, Oscar Rodriguez, Ph.D., Postdoctoral Researcher at the University of Louisville, and Ruben Martínez Barricarte, Ph.D., Assistant Professor of Medicine at Vanderbilt University, shared recent advancements in immunogenetics. This article discusses their research on genetic variation in antibody loci, antibody production processes,...
          11-06-2024, 07:24 PM
        • seqadmin
          Choosing Between NGS and qPCR
          by seqadmin



          Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...
          10-18-2024, 07:11 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, 11-08-2024, 11:09 AM
        0 responses
        211 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 11-08-2024, 06:13 AM
        0 responses
        156 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 11-01-2024, 06:09 AM
        0 responses
        80 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 10-30-2024, 05:31 AM
        0 responses
        27 views
        0 likes
        Last Post seqadmin  
        Working...
        X