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

AATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGA  >  NZ_CP009273/4058870‑4059038
                                                                                          |                                                                                
aaTTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGa                                                                                 <  1:237481/90‑1 (MQ=255)
 aTTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGaa                                                                                <  1:484716/90‑1 (MQ=255)
    tGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGa                                                                             >  2:132772/1‑90 (MQ=255)
         ttAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATc                                                                        >  2:180050/1‑90 (MQ=255)
         ttAGCCCAAATATTTTTCCGCGAGTCGAAAAAAAACCGGGGGAGAATATCGGAAAACGGGAAAAAAACCTCCCCCACAAAAAAAATaacc                                                                        <  2:430196/90‑3 (MQ=255)
             cccaaaTATTTTCCCGGGAGACGAAAAAAAACCGGGTGAGCATATGGGAAAACGGGGAAAAAAACTTCGCCACAAAAAAAATAACCCCGc                                                                    <  2:597580/86‑1 (MQ=255)
                       tttCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGc                                                          >  1:248557/1‑90 (MQ=255)
                       tttCCGCGAGTCGAAAAAAAACCCGGGGAGCATATGGGAAAACGCGATAAAAACCTTCCCCCCCAAAAAGATAACCCCGCCATAGTAAGc                                                          <  1:520829/90‑1 (MQ=255)
                       ttcccGGGAGTCGAAAAAAAACCGGGGGAGCATATGGGAAAACGCGAAAAAAACCTCCCCCACAAAAAAGAAAACCCCGCCATAGTAAGc                                                          <  2:501262/87‑1 (MQ=255)
                            gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCa                                                     >  2:403482/1‑90 (MQ=255)
                            gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCCCCACCAAAAAAATAATCCCGCCATAGTAAGCAGGCa                                                     <  1:155586/90‑1 (MQ=39)
                                         aaaCCCGGTGAGCATATGGGTAACCGGGGAAAAAAACTTCCCCACCAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATc                                        <  1:295363/90‑1 (MQ=255)
                                               ggTGAGCATATGGGAAAACGCGGAAAAAAACTCCCCCCCCAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGtata                                  <  1:180050/90‑1 (MQ=255)
                                                     cATATCGGAAAACGCGGAAAAAAACTTCCCCCCCAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAAttt                            <  1:135516/90‑1 (MQ=255)
                                                     aatatCGGAAAACGCGAAAAAAAACTTCGCCCCAAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAAttt                            <  1:379904/89‑1 (MQ=255)
                                                         ccGGAAAGCGCGGAAAAAACCTCCGCCACAAGAAAAATACCCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagg                        <  1:69487/89‑1 (MQ=255)
                                                           cggaaaCCCGGGAAAAAACCTCCCCCACAAAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagga                       <  2:340549/87‑1 (MQ=255)
                                                            gaaaaCGGGGAAAAAACCTCCGCCCCAAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga                     <  2:507494/88‑1 (MQ=255)
                                                                 cgcgGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACag                >  2:133116/1‑90 (MQ=255)
                                                                   ggggaaaaaaCTTCCCCCCCCAAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACaga               <  1:526963/86‑1 (MQ=255)
                                                                       aaaaacCTCCCCCCCCAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTAcc          <  1:139192/90‑1 (MQ=255)
                                                                             ccTTCCCCACAAAAAAGAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGc     <  1:308613/89‑1 (MQ=255)
                                                                                tCGCCACCAGAGAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCaga  >  1:365667/1‑90 (MQ=255)
                                                                                          |                                                                                
AATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGA  >  NZ_CP009273/4058870‑4059038

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