Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,634,3410CT46.2% 42.5 / 32.0 52F40F (TTC→TTT) dtpBdipeptide/tripeptide permease DtpB
Reads supporting (aligned to +/- strand):  ref base C (5/23);  new base T (24/0);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 3.00e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCT‑ATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGT  >  NZ_CP009273/3634203‑3634475
                                                                                                                                          |                                                                                                                                       
tATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTttc                                                                                                                                         >  2:64243/1‑139 (MQ=255)
     gAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTttcttttg                                                                                                                                    >  1:878088/1‑136 (MQ=255)
          ttAATAAATATGAATACAACACCACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTAATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTGTTCTTTGTTtaa                                                                                                                                >  2:529728/1‑138 (MQ=255)
          ttAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTtaa                                                                                                                                >  1:81037/1‑138 (MQ=255)
          ttAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTtaa                                                                                                                                >  1:630276/1‑138 (MQ=255)
          ttAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCTTTtaa                                                                                                                                >  2:627190/1‑138 (MQ=255)
          ttAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCGTTaaa                                                                                                                                >  2:411240/1‑138 (MQ=255)
          ttAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTaaa                                                                                                                                >  2:556762/1‑138 (MQ=255)
          ttAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTCTTCGTTaaa                                                                                                                                >  2:650626/1‑138 (MQ=255)
          ttAATAAATATGAATACAACAACACCCATGGAGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCCTTaaa                                                                                                                                >  1:122327/1‑138 (MQ=255)
                   atGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCGTTTAACAGCtttggt                                                                                                                      >  2:407113/1‑135 (MQ=255)
                   atGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTTGTTAAACAGCTTgggg                                                                                                                      >  2:872405/1‑137 (MQ=255)
                    tGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGAtt                                                                                                                     <  1:556762/139‑1 (MQ=255)
                       aTACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATtctc                                                                                                                  <  2:514398/139‑1 (MQ=255)
                          caacaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGc                                                                                                                <  2:5158/138‑1 (MQ=255)
                             caacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTTTTTAACAGCTTTGATTTTCGcca                                                                                                              >  2:884091/1‑135 (MQ=255)
                             caacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCCTTTAACCACCTTGATtctccccac                                                                                                             >  2:201292/1‑133 (MQ=255)
                                cacCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCGTTTAAAAGCTTGGGGTCTTCCAagag                                                                                                           >  2:112419/1‑137 (MQ=255)
                                     aTGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTCTTTAAAAGCTTGGGTTTTCGCAAGAGCAGGCtt                                                                                                    >  2:868962/1‑139 (MQ=255)
                                     aTGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCGTTTAACAGCTTGGGTTTTCGCGAAAGCAGGGtt                                                                                                    >  1:514268/1‑139 (MQ=255)
                                     aTGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTCTTTTAACAGCTTTGACTCTCTTCAGAGCCGGcct                                                                                                    >  1:624495/1‑137 (MQ=255)
                                       ggggATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTTTTAAACAGCTTTGGTTCTCGCCAAAGCAGGCcttt                                                                                                  >  1:877435/1‑139 (MQ=255)
                                       ggggATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTTAACAGCTTTGGTTTTTGCCAGAGCAGGCtttt                                                                                                  >  2:557569/1‑139 (MQ=255)
                                           aTGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTTAACAGCTTGGGTTTTCGCAAGAGGAGGCTTTTGTCa                                                                                              >  1:828066/1‑139 (MQ=255)
                                             gctgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTTGTTTAACAACTTGGTTTCTCGCAAGAGCAGGCTTTTTTCACt                                                                                            >  2:242866/1‑139 (MQ=255)
                                             gctgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTCGTTAAAAAGCTTGGGTTTTCGCAAGAGCAGGGTTTTTTTACt                                                                                            >  1:905440/1‑139 (MQ=255)
                                               tgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGCTTTTTTCCTTGAAAAGCTTTGTTTCTCGCAGGAGCAGGCTTTTGTTTCttt                                                                                          >  2:65347/1‑139 (MQ=255)
                                                   gcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTTTTTTAAAAACTTGGGGTCTCGCCAGAGCAGGGGTTTTTTACTTTTTGt                                                                                      >  2:207362/1‑139 (MQ=255)
                                                     aaCCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGc                                                                                    <  2:345739/139‑1 (MQ=255)
                                                     aaCCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGc                                                                                    <  2:899390/139‑1 (MQ=255)
                                                     aaCCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGc                                                                                    <  1:733991/139‑1 (MQ=255)
                                                           gCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTtg                                                                               <  2:502108/138‑1 (MQ=255)
                                                                   ttCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTg                                                                       <  1:877387/138‑1 (MQ=255)
                                                                     cATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTTAACAGCTTGGGTTTTCGCAAGAGCAGGCTTTTGTCACTTTTTGGGGTTTTTCTGCGCTGGTc                                                                    >  1:285988/1‑139 (MQ=255)
                                                                           cTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTAAACAACTTGGGTTTTCGCAAGAGCAGGCTTTTTTCACTTTTTGTGGTTTTTCTGCGCTGGTCT‑ATGGc                                                             >  1:65695/1‑139 (MQ=255)
                                                                              ttttGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCT‑ATGGc                                                             <  1:765059/136‑1 (MQ=255)
                                                                                                 cGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTTAACAGCTTGGGTTTTCGCCAGAGCAGGGTTTTGTTACTTTTTGGGCTTTTTTTTCGCTGGGCT‑ATGGCCTCATTTCCAATGGGGGCtgtg                                       >  1:28102/1‑139 (MQ=255)
                                                                                                    ttCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTttct                                                                                                                                        <  1:466587/40‑1 (MQ=255)
                                                                                                      cGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCT‑ATGGCCTCATTTCCATTGGCGGCTATGTCgg                                   <  2:285988/138‑1 (MQ=255)
                                                                                                            ctacGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCT‑ATGGCCTCATTTCCATTGGCGGCTATGTCGGCGaccac                            <  2:99439/139‑1 (MQ=255)
                                                                                                            ctacGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCT‑ATGGCCTCATTTCCATTGGCGGCTATGTCGGCGaccac                            <  1:65347/139‑1 (MQ=255)
                                                                                                              acGGCGTGCAGGGCGTACTGGCGGGTTTTTTCGTTTAACAGGTTGGGTTCTCGCAAAAAGAGGCTTTTTTCACTTTTTGTGGTTTTTCTTCGGTGGGGTCATGG‑CTCCTTTCCATTGGGGGGTATTTGGGGGGCCACct                          >  1:99439/1‑139 (MQ=255)
                                                                                                                 gCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCT‑ATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACctgct                       <  2:530294/139‑1 (MQ=255)
                                                                                                                 gCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCT‑ATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACctgct                       <  1:112419/139‑1 (MQ=255)
                                                                                                                 gCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACGTTTGGTGCTTTTGCTGCGCTGGTCT‑ATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACctgct                       <  2:122327/139‑1 (MQ=255)
                                                                                                                      cAGGGCGTACTGGCGGTTTTTTTCGTTTAACAGCTTGGATTTTCGCAAGAGCAGGCTTTTGTCACTTTTTGTGCTTTTTCTGCGCTGGGCT‑ATGGCCTCATTTCCCTTTGCGGCTATGTTGGCGACCACCCGCTggggg                  >  1:9893/1‑138 (MQ=255)
                                                                                                                        gggCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCT‑ATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGAcc                <  2:9893/139‑1 (MQ=255)
                                                                                                                           cGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCT‑ATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCaaa             <  1:64243/139‑1 (MQ=255)
                                                                                                                             tACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCT‑ATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACg           <  2:65695/139‑1 (MQ=255)
                                                                                                                             tACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCT‑ATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACg           <  1:411240/139‑1 (MQ=255)
                                                                                                                              aCTGGCGGTTTTTTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTTCTGCGCTGGTCT‑ATTGCCTCCTTTCCATTGGCGGCTATGTTGGCGGCCCCCCGCTGGGGGCCCAACGc          >  1:394564/1‑139 (MQ=255)
                                                                                                                               cTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCT‑ATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCa         <  2:630276/139‑1 (MQ=255)
                                                                                                                                  gCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCT‑ATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCa      <  1:782850/139‑1 (MQ=255)
                                                                                                                                  gCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCT‑ATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCa      <  2:878088/139‑1 (MQ=255)
                                                                                                                                      ttttCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCT‑ATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGt  <  1:489918/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                       
TATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCT‑ATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGT  >  NZ_CP009273/3634203‑3634475

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