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

TTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTA  >  NZ_CP009273/4058879‑4059044
                                                                               |                                                                                      
ttAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATc                                                                              >  1:13449/1‑90 (MQ=255)
ttAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATc                                                                              >  1:526762/1‑90 (MQ=255)
       cAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCAt                                                                       >  1:161626/1‑90 (MQ=255)
       cAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCAt                                                                       >  1:626940/1‑90 (MQ=255)
              tttCCGCGAGTCGAGAAAAAATCCGGTGAGTATATCGGTAAACGCGGAAAAAACCTTCCCCCCCAAAAAGATAACCCCGCCATAGTAAGc                                                                <  2:621846/90‑1 (MQ=255)
              ttcccGGGAGTAGAAAAAAAACCGGGGGAACATATGGGAAAACGGGAAAAAAACCTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGc                                                                <  1:439293/87‑1 (MQ=255)
              ttcccGCGAGTCGAGAAAAAACCCGGTGAGCATATCGGAAAACGCTGAAAAAAACTTCCCCCCCAAAAAGATAATCCCGCCATAGTAAGc                                                                <  2:133003/87‑1 (MQ=255)
                   gggagtcgagaAAAAACCCGGTGAGAATATCGGAAACCGCGGAAAAAACCTCCCCCACCAAAAAGATAATCCCGCCATAGTAAGCAGGCa                                                           <  2:311555/88‑1 (MQ=255)
                   gcgAGTCGAAAAAAAACCGGGTGAAAATATGGGAAAACGCGGAAAAAACCTTCCCCACCAAAAAAATAACCCCGCCATAGTAAGCAGGCa                                                           <  1:425236/90‑1 (MQ=255)
                               aaaaCCCGGTGACCATACCGGAAAACGGGGAAAAAACCTTCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCAt                                               <  1:100832/90‑1 (MQ=255)
                                aaaTCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACAAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATc                                              <  2:38589/90‑1 (MQ=255)
                                      ggTGAGCATATGGGTAAACGGGGAAAAAAACTCCGCCCCAAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGtata                                        <  2:281244/90‑1 (MQ=255)
                                            cATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAATGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAAttt                                  >  2:394088/1‑90 (MQ=255)
                                                gcGGAAAACGGGTAAAAAAACTTCCCAACAAAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagg                              <  1:280857/89‑1 (MQ=255)
                                                                     tcccccacaaAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCag         <  1:33844/86‑1 (MQ=255)
                                                                     tcccccacaaAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCag         <  1:183615/86‑1 (MQ=255)
                                                                      tCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCaga        >  1:294126/1‑90 (MQ=255)
                                                                      tCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCaga        >  1:571400/1‑90 (MQ=255)
                                                                            caaGAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTa  <  1:445851/88‑1 (MQ=255)
                                                                               |                                                                                      
TTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTA  >  NZ_CP009273/4058879‑4059044

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