Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,500,2850AC66.7% 4.5 / 12.2 15V18G (GTA→GGA) ydcOBenE family transporter YdcO
Reads supporting (aligned to +/- strand):  ref base A (4/1);  new base C (0/10);  total (4/11)
Fisher's exact test for biased strand distribution p-value = 3.66e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.27e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCAC  >  NZ_CP009273/1500197‑1500368
                                                                                        |                                                                                   
cATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATAc                                                                                    <  2:521659/90‑1 (MQ=255)
               ttGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGc                                                                     >  1:728754/1‑90 (MQ=255)
                               ccGGCGACAACCGCTGCTTGCCAGATTTTTGCCGCTGACCTGGGGTAGCCAATTAATCCCGCCAAAAACCCCGCCAGTAGCGTGGGTgga                                                     <  2:728583/90‑1 (MQ=255)
                               ccGGCGAAATTCGCTGCTTCCAAATTTTTTCCCGCTAACGGGGGGTACCAAATAAATCCCCCCAAAACCCCCGCCAGTAGCGTGGGTgga                                                     <  1:168631/90‑1 (MQ=255)
                                cGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTggag                                                    >  1:736677/1‑90 (MQ=255)
                                 ggCGAAAATCGCTGCTTGCAAAATTTTTGCCGCTGAAGGGGCGTACCAAATAAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTggagg                                                   <  1:267782/90‑1 (MQ=255)
                                               ttccccAGATTATTCCCGGTGACGGGGGGTAGCAAATAAAACCCCCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGAc                                     <  1:625540/86‑1 (MQ=255)
                                               ctcgCCAAATTATTGCCGCTAAAGGGGCGTAGCAAATAAATCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGAc                                     <  2:903195/87‑1 (MQ=255)
                                                       ttttttGCCGCTGAAGGGGGGTAGCAAATTAATCCCGCCAAAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTc                             <  2:696293/86‑1 (MQ=255)
                                                       ttttttGCCGCTAACCTGGGGTACCAAATAAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTc                             <  2:522025/86‑1 (MQ=255)
                                                                  tGAACGGGCGTAGCCAATAAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGt                  <  2:358658/90‑1 (MQ=255)
                                                                     ccgggCGTAGCAAATTAACCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTtat               <  2:831029/87‑1 (MQ=255)
                                                                             aGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACa       >  2:572972/1‑90 (MQ=255)
                                                                                aaaTAAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGc    <  2:418930/89‑1 (MQ=255)
                                                                                  aTTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCAc  >  2:465087/1‑90 (MQ=255)
                                                                                        |                                                                                   
CATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCAC  >  NZ_CP009273/1500197‑1500368

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