Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 578,568 G→A *154* (TAG→TAA rzpD → DLP12 prophage, putative murein endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913578,5680GA98.6% 228.0 / ‑6.5 72*154* (TAG→TAArzpDDLP12 prophage, putative murein endopeptidase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (35/36);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCT  >  NC_000913/578438‑578693
                                                                                                                                  |                                                                                                                             
tGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGAcc                                                                                                                       >  2:314711/1‑139 (MQ=255)
      aTGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCg                                                                                                                 <  1:392043/139‑1 (MQ=255)
            ccTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGc                                                                                                           >  1:370972/1‑139 (MQ=255)
            ccTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGc                                                                                                           >  1:1208519/1‑139 (MQ=255)
                                    cTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCaa                                                                                   >  1:285987/1‑139 (MQ=255)
                                      gAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATa                                                                                 <  1:753702/139‑1 (MQ=255)
                                      gAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATa                                                                                 <  2:877557/139‑1 (MQ=255)
                                            gATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACAcg                                                                           <  1:106177/139‑1 (MQ=255)
                                                 tttCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCtt                                                                      >  2:22788/1‑139 (MQ=255)
                                                     aCCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                  <  1:314711/139‑1 (MQ=255)
                                                       ccTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCg                                                                >  1:433229/1‑139 (MQ=255)
                                                                 agGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAAt                                                      >  2:888933/1‑139 (MQ=255)
                                                                 agGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAAt                                                      >  1:278745/1‑139 (MQ=255)
                                                                  ggCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                  <  2:198884/126‑1 (MQ=255)
                                                                  ggCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                  >  1:198884/1‑126 (MQ=255)
                                                                       aTCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATacaca                                                                             <  1:907020/110‑1 (MQ=255)
                                                                       aTCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATacaca                                                                             >  2:907020/1‑110 (MQ=255)
                                                                         cACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAg                                              <  2:433229/139‑1 (MQ=255)
                                                                             aTGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGtaat                                          >  2:1024791/1‑139 (MQ=255)
                                                                                 aaaaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATaaaa                                      <  2:285987/139‑1 (MQ=255)
                                                                                 aaaaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATAGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATaaaa                                      <  1:689699/139‑1 (MQ=255)
                                                                                   aaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCat                                                                                                                     <  2:686403/58‑1 (MQ=255)
                                                                                   aaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCat                                                                                                                     >  1:686403/1‑58 (MQ=255)
                                                                                    aacaacTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCg                                   >  1:304238/1‑139 (MQ=255)
                                                                                      caacTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCg                                                                                                                 >  2:827827/1‑59 (MQ=255)
                                                                                      caacTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCg                                                                                                                 <  1:827827/59‑1 (MQ=255)
                                                                                      caacTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAg                                 >  1:886451/1‑139 (MQ=255)
                                                                                      caacTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAg                                 <  1:888933/139‑1 (MQ=255)
                                                                                        acTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCg                                                                >  2:684917/1‑106 (MQ=255)
                                                                                        acTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCg                                                                <  1:684917/106‑1 (MQ=255)
                                                                                          tGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAAt                             <  1:528229/139‑1 (MQ=255)
                                                                                          tGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAAt                             <  1:22788/139‑1 (MQ=255)
                                                                                                     cAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGaa                  <  1:1024791/139‑1 (MQ=255)
                                                                                                      aGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAAt                 <  2:370972/139‑1 (MQ=255)
                                                                                                        aaGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATgg                                                                                                             <  2:922625/45‑1 (MQ=255)
                                                                                                        aaGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATgg                                                                                                             >  1:922625/1‑45 (MQ=255)
                                                                                                             tatTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCa                                                                                                     <  2:203876/48‑1 (MQ=255)
                                                                                                             tatTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCa                                                                                                     >  1:203876/1‑48 (MQ=255)
                                                                                                             tatTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATacacac                                                                            >  1:1267289/1‑73 (MQ=255)
                                                                                                             tatTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATacacac                                                                            <  2:1267289/73‑1 (MQ=255)
                                                                                                             tatTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAg                                              >  1:1085358/1‑103 (MQ=255)
                                                                                                             tatTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAg                                              <  2:1085358/103‑1 (MQ=255)
                                                                                                             tatTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCt          >  1:1236005/1‑139 (MQ=255)
                                                                                                              atTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTc                                                                     <  1:1014381/79‑1 (MQ=255)
                                                                                                              atTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTc                                                                     >  2:1014381/1‑79 (MQ=255)
                                                                                                              atTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCg                                                                <  2:1030276/84‑1 (MQ=255)
                                                                                                              atTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCg                                                                >  1:1030276/1‑84 (MQ=255)
                                                                                                              atTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTg         <  2:304238/139‑1 (MQ=255)
                                                                                                                  aTGAGCAGTGCAGATAAAGCTGACCATATCGATGGGc                                                                                                           <  2:135696/37‑1 (MQ=255)
                                                                                                                  aTGAGCAGTGCAGATAAAGCTGACCATATCGATGGGc                                                                                                           >  1:135696/1‑37 (MQ=255)
                                                                                                                  aTGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTAtttt                                                                                         <  1:154786/55‑1 (MQ=255)
                                                                                                                  aTGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTAtttt                                                                                         >  2:154786/1‑55 (MQ=255)
                                                                                                                  aTGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCg                                                                <  1:831059/80‑1 (MQ=255)
                                                                                                                  aTGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCAAGCg                                                                >  2:831059/1‑80 (MQ=255)
                                                                                                                     aGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTaa                                                >  1:1149906/1‑93 (MQ=255)
                                                                                                                     aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTAtttt                                                                                         >  2:1044087/1‑52 (MQ=255)
                                                                                                                     aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTAtttt                                                                                         <  1:1044087/52‑1 (MQ=255)
                                                                                                                     aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                  <  1:1133772/75‑1 (MQ=255)
                                                                                                                     aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                  >  2:1133772/1‑75 (MQ=255)
                                                                                                                     aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCgg                                                               >  2:760180/1‑78 (MQ=255)
                                                                                                                     aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCgg                                                               <  1:760180/78‑1 (MQ=255)
                                                                                                                     aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTaa                                                <  2:1149906/93‑1 (MQ=255)
                                                                                                                     aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATcc                           <  2:3421/114‑1 (MQ=255)
                                                                                                                     aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATcc                           >  1:3421/1‑114 (MQ=255)
                                                                                                                     aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCt  >  2:366554/1‑139 (MQ=255)
                                                                                                                          tGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACAcgc                                                                          <  1:608729/62‑1 (MQ=255)
                                                                                                                          tGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACAcgc                                                                          >  2:608729/1‑62 (MQ=255)
                                                                                                                          tGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                  <  1:1038733/70‑1 (MQ=255)
                                                                                                                          tGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                  >  2:1038733/1‑70 (MQ=255)
                                                                                                                            cAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCg                                                                >  2:74023/1‑70 (MQ=255)
                                                                                                                            cAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCg                                                                <  1:74023/70‑1 (MQ=255)
                                                                                                                               aTAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCg                                                                >  1:1050392/1‑67 (MQ=255)
                                                                                                                                  |                                                                                                                             
TGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCT  >  NC_000913/578438‑578693

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