Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,058,9580GA61.9% ‑6.5 / 15.8 21L57L (CTG→TTG) yihOMFS transporter
Reads supporting (aligned to +/- strand):  ref base G (8/0);  new base A (0/13);  total (8/13)
Fisher's exact test for biased strand distribution p-value = 4.91e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.99e-05
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAG  >  NZ_CP009273/4058879‑4059037
                                                                               |                                                                               
ttAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATc                                                                       >  1:534544/1‑90 (MQ=255)
            ttttcCCGGGAGTCGAGAAAAAACCCGGTGAGCATATCGGTAACCGCGGAAAAAAACTTCCCCACCAAAAAGATAATCCCGCCATAGTaa                                                           <  2:62663/90‑1 (MQ=255)
             ttttCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAg                                                          >  2:366432/1‑90 (MQ=255)
                   gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCa                                                    >  1:528058/1‑90 (MQ=255)
                   gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCa                                                    >  1:553013/1‑90 (MQ=255)
                   gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCa                                                    >  2:132483/1‑90 (MQ=255)
                           agaaaaaaCCCGGTGAGCATATCGGAAAACGCGGAAAAAAACTCCCCCACCAAAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAAt                                            <  2:5277/87‑1 (MQ=255)
                                aaaCCCGGGGAGCATATCGGAAACCGCGAAAAAAACCTCCGCCACAAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATc                                       <  2:553505/90‑1 (MQ=255)
                                      ggTGAGAATATGGGTAAACGCGAAAAAAAACTTCCCCACAAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGtata                                 <  2:528058/90‑1 (MQ=255)
                                           gCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAAtt                            >  2:551284/1‑90 (MQ=255)
                                              tatCGAAAACCGCGAAAAAAACCTCCCCCACCAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTga                         <  1:31508/90‑1 (MQ=255)
                                                   gaaaaCGCGGAAAAAAACTCCCCCACAAAAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga                    <  2:89233/88‑1 (MQ=255)
                                                   gTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCa                                                    >  1:85115/1‑58 (MQ=255)
                                                   gTAAACGCGGTAAAAAACTTCGCCACCAAAAAGATAATCCCGCCATAGTAAGCAGGCa                                                    <  2:85115/58‑1 (MQ=255)
                                                    tAAACGCGAAAAAAATCTCCCCCACAAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGAt                   <  1:289263/90‑1 (MQ=255)
                                                      aaCGCGAAAAAAACCTCCGCCACAAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATAc                 <  1:574466/90‑1 (MQ=255)
                                                           ggTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGt            >  2:273978/1‑90 (MQ=255)
                                                             aaaaaaaCTCCCCCAAAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTAc          <  1:405757/89‑1 (MQ=255)
                                                              aaaaacCTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTAcc         <  1:519211/90‑1 (MQ=255)
                                                                     tcccccacaaAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCag  <  2:258161/86‑1 (MQ=255)
                                                                     ccccccacaaAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCag  <  1:457918/86‑1 (MQ=255)
                                                                               |                                                                               
TTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAG  >  NZ_CP009273/4058879‑4059037

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