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%	206.2 / NA	63	*154* (TAG→TAA)	rzpD	DLP12 prophage, putative murein endopeptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/30);  total (33/30)

TAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAA  >  NC_000913/578442‑578700                                                                                                                               |                                                                                                                                     tAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCatat                                                                                                                          <  2:660378/139‑1 (MQ=255)              ccccGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTc                                                                                                             >  2:137466/1‑139 (MQ=255)                 cGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATg                                                                                                          <  1:1151431/139‑1 (MQ=255)                       gCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAatta                                                                                                    >  2:109301/1‑139 (MQ=255)                               cGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGc                                                                                            >  2:1429701/1‑139 (MQ=255)                                      cGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATacac                                                                                     >  1:132138/1‑139 (MQ=255)                                      cGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATacac                                                                                     >  2:195136/1‑139 (MQ=255)                                       gggATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATacaca                                                                                    <  1:137466/139‑1 (MQ=255)                                          attattTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACAcgc                                                                                 <  1:1429701/139‑1 (MQ=255)                                                  aCCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                         <  1:109301/139‑1 (MQ=255)                                                   cccTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGc                                                                        <  1:195136/139‑1 (MQ=255)                                                   cccTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGc                                                                        >  2:778797/1‑139 (MQ=255)                                                           gagagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGtata                                                                <  2:132138/139‑1 (MQ=255)                                                                gCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTAtttt                                                                                                <  2:173018/102‑1 (MQ=255)                                                                gCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTAtttt                                                                                                >  1:173018/1‑102 (MQ=255)                                                                  tGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                         >  2:247466/1‑123 (MQ=255)                                                                  tGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                         <  1:247466/123‑1 (MQ=255)                                                                    aTCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGa                                                                     >  2:467867/1‑125 (MQ=255)                                                                    aTCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGa                                                                     <  1:467867/125‑1 (MQ=255)                                                                           tGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGtaata                                                >  1:707152/1‑139 (MQ=255)                                                                           tGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGtaata                                                >  1:1369098/1‑139 (MQ=255)                                                                              aaaaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATaaaa                                             >  1:571378/1‑139 (MQ=255)                                                                                aaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGGCCATATCGATGGGCAACTCGTGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                         <  2:1114654/109‑1 (MQ=255)                                                                                aaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGGCCATATCGATGGGCAACTCGTGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                         >  1:1114654/1‑109 (MQ=255)                                                                                aaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGtata                                                                >  1:669844/1‑118 (MQ=255)                                                                                aaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGtata                                                                <  2:669844/118‑1 (MQ=255)                                                                                   caacTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTaa                                                       >  1:1224401/1‑124 (MQ=255)                                                                                   caacTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTaa                                                       <  2:1224401/124‑1 (MQ=255)                                                                                   caacTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAg                                        >  1:816738/1‑139 (MQ=255)                                                                                     acTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAAc                                                                                                               <  1:1221489/66‑1 (MQ=255)                                                                                     acTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAAc                                                                                                               >  2:1221489/1‑66 (MQ=255)                                                                                      cTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCaa                                     >  1:1357545/1‑139 (MQ=255)                                                                                       tGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAAt                                    <  2:816738/139‑1 (MQ=255)                                                                                                 ccAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGa                          >  2:506383/1‑139 (MQ=255)                                                                                                 ccAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGa                          <  1:778797/139‑1 (MQ=255)                                                                                                        tataTTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTg                   >  1:59288/1‑139 (MQ=255)                                                                                                        tataTTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTg                   <  2:707152/139‑1 (MQ=255)                                                                                                          tatTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAg                                                                                             >  2:1422917/1‑63 (MQ=255)                                                                                                          tatTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAg                                                                                             <  1:1422917/63‑1 (MQ=255)                                                                                                          tatTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCaa                                                                                          >  1:160390/1‑66 (MQ=255)                                                                                                          tatTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCaa                                                                                          <  2:160390/66‑1 (MQ=255)                                                                                                          tatTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCa                                 >  1:103114/1‑123 (MQ=255)                                                                                                          tatTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCt                 >  1:873297/1‑139 (MQ=255)                                                                                                          tatTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCt                 <  2:59288/139‑1 (MQ=255)                                                                                                          tatTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGACTAAAGTAATAAAACCGAGCAATCCa                                 <  2:103114/123‑1 (MQ=255)                                                                                                                  aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCaa                                                                                                       <  2:1380649/45‑1 (MQ=255)                                                                                                                  aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCaa                                                                                                       >  1:1380649/1‑45 (MQ=255)                                                                                                                  aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAattatt                                                                                                  >  1:879494/1‑50 (MQ=255)                                                                                                                  aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAattatt                                                                                                  <  2:879494/50‑1 (MQ=255)                                                                                                                  aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAg                                                                                             >  2:948853/1‑55 (MQ=255)                                                                                                                  aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAg                                                                                             <  1:948853/55‑1 (MQ=255)                                                                                                                  aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAAt                                    <  2:479944/112‑1 (MQ=255)                                                                                                                  aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAAt                                    >  1:479944/1‑112 (MQ=255)                                                                                                                    cAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAg                                                                                             <  2:127105/53‑1 (MQ=255)                                                                                                                    cAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAg                                                                                             >  1:127105/1‑53 (MQ=255)                                                                                                                    cAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATaaaa                                             >  2:533952/1‑101 (MQ=255)                                                                                                                    cAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATaaaa                                             <  1:533952/101‑1 (MQ=255)                                                                                                                      gTGCAGATAAAGGTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTaa                                                       <  2:1023328/89‑1 (MQ=255)                                                                                                                      gTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTaa                                                       >  1:1023328/1‑89 (MQ=255)                                                                                                                       tGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGtttt    >  2:72775/1‑139 (MQ=255)                                                                                                                       tGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGtttt    >  2:933722/1‑139 (MQ=255)                                                                                                                        gCAGATAAAGCTGACCATATCGATTGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTa   <  1:506383/139‑1 (MQ=255)                                                                                                                         cAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTaa  <  2:873297/139‑1 (MQ=255)                                                                                                                               |                                                                                                                                     TAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAA  >  NC_000913/578442‑578700

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