SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
how do i reduce misassemblies in spades? evolino Bioinformatics 0 10-30-2017 08:22 PM
SPAdes 3.9.0 on MacOS Sierra (10.12.1) osilander Bioinformatics 3 11-04-2016 09:39 AM
Complexity Time SPADES assembler mido1951 Bioinformatics 2 01-30-2016 04:12 AM
SPAdes contigs vpi Bioinformatics 3 12-02-2015 03:49 PM
SPAdes: with different read length bio_informatics Illumina/Solexa 6 06-10-2015 08:04 AM

Reply
 
Thread Tools
Old 02-08-2020, 04:45 PM   #1
bioramg
Member
 
Location: South Korea

Join Date: Apr 2014
Posts: 10
Default 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?
bioramg is offline   Reply With Quote
Old 02-09-2020, 07:28 PM   #2
SNPsaurus
Registered Vendor
 
Location: Eugene, OR

Join Date: May 2013
Posts: 521
Default

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
SNPsaurus is offline   Reply With Quote
Old 03-01-2020, 10:15 PM   #3
bioramg
Member
 
Location: South Korea

Join Date: Apr 2014
Posts: 10
Default

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
bioramg is offline   Reply With Quote
Old 03-04-2020, 12:38 PM   #4
SNPsaurus
Registered Vendor
 
Location: Eugene, OR

Join Date: May 2013
Posts: 521
Default

I would use bbtools reformat to set the number of desired reads.
__________________
Providing nextRAD genotyping and PacBio sequencing services. http://snpsaurus.com
SNPsaurus is offline   Reply With Quote
Old 03-04-2020, 11:26 PM   #5
bioramg
Member
 
Location: South Korea

Join Date: Apr 2014
Posts: 10
Default

Could you please give me with some examples. I am new to bioinformatics. How to set desired reads...
bioramg is offline   Reply With Quote
Old 03-05-2020, 11:30 AM   #6
SNPsaurus
Registered Vendor
 
Location: Eugene, OR

Join Date: May 2013
Posts: 521
Default

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
SNPsaurus is offline   Reply With Quote
Old 03-05-2020, 05:54 PM   #7
bioramg
Member
 
Location: South Korea

Join Date: Apr 2014
Posts: 10
Default

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.
bioramg is offline   Reply With Quote
Reply

Tags
bioinformactics, bioinformatic analysis, bioinformatic tools, ngs analysis platform

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 06:07 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO