Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273782,3110AC53.3% 10.9 / 16.9 15G247G (GGT→GGG) gpmA2,3‑diphosphoglycerate‑dependent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base A (6/1);  new base C (0/8);  total (6/9)
Fisher's exact test for biased strand distribution p-value = 1.40e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.15e-03
Rejected as consensus: Frequency below/above cutoff threshold.

TGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAA‑ACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATA  >  NZ_CP009273/782228‑782367
                                                                                    |                                                        
tgtgCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAA‑ACGAATGACGTTTACTTCGCTTTACCCTgg                                                    >  1:328294/1‑90 (MQ=255)
      ccATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAA‑ACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCa                                              >  2:419724/1‑90 (MQ=255)
         ttGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAA‑ACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAAcc                                           >  1:40045/1‑90 (MQ=255)
           gTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAA‑ACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAAccgc                                         >  1:14200/1‑90 (MQ=255)
                caacaaAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAA‑ACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCgctg                                    >  2:476971/1‑90 (MQ=255)
                                   tttGCAGTCGGTTTTCTCATTTAAA‑ACGAATGACGTTTCCTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTc                 <  2:536202/89‑1 (MQ=255)
                                   cTTGCAGTCGGTTTTCTCATTTTAA‑CCAAATGACGTTTCCTCCCCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTc                 <  1:291802/90‑1 (MQ=255)
                                   cTTGCAGTCGGCTTTCTCATTTTAA‑ACGAATGACTTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTc                 <  2:338699/90‑1 (MQ=255)
                                        ggTCGGCTTTCTCTTTTAAA‑ACGATTGGCTTTTACTTCGTTTTCCCCTGTTTTGAACCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCAt            <  2:375173/89‑1 (MQ=255)
                                        aGTCGGTTTTCTCATTTAAA‑ACGAATGCCGTTTCCTTCCCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCAt            <  1:555588/90‑1 (MQ=255)
                                        aGTCGGTTTTCTCATTTAAA‑ACGAATGACTTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCAt            <  2:267889/90‑1 (MQ=255)
                                        aGTCGGCTTTCTCATTTTAA‑ACGAATGACTTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCAt            <  2:519117/90‑1 (MQ=255)
                                               tttCTCCTTTTAAGACGAATTGTTTCTACTCCCCTTTACCCGCATTTTCACCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCa      <  1:748452/90‑1 (MQ=255)
                                                ttCTCATTTTAA‑ACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGa    >  2:418617/1‑90 (MQ=255)
                                                  ctcATTTTAA‑ACGAATGACTTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGata  <  2:432753/90‑1 (MQ=255)
                                                                                    |                                                        
TGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAA‑ACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATA  >  NZ_CP009273/782228‑782367

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑