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

AATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGA  >  NZ_CP009273/4058870‑4059038
                                                                                        |                                                                                
aaTTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGa                                                                                 <  2:124164/90‑1 (MQ=255)
aaTTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGa                                                                                 <  1:932064/90‑1 (MQ=255)
                cAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCAt                                                                 >  1:746361/1‑90 (MQ=255)
                            gcgAGTCGAGAAAAAACCGGGTGAGCATATCGGAAACCGGGGAAAAAACCTCCCCCCCAAAAAAAAAAACCCCGCCATAGTAAGCAGGCa                                                     <  2:152594/90‑1 (MQ=255)
                            gagagtcgagaAAAAATCCGGTGAGAATATGGGAAACCGCGGAAAAAACCTCCGCCACCAAAAAGAAAACCCCGCCATAGTAAGCAGGCa                                                     <  1:769376/88‑1 (MQ=255)
                                    aaaaaaaaCCGGGGGAGCATATCGGAAAAGGGGAAAAAAACCTCCCCCACAAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAAt                                             <  2:62322/87‑1 (MQ=255)
                                         aaaTCGGGTAAACATATCGAAAAACGCGAAAAAAACCTCCCCCGCAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATc                                        <  2:528011/90‑1 (MQ=255)
                                         aaaTCCGGTGAGCATCTCGGTAACCGCGATAAAAATCTTCGCCACAAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATc                                        <  1:880251/90‑1 (MQ=255)
                                         aaaCCCGGTGAGCATATGGGAAACCGCGGAAAAAAACTTCCCCCCCAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATc                                        <  1:97582/90‑1 (MQ=255)
                                          aaTCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCg                                       >  2:429250/1‑90 (MQ=255)
                                           accgggTGAGCATATCGGAAAACGGGGAAAAAAACTCCCCCCCAAAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCgg                                      <  2:285639/86‑1 (MQ=255)
                                               ggTGAGCATACCGGTAAACGCGAAAAAAACCTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGtata                                  <  2:108895/90‑1 (MQ=255)
                                                       tatCGGAAACCGCGGAAAAAACCTTCCCCACCAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTga                          <  1:49871/90‑1 (MQ=255)
                                                          cGGAAAACGCGGAAAAAACCTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagga                       <  2:409250/90‑1 (MQ=255)
                                                                       aaaaaaCTTCCCCACAAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTAcc          <  1:956015/90‑1 (MQ=255)
                                                                               tCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCaga  >  2:927123/1‑90 (MQ=255)
                                                                               tCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCaga  >  2:99269/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: