Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,329,8940AG45.1% 70.6 / 12.3 51D188G (GAC→GGC) sohBinner membrane protein, S49 peptidase family protein
Reads supporting (aligned to +/- strand):  ref base A (7/21);  new base G (23/0);  total (30/21)
Fisher's exact test for biased strand distribution p-value = 1.08e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTG  >  NC_000913/1329745‑1330043
                                                                                                                                                      |                                                                                                                                                               
gCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGa                                                                                                                                                                 <  1:555534/150‑1 (MQ=255)
   ccGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACa                                                                                                                                                              <  1:401863/149‑1 (MQ=255)
       cTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTTCATTATTGGCTGGGTTGCGGGCaaaa                                                                                                                                                           >  2:721835/1‑148 (MQ=255)
         ggAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGt                                                                                                                                                       <  1:95088/150‑1 (MQ=255)
            aaGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTTCAATATTGGCCGTGTGGGGGGCCAAAATtttt                                                                                                                                                     >  1:419237/1‑149 (MQ=255)
             aGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTTCATTATGGCCCGTGTGGGGGGCAAAAATttttt                                                                                                                                                    >  1:236034/1‑148 (MQ=255)
                  tggtggCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATTATGGGCTGTGTGGCGGGGAAAAATTTTTTCGCAcc                                                                                                                                             >  1:676960/1‑150 (MQ=255)
                     tggCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTATATAAAGGGCTCTGTGGGGGGGAAAAATGTTTTCCCAcccctt                                                                                                                                         >  2:10368/1‑147 (MQ=255)
                     tggCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATTATGGCCTGTGTGGGGGGCAAAAATGTTTTCGCAcc                                                                                                                                             >  1:779249/1‑147 (MQ=255)
                      ggCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTg                                                                                                                                       <  1:284274/151‑1 (MQ=255)
                        catggtgcatggtTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAACACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTTCATTATGGCCCGTGTGGGGGGCAAAAATTTTTTCTCAcc                                                                                                                                             >  1:440369/1‑144 (MQ=255)
                               catggtTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTTCATTAATGTCCCTGGGGGGGGCAAAAATGTTTCCGCCCCCGTTTgttttttgg                                                                                                                               >  2:38736/1‑143 (MQ=255)
                               catggtTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTgg                                                                                                                            <  2:647144/151‑1 (MQ=255)
                                  ggtTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGtt                                                                                                                         <  1:224688/151‑1 (MQ=255)
                                           cTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATTATGGGCTTTGTTGGGGGGAAAAATTTTTTCGCACCCCTTGCtttt                                                                                                                                 >  2:512672/1‑135 (MQ=255)
                                           cTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATTATGGCCTGTGTGGGGGGGAAAAATATTTTCCCACCCGTTTTTATTTTGGGGTTCCTTTgcg                                                                                                                 >  1:291270/1‑149 (MQ=255)
                                            tGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTTAAAGATGGCCCGTGTGGGGGGGAAAAATTTTTTCCCCCCGGTTGCTATTTTTGGGTTccttgggg                                                                                                                 >  1:655808/1‑144 (MQ=255)
                                                    tCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGgc                                                                                                       <  1:100563/149‑1 (MQ=255)
                                                      gCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATTATGGGCTGTGTGGGGGGGAAAAATTTTTTCGCACCCGTTGCTATTTTTGGGTTCCTTTggggggggggg                                                                                                         >  2:303850/1‑139 (MQ=255)
                                                       cAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCaaa                                                                                                <  1:208758/151‑1 (MQ=255)
                                                        aaCTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATTATGGGCTGTGTGGGGGGCAAAAATTTTTTCGCCCCCCTTTCTTTTTTGGGGTTCCTTTgcggggggggggcaa                                                                                                     >  1:425345/1‑136 (MQ=255)
                                                           tGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTAAATTATTGGCCGTGGGGGGGACAAAAATTTTTTCCCCCCGGTTTTTTTTGTTGGTTGCAATTGGGTGGGGGCGGaaaaacc                                                                                            >  1:196607/1‑147 (MQ=255)
                                                             cAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGGTACATGAAGGTCTCTGGGGGGGGAAAAAATTTTTTCGCACCGCTTTCTCTTTTGGGTTCCATTTGCGGGGGGGGGCgaaa                                                                                               >  2:992073/1‑145 (MQ=255)
                                                                cGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACt                                                                                      <  1:930378/151‑1 (MQ=255)
                                                                 gTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATTAATGGCTGTGGGGGGGGGAAAAATTTTTTCCCCCCCTTTTGTTTTTTGGGGTTTCTTTGCGtgggggggggaaagacc                                                                                             >  2:1040749/1‑135 (MQ=255)
                                                                 gTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGATATCTTATGGGCTTTGTGGGGGGCGAAAATTGTTTCCCACCCTTTTCTATTGTGGGGTTCTTTTGCGGGGGGGGGTGAAAGCCCCACtt                                                                                     >  2:382643/1‑151 (MQ=255)
                                                                     gcgcGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTaa                                                                                  <  2:173059/150‑1 (MQ=255)
                                                                       gcgATAAAAACATTCCTTTAACTGTGACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCAGCGTTTGCTATAGTGGCGTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTa                                                                                   <  1:489360/147‑1 (MQ=255)
                                                                          aTAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTCTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCtt                                                                            <  1:228599/151‑1 (MQ=255)
                                                                           tAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTAAATGATGGGCTGTGGGGCGGGCAAAAATGTTTTCCCCCCCGTTGCTATTTTGGGGTTCATTTGGGTGGGGGGGCGAAATCCCAACCTTTAACCCCt                                                                             >  1:71286/1‑150 (MQ=255)
                                                                            aaaaaCATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTcc                                                                          <  1:623226/151‑1 (MQ=255)
                                                                                 cATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGGCTGTGGGGGGGGCAAAAATGTTTCTCCACCCTCTGCTATTTTGGGTTTCCTTTGTGTGGGGGCGCCAAATCCCCACCTTTaccacatcctgg                                                                        >  1:345768/1‑138 (MQ=255)
                                                                                    tCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTTAATTATAGCCCGGGTGGGGGGCAAAATTGTTTTCCCCCCCTTTTCTTTTTTTGGTTTCCTTGGCGGGGGGGGGGCAAAACCCCACTTTTAACCCCTTCCGGaaa                                                                     >  1:578055/1‑149 (MQ=255)
                                                                                      cTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATTATGGGCTGTGTGGGGGGGAAAAATTTTTTCGCACCCCTTTCTATTTTGGGGTTCCTTGGCGGGGGGGGGCAAATACCCCAATTTTAACCCCtcctcgaaa                                                                     >  1:392306/1‑139 (MQ=255)
                                                                                          aaCTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGCAAAGat                                                            <  1:10369/151‑1 (MQ=255)
                                                                                              gTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATTATGGCCCGGGCGGGGGAGAAAAATTTTTTTCCCCCGGTTTTTATTGTGGGGTTCATTGGCGTGGGGGGGCAAAACCCCCACCTTTTCCCCCTCCTTAAAAACAAAGaaa                                                           >  1:1041853/1‑146 (MQ=255)
                                                                                              gTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTg                                                        <  1:394945/151‑1 (MQ=255)
                                                                                                tACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATTATTGCCTGGGTGGGGGGAAAAAATGTTTTCCCCCCGCTTTTTTTTTTTGGGTTCCTTTGCGGGGGGGGGCCAAAGGCCCACCTTTACCCCCTCCTTAAAAACAAAAAAAtt                                                         >  2:387363/1‑148 (MQ=255)
                                                                                                         cAAAGTCGCTGCCAGCGGCGGGTATATTATGGGCTGTGTGGCGGACAAAAATGTTTTCCCACCCGTTGCTTTTTTGGGGTTCATTGGGGTGGGGGGGGCAACTACCCAACTTAAACCCCTTCTTAAAAGCCAAGATATTGata                                                     >  2:722722/1‑143 (MQ=255)
                                                                                                               cGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCacac                                        <  1:465242/150‑1 (MQ=255)
                                                                                                               cGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCACACt                                       <  2:625823/151‑1 (MQ=255)
                                                                                                                 cTGCCAGCGGCGGTTACATGATGGCCTGTGTGGGGGGGAAAAATTTTTTCCCACCCCTTGCTATTTTGTGTTTCATTTGCGGGGGGGGGCAAAAACCCCACTTTTATCCCCTCCCGAAAAAAACAGAAATTTATATctca                                                >  2:796023/1‑137 (MQ=255)
                                                                                                                     cAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCACACTGCCggg                                 <  1:444462/151‑1 (MQ=255)
                                                                                                                       gcggcgGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCACACTGCCGGGCa                               <  1:1007927/151‑1 (MQ=255)
                                                                                                                        cggcggTTACATGATGGCCTGTGTGGGGGGCAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTTCATTGGCGGGGGGGGGCAAAAGCCCAACTTTTAACCCGTCCTGAAAAGAAAAAATATTGATATCGAAATGCACACTGCCggg                                 >  2:951926/1‑148 (MQ=255)
                                                                                                                           cggTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCACACTGCCGGGCAGtat                           <  2:917352/151‑1 (MQ=255)
                                                                                                                            ggTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCACACTGCCGGGCAGtata                          <  1:818598/151‑1 (MQ=255)
                                                                                                                            ggTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCACACTGCCGGGCAGtat                           >  1:994385/1‑150 (MQ=255)
                                                                                                                                       ggCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTg               >  1:44694/1‑151 (MQ=255)
                                                                                                                                                 ggCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTg      <  2:680997/151‑1 (MQ=255)
                                                                                                                                                   cGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTggg    >  1:784363/1‑151 (MQ=255)
                                                                                                                                                    ggACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTggg    >  1:810437/1‑150 (MQ=255)
                                                                                                                                                      aCAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTg  >  2:183154/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                               
GCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTG  >  NC_000913/1329745‑1330043

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