Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • SOAPdenovo2

    Hello everyone, I am new to this forum and new to using assembly tools as well. I am trying to use SOAPdenovo2 and tried to write the configuration file which is as follows :

    #maximal read length
    max_rd_len=36
    [LIB]
    #average insert size
    avg_ins=100
    #if sequence needs to be reversed
    reverse_seq=0
    #in which part(s) the reads are used
    asm_flags=3
    f=sample.fq

    I am using soapdenovo2 in steps. I tried to use the pregraph_sparse63mer version using this command :

    /pregraph_sparse63mer -s configfile -K 63 -z 1000000 -o output_graph -p 4


    The programs starts to build the graph but no nodes are allocated and it seems to parse the file but O kmer detected etc...

    I have tried with combinations of -k 23 and smaller values because my reads are very short (36 nt) in length and it is a bacterial file 88 MB

    I do not know if my parameters are correct in the config file but I am using only ONE FASTQ file.

    Thank you.

  • #2
    SOAPdenovo2

    If the file with your reads is in fastq format you should have
    q=sample.fq

    in the configuration file.

    Did you get any better results with -K 23 or smaller?
    Kmer size needs to be shorter than your reads.

    Comment


    • #3
      Hi, Thank you. Yes, it started to construct the graph and when I used the command soapdenovo2 63-mer contig -g output_graph , it produced a file named output_graph.contig

      However, the contigs are very small and in this format:

      >1 length 16 cvg_0.0_tip_0
      AAAAAAAAAAAAAAAA

      is that alright?

      I used the same pregraph_sparse command but used -K 15

      Comment


      • #4
        SOAPdenovo2

        Hi,

        A month ago, i was struggling with pregraph_sparse too.
        My problem was, as you stated "The programs starts to build the graph but no nodes are allocated and it seems to parse the file but O kmer detected etc..."

        I compiled from source in the first place. Then, i decided to give binaries a try, and it was finally working. Maybe this could solve your issue.

        Best,
        Erwan

        Comment


        • #5
          Hi thank you, I compiled from the source and the program is working but the length of the contigs are in the format that I stated above and does not seem right. Have you used Soapdenovo2 and are your contigs in the same format?

          Thank you.

          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...
            04-22-2024, 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, Yesterday, 11:49 AM
          0 responses
          15 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 04-24-2024, 08:47 AM
          0 responses
          16 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 04-11-2024, 12:08 PM
          0 responses
          61 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 04-10-2024, 10:19 PM
          0 responses
          60 views
          0 likes
          Last Post seqadmin  
          Working...
          X