euler-sr scaffolding

[tshea]

I see that mateTransformGraph has a -scaffold option:

[-scaffold] Use mate-pairs to scaffold disconnected contigs.

Is there a way to obtain fasta sequence of the resulting scaffolds? Or is there documentation to explain how one might parse any of the .graph or .edge or .path files with the goal of constructing a scaffolds fasta file?

Thanks much.

[mchaisso]

I see that mateTransformGraph has a -scaffold option:

[-scaffold] Use mate-pairs to scaffold disconnected contigs.

Is there a way to obtain fasta sequence of the resulting scaffolds? Or is there documentation to explain how one might parse any of the .graph or .edge or .path files with the goal of constructing a scaffolds fasta file?

Thanks much.

The contigs from the scaffold should be in the "reads.fasta.contig" file that is output at the end of assembly, where "reads.fasta" is the input file. If there is no reads.fasta.contig, the assembly crashed (somebody may be posting an updated scaffolder soon).

[tshea]

The contigs from the scaffold should be in the "reads.fasta.contig" file that is output at the end of assembly, where "reads.fasta" is the input file. If there is no reads.fasta.contig, the assembly crashed (somebody may be posting an updated scaffolder soon).

Thanks Mark.

The assembly did finish and there is a reads.fasta.contig file.

To clarify, is there a way to obtain a fasta file of the scaffolds themselves (linked contigs separated by the estimated number of N's in the gap). I assume that even after the mateTransformGraph stage there are still contigs which, while not overlapping, are at least linked by multiple paired reads (in my case we have 4kb pairs).

You mentioned there may be an updated scaffolder soon - is there something at present that determines and prints out scaffold fasta sequence?

Thanks again.