Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Config file for soapdenovo2

    Hello everyone,

    I am new to assembly and I have a syntax question to ask:

    This is an example contig file as given on https://github.com/aquaskyline/SOAPdenovo2

    In the [LIB] section, we already specify parameters such as avg_ins, reverse_seq, asm_flags, rank etc then why do we again define the limits outside the [LIB] at the end of the file? Are those maximum values of these parameters?

    I tried to find this in the documentation but couldnt get any information

    Code:
    #maximal read length
    max_rd_len=100
    [LIB]
    #average insert size
    avg_ins=200
    #if sequence needs to be reversed
    reverse_seq=0
    #in which part(s) the reads are used
    asm_flags=3
    #use only first 100 bps of each read
    rd_len_cutoff=100
    #in which order the reads are used while scaffolding
    rank=1
    # cutoff of pair number for a reliable connection (at least 3 for short insert size)
    pair_num_cutoff=3
    #minimum aligned length to contigs for a reliable read location (at least 32 for short insert size)
    map_len=32
    #a pair of fastq file, read 1 file should always be followed by read 2 file
    q1=/path/**LIBNAMEA**/fastq1_read_1.fq
    q2=/path/**LIBNAMEA**/fastq1_read_2.fq
    [LIB]
    avg_ins=2000
    reverse_seq=1
    asm_flags=2
    rank=2
    # cutoff of pair number for a reliable connection (at least 5 for large insert size)
    pair_num_cutoff=5
    #minimum aligned length to contigs for a reliable read location (at least 35 for large insert size)
    map_len=35
    q1=/path/**LIBNAMEB**/fastq_read_1.fq
    q2=/path/**LIBNAMEB**/fastq_read_2.fq

  • #2
    It's not outside [LIB], it's a second library that is defined! So for every library you have, you start with a line "[LIB]", and then list all the parameters (specific for each library).
    In the example there is one paired end library and one mate-pair library.

    Linnéa

    Comment


    • #3
      Thanks a lot Linnéa!

      That explains it nicely and now I also understand the 'Library rank'.
      Thank you 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
      72 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 02-05-2025, 10:34 AM
      0 responses
      113 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 02-03-2025, 09:07 AM
      0 responses
      90 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 01-31-2025, 08:31 AM
      0 responses
      49 views
      0 likes
      Last Post seqadmin  
      Working...
      X