Peregrine icon indicating copy to clipboard operation
Peregrine copied to clipboard

Published 20 hours ago • verified publisher icon cschin

shmr_map issue

Open smoretti opened this issue 5 years ago • 4 comments

I have the same issue as the one mentioned in #10 Do you know what could cause this?

smoretti avatar Oct 04 '19 11:10 smoretti

@smoretti do you have some background of the input data of your assembly run? CCS read, coverage? It looks that you don't get a draft assembly from the earlier step to me.

cschin avatar Oct 28 '19 14:10 cschin

@cschin The input data was generated using two SMRT cells v3. The mean coverage over the two SMRT cells was ~ 68x. The two .subreads.bam files were used as input to create a single fasta file using Pacbio's bam2fastx.

ghost avatar Oct 28 '19 16:10 ghost

@taneKFL it seems to me your data is PacBio's subreads rather than CCS reads. Currently, Peregrine is designed for accurate reads (>99% accuracy > 10kb length). The subreads won't have the accuracy for Peregrine to assemble.

cschin avatar Oct 29 '19 05:10 cschin

@cschin okay - so I'll try running pbccs on both my movies, convert them to a single fastq using bam2fastx, then retry running peregrine.

ghost avatar Oct 30 '19 08:10 ghost