Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,634,3410CT43.9% 35.8 / 21.4 41F40F (TTC→TTT) dtpBdipeptide/tripeptide permease DtpB
Reads supporting (aligned to +/- strand):  ref base C (6/17);  new base T (18/0);  total (24/17)
Fisher's exact test for biased strand distribution p-value = 6.89e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCA  >  NZ_CP009273/3634212‑3634468
                                                                                                                                 |                                                                                                                                
gTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTaaa                                                                                                                         <  2:397772/139‑1 (MQ=255)
gTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTaaa                                                                                                                         <  1:50719/139‑1 (MQ=255)
 ttAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTtaa                                                                                                                         >  2:622/1‑138 (MQ=255)
 ttAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTaaa                                                                                                                         >  1:167569/1‑138 (MQ=255)
 ttAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTCGTTtaa                                                                                                                         >  2:602734/1‑138 (MQ=255)
          atGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTCGTTTAAAAGCtttggt                                                                                                               >  2:487090/1‑135 (MQ=255)
           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)
                            aTGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTTGTTTAACAACTTGGGTTCTCGCACGAGGTGGC‑tc                                                                                            >  2:339460/1‑138 (MQ=255)
                            aTGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑tt                                                                                            >  1:498342/1‑139 (MQ=255)
                              ggggATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCGTTAAAAAGCTTGGGTTTTTGCCAGAACAGGC‑tttt                                                                                          >  2:500357/1‑139 (MQ=255)
                              ggggATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAACTTGGGTTTTCGCAAGAGCAGGC‑tttt                                                                                          >  2:325125/1‑139 (MQ=255)
                                    gctgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTAAAAAACTTGGGTTTTCGCAAGAGCAGGC‑TTTTGTTACt                                                                                    >  2:467629/1‑139 (MQ=255)
                                    gctgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGGTTTTTGCAAGAGCACGG‑TTTTGTCACt                                                                                    >  1:308618/1‑139 (MQ=255)
                                    gctgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTCGTTTAACAGCTTGGGTTTTCGCAAGAGCCGGG‑TTTTTTCACt                                                                                    >  2:131348/1‑139 (MQ=255)
                                      tgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTTTTAAAACGCCTCGGTTTTCGGAAGAGGAGGC‑TTTTTTTACttt                                                                                  >  1:337322/1‑139 (MQ=255)
                                         agcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTgg                                                                               <  2:520244/139‑1 (MQ=255)
                                            aaCCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGc                                                                            <  1:591917/139‑1 (MQ=255)
                                               cTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTAAAAAACTTGGGTTTTTGCAAAAGCAGGCGTTTTTTAATTTTTggggttt                                                                          >  2:591917/1‑134 (MQ=255)
                                                tCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTTAACAGCCTTGAGTCTCGCAAGAAGCGGG‑TTTTTTTACTTTTGGTGGtttt                                                                        >  1:374994/1‑139 (MQ=255)
                                                tCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTAAACAGCTTGGATTTTCGGCAGAGCAGGG‑TTTTTTCACTTTTGGTGCtttt                                                                        >  1:96412/1‑139 (MQ=255)
                                                  gCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCGTTTAACAGCTTGGGTTTTCGCAAGAGCAGGG‑TTTTTTTACTTTTTGTGCTTTtgt                                                                      >  2:408514/1‑138 (MQ=255)
                                                  gCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTtg                                                                       <  1:16119/138‑1 (MQ=255)
                                                      aTTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTgcg                                                                  <  2:498342/139‑1 (MQ=255)
                                                       tttttCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTgcgc                                                                 <  2:64085/139‑1 (MQ=255)
                                                                    tttttGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCGTTTAAAAGCTTTGGTTCTCCCAAAAGCAGGC‑TTTTGTCACTTTTTGTGGTTTTTCTGCGCTGGTCTATGGCCt                                                    >  1:201092/1‑139 (MQ=255)
                                                                     ttttGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGc                                                      <  1:408514/136‑1 (MQ=255)
                                                                              gTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCAtt                                          <  1:74985/139‑1 (MQ=255)
                                                                                       gCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTAt                                 <  1:325125/139‑1 (MQ=255)
                                                                                        cGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTTGTTAAACAGCTTGGGTTCTCGCCAGAGCAGGC‑TTTTTTCACTTTTGGTGCTTTTTCTGCGCTGGTCTATGGGCTCCTTTTCATTtgcgggcttg                                >  1:440007/1‑134 (MQ=255)
                                                                                        cGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTCGTTAAAAAGCCTGGGTTCTTGCAAGAGCGGGC‑TTTTTTTACTTTTTGTGCTTTTGCTGCGCTGGTCTATGGGCTCCTTTTCATTggcgggtctt                                >  1:145288/1‑134 (MQ=255)
                                                                                        cGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGGTTTTTTTTTGAACACCCTGGGATTTCCCCGGAGGATGC‑TTTTGTTTCTTTTTGTGATTTTTCGTCGCTGGTGTTCGGCCCCGTTTTCATTgggggggttt                                >  2:304538/1‑134 (MQ=255)
                                                                                        cGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGCTTTCTTCGGTAAACCGCTTGGGTTCTCCCAAGAGAAGGG‑TTTTGACCCCTTTTGTGGGTTTGCTGCGCTGGGCTATGGCCTCATTTTCATTtgcgggcctt                                >  1:33982/1‑134 (MQ=255)
                                                                                                   ctacGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGacc                       <  1:283861/137‑1 (MQ=255)
                                                                                                       ggCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACctgc                 <  1:487090/139‑1 (MQ=255)
                                                                                                                    tACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACg    <  1:88985/139‑1 (MQ=255)
                                                                                                                      cTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCa  <  1:131348/139‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                
GTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCA  >  NZ_CP009273/3634212‑3634468

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: ‑