Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	578,102	C→T	Q164* (CAA→TAA)	rrrD →	DLP12 prophage; putative lysozyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	578,102	0	C	T	100.0%	200.6 / NA	54	Q164* (CAA→TAA)	rrrD	DLP12 prophage; putative lysozyme
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/28);  total (26/28)

GAACATCGACGCTTCTTCGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAA  >  NC_000913/577956‑578249                                                                                                                                                   |                                                                                                                                                    gAACATCGACGCTTCTTCGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACa                                                                                                                                                 <  1:305803/151‑1 (MQ=255) gAACATCGACGCTTCTTCGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACa                                                                                                                                                 <  1:3704/151‑1 (MQ=255)    cATCGACGCTTCTTCGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGAGGGGGTAACAGTg                                                                                                                                              >  2:128620/1‑151 (MQ=255)     aTCGACGCTTCTTCGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTg                                                                                                                                              >  2:237409/1‑150 (MQ=255)       cGACGCTTCTTCGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGc                                                                                                                                           >  2:506486/1‑151 (MQ=255)          cGCTTCTTCGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGc                                                                                                                                           >  1:379099/1‑148 (MQ=255)          cGCTTCTTCGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGc                                                                                                                                           >  2:450217/1‑148 (MQ=255)              tcttcGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGt                                                                                                                                      <  2:69220/149‑1 (MQ=255)              tcttcGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTa                                                                                                                                     <  1:215751/150‑1 (MQ=255)              tcttcGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTa                                                                                                                                     >  2:45214/1‑150 (MQ=255)               cttcGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAAc                                                                                                                                   >  1:307242/1‑151 (MQ=255)                   gCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACcgcg                                                                                                                               <  1:428229/151‑1 (MQ=255)                     aaaaTAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGAt                                                                                                                             >  1:150773/1‑151 (MQ=255)                      aaaTAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGAtt                                                                                                                            >  2:599811/1‑151 (MQ=255)                      aaaTAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGAtt                                                                                                                            <  2:341785/151‑1 (MQ=255)                       aaTAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGAt                                                                                                                             <  1:468513/149‑1 (MQ=255)                               aGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGct                                                                                                                   >  2:235817/1‑151 (MQ=255)                                   cGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGa                                                                                                               <  2:253543/151‑1 (MQ=255)                                    gATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGAt                                                                                                              <  1:632162/151‑1 (MQ=255)                                         cAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCt                                                                                                         <  1:625732/151‑1 (MQ=255)                                          aaaGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTg                                                                                                        <  2:529743/151‑1 (MQ=255)                                            aGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGc                                                                                                       >  1:401957/1‑150 (MQ=255)                                             ggCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGcat                                                                                                     <  2:379092/151‑1 (MQ=255)                                                      gATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCa                                                                                             <  2:150769/150‑1 (MQ=255)                                                         cAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCt                                                                                         <  2:50012/151‑1 (MQ=255)                                                              gCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCAt                                                                                    <  2:585229/151‑1 (MQ=255)                                                              gCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCAt                                                                                    <  2:146787/151‑1 (MQ=255)                                                                gTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATgg                                                                                  >  2:563066/1‑151 (MQ=255)                                                                          aTACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGCAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCa                                                                           >  2:443244/1‑148 (MQ=255)                                                                           tACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATc                                                                        >  1:551085/1‑150 (MQ=255)                                                                                    ggTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTg                                                               >  2:249867/1‑150 (MQ=255)                                                                                      tAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGAt                                                             >  1:601390/1‑150 (MQ=255)                                                                                         gCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACg                                                         >  2:319674/1‑151 (MQ=255)                                                                                           aaTGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGc                                                        >  2:109746/1‑150 (MQ=255)                                                                                            aTGGAAAGTCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCaa                                                      <  2:287236/151‑1 (MQ=255)                                                                                                  aGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCt                                                <  2:345559/151‑1 (MQ=255)                                                                                                    gCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTAc                                              >  2:63963/1‑151 (MQ=255)                                                                                                     ccTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACa                                             <  1:254360/151‑1 (MQ=255)                                                                                                        gatgaCTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACaaa                                           <  1:599826/150‑1 (MQ=255)                                                                                                           gaCTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTcc                                       >  1:95566/1‑151 (MQ=255)                                                                                                            aCTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTcc                                       >  2:469955/1‑150 (MQ=255)                                                                                                               cgtcgtGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAgcg                                   >  2:585783/1‑151 (MQ=255)                                                                                                                  cgtGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGcgcg                                 <  1:83109/150‑1 (MQ=255)                                                                                                                     gagaTTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACa                              <  2:95566/150‑1 (MQ=255)                                                                                                                             gCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCa                     <  1:608830/151‑1 (MQ=255)                                                                                                                               gTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGcc                      >  2:61314/1‑148 (MQ=255)                                                                                                                                        gTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGaa           <  2:109637/150‑1 (MQ=255)                                                                                                                                        gTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGa            <  1:425120/149‑1 (MQ=255)                                                                                                                                        gTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGa            >  1:204587/1‑149 (MQ=255)                                                                                                                                         tttGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGc         >  2:42224/1‑151 (MQ=255)                                                                                                                                          ttGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGc         <  2:269246/150‑1 (MQ=255)                                                                                                                                           tgtgGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCt        >  2:320648/1‑150 (MQ=255)                                                                                                                                            gtgGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTa       <  1:523544/150‑1 (MQ=255)                                                                                                                                                gggTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGaa  <  1:350933/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GAACATCGACGCTTCTTCGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAA  >  NC_000913/577956‑578249

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