Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,058,9580GA52.6% ‑2.2 / 15.3 19L57L (CTG→TTG) yihOMFS transporter
Reads supporting (aligned to +/- strand):  ref base G (6/3);  new base A (0/10);  total (6/13)
Fisher's exact test for biased strand distribution p-value = 3.10e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.77e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGA  >  NZ_CP009273/4058885‑4059038
                                                                         |                                                                                
ccAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCACGCCa                                                                  >  2:613463/1‑90 (MQ=255)
  aaTATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCACGCCATa                                                                <  1:613463/90‑1 (MQ=255)
     tttttcCCGCGAGCCGAAAAAAAATCCGGTGAGCATATCGGTAAACGGGGAAAAAACCTTCCCCACCAAAAAGATAATCCCGCCATAGTa                                                             <  1:391375/89‑1 (MQ=255)
        tttCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGc                                                          >  2:102430/1‑90 (MQ=255)
             gcgaggcgagAAAAAACCCGGTGAGCATATCGGAAAACGCGGAAAAAAACTTCCCCACCAGAAAGATAACCCCGCCATAGTAAGCAGGCa                                                     <  1:629162/90‑1 (MQ=255)
             gcgAGTGGAGGAAAAACCCGGTGAGCATATCGGAAACCGCGGTAAAAACCTTCCCCCCAAAAAAGAAAACCCCGCCATAGTAAGCAGGCa                                                     <  1:271527/90‑1 (MQ=255)
             gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCa                                                     >  1:289051/1‑90 (MQ=255)
               gagtcgagtAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATc                                                   >  2:613586/1‑90 (MQ=255)
               gagtcgaaaAAAAATCCGGTGAGAATATCGGAAAACGCGGAAAAAACCTTCCCCCCAAAAAAGATAATCCCGCCATAGTAAGCAGGCATc                                                   <  1:613586/90‑1 (MQ=255)
                                          tggatAACCGCGAAAAAAAACTCCCCCACAAAAAAGAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagg                        <  2:547237/86‑1 (MQ=255)
                                          tCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagg                        >  1:465254/1‑90 (MQ=255)
                                           gggaaaaCGCGGAAAAAACCTCCCCCACAAAAAAGAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagga                       <  1:440126/86‑1 (MQ=255)
                                           cGGTAAACGCGGAAAAAACCTCCGCCACCAAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagga                       <  2:142827/90‑1 (MQ=255)
                                           cGGAAAACGCGAAAAAAACCTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagga                       <  2:168412/90‑1 (MQ=255)
                                             gTAACCGGGGAAAAAACCCCCGCCACCAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga                     <  1:377936/90‑1 (MQ=255)
                                                   gcgGAAAAAACCTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACaga               <  2:532478/90‑1 (MQ=255)
                                                     gaaaaaaacCTCCGCCACAAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGt             <  2:465254/87‑1 (MQ=255)
                                                        aaaaaaCTCCCCCACAAGAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTAcc          <  1:222657/90‑1 (MQ=255)
                                                                tCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCaga  >  1:641409/1‑90 (MQ=255)
                                                                         |                                                                                
CCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGA  >  NZ_CP009273/4058885‑4059038

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