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_000913	578,568	0	G	A	100.0%	216.6 / NA	67	*154* (TAG→TAA)	rzpD	DLP12 prophage, putative murein endopeptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (35/32);  total (35/32)

ATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAA  >  NC_000913/578444‑578703                                                                                                                             |                                                                                                                                        aTGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCg                                                                                                                           <  1:1862020/139‑1 (MQ=255)       ccTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGc                                                                                                                     >  2:648759/1‑139 (MQ=255)           cccccGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACt                                                                                                                 <  2:149211/139‑1 (MQ=255)                           acCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGa                                                                                                 >  1:586369/1‑139 (MQ=255)                             cGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGc                                                                                               >  2:449446/1‑139 (MQ=255)                               cTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAACCAACTGGAAGGAACCCAGAAGAATATTAATGAGCAGTGCAGATAAAGc                                                                                                                                      >  1:99156/1‑98 (MQ=255)                               cTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGc                                                                                                                                      <  2:99156/98‑1 (MQ=255)                                       gATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACAcg                                                                                     >  1:743971/1‑139 (MQ=255)                                            tttCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACAcgc                                                                                    >  2:1792834/1‑135 (MQ=255)                                            tttCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACAcgc                                                                                    <  1:1792834/135‑1 (MQ=255)                                               cACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCa                                                                             <  2:743971/139‑1 (MQ=255)                                                  ccTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCg                                                                          <  2:586369/139‑1 (MQ=255)                                                      agagagagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGATGTATATTAATGAGCAGTGCAGATAAAGCt                                                                                                                                     <  1:1188554/76‑1 (MQ=255)                                                      agagagagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGCGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTAtttt                                                                                                   <  2:959014/110‑1 (MQ=255)                                                      agagagagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCt                                                                                                                                     >  2:1188554/1‑76 (MQ=255)                                                      agagagagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTCTTTTGAGCAATACACACGCGCTTc                                                                               <  1:1391447/130‑1 (MQ=255)                                                      agagagagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTAtttt                                                                                                   >  1:959014/1‑110 (MQ=255)                                                      agagagagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATAATAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTAGAGCAATACACACGCGCTTc                                                                               >  2:1391447/1‑130 (MQ=255)                                                            agGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                            >  2:1900302/1‑127 (MQ=255)                                                            agGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                            <  1:1900302/127‑1 (MQ=255)                                                                  aTCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCg                                                                          <  2:285096/123‑1 (MQ=255)                                                                  aTCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCg                                                                          >  1:285096/1‑123 (MQ=255)                                                                         tGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                            <  1:321566/114‑1 (MQ=255)                                                                         tGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                            >  2:321566/1‑114 (MQ=255)                                                                            aaaaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGAcc                                                                                                                                 >  2:1639023/1‑58 (MQ=255)                                                                            aaaaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGAcc                                                                                                                                 <  1:1639023/58‑1 (MQ=255)                                                                              aaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCa                                                                                                                                <  1:388511/57‑1 (MQ=255)                                                                              aaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCa                                                                                                                                >  2:388511/1‑57 (MQ=255)                                                                              aaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGc                                                                                                                     <  2:948787/68‑1 (MQ=255)                                                                              aaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGc                                                                                                                     >  1:948787/1‑68 (MQ=255)                                                                              aaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCa                                                                                                               >  1:216364/1‑74 (MQ=255)                                                                              aaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCa                                                                                                               <  2:216364/74‑1 (MQ=255)                                                                              aaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAAcc                                              >  2:233520/1‑139 (MQ=255)                                                                                 caacTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAattatt                                                                                                     >  2:1523137/1‑81 (MQ=255)                                                                                 caacTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAattatt                                                                                                     <  1:1523137/81‑1 (MQ=255)                                                                                   acTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                            <  1:1046537/104‑1 (MQ=255)                                                                                   acTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                            >  2:1046537/1‑104 (MQ=255)                                                                                   acTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGtata                                                                   <  2:1839217/113‑1 (MQ=255)                                                                                   acTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGtata                                                                   >  1:1839217/1‑113 (MQ=255)                                                                                    cTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCaa                                        >  2:1932245/1‑139 (MQ=255)                                                                                              cccAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATACCTAAAGTAATAAAACCGAGCAATCCATTTACg                              >  1:1570572/1‑139 (MQ=255)                                                                                                 aGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCg                                                                          >  2:977418/1‑92 (MQ=255)                                                                                                 aGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCg                                                                          <  1:977418/92‑1 (MQ=255)                                                                                                      tataTTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTg                      <  1:233520/139‑1 (MQ=255)                                                                                                        tatTACTGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCGAGCCGAGTATAAATGCCTAAAGTAATAATACCGAGCAATCCATTTACGAATGTTTGCt                    <  1:449446/139‑1 (MQ=255)                                                                                                        tatTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAg                                           <  1:11836/116‑1 (MQ=255)                                                                                                        tatTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAg                                           >  2:11836/1‑116 (MQ=255)                                                                                                             aTGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTaa                                                          >  1:1219581/1‑96 (MQ=255)                                                                                                             aTGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTaa                                                          <  2:1219581/96‑1 (MQ=255)                                                                                                                aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGCGCAATACACTAGCGCTTCCAGCGGAGTATAAATGTCTTAAGTAATAAAACCGAGCAAACGATTTACGAATGATTGCTGCGTTTCt            >  1:562669/1‑139 (MQ=255)                                                                                                                aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGcc                                                             <  2:1318598/90‑1 (MQ=255)                                                                                                                aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGcc                                                             >  1:1318598/1‑90 (MQ=255)                                                                                                                aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAg                                           >  1:1435697/1‑108 (MQ=255)                                                                                                                aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAg                                           <  2:1435697/108‑1 (MQ=255)                                                                                                                aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGaa                            <  1:1959886/123‑1 (MQ=255)                                                                                                                aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGaa                            >  2:1959886/1‑123 (MQ=255)                                                                                                                aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCt            <  1:1932245/139‑1 (MQ=255)                                                                                                                aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCt            >  1:1421042/1‑139 (MQ=255)                                                                                                                  cAGTGCAGATAAAGCTGACCATATGGATGGGCAACTCATGCAATTATTTTGAGCAATACACAcgc                                                                                    >  1:1394823/1‑65 (MQ=255)                                                                                                                  cAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATacacac                                                                                      >  1:1653614/1‑63 (MQ=255)                                                                                                                  cAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATacacac                                                                                      <  2:1653614/63‑1 (MQ=255)                                                                                                                  cAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGt          >  1:1783853/1‑139 (MQ=255)                                                                                                                  cAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATGATTTTGAGCAATACACAcgc                                                                                    <  2:1394823/65‑1 (MQ=255)                                                                                                                     tGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCg                                                                          >  1:1266380/1‑72 (MQ=255)                                                                                                                     tGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCg                                                                          <  2:1266380/72‑1 (MQ=255)                                                                                                                      gCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTa      <  2:1783853/139‑1 (MQ=255)                                                                                                                          aTAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTaacaa  >  1:1333502/1‑139 (MQ=255)                                                                                                                             |                                                                                                                                        ATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAA  >  NC_000913/578444‑578703

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: ‑