This is a test version of Biostars. For the public version, visit https://www.biostars.org.
SSPACE option (-k and -a) for scaffolding

Hi,

I am trying to play with SSPACE's parameters, especially the -k and -a option.

According to the manual:

-k  Minimum number of links (read pairs) to compute scaffold (default -k 5, optional)

-a  Maximum link ratio between two best contig pairs *higher values lead to least accurate scaffolding* (default -a 0.7, optional)

But I am not sure to understand those 2 options. If someone has better and clear explanation, it will be helpful.

Thanks!

sspace
• 2,032 views • link
written by Picasa Picasa

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • Lower Assembly Quality Despite Higher Coverage?
    written by AXR • AXR
    Hi everyone, I wanted to see the effect of coverage on the assembly quality to see at which point there is a diminishing return. I …
  • How do I reduce the length of a contig?
    written by kvc0004 • kvc0004
    I have recently assembled my shotgun metagenomic reads using metaSPADES without changing any default settings. I then upload them to MG-RAST for annotation. However, some …
  • Bowtie run error
    written by mikemwanga6 • mikemwanga6
    I am running SSPACE after assembly of reads using SPADES. However SSPACE returns the following error. command: perl ~/my_scripts/SSPACE/SSPACE_Basic.pl -l libraries.txt -s contigs.fasta -k 5 …
  • AlignGraph doesn't response but i have what it needs
    written by Calvin • Calvin
    **This is the command line i was running** ``` AlignGraph --read1 W_S1_L001_R1_001.fasta --read2 W_S1_L001_R2_001.fasta --contigs 04.break.broken_assembly.fa --genome sequence.fasta.txt --distanceLow 300 --distanceHigh 1300 --extendedContig S.Pneu_extendedContigs.fa --remainingContigs …
  • SSPACE assembling Genome!
    written by fufuyou fufuyou
    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 …
  • Scaffolding with SSPACE returned scaffold the same as initial contigs set
    written by pbigbig pbigbig
    <p>Hi all,</p> <p>I have used Minia to assemble a contigs set from my paired end reads, as Minia instructed that it doesn&#39;t use pairing information …
  • choice of k-mer size for metagenomic assembly
    written by qiyunzhu qiyunzhu
    Dear community, I am doing de novo assembly of some metagenomic datasets. They are Illumina NextSeq reads (paired-end, 150 bp per read). I have tried …
  • Post Genome Assembly
    written by div • div
    I am working on Post Genome Improvement of four strains of *Staphylococcus aureus*. Here my job is to upgrade the contigs into scaffold of those …
  • Sspace And -K Parameter
    written by Ketil Ketil
    <p>I'm trying to scaffold some assemblies with SSPACE 2.1. According to the manual:</p> <pre><code>The -k minimal links and -a maximum link ratio --------------- Two parameters …
  • Abyss Speed-Up Tricks
    written by Darked89 Darked89
    <p>I am trying to assemble a 30-35Mbp diploid genome using Abyss from HighSeq Illumina runs. It takes a very long time (days) to compute assembly …