I have aligned raw 1D nanopore reads to a bacterial reference assembly using BWA –MEM and got the sam and bam file.
Is there a quick way to get the of percentage identity for each read from the sam/bam file ???
I would like to generate an histogram with reads count and the percent identity to the reference.