Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • bioramg
    Member
    • Apr 2014
    • 12

    SPAdes assembler

    Hello,
    I tried to assemble a mitochondrial genome assembly with SPAdes assembler using pair-end Illumina and Nanopore sequencing data. The size of the Illumina data is 18 GB and when I do hybrid assembly after trimming adapters, it showing that Memory is insufficient to run this program and requires at least 389 GB RAM. But, I am having 188 GB RAM.
    Then, I used BBmerge commands to sort out this issue. But when I run pair-end reads using bbmerge-auto.sh command, its saying that "try increasing the -Xmx flag and using tool specific memory related parameters".

    Please suggest me how to do hybrid assembly with both Illumina and Nanopore data?
  • SNPsaurus
    Registered Vendor
    • May 2013
    • 525

    #2
    Is your 18 Gb of Illumina data pure mtDNA or mtDNA plus nuclear? If it is pure mtDNA then you probably have many thousand-fold coverage of the mitochondrial genome and you can reduce the number of reads used as input (bbtools reformat will do that if you give the paired-end reads as in= and in2= and then use reads=10000000 or some other number). If it is not pure mtDNA then you probably still have high coverage and could cut the input by 3/4s.
    Providing nextRAD genotyping and PacBio sequencing services. http://snpsaurus.com

    Comment

    • bioramg
      Member
      • Apr 2014
      • 12

      #3
      Thank you so much for your kind reply.

      Its whole genomic data that includes, chloroplast, mitochondrial and nuclear genome data. How can I choose 3/4 read? Could you please suggest.

      Thank you so much.

      With warm regards,
      Raman. G

      Comment

      • SNPsaurus
        Registered Vendor
        • May 2013
        • 525

        #4
        I would use bbtools reformat to set the number of desired reads.
        Providing nextRAD genotyping and PacBio sequencing services. http://snpsaurus.com

        Comment

        • bioramg
          Member
          • Apr 2014
          • 12

          #5
          Could you please give me with some examples. I am new to bioinformatics. How to set desired reads...

          Comment

          • SNPsaurus
            Registered Vendor
            • May 2013
            • 525

            #6
            You used bbmerge already, so it sounds like you have bbtools installed.
            reformat.sh in=read1,fq,gz in2=read2.fq.gz out=read1_fewer.fq.gz out2=read2_fewer.fq.gz reads=10000000
            Providing nextRAD genotyping and PacBio sequencing services. http://snpsaurus.com

            Comment

            • bioramg
              Member
              • Apr 2014
              • 12

              #7
              Thank you so much.

              by the way, I have used this command to estimate genome size using the following command:
              java -ea -Xmx188g -Xms188g -cp /home/pmslab/Desktop/Raman/bin/bin/bbmap/current/ jgi.KmerCountExact in=trim_CKEMJ1.fastq in2=trim_CKEMJ2.fastq khist=khist_trim_CKEMJ.txt peaks=peaks_trim_CKEMJ.txt -Xmx188g

              But i got an error message:
              ./../bin/bin/bbmap/kmercountexact.sh: line 138: 30199 Killed

              What could be the problem?
              Thank you.

              Comment

              Latest Articles

              Collapse

              • seqadmin
                New Genomics Tools and Methods Shared at AGBT 2025
                by seqadmin


                This year’s Advances in Genome Biology and Technology (AGBT) General Meeting commemorated the 25th anniversary of the event at its original venue on Marco Island, Florida. While this year’s event didn’t include high-profile musical performances, the industry announcements and cutting-edge research still drew the attention of leading scientists.

                The Headliner
                The biggest announcement was Roche stepping back into the sequencing platform market. In the years since...
                03-03-2025, 01:39 PM
              • seqadmin
                Investigating the Gut Microbiome Through Diet and Spatial Biology
                by seqadmin




                The human gut contains trillions of microorganisms that impact digestion, immune functions, and overall health1. Despite major breakthroughs, we’re only beginning to understand the full extent of the microbiome’s influence on health and disease. Advances in next-generation sequencing and spatial biology have opened new windows into this complex environment, yet many questions remain. This article highlights two recent studies exploring how diet influences microbial...
                02-24-2025, 06:31 AM

              ad_right_rmr

              Collapse

              News

              Collapse

              Topics Statistics Last Post
              Started by seqadmin, 03-20-2025, 05:03 AM
              0 responses
              17 views
              0 reactions
              Last Post seqadmin  
              Started by seqadmin, 03-19-2025, 07:27 AM
              0 responses
              18 views
              0 reactions
              Last Post seqadmin  
              Started by seqadmin, 03-18-2025, 12:50 PM
              0 responses
              19 views
              0 reactions
              Last Post seqadmin  
              Started by seqadmin, 03-03-2025, 01:15 PM
              0 responses
              185 views
              0 reactions
              Last Post seqadmin  
              Working...