Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,086,4990CA47.4% 6.4 / 14.0 19V100V (GTG→GTTrhaML‑rhamnose mutarotase
Reads supporting (aligned to +/- strand):  ref base C (4/6);  new base A (9/0);  total (13/6)
Fisher's exact test for biased strand distribution p-value = 1.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCA  >  NZ_CP009273/4086419‑4086582
                                                                                |                                                                                   
atGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGc                                                                            <  1:602695/90‑1 (MQ=255)
   tttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACTTCTTTCAGc                                                                         >  2:413677/1‑90 (MQ=255)
         ttttCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGGAAAAACTTCTTTGAGCTCGCTa                                                                   >  1:4812/1‑90 (MQ=255)
          tttCAAATCGAAATTTAAAATATTGTGCCGGAGGCGTCTCTGGCACATTGGGCAATTACGGCAGGGAAAACACCTCTTGCAGCTCGCCAc                                                                  >  2:290904/1‑90 (MQ=255)
                                    gCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCg                                        <  2:631635/90‑1 (MQ=255)
                                    gCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCg                                        <  1:413625/90‑1 (MQ=255)
                                            cATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACTTTTTTCAGCTCGCTCCTCCCCCGGCTGTTTTTCGGGGTTGCGGGCATaaa                                >  1:631635/1‑89 (MQ=255)
                                               ctctGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATc                             <  2:202344/90‑1 (MQ=255)
                                                         ttGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCatatat                   >  2:255677/1‑90 (MQ=255)
                                                         ttGGGCAATTACGGCAGGGAAAAAACTTCTTTCCGCTCCCCACTCACCCGGCTGTTTTTCGGGGTAGCGGGCATAACATCGGTCatatat                   >  1:323099/1‑90 (MQ=255)
                                                         ttGGGCAATTACGGCAGGGAAAAAACTTCTTTCAGCTCGCTACTCCCCCGGCCGTTTTCCCGGGTAGCGGGGATAAAATCGGGCatatat                   >  1:66728/1‑90 (MQ=255)
                                                         ttGGGCAATTACGGCAGGGAAAAAACCTCTTTTAGCTCGCCACTCACCCGGCTGTTTTTCCGGGTAGCGGGCATAAAAACCGTCatatat                   >  1:410029/1‑90 (MQ=255)
                                                         ttGGGCAATTACGGCAGGGAAAAAAATTTTTTCAGCTCCCTACTCACCCGGCCGTTTTCCCGGGTAGGGGGGATAAAATCGGGCatatat                   >  1:332761/1‑90 (MQ=255)
                                                          tGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATAtt                  <  1:334075/90‑1 (MQ=255)
                                                           gggCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATAttt                 >  2:2505/1‑90 (MQ=255)
                                                           gggCAATTACGGCAGGGAAAAAACTTCTTTCAGCCCGCTACTCACCCGGCTGTTATCCGGGGTAGCGGGCATAAAATCGGGCATATAttt                 >  1:61281/1‑90 (MQ=255)
                                                                 ttACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAAAATCGGTCATATATTTccccca           >  2:145514/1‑90 (MQ=255)
                                                                     ggCAGGTAAAAAACTTCTTTTAGCTCGCTACTCACCCGGCGGTTATTCGGGGTAGCGGGGATAAAAACGGTCATATTTTTTCACCCACGt       >  1:542416/1‑90 (MQ=255)
                                                                          gTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCa  <  1:413677/90‑1 (MQ=255)
                                                                                |                                                                                   
ATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCA  >  NZ_CP009273/4086419‑4086582

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