DNA sequence of Streptomyces sp. RKAG293 from Arctic marine sediment
Contigs assembly was done using what is referred as theHGAPworkflow(Chin et al. 2013). Briefly, raw subreads where generated from raw .bax.h5 PacBio data files. A subread length cutoff value was extracted from subreads and used into the preassembly (BLASR(Chaisson and Tesler 2012)) step which consists of aligning short subreads on long subreads. Since errors in PacBio reads is random, the alignment of multiple short reads on longer reads allows to correct sequencing error on long reads. These long corrected reads are then used as seeds into assembly (Celera Assembler(Myers et al. 2000)) which gives contigs. These contigs are then 'polished' by aligning raw reads on contigs (BLASR) that are then processed through a variant calling algorithm (Quiver) that generates high quality consensus sequences using local realignments and PacBio quality scores.