Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,412,000	C→T	*412* (TAG→TAA)	intR ←	Rac prophage, integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,412,000	0	C	T	100.0%	229.4 / NA	70	*412* (TAG→TAA)	intR	Rac prophage, integrase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (40/30);  total (40/30)

AGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATGCCGGGTTCGCTC  >  NC_000913/1411884‑1412130                                                                                                                     |                                                                                                                                   aGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACa                                                                                                              >  2:588448/1‑139 (MQ=255) aGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACa                                                                                                              >  2:1608977/1‑139 (MQ=255)             aaTTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGtt                                                                                                  <  1:588448/139‑1 (MQ=255)                                 ttctttctTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACcg                                                                                                                         >  1:1493382/1‑96 (MQ=255)                                 ttctttctTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACcg                                                                                                                         <  2:1493382/96‑1 (MQ=255)                                 ttctttctTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTtcg                                                                                                <  1:1003980/121‑1 (MQ=255)                                 ttctttctTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTtcg                                                                                                >  2:1003980/1‑121 (MQ=255)                                  tctttcttTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATg                                                                                       <  1:1399497/129‑1 (MQ=255)                                  tctttcttTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATg                                                                                       >  2:1399497/1‑129 (MQ=255)                                    tttctttGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCcac                                                                           >  2:952831/1‑139 (MQ=255)                                      tctttGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACAt                                                                         >  1:1154669/1‑139 (MQ=255)                                      tctttGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACAt                                                                         >  1:993414/1‑139 (MQ=255)                                          tGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCg                                                                     >  1:1245035/1‑139 (MQ=255)                                             aaTTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTaa                                                                  <  1:1439919/139‑1 (MQ=255)                                               tttGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTaa                                                                  >  1:1586759/1‑137 (MQ=255)                                               tttGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTaa                                                                  <  2:1586759/137‑1 (MQ=255)                                               tttGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAAt                                                                <  2:993414/139‑1 (MQ=255)                                                           cTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCAt                                                    <  1:917378/139‑1 (MQ=255)                                                             tGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCt                                                  >  1:543323/1‑139 (MQ=255)                                                                attattTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTgcg                                               <  2:1154669/139‑1 (MQ=255)                                                                    tttCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAAc                                                                                                               <  2:215375/71‑1 (MQ=255)                                                                    tttCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAAc                                                                                                               >  1:215375/1‑71 (MQ=255)                                                                         gTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACcg                                                                                                                         >  1:1159998/1‑56 (MQ=255)                                                                         gTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACcg                                                                                                                         <  2:1159998/56‑1 (MQ=255)                                                                          ttCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCAt                                                    <  1:303689/124‑1 (MQ=255)                                                                          ttCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCAt                                                    >  2:303689/1‑124 (MQ=255)                                                                          ttCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATg                                     <  1:1608977/139‑1 (MQ=255)                                                                              ctGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTAt                                                                                                       <  2:938638/69‑1 (MQ=255)                                                                              ctGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTAt                                                                                                       >  1:938638/1‑69 (MQ=255)                                                                              ctGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGtt                                                                                           <  1:334755/81‑1 (MQ=255)                                                                              ctGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGtt                                                                                           >  2:334755/1‑81 (MQ=255)                                                                               tGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCa                                >  1:1358483/1‑139 (MQ=255)                                                                               tGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCa                                >  2:1570626/1‑139 (MQ=255)                                                                               tGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCa                                >  2:1337918/1‑139 (MQ=255)                                                                                  tACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCAttt                             >  1:700320/1‑139 (MQ=255)                                                                                    cTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGtt                                                                                                  <  2:175667/68‑1 (MQ=255)                                                                                    cTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGtt                                                                                                  >  1:175667/1‑68 (MQ=255)                                                                                    cTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGTCCCATTtgg                           <  2:543323/139‑1 (MQ=255)                                                                                      aaaTGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtggct                         >  1:1648847/1‑139 (MQ=255)                                                                                       aaTGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCAtt                                                                                                                   <  2:634028/48‑1 (MQ=255)                                                                                       aaTGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCAtt                                                                                                                   >  1:634028/1‑48 (MQ=255)                                                                                       aaTGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtggctg                        >  1:847142/1‑139 (MQ=255)                                                                                       aaTGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtggctg                        >  1:646391/1‑139 (MQ=255)                                                                                         tGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtcatca                                                                                 >  1:881824/1‑80 (MQ=255)                                                                                         tGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtcatca                                                                                 <  2:881824/80‑1 (MQ=255)                                                                                           gggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCta                    >  1:52208/1‑139 (MQ=255)                                                                                            ggCAAATTGGGGGCAAACTTTGCAATTACGATAACcg                                                                                                                         <  2:1331240/37‑1 (MQ=255)                                                                                            ggCAAATTGGGGGCAAACTTTGCAATTACGATAACcg                                                                                                                         >  1:1331240/1‑37 (MQ=255)                                                                                            ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtcatc                                                                                  <  1:566060/76‑1 (MQ=255)                                                                                            ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtcatc                                                                                  >  2:566060/1‑76 (MQ=255)                                                                                            ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtcat                                                                                   >  2:737635/1‑75 (MQ=255)                                                                                            ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtcat                                                                                   <  1:737635/75‑1 (MQ=255)                                                                                            ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtggctgg                       <  1:44836/135‑1 (MQ=255)                                                                                            ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtggctgg                       >  2:44836/1‑135 (MQ=255)                                                                                            ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtgg                           <  2:1454676/131‑1 (MQ=255)                                                                                            ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtgg                           >  1:1454676/1‑131 (MQ=255)                                                                                            ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCtat                   >  2:1170150/1‑139 (MQ=255)                                                                                            ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCtat                   >  1:1655928/1‑139 (MQ=255)                                                                                            ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCtat                   >  2:512811/1‑139 (MQ=255)                                                                                             gCAAATTGGGGGCAAACTTTGCAATTACGATAACCgcgc                                                                                                                      >  1:379364/1‑39 (MQ=255)                                                                                             gCAAATTGGGGGCAAACTTTGCAATTACGATAACCgcgc                                                                                                                      <  2:379364/39‑1 (MQ=255)                                                                                                 aTTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtca                                                                                    >  2:446906/1‑69 (MQ=255)                                                                                                 aTTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtca                                                                                    <  1:446906/69‑1 (MQ=255)                                                                                                  ttGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATGc             <  2:1655928/139‑1 (MQ=255)                                                                                                      gggCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtgg                           <  2:819264/121‑1 (MQ=255)                                                                                                      gggCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtgg                           >  1:819264/1‑121 (MQ=255)                                                                                                           aaCTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtca                                                                                    >  2:1317061/1‑59 (MQ=255)                                                                                                           aaCTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtca                                                                                    <  1:1317061/59‑1 (MQ=255)                                                                                                               ttGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATGCCGGGTTCGCTc  >  1:960313/1‑137 (MQ=255)                                                                                                               ttGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATGCCGGGTTCGCTc  <  2:960313/137‑1 (MQ=255)                                                                                                                     |                                                                                                                                   AGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATGCCGGGTTCGCTC  >  NC_000913/1411884‑1412130

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