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.5% 218.3 / NA 68*154* (TAG→TAArzpDDLP12 prophage, putative murein endopeptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (33/34);  minor base T (1/0);  total (34/34)
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

CCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAAC  >  NC_000913/578450‑578701
                                                                                                                      |                                                                                                                                     
ccTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGc                                                                                                                   >  1:225086/1‑139 (MQ=255)
     ccccGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGa                                                                                                                                 >  2:864990/1‑120 (MQ=255)
     ccccGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGa                                                                                                                                 <  1:864990/120‑1 (MQ=255)
      cccGACTGGCAGACACCGCTGAACGGGATTATTTCACCCACAGAGAGTGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCa                                                                                                             <  2:1783417/139‑1 (MQ=255)
       ccGACTGGCAGACACCGCTGCACGGGATTCTTCCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTCTATTAATGAGCAGTGCAGATATAGCTGACAATATCGATGGGCAACTCAt                                                                                                            >  1:351652/1‑139 (MQ=255)
                        cTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCa                                                                                                             >  2:355488/1‑121 (MQ=255)
                        cTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCa                                                                                                             <  1:355488/121‑1 (MQ=255)
                                 attattTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACAcgc                                                                                  <  1:1463808/139‑1 (MQ=255)
                                     tttCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCt                                                                                                                                   <  2:1656873/86‑1 (MQ=255)
                                     tttCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCt                                                                                                                                   >  1:1656873/1‑86 (MQ=255)
                                         aCCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                          <  2:1124179/139‑1 (MQ=255)
                                          cccTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGc                                                                         >  2:512202/1‑139 (MQ=255)
                                             tCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGa                                                                      <  2:757575/139‑1 (MQ=255)
                                               agagagagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTc                                                                             <  2:1733209/130‑1 (MQ=255)
                                               agagagagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTc                                                                             >  1:1733209/1‑130 (MQ=255)
                                               agagagagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGt                                                                    >  1:1468149/1‑139 (MQ=255)
                                               agagagagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGt                                                                    >  2:807251/1‑139 (MQ=255)
                                               agagagagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGt                                                                    >  2:381497/1‑139 (MQ=255)
                                                gagagagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTg                                                                                                <  2:1164857/110‑1 (MQ=255)
                                                gagagagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTg                                                                                                >  1:1164857/1‑110 (MQ=255)
                                                gagagagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACAcg                                                                                   >  2:1736276/1‑123 (MQ=255)
                                                gagagagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACAcg                                                                                   <  1:1736276/123‑1 (MQ=255)
                                                  gagagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGtata                                                                 <  2:742758/139‑1 (MQ=255)
                                                     agGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACAcg                                                                                   <  2:743045/118‑1 (MQ=255)
                                                     agGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACAcg                                                                                   >  1:743045/1‑118 (MQ=255)
                                                        cTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGcc                                                           <  1:512202/139‑1 (MQ=255)
                                                           aTCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTg                                                                                                >  1:1480596/1‑99 (MQ=255)
                                                           aTCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTg                                                                                                <  2:1480596/99‑1 (MQ=255)
                                                           aTCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCa                                                                                            <  2:1395644/103‑1 (MQ=255)
                                                           aTCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCa                                                                                            >  1:1395644/1‑103 (MQ=255)
                                                                 aTGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGtaat                                                  >  2:877475/1‑139 (MQ=255)
                                                                 aTGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGAACATATCGATGGGCAACTAATGCAATTATGTTGGGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGtaaa                                                  >  1:1695026/1‑138 (MQ=255)
                                                                     aaaaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATaaaa                                              <  1:807251/139‑1 (MQ=255)
                                                                          caacTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAg                                         <  2:1695026/139‑1 (MQ=255)
                                                                            acTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCt                                                                                                                                   <  2:1621392/47‑1 (MQ=255)
                                                                            acTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCt                                                                                                                                   >  1:1621392/1‑47 (MQ=255)
                                                                            acTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTAtttt                                                                                                 <  1:557238/81‑1 (MQ=255)
                                                                            acTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTAtttt                                                                                                 >  2:557238/1‑81 (MQ=255)
                                                                              tGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAAt                                     <  2:1468149/139‑1 (MQ=255)
                                                                                          aGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                          >  2:1418257/1‑90 (MQ=255)
                                                                                          aGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGAAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                          <  1:1418257/90‑1 (MQ=255)
                                                                                                 tatTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCt                  >  2:778567/1‑139 (MQ=255)
                                                                                                    tAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGtaa                                                   >  2:578924/1‑103 (MQ=255)
                                                                                                    tAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGtaa                                                   <  1:578924/103‑1 (MQ=255)
                                                                                                      aTGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCa                                                                                                             >  1:506258/1‑43 (MQ=255)
                                                                                                      aTGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCa                                                                                                             <  2:506258/43‑1 (MQ=255)
                                                                                                      aTGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAg                                                                                              >  2:757148/1‑58 (MQ=255)
                                                                                                      aTGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAg                                                                                              <  1:757148/58‑1 (MQ=255)
                                                                                                         aGCAGTGCAGATAAAGCTGATCATATCGATGGGCAACTCATGCAATTCTTTTGAGCAATACACACGCGCTTCCAGCGGAGtata                                                                 <  2:1197838/84‑1 (MQ=255)
                                                                                                         aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTAtttt                                                                                                 >  2:1271008/1‑52 (MQ=255)
                                                                                                         aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTAtttt                                                                                                 <  1:1271008/52‑1 (MQ=255)
                                                                                                         aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGtata                                                                 >  1:1197838/1‑84 (MQ=255)
                                                                                                         aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTg                 >  2:485806/1‑132 (MQ=255)
                                                                                                         aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTg                 <  1:485806/132‑1 (MQ=255)
                                                                                                         aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCt          >  2:16633/1‑139 (MQ=255)
                                                                                                           cAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCa                                       <  2:321225/108‑1 (MQ=255)
                                                                                                           cAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCa                                       >  1:321225/1‑108 (MQ=255)
                                                                                                           cAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGt        >  1:539918/1‑139 (MQ=255)
                                                                                                              tGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                          <  2:496681/70‑1 (MQ=255)
                                                                                                              tGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                          <  2:18443/70‑1 (MQ=255)
                                                                                                              tGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                          >  1:496681/1‑70 (MQ=255)
                                                                                                              tGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                          >  1:18443/1‑70 (MQ=255)
                                                                                                                cAGATAAAGCTGACCATATCGATGGGCAACTCATGCaa                                                                                                        >  1:994346/1‑38 (MQ=255)
                                                                                                                cAGATAAAGCTGACCATATCGATGGGCAACTCATGCaa                                                                                                        <  2:994346/38‑1 (MQ=255)
                                                                                                                 aGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTaac  <  1:16633/139‑1 (MQ=255)
                                                                                                                   aTAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTaa                                                        <  2:261091/83‑1 (MQ=255)
                                                                                                                   aTAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTaa                                                        >  1:261091/1‑83 (MQ=255)
                                                                                                                   aTAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCa                                       <  1:213159/100‑1 (MQ=255)
                                                                                                                      |                                                                                                                                     
CCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAAC  >  NC_000913/578450‑578701

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