Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • gsAssembly (Newbler) de novo behaviour, inputs and outputs

    Hello there,
    It appears that the content of the ACE file is quite inconsistent with the fasta files content and the 454NewblerMetrics file:
    Sometimes, the fna files are empty, 454NewblerMetrics has no contigs nor isotigs but...
    The corresponding ACE will have a hundred contigs/isotigs with a mean depth of 10 000 and minimum size of 120.
    Looking into the ace file show that these sequences are fulled with "*".

    So I wonder: is the ace file a temporary (unfiltered) one ?

    Another point: it seems that running a de novo assembly starting from fasta and not sff improve the assembly quality (less "*" in the sequences)

    Anyone has seen that ? Is there some Newbler options to use for fixing this ?
    Cheers

  • #2
    Based on the fact that you state that the 454NewblerMetrics file reports no contigs, I would assume there is something that prevents assembly (or is this mapping?). Any contigs in the ace file are probably artifacts.

    What kind of project/sample/reads do you have?

    The '*' symbols represent gaps introduced to optimize the alignment, and 454 ace files have tons of them due to the homopolymer errors (or rather, variation in homopolymer length between reads). I would not deem the amount of these '*' a measure of quality. If you want to know the effect of using fasta input over sff, I would take reads from a known genome and check the correctness of the contigs relative to the reference sequence. Although I expect sff input to give better results, you never know...

    If somebody has done this already, let me know :-)

    Comment


    • #3
      Hi flxlex,
      I'm considering de novo assemblies only, these observations are based on the study of dozens of assemblies done with Newbler 2.3 and this problem occurs:
      - on gDNA from bacteria
      - on gDNA from eukarya
      - on cDNA (option -cdna) from plantae
      But I all cases have some example cons.

      Using the same raw datas we obtained different behaviour from newbler.
      If we provide a sff file, contigs will be mainly composed by '*' what is not true if we basicly convert this file into fasta and assmble it with newbler. If doing so the asmmbly will look just fine.
      Where could the behaviour difference could come from ?

      Still investigate...
      Last edited by nicolallias; 10-15-2010, 12:02 AM.

      Comment


      • #4
        Originally posted by nicolallias View Post
        Hi flxlex,
        I'm considering de novo assemblies only, these observations are based on the study of dozens of assemblies done with Newbler 2.3 and this problem occurs:
        - on gDNA from bacteria
        - on gDNA from eukarya
        - on cDNA (option -cdna) from plantae
        But I all cases have some example cons.
        I don't understand, you get empty fna files and ace files with many '*' contigs for all these assemblies?

        Using the same raw datas we obtained different behaviour from newbler.
        If we provide a sff file, contigs will be mainly composed by '*' what is not true if we basicly convert this file into fasta and assmble it with newbler. If doing so the asmmbly will look just fine.
        Where could the behaviour difference could come from ?
        I am not sure, but what I would be more interested in is how correct the contigs are between assemblies using sff and fasta files. if you do this for a known genome (e.g. E. coli) what do you find?

        Comment


        • #5
          I don't understand, you get empty fna files and ace files with many '*' contigs for all these assemblies?
          Yes
          I am not sure, but what I would be more interested in is how correct the contigs are between assemblies using sff and fasta files. if you do this for a known genome (e.g. E. coli) what do you find?
          I will do that.

          Edit: Both alignments on sff and fasta have been proceeded on E.coli, both gave ace and fna files full of contigs. But I won't look after the best method (fasta or sff) for alignment : the main question here is why ?
          Sometimes we have nothing (ace fulled of "*", fna files empty) with the sff while we have something (ace and fna files have contigs) when starting with a fasta ?

          We're getting contact with Roche...
          Last edited by nicolallias; 10-19-2010, 12:42 AM.

          Comment


          • #6
            The problem seems to be the sff file, software developers at 454 are working on this.
            Hope some good news soon ;-)

            Comment


            • #7
              Some news : the matter seems to be located on the sff file generation...

              Comment

              Latest Articles

              Collapse

              • seqadmin
                Exploring the Dynamics of the Tumor Microenvironment
                by seqadmin




                The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
                07-08-2024, 03:19 PM
              • seqadmin
                Exploring Human Diversity Through Large-Scale Omics
                by seqadmin


                In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
                06-25-2024, 06:43 AM

              ad_right_rmr

              Collapse

              News

              Collapse

              Topics Statistics Last Post
              Started by seqadmin, Today, 11:09 AM
              0 responses
              14 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 07-19-2024, 07:20 AM
              0 responses
              146 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 07-16-2024, 05:49 AM
              0 responses
              120 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 07-15-2024, 06:53 AM
              0 responses
              111 views
              0 likes
              Last Post seqadmin  
              Working...
              X