Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Summary Statistics in Newbler?

    The Mapper output comes with handy xls file with a lot of great statistics in it. The Newbler output doesn't seem to have anything like this? I am particularly looking for statistics about the coverage in the contigs (avg, min, max, etc.)

    Anyone know where this info is or a good program to use to find it?

  • #2
    Hi,

    you might want to have a look at 454NewblerMetrics.txt (or csv) for assembly statistics and 454ReadStatus.txt to see what happened to your reads ...

    A detailed contig description (for parsing) gives you 454AlignmentInfo.tsv .

    hth,
    Sven

    Comment


    • #3
      And/or have newbler produce the 454ContigGraph.txt file which gives you average coverage per contig (not min/max). This can be done by setting the -g flag (note: this sets the minimum contig length to report from 100 bp to 1 bp...).

      Comment


      • #4
        Thanks for your reply, Flxlex. Couple questions: How do you set the g-flag? Is there no way to get min/max coverage?

        Thanks.

        Comment


        • #5
          Thanks, Sklages. The metrics files don't seem to have info about coverage, right?

          The 454AlignmentInfo.tsz files only gives me partial info on one contig. Is that right? It can't be. Any ideas?

          Thanks so much.

          Comment


          • #6
            The 454NewblerMetrics.txt reports only common stats about the assembly.

            Well, and the tsv file covers all contigs:

            Code:
            $ fgrep -c \>contig 454AlignmentInfo.tsv
            430
            $ fgrep -c \> 454AllContigs.fna
            430
            cheers,
            Sven

            Comment


            • #7
              gaster,

              Are you using the Roche GUI (gsAssembler) to view the 454AlignmentInfo.tsv? If so this viewer truncates the display at 50,000 lines. The file is located in the assembly subdirectory of your main project directory. You can view the entire file with a standard text editor. (Careful, if your assembly is large this file will get huge.) You will need to process the file using some sort of script to extract and calculate the coverage information you asked about (ave, min, max, etc.).

              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, 07-19-2024, 07:20 AM
              0 responses
              28 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 07-16-2024, 05:49 AM
              0 responses
              41 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 07-15-2024, 06:53 AM
              0 responses
              46 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 07-10-2024, 07:30 AM
              0 responses
              43 views
              0 likes
              Last Post seqadmin  
              Working...
              X