Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,457,1050AG46.4% 61.6 / 19.7 56I441V (ATC→GTC) mplUDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base A (1/29);  new base G (26/0);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 4.06e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGA  >  NC_000913/4456965‑4457249
                                                                                                                                            |                                                                                                                                                
ggcaggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGaa                                                                                                                                        <  1:361028/151‑1 (MQ=255)
       ggcagaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAActgctg                                                                                                                                 <  2:577419/151‑1 (MQ=255)
       ggcagaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAActgctg                                                                                                                                 <  2:408334/151‑1 (MQ=255)
        gcagaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATACTCTGGTGATGAGCAACGGCGGGTTTGGTGGGGTCCATCAGAAACTGCTgg                                                                                                                                >  2:423748/1‑151 (MQ=255)
          agaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGAt                                                                                                                              <  2:92432/151‑1 (MQ=255)
                ggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTTCATCAGAAACTGCTGGAGGGGGTg                                                                                                                        >  1:789283/1‑151 (MQ=255)
                    gaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCga                                                                                                                    <  2:425751/151‑1 (MQ=255)
                      agCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGTGGGGTCCCGCAGAAACTGCTGggtggggtgg                                                                                                                       >  1:152647/1‑146 (MQ=255)
                      agCCTGCGTTCAGCCTGCACACTGGAGTGGCGATATGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGAACATATTCTGGTGATGAGCAACGGCGGATTTTGTGGGGTGCATACAAAAATGGCGGGTGGGCTGGCg                                                                                                                     >  1:457022/1‑148 (MQ=255)
                          tGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGGTGGGCTGGCGAAGAAgg                                                                                                              >  2:16602/1‑151 (MQ=255)
                           gCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATGCGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGAGATGAGCAACGGCGGTTTTGGTGCGATCCATCAGAAACTGCTGGATGGGCTGACgaagaatg                                                                                                              <  1:535815/150‑3 (MQ=255)
                             gTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAACTGCTGGGTGGGCTGGCGGAGAAGgggg                                                                                                           >  2:5829/1‑151 (MQ=255)
                                cAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGaa                                                                                                         <  2:763237/150‑1 (MQ=255)
                                 aGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAga                                                                                                       <  2:175080/151‑2 (MQ=255)
                                  gCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCCTCAGAAAATGCTGGATGGGCTGGCGAAGAAGGCg                                                                                                            >  2:661517/1‑145 (MQ=255)
                                   ccggCACACTGGAGTGGCGATGTGGATTCGCTGGCAGATATGGTGGTGAAAACCGCTTAGCCTGGCGACCATACTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTTGCTGGGATGGCGAAGAAGGCGGAAGcg                                                                                                      <  1:576892/147‑2 (MQ=255)
                                     tGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGGGGGGTCCCTCAGGAAATGCTGGATGGGGTGGGGAAGAAGgggg                                                                                                           >  1:706324/1‑143 (MQ=255)
                                          acTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCCTCCGAAACTGCTGGGTGGGGTGGCGAAGAAGGGGGGAGcc                                                                                                      >  1:359109/1‑143 (MQ=255)
                                           cTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTGGGTGGGATCCATCAGAAACTGCTGGAGGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTa                                                                                              <  1:49425/150‑1 (MQ=255)
                                            tGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCCTCAGAAACTGCTGGATGGGCTGGCGAAGAAAGCGGGAGCCGCGCCGTaa                                                                                             >  2:207536/1‑150 (MQ=255)
                                             ggAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCTGGGTGGGGTGGCGAAGGAAGCGGAAGCCGCCCCGGAAtt                                                                                           >  1:665310/1‑151 (MQ=255)
                                              gAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTc                                                                                          <  1:597794/151‑1 (MQ=255)
                                               aGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCg                                                                                         <  1:446471/151‑1 (MQ=255)
                                                        gtgGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGcc                                                                                <  1:700801/151‑1 (MQ=255)
                                                             tACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAGAAAATGCTGGATGGGCTGGCGAAGAAGGCGGGAGGCGCGCAGGAATTCCGCCTTAGGCAggt                                                                            >  2:172882/1‑148 (MQ=255)
                                                                 cTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAAATGCTGGATGGGCTGGCGAAGAAAGCGGAAGGCGCGCAGGAATTCGGCCTTAGCCAGATATAg                                                                        >  2:552616/1‑150 (MQ=255)
                                                                  tGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCCTCCGAAAATGCTGGATGGGCTGGCGGAGAAGGCGGAAGCCGCGCCGTAATTTGGCCTTAGCCCGGtata                                                                         >  2:380826/1‑148 (MQ=255)
                                                                   ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCCTCAGAAACTGCTGGATGGGCTGGCGGAGAAGGCGGGAGCCGCGCCGTAATTTGGCCTTAGCCAGATATAGt                                                                       >  1:267172/1‑149 (MQ=255)
                                                                   ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTa                                                                      <  2:471108/150‑1 (MQ=255)
                                                                   ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTAt                                                                     <  2:159451/151‑1 (MQ=255)
                                                                     cAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGc                                                                   <  1:210882/151‑1 (MQ=255)
                                                                         tatGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCCTCAGAAACTGCTGGATGGGCTGGGGAAGAAGGCGGAAGCCGCGCCGTAATTCGGCCTTTGCCCGATATAGTATGCg                                                                  >  1:243260/1‑148 (MQ=255)
                                                                          atGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGATGGGCTGGCGAAGAAGGCGGGAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAgg                                                              >  1:779227/1‑151 (MQ=255)
                                                                          atGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAgg                                                              <  2:83790/151‑1 (MQ=255)
                                                                               ggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCgg                                                          <  1:14767/150‑1 (MQ=255)
                                                                                  gAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATa                                                       <  1:657601/150‑1 (MQ=255)
                                                                                     aaCCGCTCAGCCGGGCGACCATATTCTGGTGATGAGCAACGGCGGCTTTGGTGGGGGCCATTAGAAAATTCTGGTTGGGCCGGCGGAGAAGGCGGAAACCGCGGTGTAATTTGGCCTTTAATAAAAAAAATATGCGTAGGGCGGGTaa                                                      >  1:305078/1‑148 (MQ=255)
                                                                                           tCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCCTCAGAAACTGCTGGGTGGGCTGGCGGAGAAGGCGGGAGCCGCGCCGTAATTCGGGCCTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGGTCa                                              >  2:616416/1‑150 (MQ=255)
                                                                                            cAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAGAAAATGCTGGATGGGGTGGCGGAGAAGGCGGGAGGCGCGCCGTAATTCGGCCTTAGCCCGATATAGTATGGGTAGGCCGGGAAAGGCGTTCa                                              >  1:669007/1‑149 (MQ=255)
                                                                                             aGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAACTGCTGGGTGGGCTGGCGGAGAAAGCGGGAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGTTCAcgc                                           >  1:92280/1‑151 (MQ=255)
                                                                                               ccTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCCTCAGAAACTGCTGGATGGGCTGGCGAAGAAAGCGGGAGGCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAcgcc                                          >  2:740136/1‑150 (MQ=255)
                                                                                                cTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCTGGATGGGGGGGCGGAGAAAGCGGGAGGCGCGCAGTAAATCGGCCTTAGCCAGGTATAGTATGCGTAGGCCGGGTAAGGCGTTCAcgccgc                                        >  1:219883/1‑151 (MQ=255)
                                                                                                  ggCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCAt                                      <  2:676918/151‑1 (MQ=255)
                                                                                                         atatTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAGAAACTGCTGGGTGGGGTGGCGGAGAAAGGGGGAGCCGCGCCGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGTTCACGCCGCATcc                                    >  1:406989/1‑146 (MQ=255)
                                                                                                           atTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAtt                              <  1:160132/150‑1 (MQ=255)
                                                                                                           atTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAtt                              <  2:700762/150‑1 (MQ=255)
                                                                                                            ttCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTTCTGGATGGGCTGGCGGAGAAGGGGGGAGCCGCGCAGGAATTCGGCCTTAGCCAGATATAATATGCGGAGGCCGGGTAAGGCGGTCACGCCCCATCCGGCAttt                             >  1:432280/1‑150 (MQ=255)
                                                                                                            ttCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAttt                             <  1:8818/150‑1 (MQ=255)
                                                                                                                ggTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAg                          <  2:675810/149‑1 (MQ=255)
                                                                                                                  tgatgaGCAACGGCGGTTTTGGTGGGGTCCAACAGAAAATTCTGGATGGGGTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATTCGTAGGCCGGATAAGGCGTTCACGCCGCCTCCGGCATTTGAggct                       >  1:349418/1‑147 (MQ=38)
                                                                                                                   gatgaGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAg                     <  2:703523/151‑1 (MQ=255)
                                                                                                                    atgaGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGt                    <  1:109722/151‑1 (MQ=255)
                                                                                                                     tgaGCAACGGCGGTTTTTGGGGGGTCCATCAGAAACTGCTGGATGGGGTGGCGGAGAAGGCGGGAGGCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGTTCACGCCGCATCCGGCATTTTAGCATAGTg                   >  1:183913/1‑151 (MQ=255)
                                                                                                                              gcggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCg          <  2:780176/151‑1 (MQ=255)
                                                                                                                                    ttGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCt     <  2:789209/150‑1 (MQ=255)
                                                                                                                                      ggtggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGa  <  1:742313/151‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                
GGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGA  >  NC_000913/4456965‑4457249

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