Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,663,7650TA62.5% 3.8 / 15.9 16I29F (ATT→TTT) ynfLLysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (0/6);  new base A (10/0);  total (10/6)
Fisher's exact test for biased strand distribution p-value = 1.25e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.53e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCA  >  NZ_CP009273/1663683‑1663848
                                                                                  |                                                                                   
cTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGGAAAATTTAgg                                                                              >  1:740225/1‑90 (MQ=255)
cTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGGTGGGGAAAATTTAgg                                                                              >  1:302225/1‑90 (MQ=255)
cTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGGTGGGAAAAAtttaag                                                                              >  1:18802/1‑86 (MQ=255)
  gCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGGTGCGGAAAATTTAAGCg                                                                            >  2:761663/1‑90 (MQ=255)
                   gcagTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAAAATTTAGGCGGGGGGCGGCGCGCCCGa                                                           >  2:300098/1‑90 (MQ=255)
                   gcagTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGGTGGGGAAAATTTAGGGGGGGGGGGGGGCGccccg                                                           >  1:34179/1‑88 (MQ=255)
                      gTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAAAATTCAGGGGGGGAGCGGCGCGCCCGcaaa                                                        >  2:70816/1‑89 (MQ=255)
                                            cAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCgcaacagcaac                                  <  1:719798/90‑1 (MQ=255)
                                              gcgcCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGcaacaa                                <  1:423823/90‑1 (MQ=255)
                                                          gctgACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGAt                    <  1:70816/90‑1 (MQ=255)
                                                          gctgACTTAGCGGCGGTTGCGAAAAATTTAGGGGGGCGGCGGCGCGCCCCAAAAGCAGCCCTTCCCCAACAACAACAAAAAAAAGCAGaa                    >  1:32421/1‑89 (MQ=255)
                                                          gctgACTTAGCGGCGGTTGCGAAAAATTCAGGCGGGCAGCGGCGCGCCCCAAAAGCAGCTCTTCCCCAACAGCAAAAAAGTAAAGCAGaa                    >  2:12318/1‑89 (MQ=255)
                                                               cTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGa               <  1:761663/90‑1 (MQ=255)
                                                                ttAGCGGCGGTTGCGAAAAATTCAGGCGGGCAGCGGCGCGCCCCAAAAGCAGCTCTTCCCCAACAACAACAAAGAAAAGCAGATGACGaa              >  2:281638/1‑90 (MQ=255)
                                                                    cggcggTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTc          <  1:550406/90‑1 (MQ=255)
                                                                            gCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCa  <  2:34179/90‑1 (MQ=255)
                                                                                  |                                                                                   
CTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCA  >  NZ_CP009273/1663683‑1663848

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