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
          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
        • seqadmin
          Non-Coding RNA Research and Technologies
          by seqadmin




          Non-coding RNAs (ncRNAs) do not code for proteins but play important roles in numerous cellular processes including gene silencing, developmental pathways, and more. There are numerous types including microRNA (miRNA), long ncRNA (lncRNA), circular RNA (circRNA), and more. In this article, we discuss innovative ncRNA research and explore recent technological advancements that improve the study of ncRNAs.

          Nobel Prize for MicroRNA Discovery
          This week,...
          10-07-2024, 08:07 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, 11-01-2024, 06:09 AM
        0 responses
        11 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 10-30-2024, 05:31 AM
        0 responses
        14 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 10-24-2024, 06:58 AM
        0 responses
        24 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 10-23-2024, 08:43 AM
        0 responses
        52 views
        0 likes
        Last Post seqadmin  
        Working...
        X