Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,634,3590AG42.1% 52.0 / 13.4 38G46G (GGA→GGGdtpBdipeptide/tripeptide permease DtpB
Reads supporting (aligned to +/- strand):  ref base A (5/17);  new base G (16/0);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 1.04e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGC  >  NZ_CP009273/3634223‑3634496
                                                                                                                                        |                                                                                                                                          
tGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTTAACAGCTTGGGtt                                                                                                                                          >  2:342394/1‑139 (MQ=255)
tGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGAtt                                                                                                                                          <  1:622/139‑1 (MQ=255)
         caacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAaga                                                                                                                                 <  1:419369/139‑1 (MQ=255)
         caacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTCGGTGAAAAAGTTGGGTTTTCGCaa                                                                                                                                   >  1:180474/1‑137 (MQ=255)
              cccATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCTTTTAACAGCTTGGATTTTCGCAAGGGCAgg                                                                                                                            >  2:411982/1‑139 (MQ=255)
              cccATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAgg                                                                                                                            <  1:194451/139‑1 (MQ=255)
                 aTGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTTGTTTAACAACTTGGGTTCTCGCACGAGGTGGCtc                                                                                                                        >  2:339460/1‑138 (MQ=255)
                 aTGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCtt                                                                                                                        >  1:498342/1‑139 (MQ=255)
                   ggggATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCGTTAAAAAGCTTGGGTTTTTGCCAGAACAGGCtttt                                                                                                                      >  2:500357/1‑139 (MQ=255)
                   ggggATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAACTTGGGTTTTCGCAAGAGCAGGCtttt                                                                                                                      >  2:325125/1‑139 (MQ=255)
                         gctgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTAAAAAACTTGGGTTTTCGCAAGAGCAGGCTTTTGTTACt                                                                                                                >  2:467629/1‑139 (MQ=255)
                         gctgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGGTTTTTGCAAGAGCACGGTTTTGTCACt                                                                                                                >  1:308618/1‑139 (MQ=255)
                         gctgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTCGTTTAACAGCTTGGGTTTTCGCAAGAGCCGGGTTTTTTCACt                                                                                                                >  2:131348/1‑139 (MQ=255)
                           tgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTTTTAAAACGCCTCGGTTTTCGGAAGAGGAGGCTTTTTTTACttt                                                                                                              >  1:337322/1‑139 (MQ=255)
                              agcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTgg                                                                                                           <  2:520244/139‑1 (MQ=255)
                                 aaCCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGc                                                                                                        <  1:591917/139‑1 (MQ=255)
                                    cTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTAAAAAACTTGGGTTTTTGCAAAAGCAGGCGTTTTTTAATTTTTggggttt                                                                                                      >  2:591917/1‑134 (MQ=255)
                                     tCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTTAACAGCCTTGAGTCTCGCAAGAAGCGGGTTTTTTTACTTTTGGTGGtttt                                                                                                    >  1:374994/1‑139 (MQ=255)
                                     tCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTAAACAGCTTGGATTTTCGGCAGAGCAGGGTTTTTTCACTTTTGGTGCtttt                                                                                                    >  1:96412/1‑139 (MQ=255)
                                       gCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCGTTTAACAGCTTGGGTTTTCGCAAGAGCAGGGTTTTTTTACTTTTTGTGCTTTtgt                                                                                                  >  2:408514/1‑138 (MQ=255)
                                       gCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTtg                                                                                                   <  1:16119/138‑1 (MQ=255)
                                           aTTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTgcg                                                                                              <  2:498342/139‑1 (MQ=255)
                                            tttttCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTgcgc                                                                                             <  2:64085/139‑1 (MQ=255)
                                                         tttttGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCGTTTAAAAGCTTTGGTTCTCCCAAAAGCAGGCTTTTGTCACTTTTTGTGGTTTTTCTGCGCTGGTCTATGGCCt                                                                                >  1:201092/1‑139 (MQ=255)
                                                          ttttGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGc                                                                                  <  1:408514/136‑1 (MQ=255)
                                                                   gTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCAtt                                                                      <  1:74985/139‑1 (MQ=255)
                                                                            gCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTAt                                                             <  1:325125/139‑1 (MQ=255)
                                                                             cGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTTGTTAAACAGCTTGGGTTCTCGCCAGAGCAGGCTTTTTTCACTTTTGGTGCTTTTTCTGCGCTGGTCTATGGGCTCCTTTTCATTtgcgggcttg                                                            >  1:440007/1‑134 (MQ=255)
                                                                             cGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTCGTTAAAAAGCCTGGGTTCTTGCAAGAGCGGGCTTTTTTTACTTTTTGTGCTTTTGCTGCGCTGGTCTATGGGCTCCTTTTCATTggcgggtctt                                                            >  1:145288/1‑134 (MQ=255)
                                                                             cGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGGTTTTTTTTTGAACACCCTGGGATTTCCCCGGAGGATGCTTTTGTTTCTTTTTGTGATTTTTCGTCGCTGGTGTTCGGCCCCGTTTTCATTgggggggttt                                                            >  2:304538/1‑134 (MQ=255)
                                                                             cGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGCTTTCTTCGGTAAACCGCTTGGGTTCTCCCAAGAGAAGGGTTTTGACCCCTTTTGTGGGTTTGCTGCGCTGGGCTATGGCCTCATTTTCATTtgcgggcctt                                                            >  1:33982/1‑134 (MQ=255)
                                                                                        ctacGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGacc                                                   <  1:283861/137‑1 (MQ=255)
                                                                                            ggCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACctgc                                             <  1:487090/139‑1 (MQ=255)
                                                                                                         tACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACg                                <  1:88985/139‑1 (MQ=255)
                                                                                                           cTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCa                              <  1:131348/139‑1 (MQ=255)
                                                                                                                          gTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAg               <  1:467629/139‑1 (MQ=255)
                                                                                                                                  gCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTtgtg       <  1:258773/139‑1 (MQ=255)
                                                                                                                                   cTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGc      >  2:97126/1‑139 (MQ=255)
                                                                                                                                        aTTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGc  >  2:38394/1‑138 (MQ=255)
                                                                                                                                        |                                                                                                                                          
TGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGC  >  NZ_CP009273/3634223‑3634496

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: