Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Joining Contigs

    Hi all,

    after comparing a newbler assembly with a reference genome I have identified two contigs which contain the beginning and end (respectively) of the same gene. A PCR was run to bridge the gap and the product was sequenced yielding a 790bp fragment which overlaps both contigs
    I am curious to know which tool is best to close the gap between the contigs?

    regards

    Brian

  • #2
    Did you ask newbler to produce .ace files? -- I'd look at the ends of the contigs to see if there's any noticeable increase in read disagreement, or decrease in depth. It's possible that you simply had too few reads, and by chance did not span the gap ... or that other parts of your genome were similar enough to break the contig at that point. The .ace file should tell you the former, but you may need to run some blasts (or maybe pore through the alignment result files??) for the latter.

    Comment


    • #3
      If you had told newbler to produce consed output, I'd use consed [1] for that
      purpose. At least, that's what we are doing :-)
      Most other software (except Staden's gap4 [2]) has very limited editing capabilities AFAIK.

      cheers,
      Sven

      [1] = http://bozeman.mbt.washington.edu/consed/consed.html
      [2] = http://sourceforge.net/projects/staden

      Comment


      • #4
        Gap4 works ok, but you need external tools (eg caftools from Sanger) to convert ace to gap4 as gap4's input format handling is shockingly bad still.

        Gap5 is getting there. It's not really production ready, but a version exists (also at source forge) and I plan on making a new release again soon. The CVS tree (and hence next release) has full support for reading ACE files.

        Comment


        • #5
          thanks everyone,

          I did it in the end using Phrap. Although I had to make a quality file, on a temp basis, for the PCR product (quality of 20).

          thanks again

          Comment

          Latest Articles

          Collapse

          • seqadmin
            Quality Control Essentials for Next-Generation Sequencing Workflows
            by seqadmin




            Like all molecular biology applications, next-generation sequencing (NGS) workflows require diligent quality control (QC) measures to ensure accurate and reproducible results. Proper QC begins at nucleic acid extraction and continues all the way through to data analysis. This article outlines the key QC steps in an NGS workflow, along with the commonly used tools and techniques.

            Nucleic Acid Quality Control
            Preparing for NGS starts with isolating the...
            02-10-2025, 01:58 PM
          • seqadmin
            An Introduction to the Technologies Transforming Precision Medicine
            by seqadmin


            In recent years, precision medicine has become a major focus for researchers and healthcare professionals. This approach offers personalized treatment and wellness plans by utilizing insights from each person's unique biology and lifestyle to deliver more effective care. Its advancement relies on innovative technologies that enable a deeper understanding of individual variability. In a joint documentary with our colleagues at Biocompare, we examined the foundational principles of precision...
            01-27-2025, 07:46 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, 02-07-2025, 09:30 AM
          0 responses
          65 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 02-05-2025, 10:34 AM
          0 responses
          101 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 02-03-2025, 09:07 AM
          0 responses
          81 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 01-31-2025, 08:31 AM
          0 responses
          45 views
          0 likes
          Last Post seqadmin  
          Working...
          X