pacbio data SNPs

Hello,everyone,
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

Answers

Sign In or Register to comment.