Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,500,2850AC56.0% 24.2 / 32.0 25V18G (GTA→GGA) ydcOBenE family transporter YdcO
Reads supporting (aligned to +/- strand):  ref base A (9/2);  new base C (0/14);  total (9/16)
Fisher's exact test for biased strand distribution p-value = 2.69e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.44e-05
Rejected as consensus: Frequency below/above cutoff threshold.

CATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACC  >  NZ_CP009273/1500197‑1500369
                                                                                        |                                                                                    
cATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATAc                                                                                     <  1:1175282/90‑1 (MQ=255)
  tCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCg                                                                                   <  1:1312191/90‑1 (MQ=255)
         agagTTTTGGGCTGTGGGGGCCCCGGCGAAAATCGCTGCTGGCAAGATTATTGCCGCTGAATTGGCGAACCAAATAAATCCCGCCAGaac                                                                            <  1:949632/90‑2 (MQ=255)
                        ggtggCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCgg                                                             >  1:1418432/1‑89 (MQ=255)
                           ggCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCgtggg                                                          >  2:3599/1‑90 (MQ=255)
                           ggCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCgtggg                                                          >  1:1303338/1‑90 (MQ=255)
                                cGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTggag                                                     >  1:1095640/1‑90 (MQ=255)
                                cGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTggag                                                     >  1:837009/1‑90 (MQ=255)
                                cGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTggag                                                     >  1:181295/1‑90 (MQ=255)
                                 ggCAACAACCGCTGCTGGCCATTTTTTGGCCGCGGACGGGGGGTAGAAAATAAACCCCGCAAAAACCCCCGCCAGTAGCGTGGGTggagg                                                    <  1:703359/90‑1 (MQ=255)
                                               ctTGCCAGATTATTGCCGCTGAACGGGGGTACCCAATTAATCCCCCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGAc                                      <  2:1384280/90‑1 (MQ=255)
                                               ctTCCCAGATTATTGCCGCTGAACGGGCGTACCAAATAAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGAc                                      <  2:610592/90‑1 (MQ=255)
                                                 tcccAGATTATCGCCGCTGACCGGGCGTCACAAATAAATCCCCCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGc                                    <  2:1392151/88‑1 (MQ=255)
                                                 tcccAAATTATTGCCGCTGAACGGGGGTAGCCAATTAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGc                                    <  2:537899/88‑1 (MQ=255)
                                                  gCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCa                                   >  1:1052227/1‑90 (MQ=255)
                                                   ccAGATTATTGCCGCTGAACGGGCGTAGCCAATTAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCAt                                  <  2:486442/90‑1 (MQ=255)
                                                           ttGCCGCAGAACGGGGGAAGCAAATAAACCCCCCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCt                          <  2:1042858/90‑1 (MQ=255)
                                                             cccGCTGACCGGGCGTACCCAATTAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCg                        <  2:8494/89‑1 (MQ=255)
                                                                gCTGAACGGGAGTACCAAATTAATCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGc                     <  2:1095640/90‑1 (MQ=255)
                                                                      cgggCGTACCCAATAAATCCCCCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTtata               <  2:1418432/88‑1 (MQ=255)
                                                                      cgggCGTACCAAATAAACCCCCCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTtata               <  2:251902/88‑1 (MQ=255)
                                                                        ggCGTAGCAAATTAATCCCCCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGc             <  1:3599/90‑1 (MQ=255)
                                                                             aGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACa        >  1:356416/1‑90 (MQ=255)
                                                                                  aTTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCAc   >  1:586056/1‑90 (MQ=255)
                                                                                   ttAATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAGCAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCAcc  <  2:512697/90‑1 (MQ=255)
                                                                                        |                                                                                    
CATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACC  >  NZ_CP009273/1500197‑1500369

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