Question: Using SSPACE: how to calculate 'expected/observed insert size between paired reads and the minimum allowed error?'
0
gravatar for freuv
3 days ago by
freuv0
freuv0 wrote:

Hi,

I am using SSPACE for scaffolding of AbySS contigs -- the tool requires a text library file with some info re: the reads.

It is unclear what and how the following required numbers for this text library file are calculated.

  • expected/observed insert size between paired reads
  • minimum allowed error

Any guidance would be appreciated. Thanks!

F

ADD COMMENTlink written 3 days ago by freuv0

You can calculate the insert sizes using BBMap (C: Target fragment size versus final insert size ). Adjust value of k= depending size of your reads.

ADD REPLYlink written 3 days ago by genomax39k

If you don't mind, can you please explain what these two values mean?

  • expected/observed insert size between paired reads
  • minimum allowed error

My raw reads are 150bp -- what value of 'k' would I use?

Thanks for your help!

ADD REPLYlink written 2 days ago by freuv0

I don't know the full context of this question so here goes. Expected insert size is something you will have determined during your library prep/QC (using e.g. bioanalyzer etc). Observed insert size is what you will get from bbmap.sh or bbmerge.sh method. This would be the reality present in your library. In general, smaller inserts cluster more efficiently on Illumina flowcells, so unless you are very careful about restricting them you will see insert sizes smaller than you expect to see. Take a look at the bbmerge guide here.

ADD REPLYlink modified 2 days ago • written 2 days ago by genomax39k
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: 1492 users visited in the last hour