Sean had previously attempted to run PB Jelly, but ran into some issues running on Hyak, so I decided to try this on Emu.
Here’s a brief rundown of how this was run:
- Default PB Jelly settings (including default settings for blasr).
- Illumina reference FASTA: Sean’s Platanus kmer=22 assembly
- PacBio reads for mapping
- Protocol.xml file needed for PB Jelly to run
See the Jupyter Notebook for full details of run (see Results section below).
This completed very quickly (like, just a couple of hours). I also didn’t experience the woes of multimillion temp file production that killed Sean’s attempt at running this on Mox (Hyak).
However, it doesn’t seem to have produced an assembly!
Looking through the output, it seems as though it didn’t produce an assembly because there weren’t any gaps to fill in the reference. This makes sense (in regards to the lack of gaps in the reference Illumina assembly) because I used the Platanus contig FASTA file (i.e. not a scaffolds file). I didn’t realize PB Jelly was just designed for gap filling. Guess I’ll give this another go using the BGI scaffold FASTA file and see what we get.
Jupyter Notebook (GitHub): 20171113_emu_pbjelly_22mer_plat.ipynb