Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,291,0760AC52.4% 11.4 / 26.3 21V141G (GTG→GGG) yraQpermease
Reads supporting (aligned to +/- strand):  ref base A (8/2);  new base C (0/11);  total (8/13)
Fisher's exact test for biased strand distribution p-value = 2.21e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.06e-04
Rejected as consensus: Frequency below/above cutoff threshold.

CCAGCCGAGGACAAAGCCCATAAACACCAGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGCCCGGCAACGAAAACAGCGTTCCCAG  >  NZ_CP009273/3290988‑3291158
                                                                                        |                                                                                  
ccAGCCGAGGACAAAGCCCATAAACACCAGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGacac                                                                                   >  1:890658/1‑90 (MQ=255)
      gAGGACAAAGCCCATAAACACCAGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTg                                                                             >  2:392787/1‑90 (MQ=255)
                   aTAAACACCAGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTGGCGACGCATTccc                                                                >  2:25172/1‑90 (MQ=255)
                            aGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTGGCGACGCATTCCCGCCGCGAcc                                                       >  2:354660/1‑90 (MQ=255)
                               gTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTGGCGACGCATTCCCGCCGCGACCGGa                                                    >  1:349930/1‑90 (MQ=255)
                                 cggcGGGTTTAACCCGGGGTCGCCCATCCATAATGACGGCGCACCGCCCATCGCCCCCTGTTGGCGACGCTTTCCCGCCGCGACCGGAGc                                                  <  2:901065/87‑1 (MQ=255)
                                 cGCCGGGTTTAACACCGGATTGCCCACCCAGAACGCCAGCGCCCCCCCCACCGACACCTGTTGGCGACGCATTCCCGCCGCGACCGGAGc                                                  <  1:358238/90‑1 (MQ=255)
                                       gTTTAACCCCGGATTGCCCATCCAGAAGGCCAGCGCCCCGCCCATCGACCCCTGTTGGCGACGCATGCCCGCCGCGACCGGAGCCGcaca                                            <  2:590151/90‑1 (MQ=255)
                                       gTTAAACCCCGGATTGCCCACCCAGAATGCCAGCGCCCCGCCCACCGCCCCCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGcaca                                            <  1:666635/90‑1 (MQ=255)
                                         ttAACCCCGGATTGCCCATCCAAAATGCCAGCGCACCGCCCACCGACCCCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCAcagc                                          <  2:147257/90‑1 (MQ=255)
                                         ttAACCCCGGATTCCCCACCAAGAATGCCAGCCCCCCGCCCATCGACCCCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCAcagc                                          <  2:737671/90‑1 (MQ=255)
                                         gagacagCGGATTCCCCACCCAGAACGCCCGCGCCCCCCCCATCGACCCCTGTTGGCGACGCATCCCCGCCGCGACCGGAGCCGCAcagc                                          <  2:716173/87‑1 (MQ=255)
                                               ccGGATTGCCCACCCAAAATGCCAGCGCCCCGCCCATCGACCCCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTac                                    <  2:786224/90‑1 (MQ=255)
                                                           tCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGCCCg                        >  1:848924/1‑90 (MQ=255)
                                                            ccAGAATGCCAGCCCCCCCCCCATCGCCCCCTTTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGCCCgg                       <  2:89375/90‑1 (MQ=255)
                                                                aacgCCAGCGCCCCCCCCATCGACCCCTGTGGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGCCCGGCAAc                   <  1:93681/87‑1 (MQ=255)
                                                                         gcACCGCCCATCGACACCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGCCCGGCAACGAAAACAGc          >  1:919338/1‑90 (MQ=255)
                                                                          ccccGCCCATCGACCCCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGCCCGGCAACGAAAACAGCg         <  2:408475/88‑1 (MQ=255)
                                                                          ccccGCCCACCGACCCCTGTGGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGCCCGGCAACGAAAACAGCg         <  2:919338/88‑1 (MQ=255)
                                                                          cACCGCCCATCGCCACCTGTTGGCGACGCTTTCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGCCCGGCAACGAAAACAGCg         <  2:123221/90‑1 (MQ=255)
                                                                             cGCCCATCGACACCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGCCCGGCAACGAAAACAGCGTTc      >  1:1121868/1‑90 (MQ=255)
                                                                                 cATCGACACCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGCCCGGCAACGAAAACAGCGTTCCcag  >  1:942213/1‑90 (MQ=255)
                                                                                        |                                                                                  
CCAGCCGAGGACAAAGCCCATAAACACCAGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTACACATCATGCCCGGCAACGAAAACAGCGTTCCCAG  >  NZ_CP009273/3290988‑3291158

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