Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,628,6250AC23.3% 91.9 / 12.0 44V105G (GTC→GGC) pdeFcyclic di‑GMP phosphodiesterase PdeF
Reads supporting (aligned to +/- strand):  ref base A (30/3);  new base C (0/10);  total (30/14)
Fisher's exact test for biased strand distribution p-value = 1.49e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAA  >  NC_000913/2628490‑2628764
                                                                                                                                       |                                                                                                                                           
gCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCt                                                                                                                                         <  1:26977/140‑1 (MQ=255)
     tAATCCTAAAAATGCAGAAAACGGGAAAAGTAAAAGAAACAGCGTGGCAGAAAACACAATCTGCCAGAATATACGCGGGAAAATAAAGCGGGTATCCCCATGTGAAACGTTGTTGCGCGGGGGAGTAAAGGCCCTGTAAc                                                                                                                                    <  2:162674/140‑1 (MQ=255)
     aaaTCCTACAAATGCAGCAAACTGAAAAAGTATCAGAAACGGGGTTGGAGAAAACAAAATCGGCAAGAATATCCGCTGGAAAATTAAACGGTTACCACCATGTAAAGCGTTGTTACGCGGGGGAGTAAAGCCCCTGTAAc                                                                                                                                    <  2:29576/139‑1 (MQ=255)
          cTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGCATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGccc                                                                                                                               >  1:318211/1‑140 (MQ=255)
            aCAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCa                                                                                                                             >  2:523744/1‑140 (MQ=255)
            aCAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCa                                                                                                                             >  1:421391/1‑140 (MQ=255)
              aaaTGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGCAGTAAAGACCCTGTAACCCCCCCAAc                                                                                                                           >  2:1340463/1‑140 (MQ=255)
               aaTGCAGCAAACTGAAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCAAGAATATACGCTGGAAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCGGGGGAGTAAAGCCCCTGTAACCGCCCCAACa                                                                                                                          <  1:170541/140‑1 (MQ=255)
                 tGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGCCCCTGTAACCGCCCCAACata                                                                                                                        <  2:858031/140‑3 (MQ=255)
                 ggcagcaAGCGGGAAAAGTATCAGAAACGGCGTTGCAGGAAACACAATCTCCCAGAATATCCGCGGGGAAATAAAGGGGTTATCACCATGTGAAGCGTTTTTCCCCGGGGGAGTAAAGCCCCTGTAACCGCCCCAACaga                                                                                                                        <  2:672732/139‑1 (MQ=255)
                   aagaaaaCTGAAAAAAAATCAGAAACAGCGTGGCAGGAAACCCAATCTCCAAAAATATACGCTGGAAAATTAACCGGGTACCCCCATGTGAAACTTTGTTACGCGGGGGAGTAAAGACCCTGTAACCGCCCCAACAgaga                                                                                                                      <  2:812360/136‑1 (MQ=255)
                        aaCTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaatgat                                                                                                                 >  1:1442288/1‑140 (MQ=255)
                                     acaaaaaCAGCGTTGCAGAAAACACAATCTGCAAGAATATACGCGGGAAAATAAAACGGGTATCCCCATGTGAAACGTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATaaaa                                                                                                    <  2:149705/136‑1 (MQ=255)
                                            cAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCg                                                                                             >  1:924796/1‑140 (MQ=255)
                                               cGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTa                                                                                          >  1:606011/1‑140 (MQ=255)
                                                                 cTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaaa                                                                        >  1:263694/1‑140 (MQ=255)
                                                                     aaGAATATACGCTGGGAATTTAAACGGGTATCACCATGTGAAACGTTGTTCCGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAcc                                                                    <  1:1043119/139‑1 (MQ=255)
                                                                        aaTATACGCGGGGAAATTAAACGGGTATCCCCATGTAAAACGTTGTTCCGCGGGGGAAAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg                                                                 <  2:421391/140‑1 (MQ=255)
                                                                          tataCGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct                                                               >  2:1331083/1‑140 (MQ=255)
                                                                          tataCGCTGGAAAATTAACGGGTTATCCCCATGTAAAACGTTGTTCCGCGGGGGAGAAAAGCCCCTTTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct                                                               <  2:606011/140‑1 (MQ=255)
                                                                           ataCGCTGGAAAATAAACGGGGTATCCCAAGGTAAAATTTTTTTACGCGGGGGAGAAAAGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGctc                                                              <  2:1119915/140‑1 (MQ=255)
                                                                              cGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGc                                                           >  1:1467796/1‑140 (MQ=255)
                                                                              cGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGc                                                           >  2:846022/1‑140 (MQ=255)
                                                                              cGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTCTCGTAAAAACCCAGCTCTGc                                                           >  1:612370/1‑140 (MQ=255)
                                                                                 tGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCata                                                        >  1:467063/1‑140 (MQ=255)
                                                                                  gggAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCatat                                                       >  1:1262745/1‑140 (MQ=255)
                                                                                  gcgAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCatat                                                       >  2:436627/3‑140 (MQ=255)
                                                                                   ggAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAtt                                                      >  1:93397/1‑140 (MQ=255)
                                                                                   ggAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAtt                                                      >  2:661403/1‑140 (MQ=255)
                                                                                   ggAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAtt                                                      >  1:646431/1‑140 (MQ=255)
                                                                                   ggAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAtt                                                      >  2:1458750/1‑140 (MQ=255)
                                                                                   ggAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAtt                                                      >  1:1119915/1‑140 (MQ=255)
                                                                                             cGGGTATCCCCATGGGAAACGTTTTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaa                                            <  1:1338999/140‑1 (MQ=255)
                                                                                                        ttgtaAAACGTTGTTCCGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAgg                                 <  2:924796/139‑1 (MQ=255)
                                                                                                          gtgAAACTTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg                               <  2:70544/140‑1 (MQ=255)
                                                                                                          gtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg                               >  2:1027773/1‑140 (MQ=255)
                                                                                                          gtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg                               >  1:865934/1‑140 (MQ=255)
                                                                                                           tgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGa                              >  2:769487/1‑140 (MQ=255)
                                                                                                             aaaCGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaaa                                           >  2:343956/1‑125 (MQ=255)
                                                                                                             aaaCGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaaa                                           <  1:343956/125‑1 (MQ=255)
                                                                                                                 gttgttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGATAAAAATCCCTAAGTCGATCCCTg                        >  2:1444690/1‑140 (MQ=255)
                                                                                                                 gttgttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTg                        >  2:605024/1‑140 (MQ=255)
                                                                                                                 gttgttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTg                        >  1:70544/1‑140 (MQ=255)
                                                                                                                  ttgttACGCCGGGGAGTAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTgg                       <  1:769487/140‑1 (MQ=255)
                                                                                                                     ttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCaa                    >  1:812812/1‑140 (MQ=255)
                                                                                                                            ggggAGTAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATgccgcc             <  1:1208628/140‑1 (MQ=255)
                                                                                                                             gggAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGccc            >  1:1131499/1‑140 (MQ=255)
                                                                                                                               gAGTAAAGCCCTTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAg          <  1:979442/140‑1 (MQ=255)
                                                                                                                                agaaaaGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGc         <  1:1277801/137‑1 (MQ=255)
                                                                                                                                  tAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTg       >  2:562176/1‑140 (MQ=255)
                                                                                                                                       aCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaaaaa  >  1:977857/1‑140 (MQ=255)
                                                                                                                                       aCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaaaaa  <  2:263694/140‑1 (MQ=255)
                                                                                                                                       aCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaaaaa  >  2:193957/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                           
GCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAA  >  NC_000913/2628490‑2628764

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