Question: SSPACE assembling Genome!
1
gravatar for fufuyou
3.9 years ago by
fufuyou110
United States
fufuyou110 wrote:

Hi ,

I am assembling a genome using SSPACE. I met a problem as:

perl /scratch/snyder/f/fu115/Genome_assembly/SSPACE/41SSPACE-STANDARD-3.0_linux-x86_64/SSPACE-STANDARD-3.0_linux-x86_64/SSPACE_Standard_v3.0.pl -l libraries.txt -s contigs_abyss.fasta -k 5 -a 0.7 -x 0 -T 20 -b ecoli_scaffolds_no_extension
Your inserted inputs on [SSPACE_Standard_v3.0_linux] at Thu Nov 19 16:55:24 2015:
Required inputs: 
        -l = libraries.txt
                Number of paired files = 1
        -s = contigs_abyss.fasta
        -b = ecoli_scaffolds_no_extension

Optional inputs:
        -x = 0
        -z = 0
        -k = 5
        -g = 0
        -a = 0.7
        -n = 15
        -T = 20
        -p = 0


=>Thu Nov 19 16:55:24 2015: Reading, filtering and converting input sequences of library file initiated
                                                
------------------------------------------------------------

=>Thu Nov 19 16:55:50 2015: Storing contigs to format for scaffolding
                
LIBRARY lib1
------------------------------------------------------------

=>Thu Nov 19 16:55:51 2015: Reading contig file
                
=>Thu Nov 19 16:55:51 2015: Building BWA index for contigs
**************************************************

Process 'mapping reads' failed on Thu Nov 19 16:55:51 2015

I do not know what it is happen.

Thanks,

Fuyou

myposts assembly • 2.3k views
ADD COMMENTlink modified 2.2 years ago by Biostar ♦♦ 20 • written 3.9 years ago by fufuyou110
1
gravatar for tstrobaugh
3.8 years ago by
tstrobaugh10
USA/Phila/USDA
tstrobaugh10 wrote:

-T=20? Do you really have 20 cores or are you using hyperthreading? Back down to 8 or so. I'm only answering this question because I am having the same problem. Anyway, I believe this to be a "chmod" problem somewhere with BWA. I got it to work by replacing "BWA" with "bowtie" in the library file. Still trying to get bwa to work, problem is, it worked earlier and I didn't change anything.

ADD COMMENTlink written 3.8 years ago by tstrobaugh10

It is the problem about RAM? I found if I change to -T=1, it worked well. But it is very slow. How about your work now?

Thanks,

Fuyou

ADD REPLYlink written 3.8 years ago by fufuyou110
1
gravatar for sutturka
2.3 years ago by
sutturka150
USA
sutturka150 wrote:

I am sharing my experience with SSPACE as it may be useful for someone.

Data: We have multiple Paired-end and mate-pair libraries from a large eukaryotic genome. We performed assembly through ABySS and ran SSPACE on ABySS contigs to perform scaffolding and contig extension.

Issue: SSPACE scaffolding finished without any problems but contig extension (x = 1) was consistently failing.

Solution: Later, I found SSPACE uses unmapped reads from given library for extension and stores the reads in the memory. Therefore, the higher number of reads being stored in the memory was the issue for failing the SSPACE run. A simple solution is to divide your library into multiple files (I divided PE library file into 4 parts) and sequentially ran contig extension step to finish it successfully. Ultimately, a number of unmapped reads from the library (which are stored in the virtual memory) dictates the SSPACE run and run can be finished on low memory machines by splitting the files (in the case of MP library, I split my file into 14 parts). Hope this helps.

ADD COMMENTlink modified 2.3 years ago • written 2.3 years ago by sutturka150

can you give an example of what your library.txt file looked like?

ADD REPLYlink written 5 months ago by AGE20
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2631 users visited in the last hour