Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,125,8250TC55.6% 13.7 / 20.2 18G165G (GGA→GGGflgAflagellar basal body P‑ring formation chaperone FlgA
Reads supporting (aligned to +/- strand):  ref base T (6/2);  new base C (0/10);  total (6/12)
Fisher's exact test for biased strand distribution p-value = 1.51e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.11e-03
Rejected as consensus: Frequency below/above cutoff threshold.

GCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGC  >  NZ_CP009273/1125745‑1125894
                                                                                |                                                                     
gCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCtttt                                                              >  2:203200/1‑90 (MQ=255)
                     gcCTGCCCTTCTGCGTTGGCGCTAACCCCATCCCCGCTGGCGATCACATTGACGCGTTGCCCCGCTTTTCCCCGCCATGCCTGGCGAAAc                                         <  2:216355/90‑1 (MQ=255)
                        tGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTgg                                      >  1:270799/1‑90 (MQ=255)
                         ccccTTTTGCTTTGGCGCAAACCCCCCCCCCGCGGGCGATCACTTGGCCGCGTTGCCCCCGTTTTACCCGCCATGCCTGGCGAAACTggg                                     <  1:803086/88‑1 (MQ=255)
                             ttCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTaa                                 >  1:753362/1‑90 (MQ=255)
                              tCTGGGTGGGCGCAAACCCCACCCCCGCGGGCGATCACTTTGCCGCGTTGCCCCTTTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAAc                                <  2:650999/90‑1 (MQ=255)
                              tCTGCGTTGGCGCTAACCCCACCCCCGCGGGCGATCACATTGACGCGTTGCCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAAc                                <  2:681228/90‑1 (MQ=255)
                               cTGCGTTGGCGCAAACCCCACCCCCGCTGGCGATCACATTGACGCGTTGCCCCGTTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACt                               <  2:736279/90‑1 (MQ=255)
                                tGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTg                              >  1:759015/1‑90 (MQ=255)
                                    ttGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATa                          >  2:669864/1‑90 (MQ=255)
                                           tAACCCCACCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGAt                   <  2:269920/90‑1 (MQ=255)
                                           tAAACCCACCCCCGCTGGCGATCACATTGCCGCGTTGCCCCGTTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGAt                   <  2:440307/90‑1 (MQ=255)
                                             accccATCCCCGCTGGCGATCACATTGCCGCGTTGCCCCGTTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCg                 <  2:813645/88‑1 (MQ=255)
                                             accccATCCCCGCTGGCGATCACATTGACGCTTTGTCCCGTTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCg                 <  2:91392/88‑1 (MQ=255)
                                             accccACCCCCGCTGGCGATCACATTGCCGCGTTGCCCCGTTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCg                 <  2:563486/88‑1 (MQ=255)
                                             accccACCCCCGCGGGCGATCACTTTGCCGGTTTGCCCCTTTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCg                 <  1:669864/88‑1 (MQ=255)
                                                        gCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATc      >  1:20051/1‑90 (MQ=255)
                                                            gCGATCACATTGACGCGTTGCCCCGCTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCgcgc  <  2:297829/90‑1 (MQ=255)
                                                                                |                                                                     
GCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGC  >  NZ_CP009273/1125745‑1125894

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