Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273117,2520TG52.9% 9.5 / 15.4 17T263P (ACC→CCC) aroParomatic amino acid transporter AroP
Reads supporting (aligned to +/- strand):  ref base T (3/5);  new base G (9/0);  total (12/5)
Fisher's exact test for biased strand distribution p-value = 9.05e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.66e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CGATGTTCAGCGCATTCGCCACAAAGGTATCGCCTAACTCGTGGAAGATCAGCACAAACGGACTGGTATCGGCGGTAACGCGGGTCCACGGCATCAGTGAGAGCAGAACGGCTAACGAACCAATATAGAAAATCAGGATGCGGTAGATAACCTGGTTAGTTGCTTT  >  NZ_CP009273/117168‑117333
                                                                                    |                                                                                 
cGATGTTCAGCGCATTCGCCACAAAGGTATCGCCTAACTCGTGGAAGATCAGCACAAACGGACTGGTATCGGCGGTAACGCGGGGCCACg                                                                              >  2:308448/1‑90 (MQ=255)
            cATTCGCCACAAAGGTATCGCCTAACTCGTGGAAGATCAGCACAAACGGACTGGTATCGGCGGTAACGCGGGGGCCCCGCATTAGGgaga                                                                  >  2:213713/1‑90 (MQ=255)
              ttCGCCACAAAGGTATCGCCTAACTCGTGGAAGATCAGCACAAACGGACTGGTATCGGCGGGAAAGCGGGGCCACGGCATCAGGGagagg                                                                >  1:289886/1‑89 (MQ=255)
                          gTATCGCCTAACTCGTGGAAGATCAGCACAAACGGACTGGTATCGGCGGTAACGCGGGGCCACGGCATCAGTGAGAGCAGAACGGCTAAc                                                    >  2:382304/1‑90 (MQ=255)
                            aTCGCCTAACTCGTGGAAGATCAGCACAAACGGACTGGTATCGGCGGGAAACCGGGGCCCCGGCATCAGTGGGAGGAGAACGGGTAAAGa                                                  >  1:377745/1‑90 (MQ=255)
                                        gTGGAAGATCAGCACAAACGGACTGGTATCGGCGGTAACGCGGGGCCCCGGCATCAGTGAGAGCAGAACGGCTAACGAACCAATATAGaa                                      >  2:264595/1‑90 (MQ=255)
                                        gTGGAAGATCAGCACAAACGGACTGGTATCGGCGGGAACGCGGGGGCACGGGGTCAGGGAGGGGAGAACGGCGTAAGCAACAATATAgga                                      >  2:380802/1‑88 (MQ=255)
                                                    cacaAACGGACTGGTATCGGCGGTAACGCGGGTCCACGGCATCAGTGAGGGCAGAACGGCTAACGAACCAATATAGAAAATCAGGATGCg                          <  2:351635/90‑1 (MQ=255)
                                                    cacaAACGGACTGGTATCGGCGGTAACGCGGGTCCACGGCATCAGTGAGAGCAGAACGGCTAACGAACCAATATAGAAAATCAGGATGCg                          <  1:308448/90‑1 (MQ=255)
                                                             aCTGGTATCGGCGGTAACCCGGGTCCACGGCATCAGTGAGAGCAGAACGGCTAACGAACCAATATAGAAAAACAGGATGCGGTAGATAAc                 >  2:225123/1‑90 (MQ=255)
                                                                ggTATCGGCGGTAAAGCGGGGCCACGGCATCAGTGAGAGCAGAAAGGCTAACGAAACAAAATAGAAAAACAGGAAGCGGGAGATAAACTg              >  1:435828/1‑90 (MQ=255)
                                                                      ggcggTAACGCGGGTCCACGGCATCAGTGAGAGCAGAACGGCTAACGAACCAATATAGAAAATCAGGATGCGGTAGATAACCTGGTTAGt        <  2:71948/90‑1 (MQ=255)
                                                                      ggcggTAACGCGGGTCCACGGCATCAGGGAGAGCAGAACGGCTAACGAACCAATATAGAAAATCAGGATGCGGGAGATAACCTGGGTAGt        >  2:318942/1‑90 (MQ=255)
                                                                          gTAACGCGGGTCCACGGCATCAGTGAGAGCAGAACGGCTAACGAACCAATATAGAAAATCAGGATGCGGTAGATAACCTGGTTAGTTGCt    <  2:123194/90‑1 (MQ=255)
                                                                          gTAACGCGGGTCCACGGCATCAGTGAGAGCAGAACGGCTAACGAACCAATATAGAAAATCAGGATGCGGTAGATAACCTGGTTAGTTGCt    <  2:435484/90‑1 (MQ=255)
                                                                          gTAACGCGGGTCCACGGCATCAGGGAGAGCAGAACGGCTAACGAACCAATATAAAAAAACAAGAGGGGGTAGATAAACTGGTGATTTGCt    >  2:229393/1‑90 (MQ=255)
                                                                            aaCGCGGGGCCACGGCATCAGTGAGAGCAGAAAGGCTAACGAACCAATATAGAAAAACAGGATGCGGGAGATAACCTGGGTAGTTGCttt  >  1:159163/1‑90 (MQ=255)
                                                                                    |                                                                                 
CGATGTTCAGCGCATTCGCCACAAAGGTATCGCCTAACTCGTGGAAGATCAGCACAAACGGACTGGTATCGGCGGTAACGCGGGTCCACGGCATCAGTGAGAGCAGAACGGCTAACGAACCAATATAGAAAATCAGGATGCGGTAGATAACCTGGTTAGTTGCTTT  >  NZ_CP009273/117168‑117333

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