pacbio data SNPs

I am here to ask for your help.I have run a bacteria with Pacbio,and I have finished the genome.I want to find the SNPs between the fitered-subreads from sequencing and my finished_genome.I mean to find at the same postion of the reference,how many reads support the same base with the reference,and how many reads support another base,I want the propotion.Could I use the SMRT .bam file and GATK to generate the result which I want??Please instruct me.

Best Answer


Sign In or Register to comment.