Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Help with interpreting IDBA_UD assembly output

    Hello,

    I was trying the IDBA_UD assembler for my metagenomes for the first time. I have been trying to find a documentation on interpreting the outputs but haven't found any. Could anyone help me with the output? I also went through old postings here but was not able to find the answer for my question. Any advice would be greatly appreciated!

    An output directory containing multiple fasta files was generated after running the IDBA_UD. I assume that I have to use "contig.fa" for the further analysis. Am I correct? Could you explain differences between them?

    I have hard time to interpret the "log" file. Below is the one I got. Which one is the stats for the final "contig.fa"? Why are there two "kmer of 60"?

    number of threads 1
    reads 46326304
    long reads 0
    extra reads 0
    read_length 100
    kmer 60
    kmers 77958510 76121085
    merge bubble 1215
    contigs: 75437 n50: 119 max: 5289 mean: 118 total length: 8910129 n80: 63
    aligned 1835482 reads
    confirmed bases: 2768028 correct reads: 410078 bases: 343009
    kmer 20
    kmers 209804973 209944022
    merge bubble 142969
    contigs: 1324264 n50: 72 max: 6189 mean: 51 total length: 68617304 n80: 28
    aligned 936011 reads
    confirmed bases: 5037462 correct reads: 159618 bases: 20497
    distance mean 271.123 sd 122.238
    seed contigs 51438 local contigs 2648528
    kmer 40
    kmers 87229157 87354041
    merge bubble 9179
    contigs: 283618 n50: 261 max: 8023 mean: 127 total length: 36149932 n80: 71
    aligned 1622279 reads
    confirmed bases: 8349992 correct reads: 299441 bases: 27070
    distance mean 280.19 sd 148.882
    seed contigs 53979 local contigs 567236
    kmer 60
    kmers 52067154 52217769
    merge bubble 2684
    contigs: 136775 n50: 361 max: 12385 mean: 212 total length: 29072718 n80: 118
    aligned 3012828 reads
    confirmed bases: 9798172 correct reads: 417627 bases: 31962
    distance mean 281.18 sd 147.978
    seed contigs 47029 local contigs 273550
    kmer 80
    kmers 33745199 33462753
    merge bubble 1343
    contigs: 69296 n50: 478 max: 8256 mean: 343 total length: 23776057 n80: 301
    aligned 9589068 reads
    confirmed bases: 10063129 correct reads: 620205 bases: 192024
    distance mean 261.349 sd 149.598
    seed contigs 39512 local contigs 138592
    kmer 100
    kmers 20137243 19250177
    merge bubble 959
    contigs: 34622 n50: 594 max: 10832 mean: 569 total length: 19702800 n80: 377
    reads 46326304
    aligned 5056106 reads
    distance mean 271.986 sd 147.421
    expected coverage 0.310699
    edgs 137
    contigs: 34555 n50: 593 max: 16184 mean: 569 total length: 19690114 n80: 377

  • #2
    I have se same question!

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Essential Discoveries and Tools in Epitranscriptomics
      by seqadmin


      The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist on Modified Bases...
      Yesterday, 07:01 AM
    • seqadmin
      Current Approaches to Protein Sequencing
      by seqadmin


      Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
      04-04-2024, 04:25 PM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 04-11-2024, 12:08 PM
    0 responses
    37 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-10-2024, 10:19 PM
    0 responses
    41 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-10-2024, 09:21 AM
    0 responses
    35 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-04-2024, 09:00 AM
    0 responses
    55 views
    0 likes
    Last Post seqadmin  
    Working...
    X