Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,562,0430AC51.8% 22.6 / 35.4 27G276G (GGT→GGG) eutGethanolamine utilization ethanol dehydrogenase EutG
Reads supporting (aligned to +/- strand):  ref base A (10/3);  new base C (0/14);  total (10/17)
Fisher's exact test for biased strand distribution p-value = 3.39e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.43e-05
Rejected as consensus: Frequency below/above cutoff threshold.

GTTGGCAGCAACATGGCGTTCGCGAGACCGTGCGGAATATGCAGCGCCGCGCCCGGCTGATGCGCCATCGCGTGGCACAACCCAAGACCCGCACTGGAAAACGCCATTCCCGCCATACATGAAGCCAGCAACATGCTCTCGCGCGCGGCAAGGTCGTGACCGTAGCCCACCGCT  >  NZ_CP009273/2561957‑2562130
                                                                                      |                                                                                       
gTTGGCAGCAACATGGCGTTCGCGAGACCGTGCGGAATATGCAGCGCCGCGCCCGGCTGATGCGCCATCGCGTGGCACAACCCAAGAccc                                                                                      <  2:681224/90‑1 (MQ=255)
gTTGGCAGCAACATGGCGTTCGCGAGACCGTGCGGAATATGCAGCGCCGCGCCCGGCTGATGCGCCATCGCGTGGCACAACCCAAGAccc                                                                                      >  2:380096/1‑90 (MQ=255)
  tGGCAGCAACATGGCGTTCGCGAGACCGTGCGGAATATGCAGCGCCGCGCCCGGCTGATGCGCCATCGCGTGGCACAACCCAAGACCCGc                                                                                    >  2:536776/1‑90 (MQ=255)
        caACATGGCGTTCGCGAGACCGTGCGGAATATGCAGCGCCGCGCCCGGCTGATGCGCCATCGCGTGGCACAACCCAAGACCCGCACTGGa                                                                              >  2:1070183/1‑90 (MQ=255)
                  ttCGCGAGACCGTGCGGAATATGCAGCGCCGCGCCCGGCTGATGCGCCATCGCGTGGCACAACCCAAGACCCGCACTGGAAAACGCCAtt                                                                    >  2:1453518/1‑90 (MQ=255)
                    cgcgAGCCCGTGCGGAATATGCAGCCCCGCCCCCGGCTGATGCGCCACCGCGTGGCACACCCCAAGCCCCGCACTGAAAACCGCCATTcc                                                                  <  1:1415353/90‑1 (MQ=255)
                        agaCCGTGCGAAATATGCGGCGCCGCCCCCGCCTGAGGCGCCCGCGCGGGGCACAACCCAAGCCCCGCACTGGAAAACGCCATTCCCGcc                                                              <  1:1099867/90‑1 (MQ=255)
                          cccGGGCGAAATATGCACCGCCGCCCCCGGCTGATGCCCCACCGCGTGGCCCACCCCAACCCCCGCCCTGGAAACCGCCATTCCCGCCAt                                                            <  1:1302951/89‑1 (MQ=255)
                             gTGCGGAATATGCAGCGCCGCGCCCGGCTGATGCGCCATCGCGTGGCACAACCCAAGACCCGCACTGGAAAACGCCATTc                                                                   <  2:1064868/80‑1 (MQ=255)
                             gTGCGGAATATGCAGCGCCGCGCCCGGCTGATGCGCCATCGCGTGGCACAACCCAAGACCCGCACTGGAAAACGCCATTc                                                                   >  1:1064868/1‑80 (MQ=255)
                                cgaaaTATGCAGCGCCGCGCCCGGCTGATGCCCCATCGCGTGGCACAACCCAAGCCCCCCACTGAAAAACCCCATTCCCGCCATACATGa                                                      <  2:1570109/87‑1 (MQ=255)
                                cGGAATATGCGGCCCCGCGCCCGGCTGACGCCCCACCGCGTGCCACACCCAAAGCCCCGCACGGAAAACCGCCATTCCCGCCATACATGa                                                      <  1:203341/90‑1 (MQ=255)
                                 ggAATATGCAGCGCCGCGCCCGGCTGACGCCCCCCCGCGTGGCACAACCCAAGCCCCGCACTGAAAAACGCCATTCCCGCCATACATGaa                                                     <  1:1226775/90‑1 (MQ=255)
                                  gAATATGCCGCCCCGCGCCCGGGGGAGGCCCCACCGGGGGGCACACCCCAAGCCCCCCATTGAAAAACCCCATTCCCGCCATACATGAAg                                                    <  1:536776/90‑1 (MQ=255)
                                    atacGCAGCCCCGCCCCCGGCTGATGCCCCACCGCGGGGCACACCCCAAGCCCCGCACGGAAAACCGCCATCCCCGCCATACATGAAGcc                                                  <  1:946204/90‑1 (MQ=255)
                                       tGCAGCGCCGCGCCCGGCTGATGCGCCATCGCGTGGCACAACCCAAGACCCGCACTGGAAAACGCCATTCCCGCCATACATGAAGCCAGc                                               >  2:1401052/1‑90 (MQ=255)
                                       tGCAGCGCCGCGCCCGGCTGATGCGCCATCGCGTGGCACAACCCAAGACCCGCACTGGAAAACGCCATTCCCGCCATACATGAAGCCAGc                                               >  2:1562831/1‑90 (MQ=255)
                                       tGCAGCGCCGCGCCCGGCTGATGCCCCATCGCGTGGCACAACCCAACCCCCGCACTGAAAACCCCCATTCCCGCCATACATGAAGCCAGc                                               <  1:1453518/90‑1 (MQ=255)
                                       tGCAGCCCCGCGCCCGGCTGATGCCCCACCGCGTGCCACAACCCAAGCCCCGCACTGAAAAACCCCATTCCCGCCATACATGAAGCCAGc                                               <  1:173166/90‑1 (MQ=255)
                                        gCAGCCCCGCCCCGGGCTGAGGCCCCCCCGCGGGGCACAACCCAAGCCCCGCCCGGAAAAACGCCATCCCCGCCATACATGAAGCCAGCa                                              <  1:1293120/90‑1 (MQ=255)
                                           gcgcccccccCGGCTGATGCCCCATCCCGTTGCAAAACCCAAGCCCCCCACTGAAAACCGCCATTCCCGCCATACATGAAGCCAGCAACa                                           <  2:1477447/90‑1 (MQ=255)
                                            cgccgcgccCGGCTGATGCGCCATCGCGTGGCACAACCCAAGACCCGCACTGGAAAACGCCATTCCCGCCATACATGAAGCCAGCAACAt                                          >  1:234975/1‑90 (MQ=255)
                                                      ggCTGATGCGCCACCGCGTGCCACACCCCAAGCCCCCCACTGAAAACCGCCATTCCCGCCATACATGAAGCCAGCAACATGCTCTcgcgc                                <  1:1453339/90‑1 (MQ=255)
                                                      ggCTGAGGCCCCATCGCGGGCCACACCCCAAGACCCGCACTGGAAACCGCCATTCCCGCCATACATGAAGCCAGCAACATGCTCTcgcgc                                <  2:978543/90‑1 (MQ=255)
                                                             gcgcCGCCGCGTGGCACACCCCAAGCCCCGCACGGAAAAACGCCATTCCCGCCATACATGAAGCCAGCAACATGCTCTCGCGCGCGGCaa                         <  1:1312253/90‑1 (MQ=255)
                                                                                    aGACCCGCACTGGAAAACGCCATTCCCGCCATACATGAAGCCAGCAACATGCTCTCGCGCGCGGCAAGGTCGTGACCGTAGCCCACCGCt  >  1:404031/1‑90 (MQ=255)
                                                                                    aGACCCGCACTGGAAAACGCCATTCCCGCCATACATGAAGCCAGCAACATGCTCTCGCGCGCGGCAAGGTCGTGACCGTAGCCCACCGCt  >  2:1575505/1‑90 (MQ=255)
                                                                                      |                                                                                       
GTTGGCAGCAACATGGCGTTCGCGAGACCGTGCGGAATATGCAGCGCCGCGCCCGGCTGATGCGCCATCGCGTGGCACAACCCAAGACCCGCACTGGAAAACGCCATTCCCGCCATACATGAAGCCAGCAACATGCTCTCGCGCGCGGCAAGGTCGTGACCGTAGCCCACCGCT  >  NZ_CP009273/2561957‑2562130

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