Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • the N50 is so low from soapdenovo

    Hi, i used soapdenovo for my rice genome assembly,i tried many times,but the N50 is only 920(not an encouraged result it seems ) i also tried velvet, but it keep running for two week.
    raw data: PE of reads length 90bp,20*coverage, insert size 500bp(only one library).
    /soapdenovo63kmer all -s rice.config -o rice -k (31..59) -R -d -F -L 100

    is there something wrong with my parameter set or is there any other tool for large genome assembly?

    all reply is appreciated!

    heiya

  • #2
    With only 20 x coverage of 90 bp PE reads, this is about the best assembly you can expect. You would need at the very least to double the coverage, and preferably add Mate Pair reads as well (to get scaffolds)

    Comment


    • #3
      Originally posted by flxlex View Post
      With only 20 x coverage of 90 bp PE reads, this is about the best assembly you can expect. You would need at the very least to double the coverage, and preferably add Mate Pair reads as well (to get scaffolds)
      thanks, flxlex!

      Comment


      • #4
        Originally posted by heiya View Post
        Hi, i used soapdenovo for my rice genome assembly,i tried many times,but the N50 is only 920(not an encouraged result it seems ) i also tried velvet, but it keep running for two week.
        raw data: PE of reads length 90bp,20*coverage, insert size 500bp(only one library).
        /soapdenovo63kmer all -s rice.config -o rice -k (31..59) -R -d -F -L 100

        is there something wrong with my parameter set or is there any other tool for large genome assembly?

        all reply is appreciated!

        heiya
        Hi, heiya,
        Have you find the way to improve your assembly's N50 ?
        I'm dealing with the same illumina 101 PE data looks the same to yours.
        Before assembly, quality filtering, trimming and errors correction would be essential.
        I'd preprocessing and got assembly contig N50= 2.5Kb and scaffold N50=15Kb with Soapdenovo2 (Ave_ins 232, -K50).
        But I'm wondering if the assembly can be done better.
        Additionally, by mapping reads to contig sequence, the estimated average insert size was 232 +/-63 bp, looks different from the expected size (~ 360 bp)
        Thankx
        Last edited by lovenlong; 05-31-2013, 05:17 PM.

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Genetic Variation in Immunogenetics and Antibody Diversity
          by seqadmin



          The field of immunogenetics explores how genetic variations influence immune responses and susceptibility to disease. In a recent SEQanswers webinar, Oscar Rodriguez, Ph.D., Postdoctoral Researcher at the University of Louisville, and Ruben Martínez Barricarte, Ph.D., Assistant Professor of Medicine at Vanderbilt University, shared recent advancements in immunogenetics. This article discusses their research on genetic variation in antibody loci, antibody production processes,...
          Yesterday, 07:24 PM
        • 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

        ad_right_rmr

        Collapse

        News

        Collapse

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