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%	212.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 (32/36);  total (32/36)

TAGTGCTAACTCCGGTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATTGT  >  NC_000913/1432141‑1432398                                                                                                                                     |                                                                                                                              tAGTGCTAACTCCGGTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACt                                                                                                                         <  1:550064/139‑1 (MQ=255)        aaCTCCGGTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTa                                                                                                                  >  2:113159/1‑139 (MQ=255)                     gTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATAc                                                                                                     >  2:1412028/1‑139 (MQ=255)                             tCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCg                                                                                             >  2:1331272/1‑139 (MQ=255)                                acaaGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTg                                                                                          >  1:451936/1‑139 (MQ=255)                                       tGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTgg                                                                                   <  1:785795/139‑1 (MQ=255)                                       tGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTgg                                                                                   <  2:1322759/139‑1 (MQ=255)                                         tCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGtt                                                                                 >  2:1356865/1‑139 (MQ=255)                                            gttgtGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTcac                                                                              <  2:140073/139‑1 (MQ=255)                                             ttgtGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTcaca                                                                             <  2:1092131/139‑1 (MQ=255)                                              tgtgCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACAt                                                                            >  2:957334/1‑139 (MQ=255)                                                 gCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGa                                                                         <  2:164650/139‑1 (MQ=255)                                                       tCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATa                                                                                                                           >  1:136007/1‑83 (MQ=255)                                                       tCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATa                                                                                                                           <  2:136007/83‑1 (MQ=255)                                                         aaaaCTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCAt                                                                 <  2:3050/139‑1 (MQ=255)                                                          aaaCTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACatcatc                                                                <  1:957334/139‑1 (MQ=255)                                                          aaaCTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATc                                                                <  1:113159/139‑1 (MQ=255)                                                            aCTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCAc                                                              >  1:170220/1‑139 (MQ=255)                                                               aTGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCt                                                                                                          >  1:186872/1‑92 (MQ=255)                                                               aTGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCt                                                                                                          <  2:186872/92‑1 (MQ=255)                                                               aTGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGt                                                           >  1:308229/1‑139 (MQ=255)                                                                       tcatcTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACgctg                                                   >  2:1427491/1‑139 (MQ=255)                                                                        catcTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACgctgc                                                  >  1:913574/1‑139 (MQ=255)                                                                          tcTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGctgctg                                                >  2:1347469/1‑139 (MQ=255)                                                                                 gcgcACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATa                                                                                                                           >  2:163931/1‑57 (MQ=255)                                                                                 gcgcACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATa                                                                                                                           <  1:163931/57‑1 (MQ=255)                                                                                 gcgcACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAAcgc                                         >  1:262043/1‑139 (MQ=255)                                                                                    cacaTACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGa                                      >  2:503978/1‑139 (MQ=255)                                                                                         aCCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAaca                                 <  1:1427491/139‑1 (MQ=255)                                                                                         aCCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAaca                                 <  1:1412028/139‑1 (MQ=255)                                                                                                cactGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCaa                          >  1:1132902/1‑139 (MQ=37)                                                                                                 actGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTgg                                                                                   >  2:1098859/1‑81 (MQ=255)                                                                                                 actGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTgg                                                                                   <  1:1098859/81‑1 (MQ=255)                                                                                                 actGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCaaa                         >  2:967328/1‑139 (MQ=25)                                                                                                   tGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTcac                                                                              <  2:765841/84‑1 (MQ=255)                                                                                                   tGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTcac                                                                              >  1:765841/1‑84 (MQ=255)                                                                                                   tGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAAc                                           <  1:1307594/119‑1 (MQ=25)                                                                                                   tGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAAc                                           >  2:1307594/1‑119 (MQ=25)                                                                                                   tGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTACGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCaaaaa                       >  1:842281/1‑139 (MQ=21)                                                                                                    gTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAAc                      >  1:173528/1‑139 (MQ=25)                                                                                                     tCCGGTACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACa                     >  2:346280/1‑139 (MQ=19)                                                                                                     tCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACa                     >  2:1288685/1‑139 (MQ=25)                                                                                                     tCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACa                     <  2:262043/139‑1 (MQ=25)                                                                                                      ccGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACAt                    <  2:1443864/139‑1 (MQ=37)                                                                                                      ccGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACAt                    <  2:173528/139‑1 (MQ=25)                                                                                                       cGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTCGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATc                   >  2:1291827/1‑139 (MQ=255)                                                                                                         gCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTgg                                                                                   >  2:988917/1‑73 (MQ=255)                                                                                                         gCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTgg                                                                                   <  1:988917/73‑1 (MQ=255)                                                                                                          cACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCa                <  2:1132902/139‑1 (MQ=37)                                                                                                          cACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCa                <  2:842281/139‑1 (MQ=21)                                                                                                          cACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCa                <  2:451936/139‑1 (MQ=21)                                                                                                               ctgcAAGCGCAGGTGTACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACTGTCAAAAACATCGCa                <  1:984359/134‑1 (MQ=255)                                                                                                               ctgcAAGCGCAGGTGTACACGCACATACTGTCGGTATAGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTTCTGGTAACGCGGAAAACACCGTCAAAAACATCGCa                >  2:984359/1‑134 (MQ=255)                                                                                                               ctgcAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCa                <  1:44206/134‑1 (MQ=37)                                                                                                               ctgcAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCa                >  2:44206/1‑134 (MQ=25)                                                                                                               ctgcAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCAtt              <  1:1193342/136‑1 (MQ=25)                                                                                                               ctgcAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCAtt              >  2:915163/1‑136 (MQ=25)                                                                                                               ctgcAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCAtt              <  1:915163/136‑1 (MQ=25)                                                                                                               ctgcAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCAtt              >  2:1193342/1‑136 (MQ=37)                                                                                                               ctgcAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTaa           <  1:346280/139‑1 (MQ=25)                                                                                                               ctgcAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTaa           <  2:1020650/139‑1 (MQ=25)                                                                                                               ctgcAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTaa           <  1:1356865/139‑1 (MQ=37)                                                                                                               ctgcAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTaa           <  2:1308982/139‑1 (MQ=25)                                                                                                               ctgcAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTaa           <  2:1439750/139‑1 (MQ=21)                                                                                                                        cAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATtgt  <  1:503978/139‑1 (MQ=21)                                                                                                                        cAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATtgt  <  2:471794/139‑1 (MQ=25)                                                                                                                         aGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCa                >  2:280057/1‑124 (MQ=25)                                                                                                                         aGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCa                <  1:280057/124‑1 (MQ=25)                                                                                                                                     |                                                                                                                              TAGTGCTAACTCCGGTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATTGT  >  NC_000913/1432141‑1432398

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