Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,622,8160GA44.6% 50.2 / 42.6 56intergenic (+75/+144)yneM/dgcZinner membrane‑associated protein/diguanylate cyclase, zinc‑sensing
Reads supporting (aligned to +/- strand):  ref base G (27/4);  new base A (0/25);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 4.32e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACT  >  NC_000913/1622668‑1622966
                                                                                                                                                        |                                                                                                                                                      
ttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGa                                                                                                                                                       >  2:321201/1‑150 (MQ=255)
 tATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaa                                                                                                                                                      <  2:934799/150‑1 (MQ=255)
   tGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaaaa                                                                                                                                                    >  1:880002/1‑150 (MQ=255)
    ggCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaaaaa                                                                                                                                                   >  2:439413/1‑150 (MQ=255)
          aCTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGc                                                                                                                                            >  1:795101/1‑151 (MQ=255)
          aCTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGc                                                                                                                                            <  2:241339/151‑1 (MQ=255)
            tGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGcc                                                                                                                                           >  1:914384/1‑150 (MQ=255)
               gAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGt                                                                                                                                        >  1:14085/1‑150 (MQ=255)
                   aatTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAAcc                                                                                                                                   >  2:902765/1‑151 (MQ=255)
                     ttttattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcgg                                                                                                                                 >  1:571449/1‑151 (MQ=255)
                     ttttattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcgg                                                                                                                                 >  2:430642/1‑151 (MQ=255)
                             tctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGaa                                                                                                                          >  1:6318/1‑150 (MQ=255)
                              ctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATg                                                                                                                        >  1:431990/1‑151 (MQ=255)
                                ggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCt                                                                                                                      >  1:164161/1‑151 (MQ=255)
                                 tttttcgggCCGCTTATTCCAGCAAAAAGGGGGATTATTAATGAGCGAAAATAACCCCCCACCAACCCGCCCCTTTGTTACATTTGTGAAAAGGGGGCTTGTTTTTTAGGGCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCt                                                                                                                      <  1:607493/149‑1 (MQ=255)
                                    ttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                                   >  2:74998/1‑150 (MQ=255)
                                        ggCCGCTTTTTCCACCCAAAAAGGGGATGACTAATGAAGGGAAATAACCCCTCCCCAACCCGCCCCTTTTTTACATTTGTGTAAAGGGGGCTTGTTTTTTATGACGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCa                                                                                                                  <  1:988182/147‑1 (MQ=255)
                                        gcccgattattacgcccACAAATGGTATTAATAATGAGGGGAAAAATCCCCTCACCAACCCGGCCCTTTTTTACATTTGTGAAAAGGGGGCGTTTTTTTTATGAGGGGAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATgg                                                                                                               <  2:478673/135‑1 (MQ=255)
                                              gTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAtt                                                                                                        >  1:105372/1‑151 (MQ=255)
                                               tATTTCACCCAAAAAGGGGATGACTAATGAAGGGAAATAACCCCTCCCCTACCGGCCCCCTTTTTAAGTTTGTGTAAAGGGGGCCTGTTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAtt                                                                                                        <  1:134007/150‑1 (MQ=255)
                                                          aaaaGGGGATAACTAATGAAGGGAAATACCCCCCCCCCAACCCGGCCCTTTGTTACGGTTGTGTAAAGGGGGCCTGTTTTTTAGGCGGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtt                                                                                            <  1:203570/150‑1 (MQ=255)
                                                           aaaTGGGAGGATAAATGAACGAAGATATCCCCCCCCCTACCGGCCCCTTTTTTACATTTGTGAACAGGGGCCTTGTTTTTTATAAGGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCttt                                                                                           <  1:433603/151‑1 (MQ=255)
                                                               gggATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGTTTTTTATGACGGCGAAAAAAAGCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACg                                                                                       <  2:224697/151‑1 (MQ=255)
                                                               gggATGACTAAGGACGGGAAATAACCCCTCACCTCCCCGCCCCTTTTTTACATTTTTAAACAAGGGGTTTTTTTTTTGTGGCGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACg                                                                                       <  2:940419/151‑1 (MQ=255)
                                                                   tgaaaaaTGAAGGGAAATACTCCCTCCCAAACCCGCCCCTTGTTTATAGTGGTGCAGAGGGGGCCTATTTTTTATGGCGGAGAAAAAAAGCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCt                                                                                     <  1:813292/144‑1 (MQ=255)
                                                                     aCTAATGAGGGGAAATATCCCCCCCCATACCCGCCCCTTTTTTACTTTTGTCAAAAGGGGGCTTGTTTTTTATGGCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACa                                                                                 <  2:819459/151‑1 (MQ=255)
                                                                     aCTAATGAAGGGAAAAAACCCCTCCCCTACCCGGCCCTTTGTTACATTTGTGTACAGGGGGCTTGTTTTTTATGCGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACa                                                                                 <  1:826876/151‑1 (MQ=255)
                                                                        aaTGAGGGGAGATACCCCCTCCCCAACCGGCCCTCTTTTTAGATTTGTGAACAGGGGGCTTTTTTTTTATGCCGGGAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAg                                                                                <  1:821586/149‑1 (MQ=255)
                                                                        aaTAAACGGAGATACCCCCTCCCCAACCCGCCCCTTTGTTACTGTTGTGTAAAAGGGGCTTGTTTTTTAGGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAg                                                                                <  1:522211/149‑1 (MQ=255)
                                                                             aCGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAtttt                                                                         >  1:386770/1‑151 (MQ=255)
                                                                                  gaTAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGc                                                                    >  2:282298/1‑151 (MQ=255)
                                                                                        ccccccccaacccGCCCCTTTGTTAAATTTGTTAAAAGGGGGCTTTTTTTTTAGGGCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCAc                                                                  <  1:530029/136‑1 (MQ=255)
                                                                                        cccTCCCCTACCCGCCCCTTTTTTACATTTGTGTACAGGGGGCTTTTTTTTTATGCCGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACt                                                                 <  1:36973/148‑1 (MQ=255)
                                                                                         ccTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGGGAATGGCTGCATGGGTAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTTCTGCCCt                                                                 >  2:384888/1‑147 (MQ=255)
                                                                                            ccccTACCCGGCCCTTTGTTACATTTGTGTACAGGGGGCTTGTTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATaa                                                            <  1:201267/147‑1 (MQ=255)
                                                                                             aCCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACg                                                          >  1:945470/1‑150 (MQ=255)
                                                                                                    cGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTg                                                  >  1:515935/1‑151 (MQ=255)
                                                                                                       ccctTTGTTACAGTTGTGTACAGGGGGCTTTTTTTTTATGACGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGa                                                 <  1:141585/149‑1 (MQ=255)
                                                                                                          cccttggttaagtttGTGTAAAGGGGCCTTGTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTg                                                  <  2:838026/135‑1 (MQ=255)
                                                                                                                 ttacttttGTGAAAAGGGGGCTTTTTTTTTAGGCCGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTa                                          <  1:844999/143‑1 (MQ=255)
                                                                                                                 ttacgtttGTGTAAAGGGGGCTTGTTTTTTAGGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAg                                         <  1:846791/144‑1 (MQ=255)
                                                                                                                 ttacatgtgggaaCAGGGGGCCTTTTTTTTATGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTa                                          <  1:539115/137‑1 (MQ=255)
                                                                                                                 taactgTTGTGTAAAGGGGCCTTGTTTTTTAAGACGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAg                                         <  1:902875/145‑1 (MQ=255)
                                                                                                                 atACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGCGAGATTTGGGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTa                                          >  2:185462/2‑149 (MQ=255)
                                                                                                                   aCAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCAt                                      >  1:298939/1‑151 (MQ=255)
                                                                                                                   aCAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCAt                                      >  2:132591/1‑151 (MQ=255)
                                                                                                                   aCAGTTGTGAACAGGGGGCTTTTTTTTTATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCAt                                      <  1:255410/151‑1 (MQ=255)
                                                                                                                       ttGTGAAAAGGGGGCCTTTTTTTTATGGCGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGc                                        <  1:264236/145‑1 (MQ=255)
                                                                                                                          tgtACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGtttt                               >  2:868784/1‑151 (MQ=255)
                                                                                                                            tACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtct                             >  1:477920/1‑151 (MQ=255)
                                                                                                                                  gggccttgttttttATGAGGAAAAAAAAAATCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGt                                  <  2:164148/132‑1 (MQ=255)
                                                                                                                                        tGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGtt                  >  1:850904/1‑151 (MQ=255)
                                                                                                                                             tttgtgcgggCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAAc              <  2:641435/142‑1 (MQ=255)
                                                                                                                                                  gACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCt        >  2:709833/1‑151 (MQ=255)
                                                                                                                                                     ggCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAg          <  1:731359/146‑1 (MQ=255)
                                                                                                                                                       cGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAAc   >  1:813681/1‑151 (MQ=255)
                                                                                                                                                        gAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACt  >  2:564785/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                      
TTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACT  >  NC_000913/1622668‑1622966

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