Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,432,273	G→A	A1075A (GCG→GCA)	stfR →	Rac prophage, putative tail fiber protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,432,273	0	G	A	100.0%	207.1 / NA	68	A1075A (GCG→GCA)	stfR	Rac prophage, putative tail fiber protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/35);  total (33/35)

CTAGTGCTAACTCCGGTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATTGTGAGGCTT  >  NC_000913/1432140‑1432405                                                                                                                                      |                                                                                                                                     cTAGTGCTAACTCCGGTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATAc                                                                                                                                 >  1:1409421/1‑139 (MQ=255)  tAGTGCTAACTCCGGTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACt                                                                                                                                <  2:172894/139‑1 (MQ=255)  tAGTGCTAACTCCGGTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACt                                                                                                                                <  2:526390/139‑1 (MQ=255)               ggTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTg                                                                                                                   >  2:196409/1‑139 (MQ=255)                 tGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCt                                                                                                                 >  1:968217/1‑139 (MQ=255)                 tGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCt                                                                                                                 <  2:245631/139‑1 (MQ=255)                  gCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTc                                                                                                                >  1:1248600/1‑139 (MQ=255)                      gTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTc                                                                                                                             >  2:1042904/1‑122 (MQ=255)                      gTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTc                                                                                                                             <  1:1042904/122‑1 (MQ=255)                        aGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGc                                                                                                          <  1:1112715/139‑1 (MQ=255)                           gCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACt                                                                                                       <  1:196409/139‑1 (MQ=255)                            cATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCg                                                                                                                            <  1:138632/117‑1 (MQ=255)                            cATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCg                                                                                                                            >  2:138632/1‑117 (MQ=255)                            cATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTc                                                                                                      >  2:431132/1‑139 (MQ=255)                              tCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGTCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCg                                                                                                    >  2:879618/1‑139 (MQ=255)                                        tGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTgg                                                                                          <  2:9630/139‑1 (MQ=255)                                             gttgtGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTcac                                                                                     <  1:431132/139‑1 (MQ=255)                                               tgtgCATAATCAAAACTATGCAACATCGTCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACAt                                                                                   >  1:1337370/1‑139 (MQ=255)                                                  gCATAATCAAAACTATGTAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTTTTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATggc                                                                                <  1:879618/139‑2 (MQ=255)                                                  gCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGa                                                                                <  2:647154/139‑1 (MQ=255)                                                   cATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGac                                                                               >  2:137581/1‑139 (MQ=255)                                                   cATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGac                                                                               >  2:1224791/1‑139 (MQ=255)                                                        tCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATa                                                                                                                                  >  1:525450/1‑83 (MQ=255)                                                        tCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATa                                                                                                                                  <  2:525450/83‑1 (MQ=255)                                                        tCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACAcc                                                                          >  1:668923/1‑139 (MQ=255)                                                         cAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCa                                                                         <  1:246101/139‑1 (MQ=255)                                                          aaaaCTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCAt                                                                        <  1:1224791/139‑1 (MQ=255)                                                           aaaCTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATc                                                                       >  2:1169219/1‑139 (MQ=255)                                                             aCTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTACGATTGGTTCACATGGACACACCATCAc                                                                     >  1:857893/1‑139 (MQ=255)                                                                aTGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGt                                                                  >  2:733239/1‑139 (MQ=255)                                                                  gCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTa                                                                >  2:407403/1‑139 (MQ=255)                                                                     aCATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACg                                                             <  2:1409421/139‑1 (MQ=255)                                                                         catcTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACgctgc                                                         >  2:894890/1‑139 (MQ=255)                                                                         catcTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACgctgc                                                         >  1:273949/1‑139 (MQ=255)                                                                         catcTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACgctgc                                                         >  2:497453/1‑139 (MQ=255)                                                                         catcTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGCGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACgctgc                                                         <  2:215648/139‑1 (MQ=255)                                                                           tcTGCTGGCGCACATACCCGCTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGtt                                                                                        <  2:1210467/106‑1 (MQ=255)                                                                           tcTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGtt                                                                                        >  1:1210467/1‑106 (MQ=255)                                                                              gctgGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTa                                                    <  1:1177934/139‑1 (MQ=255)                                                                                 ggCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCTGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTaagg                                                 <  2:857893/139‑3 (MQ=255)                                                                                  gcgcACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAAcgc                                                >  2:423827/1‑139 (MQ=255)                                                                                          aCCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAaca                                        <  2:548098/139‑1 (MQ=255)                                                                                          aCCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAaca                                        <  1:733239/139‑1 (MQ=255)                                                                                                    tGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAaca                                        <  1:275824/129‑1 (MQ=25)                                                                                                    tGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAaca                                        >  2:275824/1‑129 (MQ=25)                                                                                                    tGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCaaaaa                              >  2:972074/1‑139 (MQ=25)                                                                                                      tCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACa                            >  1:899408/1‑139 (MQ=21)                                                                                                      tCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACa                            >  1:1270929/1‑139 (MQ=37)                                                                                                      tCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACa                            >  2:1379117/1‑139 (MQ=25)                                                                                                       ccGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACAt                           <  2:968217/139‑1 (MQ=25)                                                                                                        cGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACgc                                                            >  1:1115847/1‑105 (MQ=255)                                                                                                        cGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACgc                                                            <  2:1115847/105‑1 (MQ=25)                                                                                                        cGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATc                          >  1:1190236/1‑139 (MQ=21)                                                                                                         ggCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGACTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCg                         <  2:362135/139‑1 (MQ=255)                                                                                                          gCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATAc                                                                                                            <  2:113975/55‑1 (MQ=255)                                                                                                          gCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATAc                                                                                                            >  1:113975/1‑55 (MQ=255)                                                                                                           cACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCa                       <  1:894890/139‑1 (MQ=18)                                                                                                           cACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCa                       <  1:407403/139‑1 (MQ=18)                                                                                                           cACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCa                       <  2:899408/139‑1 (MQ=18)                                                                                                             ctgctgCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCAtt                     <  1:1379117/139‑1 (MQ=37)                                                                                                                ctgcAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGa                                                                                >  2:569828/1‑77 (MQ=255)                                                                                                                ctgcAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGa                                                                                <  1:569828/77‑1 (MQ=255)                                                                                                                    aaGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACtat              >  2:643054/1‑139 (MQ=25)                                                                                                                         cAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATtgt         <  2:390004/139‑1 (MQ=37)                                                                                                                          aGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATtgtg        <  1:198083/139‑1 (MQ=21)                                                                                                                              gCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATTGTGAGGc    <  2:668923/139‑1 (MQ=17)                                                                                                                              gCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATTGTGAGGc    <  1:1169219/139‑1 (MQ=14)                                                                                                                                acacGCACATACTGTCGGTTTTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATTGTGAGGCtt  >  2:153294/1‑139 (MQ=0)                                                                                                                                      |                                                                                                                                     CTAGTGCTAACTCCGGTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATTGTGAGGCTT  >  NC_000913/1432140‑1432405

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