I've downloaded this and started using it. I have to say, I am either doing something wrong or it's not nearly as useful as I was hoping. Hopefully I'm doing something wrong here.
This is what I have tried.
Took 4 SMRT cells of long read data. Isolated out reads that are >5000. Imported these >5K reads into SMRT Portal. Ran RS_Mapping_QVs to get the new cmp.h5 file.
Finally, ran Quiver and output the reads as fastq.
Was left with about 12X coverage of my organism of interest. Tried Celera to assemble, and it essentially failed...most of the reads aren't even used.
I've taken the same data, used PacBioToCA to correct with Illumina, and assembled using about the same coverage and got less than 100 contigs...
Ideas?
This is what I have tried.
Took 4 SMRT cells of long read data. Isolated out reads that are >5000. Imported these >5K reads into SMRT Portal. Ran RS_Mapping_QVs to get the new cmp.h5 file.
Finally, ran Quiver and output the reads as fastq.
Was left with about 12X coverage of my organism of interest. Tried Celera to assemble, and it essentially failed...most of the reads aren't even used.
I've taken the same data, used PacBioToCA to correct with Illumina, and assembled using about the same coverage and got less than 100 contigs...
Ideas?
Comment