Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,875,3310AG65.2% 46.4 / 41.5 46W246R (TGG→CGG) dmlRDNA‑binding transcriptional regulator DmlR
Reads supporting (aligned to +/- strand):  ref base A (4/12);  new base G (12/18);  total (16/30)
Fisher's exact test for biased strand distribution p-value = 3.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCA  >  NZ_CP009273/1875194‑1875454
                                                                                                                                         |                                                                                                                           
tAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCccac                                                                                                                            <  2:1125115/139‑1 (MQ=255)
   aaaTTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAg                                                                                                                         <  1:896452/139‑1 (MQ=255)
           ccGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTCCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAAGGTTTTCGCTAACACCCCGCCAGGAACGCa                                                                                                                  <  2:838009/138‑1 (MQ=255)
              ttATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGc                                                                                                              <  2:857266/139‑1 (MQ=255)
                aTCCGCTCTTTCGCTGGCGTCTCCAGCCTTGAAACATAAACGGCCCAGACGTGCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCTGGAACGCAGAGCa                                                                                                             <  2:808695/138‑1 (MQ=255)
                       cTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCGCCAGGAACGCAGAGCAATCCCCTg                                                                                                     <  2:671543/139‑1 (MQ=255)
                          tCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGcc                                                                                                   >  1:1008474/1‑138 (MQ=255)
                           cGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCGCCAGGAACGCAGAGCAATCCCCTGccc                                                                                                  <  2:896392/138‑1 (MQ=255)
                            gCTGACGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCAt                                                                                                <  1:6505/139‑1 (MQ=255)
                                aCGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCGCCAGGAACGCAGAGCAATCCCCTGCCCATc                                                                                               >  2:84764/1‑136 (MQ=255)
                                aCGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCga                                                                                             <  1:178693/138‑1 (MQ=255)
                                 cGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCGCCAGGAACGCAGAGCAATCCCCTGCCCATCgaga                                                                                           <  2:266643/139‑1 (MQ=255)
                                 cGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCgaga                                                                                           <  1:1214751/139‑1 (MQ=255)
                                     gCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGacac                                                                                        <  2:769346/138‑1 (MQ=255)
                                              gAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATT‑CGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGa                                                                              <  1:503005/138‑1 (MQ=255)
                                               aaaCATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCGCCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGa                                                                              >  2:390040/1‑138 (MQ=255)
                                                        cGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCGCCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTcc                                                                    <  2:1018504/139‑1 (MQ=255)
                                                            ccAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCGCCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGa                                                                <  1:398593/139‑1 (MQ=255)
                                                            ccAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGa                                                                <  2:347841/139‑1 (MQ=255)
                                                             cAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGAt                                                               >  2:130358/1‑139 (MQ=255)
                                                                 cGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCGCCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAga                                                           <  1:226049/139‑1 (MQ=255)
                                                                     cGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCGCCAGGAACGCagag                                                                                                               <  1:168729/83‑1 (MQ=255)
                                                                     cGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCGCCAGGAACGCagag                                                                                                               >  2:168729/1‑83 (MQ=255)
                                                                     cGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCGCCAGGAACGCAGAGCAATCCCCTGCCCATCga                                                                                             <  1:942297/101‑1 (MQ=255)
                                                                     cGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCGCCAGGAACGCAGAGCAATCCCCTGCCCATCga                                                                                             >  2:942297/1‑101 (MQ=255)
                                                                     cGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGCTGGACAATTTCTCCATGATTAga                                                           >  2:777726/1‑135 (MQ=255)
                                                                     cGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAga                                                           <  1:777726/135‑1 (MQ=255)
                                                                              gTAATATTCCGGTAGCACTTGCACTAAATGGCCACTGGCAATGTAGTCGCTAACATCCCGCCAGGAACGCAGAGCAATCCCCAGCCCATCGAGACACCACTGATGGACAATTTCGCCATGATTAGACGACAGTGGGcag                                              <  2:1274957/139‑3 (MQ=255)
                                                                              gTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCGCCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCg                                              <  1:300570/139‑1 (MQ=255)
                                                                                         gTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCGCCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAAt                                   >  2:420834/1‑139 (MQ=255)
                                                                                         gTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCGCCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAAt                                   <  1:313907/139‑1 (MQ=255)
                                                                                         gTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCGCCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAAt                                   >  2:961587/1‑139 (MQ=255)
                                                                                         gTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAAt                                   <  1:22833/139‑1 (MQ=255)
                                                                                               cTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCGCCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATg                             <  1:84764/139‑1 (MQ=255)
                                                                                               cTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATg                             >  1:574977/1‑139 (MQ=255)
                                                                                                  gCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGCCAGTGGGCCGGTGACTTTAATCGCATgg                            <  1:961592/137‑1 (MQ=255)
                                                                                                   cACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCGCCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATgggg                          <  1:420834/138‑1 (MQ=255)
                                                                                                        aaTGACCACTGGCAATGTTTTCGCTAACATCCCGCCAGGAACGCAGAGCa                                                                                                             <  2:3472/50‑1 (MQ=255)
                                                                                                        aaTGACCACTGGCAATGTTTTCGCTAACATCCCGCCAGGAACGCAGAGCa                                                                                                             >  1:3472/1‑50 (MQ=255)
                                                                                                        aaTGACCACTGGCAATGTTTTCGCTAACATCCCGCCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGCCTTTAATCGCATGGGGGCcttc                    >  2:341795/1‑139 (MQ=255)
                                                                                                          tGACCACTGGCAATGTTTTCGCTAACATCCCGCCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAga                                                           >  2:644841/1‑98 (MQ=255)
                                                                                                          tGACCACTGGCAATGTTTTCGCTAACATCCCGCCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAga                                                           <  1:644841/98‑1 (MQ=255)
                                                                                                             ccACTGGCAATGTTTTCGCTAACATCCCGCCAGGAACGCagag                                                                                                               >  2:1190846/1‑43 (MQ=255)
                                                                                                             ccACTGGCAATGTTTTCGCTAACATCCCGCCAGGAACGCagag                                                                                                               <  1:1190846/43‑1 (MQ=255)
                                                                                                              cACTGGCAATGTTTTCGCTAACATCCCGCCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTAtt              >  1:30352/1‑139 (MQ=255)
                                                                                                                  ggCAATGTTTTCGCTAACATCCCGCCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACAACACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCg           >  1:1079512/1‑138 (MQ=255)
                                                                                                                           ttCGCTAACATCCCGCCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCa  <  1:1273124/138‑1 (MQ=255)
                                                                                                                                         |                                                                                                                           
TAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCA  >  NZ_CP009273/1875194‑1875454

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