Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,058,9580GA63.2% ‑5.5 / 12.0 19L57L (CTG→TTG) yihOMFS transporter
Reads supporting (aligned to +/- strand):  ref base G (7/0);  new base A (0/12);  total (7/12)
Fisher's exact test for biased strand distribution p-value = 1.98e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.93e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGA  >  NZ_CP009273/4058874‑4059038
                                                                                    |                                                                                
tGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGCAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGa                                                                             >  2:318002/1‑90 (MQ=255)
        ggCCAAATTTTTTCCCGGGAGTCGAAAAAAAACCGGGGGAGCATATCGGAAAACGCGGAAAAAAACTTCGCCACAAAAAAAATAACCCCg                                                                     <  1:32090/90‑1 (MQ=255)
                ttttttCCGCGAGTCGAGAAAAAATCCGGTGAGCATATCGGAAACCGCGGAAAAAAACTCCCCCCCCAAAAAGATAACCCCGCCATAGTa                                                             <  2:46088/89‑1 (MQ=255)
                        gggagtcgaaaAAAAACCGGGGGAGCATATGGGAAAACGCGGAAAAAACCTCCCCCACAAAAAAAATAATCCCGCCATAGTAAGCAGGCa                                                     <  1:61814/88‑1 (MQ=255)
                        gcgGGTCGAGTAAAAACCGGGTGAGCATATCGGTAAACGGGTTAAAAAACTCCCCCCCCAAAAAGAAAATCCCGCCATAGTAAGCAGGCa                                                     <  2:336393/90‑1 (MQ=255)
                        gcgAGTCGAAAAAAAACCGGGGGAGAATATGGGAAAACGCGAAAAAAACCTCCCCCCCAAAAAAAATAACCCCGCCATAGTAAGCAGGCa                                                     <  1:467181/90‑1 (MQ=255)
                                    aaaaCCGGGTGAGCATATCGAAAAACGGGGAAAAAACCTTCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCAt                                         <  1:107345/90‑1 (MQ=255)
                                      aaTCCGGGGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGGAAGCAGGCATCCCTAAGTCATCg                                       >  1:274938/1‑90 (MQ=38)
                                          cGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGtat                                   >  1:465119/1‑90 (MQ=255)
                                                     tCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagg                        >  1:34428/1‑90 (MQ=255)
                                                     tCGGAAACCGCGAAAAAAAACTTCGCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagg                        <  2:418043/90‑1 (MQ=255)
                                                        gaaaaCGCGTAAAAAAACTCCCCCCCCAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga                     <  2:507743/88‑1 (MQ=255)
                                                        gTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga                     >  1:6021/1‑90 (MQ=255)
                                                        gTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga                     >  1:32843/1‑90 (MQ=255)
                                                           accgcgGAAAAAACCTCCCCCACAAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATAc                  <  2:377079/88‑1 (MQ=255)
                                                           accgcgGAAAAAAACTCCCCCCCAAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATAc                  <  2:115171/88‑1 (MQ=255)
                                                                ggaaaaaacCTCCCCCCCCAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGt             <  2:464238/87‑1 (MQ=255)
                                                                          tcccccacaaAAAAAATACTCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCag   <  2:465119/86‑1 (MQ=255)
                                                                           tCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCaga  >  1:339973/1‑90 (MQ=255)
                                                                                    |                                                                                
TGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGA  >  NZ_CP009273/4058874‑4059038

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: