Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,875,4270G.73.5% 62.7 / 37.1 34coding (640/924 nt)dmlRDNA‑binding transcriptional regulator DmlR
Reads supporting (aligned to +/- strand):  ref base G (5/4);  new base . (13/12);  total (18/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCA  >  NZ_CP009273/1875294‑1875562
                                                                                                                                     |                                                                                                                                       
aCTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGc                                                                                                                                    <  1:515522/139‑1 (MQ=255)
    aaTGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGCcttct                                                                                                                               >  1:593539/1‑139 (MQ=255)
         ccACTGGCAATGTTTTCGCTAACATCACACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGCATTCTGTAtt                                                                                                                          >  2:614484/1‑139 (MQ=255)
          cACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGCCTTCTTTATTg                                                                                                                         >  1:107079/1‑139 (MQ=255)
              ggCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCg                                                                                                                       >  1:22523/1‑138 (MQ=255)
              ggCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGCCTTCTTTATTGCGt                                                                                                                      >  1:749017/1‑138 (MQ=255)
                        tCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGCCTTCTTTATTGCGTAATTGCCAGa                                                                                                            <  2:225097/138‑1 (MQ=255)
                        tCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGCCTTCTTTATTGCGTAATTGCCAGa                                                                                                            <  1:670337/138‑1 (MQ=255)
                              aCATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACAcc                                                                                                        <  2:22523/137‑1 (MQ=255)
                              aCATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCaa                                                                                                      <  1:619819/139‑1 (MQ=255)
                                          gAACGCAGCGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATcgcg                                                                                         >  1:400033/1‑139 (MQ=255)
                                          gAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATcgcg                                                                                         <  1:150921/139‑1 (MQ=255)
                                          gAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATcgcg                                                                                         <  2:107079/139‑1 (MQ=255)
                                            aCGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCgcgc                                                                                        <  2:637536/139‑1 (MQ=255)
                                            aCGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGct                                                                                       <  2:72386/139‑1 (MQ=255)
                                                                        caccacTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCAGGCACGGCAATg                                                               >  1:588728/1‑136 (MQ=255)
                                                                                aTGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTa                                                                                  <  1:425260/108‑1 (MQ=255)
                                                                                aTGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTa                                                                                  >  2:425260/1‑108 (MQ=255)
                                                                                   gACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCAGGCACgg                                                                    >  1:63413/1‑119 (MQ=255)
                                                                                   gACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCAGGCACgg                                                                    <  2:63413/119‑1 (MQ=255)
                                                                                              cATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCg                                       <  1:512784/137‑1 (MQ=255)
                                                                                                  aTTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGgcgtg                                  <  2:400033/138‑3 (MQ=255)
                                                                                                  aTTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGgcgcgc                                 >  1:676135/1‑139 (MQ=255)
                                                                                                            aGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTgcg                        >  2:443065/1‑139 (MQ=255)
                                                                                                            aGTGGGCCGGTGACTTTAATCGCATGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCg                                       >  2:771964/1‑123 (MQ=255)
                                                                                                            aGTGGGCCGGTGACTTTAATCGCATGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCg                                       <  1:771964/123‑1 (MQ=255)
                                                                                                                ggCCGGTGACTTTAATCGCATGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAAt                    >  1:783280/1‑138 (MQ=255)
                                                                                                                   cGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATaaa                 >  1:426570/1‑139 (MQ=255)
                                                                                                                     gTGACTTTAATCGCATGGGGCCTTCTTTATTGCGTAATTGCCAGACAcc                                                                                                        >  1:726784/1‑49 (MQ=255)
                                                                                                                     gTGACTTTAATCGCATGGGGCCTTCTTTATTGCGTAATTGCCAGACAcc                                                                                                        <  2:726784/49‑1 (MQ=255)
                                                                                                                     gTGACTTTAATCGCATGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGcc                                <  1:92356/121‑1 (MQ=255)
                                                                                                                     gTGACTTTAATCGCATGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGcc                                >  2:92356/1‑121 (MQ=255)
                                                                                                                     gTGACTTTAATCGCATGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTc              >  1:473964/1‑139 (MQ=255)
                                                                                                                                    tGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCa  >  2:425288/1‑137 (MQ=255)
                                                                                                                                     |                                                                                                                                       
ACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCA  >  NZ_CP009273/1875294‑1875562

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: