Hello

I have sequenced a BAC clone with PacBio RSII To make the assembly I used Facon through pbbioconda and for polishing I used quiver

To have an estimation of the consensus quality I remap the original bam reads file against the consensus fasta file

How to estimate a mean quality value, in other world, a consensus Phred score for the base calls of the consensus ... :-)

Thank you very much in advance Philippe