Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,125,8300CT51.3% 24.8 / 50.4 37A164T (GCG→ACG) flgAflagellar basal body P‑ring formation chaperone FlgA
Reads supporting (aligned to +/- strand):  ref base C (15/3);  new base T (0/19);  total (15/22)
Fisher's exact test for biased strand distribution p-value = 8.71e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.22e-06
Rejected as consensus: Frequency below/above cutoff threshold.

CTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGAT  >  NZ_CP009273/1125741‑1125919
                                                                                         |                                                                                         
cTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGc                                                                                           >  1:1421333/1‑90 (MQ=255)
      gACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTAc                                                                                     >  2:46803/1‑90 (MQ=255)
            tGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCa                                                                               >  2:1861610/1‑90 (MQ=255)
               tttGTTCAGCGCCTGCCTTTCTGCGTTGGCGCAAACCCCACCCCCGCGGGCGATCACATTGCCGCGTTGTCCCGTTTTTCCCCGCCATGc                                                                            <  1:1426506/89‑1 (MQ=255)
               aTTGTTCAGCGCCTGCCCTTCTGGGTTGGCGCTAAACCCATCCCCGCGGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGc                                                                            <  2:1281822/90‑1 (MQ=255)
                  gTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTg                                                                         >  2:92202/1‑90 (MQ=255)
                  gTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTg                                                                         >  2:143259/1‑90 (MQ=255)
                       gcgcCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGaa                                                                    >  2:1426507/1‑90 (MQ=255)
                       gcgcCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGaa                                                                    >  1:320715/1‑90 (MQ=255)
                            tGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTgg                                                               >  1:69786/1‑90 (MQ=255)
                                 ttCTGGGTTGGCGCAAAACCCATCCCCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTaa                                                          <  2:1421334/90‑1 (MQ=255)
                                 ttCTGCTTGGGCCCAAACCCCACCCCCGCGGGCGACCACATTGACGCTTTGCCCCGTTTTTCCCCGCCATGCCTGGCGAAACTGGGTTaa                                                          <  1:192/90‑1 (MQ=255)
                                  tcggCGTTGGCGCAAACCCCACCCCCGCGGGCGATCACATTGACGCGTTGCCCCGTTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAAc                                                         <  1:430936/87‑1 (MQ=255)
                                        tggccgcAAACCCCATCCCCGCGGGGAATCATTTTGGCGGGTTCCCCCGTTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATa                                                   <  1:111949/86‑1 (MQ=255)
                                         tGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAg                                                  >  2:699218/1‑90 (MQ=255)
                                             gcaaaaCCCACCCCCGCTGGCGATCACATTGACGGGTTGTCCCGTTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTg                                              <  2:758286/87‑1 (MQ=255)
                                               tAACCCCATCCCCGCTGGCGATCACATTGCCGGTTTGCCCCGTTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGAt                                            <  2:69786/90‑1 (MQ=255)
                                               tAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGAt                                            >  1:776477/1‑90 (MQ=255)
                                               aaaccccACCCCCGCTGGGGATCACATTGCCGCGTTGCCCCGTTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGAt                                            <  2:531045/86‑1 (MQ=255)
                                                 accccATCCCCGCTGGGGATCACATTGACGCGTTGCCCCGTTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCg                                          <  2:975886/88‑1 (MQ=255)
                                                 accccATCCCCGCTGGCGATCACATTGACGGGTTGTCCCGTTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCg                                          <  2:715969/88‑1 (MQ=255)
                                                 aaCCCATCCCCGCGGGCGATCACATTGACGTTTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCg                                          <  2:643290/90‑1 (MQ=255)
                                                    ccATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGt                                       >  2:1509705/1‑90 (MQ=255)
                                                        ccccGCGGGCGACCACATTGACGGTTTGCCCCTTTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATa                                   <  1:699218/88‑1 (MQ=255)
                                                            gCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATc                               >  2:1588465/1‑90 (MQ=255)
                                                            gCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATc                               >  1:186181/1‑90 (MQ=255)
                                                            gCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATc                               >  1:439757/1‑90 (MQ=255)
                                                                gCGATCAGATGGACGCTTTGCCCTGTTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCgcgc                           <  1:414159/90‑1 (MQ=255)
                                                                gCGATCACATTGCCGCGTTGTCCCGTTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCgcgc                           <  2:186181/90‑1 (MQ=255)
                                                                 cGATCACATTGACGTGTTGCCCCGTTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCa                          <  2:1813695/90‑1 (MQ=255)
                                                                        aTTGCCGCGTTGCCCCGTTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAAt                   <  1:1507070/90‑1 (MQ=255)
                                                                         ttGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGAGTGATAGATCGCGCAGGCTAATg                  >  1:370207/1‑90 (MQ=255)
                                                                                 ttGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCgcgc                           >  2:1567742/1‑73 (MQ=255)
                                                                                 ttGCCCCGTTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCgcgc                           <  1:1567741/73‑1 (MQ=255)
                                                                                  tcccccGTTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATcaaca         <  2:62171/87‑1 (MQ=255)
                                                                                  tGTCCCTTTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATcaaca         <  1:1063653/90‑1 (MQ=255)
                                                                                  tGTCCCGTTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATcaaca         <  2:265405/90‑1 (MQ=255)
                                                                                  tGTCCCGTTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATcaaca         <  1:202418/90‑1 (MQ=255)
                                                                                         cTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGAGCGCGCAGGCTAATGGCATCAACAAGttgat  >  2:1777867/1‑90 (MQ=255)
                                                                                         |                                                                                         
CTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGAT  >  NZ_CP009273/1125741‑1125919

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