Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,776,2790TA40.5% 27.0 / 11.4 38D248V (GAT→GTT) gpsANAD(P)H‑dependent glycerol‑3‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (11/11);  new base A (15/0);  total (26/12)
Fisher's exact test for biased strand distribution p-value = 9.04e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTG  >  NZ_CP009273/3776193‑3776365
                                                                                      |                                                                                      
gcgcgcTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGGTGTCGGGACAGGGAAGGACCAAATCGc                                                                                     >  1:442743/1‑90 (MQ=255)
  gcgcTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTTCCCGACTGGTTTTCGGGACCGGTAAACACCAAAACGCCa                                                                                   >  2:86483/1‑90 (MQ=255)
              tCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTTCCCCCCCTGGTGTGGGGAAAGGTAAACACCAAAACGCCAAACccccccccg                                                                       >  2:631402/1‑87 (MQ=255)
              tCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATg                                                                       <  1:405282/90‑1 (MQ=255)
              tCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATg                                                                       <  1:494639/90‑1 (MQ=255)
              accgtgccctgcccGAGCACCATGCCAAAACGGCGGTTACGCGATTGGTTGTCGGTAGAGGTAAGCCCAAAACCGCCAAGCCCCGCCATg                                                                       <  2:606612/86‑1 (MQ=255)
               ccATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACCTGGTGTCGGGACAGGGAAACACCCAAACGCCAAGCcccccccccc                                                                      >  2:405282/1‑86 (MQ=255)
                 aTGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGGCTGGTTGTCGGGACAGGTAAGCACCAAAACGCCAAGCCCCCCCCTGccc                                                                    >  1:412177/1‑90 (MQ=255)
                 aTGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGGCTGGTTGTCGGGACAGGGAAGCACCAGATCGCCAAGCCCCCCCCTGccc                                                                    >  2:691826/1‑90 (MQ=255)
                 aTGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATg                                                                       <  2:249020/87‑1 (MQ=255)
                 aTGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGGACAGGTAAGCACCAGATCGCCAAGCCCCGCCATg                                                                       >  1:249020/1‑87 (MQ=255)
                 aTGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGGACAGGGAAACACCAAAACGCCAAACCCCCCCCTGccc                                                                    >  1:540144/1‑90 (MQ=255)
                   gCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGGAAGCACCAGATCGCCAAGCCCCGCCATGCCCAt                                                                  >  1:579354/1‑90 (MQ=255)
                        gACCGAGCATCATGCCAAAACGGCGGTTACGCGGCTGGGTGTCGGTACAGGGAAGCACCAGAACGCCAAACCCCCCCCTGCCCCCAAAgg                                                             >  2:128993/1‑90 (MQ=255)
                        gACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGGACAGGGAAGCACCAGAACGCCAAGCCCCCCCCTGCCCCTAAAgg                                                             >  2:111218/1‑90 (MQ=255)
                        gACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGGTGTCGGGACAGGGAAGCACCAAAACGCCCAACCCCCCCCTGCCCCTaaaag                                                             >  1:447700/1‑88 (MQ=255)
                        gACCGAGCATCATGCCAAAACGGCGGGTACGCGACTGGTTGTCGGGACAGGGAAGCACCAGAACGCCCAACCCCCCCACGCCCATaaaag                                                             >  2:494639/1‑88 (MQ=255)
                            gAGCATCATGCCAAAACGGCGGTTACGCGGCTGGTTGTCGGGACAGGTAAGCACCAAAACGCCAAGCCCCCCCCTGCCCCTAAAAGGGGc                                                         >  2:647029/1‑90 (MQ=255)
                             aGCATCATCCCAAAACGGCGTTTCCGCGACTGGTTGTCGGTACGGGTAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCa                                                        <  2:442743/90‑1 (MQ=255)
                              gCATCATGCCAAAACGGCGGTTACGCGGCTGGTTGTCGGCGCAGGGAAACACCAAAAACCCCACCCCCCCCCTGCCCCTAAAAGGGGCAg                                                       >  2:111194/1‑90 (MQ=255)
                               catcatGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAgg                                                      <  1:691826/90‑1 (MQ=255)
                                   atGCCAAAACGGCGGTTACGCGACTGGTTGTCGGGACAGGGAAACACCCAAACGCCAAGCCCCCCCCTGCCCCTAAAGGGGGGAGGGGCg                                                  >  2:343009/1‑90 (MQ=255)
                                   atGCCAAAACGGCGGTTACGCGACTGGGTGTCGGGACAGGGAAGCACCAAAACGCCAAACCCCCCCCTCCCCCTAAAAGGGGGAGGGGCg                                                  >  2:397727/1‑90 (MQ=255)
                                    tGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCgg                                                 <  1:343009/90‑1 (MQ=255)
                                    tGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCgg                                                 <  1:86483/90‑1 (MQ=255)
                                         aaaCGGCGGTTACGCGACTGGGTGGCGGGACCGGGAAGCACCAAAACCCCAACCCCCGCCATGCCCCAAAAGGAGGCAGGGGCGGCAccc                                            >  1:565400/1‑90 (MQ=255)
                                              gcggTTCCGCGCCTGGTTGTGGGTACGGGAAAGCCCCAGATCGCCAACCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAgcgc                                       <  2:540144/90‑1 (MQ=255)
                                              gcggTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAgcgc                                       <  1:128993/90‑1 (MQ=255)
                                              gcggTTACGCGACTGGGTGTTGGTGCAGGGAATCACCAAAACGCCCACCCCCCCCCTGCCCCCAAAGGGGGGGGCGGCGGGACCCatacg                                       >  2:409260/1‑86 (MQ=255)
                                                           tGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGAcgc                          >  2:162440/1‑90 (MQ=255)
                                                                       cAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTcagcc              >  2:247155/1‑90 (MQ=255)
                                                                         ggTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCTGCACCCAGCGCCGCACCAAGACGCGACATTTcagccag            >  1:512570/1‑90 (MQ=255)
                                                                            aaGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGccc         >  2:642800/1‑90 (MQ=255)
                                                                              gCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCAc       >  1:431581/1‑90 (MQ=255)
                                                                               cACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACg      <  1:542767/90‑1 (MQ=255)
                                                                                 caaGATCCCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACggg    <  2:107508/88‑1 (MQ=255)
                                                                                   aGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCGCGGGTg  >  1:44441/1‑90 (MQ=255)
                                                                                   aGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTg  >  1:682754/1‑90 (MQ=255)
                                                                                      |                                                                                      
GCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTG  >  NZ_CP009273/3776193‑3776365

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