Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,058,9580GA53.8% 10.9 / 30.3 26L57L (CTG→TTG) yihOMFS transporter
Reads supporting (aligned to +/- strand):  ref base G (6/6);  new base A (0/14);  total (6/20)
Fisher's exact test for biased strand distribution p-value = 4.01e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.92e-05
Rejected as consensus: Frequency below/above cutoff threshold.

AAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAG  >  NZ_CP009273/4058869‑4059039
                                                                                          |                                                                                 
aaaTTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAg                                                                                   <  1:1677/90‑1 (MQ=255)
 aaTTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGa                                                                                  <  2:466047/90‑1 (MQ=255)
 aaTTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGa                                                                                  <  1:97982/90‑1 (MQ=255)
     tGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGa                                                                              <  2:549765/90‑1 (MQ=255)
       ccTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATa                                                                            <  2:902415/90‑1 (MQ=255)
              gCCAAATTTTTTTCGGCGAGCCGAGAAAAAATCCGGTGAGAATATCGGAAAACGCGGAAAAAACCTCCCCCACAAAAAAGATACCCCCGc                                                                     <  2:890406/90‑1 (MQ=255)
                ccAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCa                                                                   >  1:564428/1‑90 (MQ=255)
                       tttcCCGGGAGCCGAGAAAAAACCCGGTGAGAATATCGGAAAACGCGGAAAAAACCTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAg                                                            <  1:91523/90‑1 (MQ=255)
                        tttCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGaa                                                                                 <  2:868079/68‑1 (MQ=255)
                        tttCCGCGAGTCGAGTAAAAATCCGGCGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGaa                                                                                 >  1:868079/1‑68 (MQ=255)
                        ttcccGCGAGCCGAAAAAAAATCCGGTGAAAATATCGAAAAACGCGGAAAAAACCTCCCCCAAAAAAAAAATAACCCCGCCATAGTAAGc                                                           <  2:575813/87‑1 (MQ=255)
                            cgcgAGTCGAGAAAAAACCCGGTGAGCATATGGGAAAACGCGGAAAAAAACTTCCCCACCAAAAAGATAACCCCGCCATAGTAAGCAGGc                                                       <  1:837101/90‑1 (MQ=255)
                             gcgGGTCGAAAAAAAATCGGGTGAGAATAGGGGAAAACGGGAAAAAAACCTCCCCCACAAAAAAAATAATCCCGCCATAGTAAGCAGGCa                                                      <  1:443219/90‑1 (MQ=255)
                             gcgAGTGGAGAAAAAATCCGGTGAGCATATCGGAAAACGCGGAAAAAAACTCCCCCACAAAAAAGATAACCCCGCCATAGTAAGCAGGCa                                                      <  1:54504/90‑1 (MQ=255)
                             gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCa                                                      >  2:151197/1‑90 (MQ=255)
                             gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCa                                                      >  2:837101/1‑90 (MQ=255)
                             gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCa                                                      >  2:659376/1‑90 (MQ=255)
                             gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGAAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCa                                                      <  1:552918/90‑1 (MQ=255)
                                         aaacTCGGGTGAGCATATGGGAAAACGCGGTAAAAACCTTCCCCCCCAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCAt                                          <  1:293905/90‑1 (MQ=255)
                                         aaacTCCGGTGAAAATACGGAAAACCGGGAAAAAAACCTCCCCCCCAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCAt                                          <  2:188773/90‑1 (MQ=255)
                                         aaaaTCGGGTGACCATACCGTAAAACGCGGAAAAAAACTCCCCCCCAAAAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCAt                                          <  1:155920/90‑1 (MQ=255)
                                                ggTGAGCATATAGAAAAACGCGAAAAAAACCTTCCCCACAAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGtata                                   <  2:879104/90‑1 (MQ=255)
                                                        tatCGAAAAACGGGGAAAAAACCTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTga                           <  2:465674/90‑1 (MQ=255)
                                                           cgtaaaaCGCGAAAAAAACCCCCGCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagga                        <  2:326637/86‑1 (MQ=255)
                                                                     gaaaaaaacCTCCGCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGt              <  2:537137/87‑1 (MQ=255)
                                                                              cTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCa     >  1:541441/1‑90 (MQ=255)
                                                                                 ccccacaaAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCagag  <  2:541441/88‑1 (MQ=255)
                                                                                          |                                                                                 
AAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAG  >  NZ_CP009273/4058869‑4059039

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: