Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,432,289	G→A	G1081S (GGT→AGT)	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,289	0	G	A	100.0%	220.7 / NA	63	G1081S (GGT→AGT)	stfR	Rac prophage, putative tail fiber protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/32);  total (31/32)

GTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATTGTGAGGCTTGCA  >  NC_000913/1432155‑1432408                                                                                                                                       |                                                                                                                        gTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGc                                                                                                                     >  1:568563/1‑139 (MQ=255)  tGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCt                                                                                                                    <  1:469855/139‑1 (MQ=255)          gTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCa                                                                                                            <  2:568563/139‑1 (MQ=255)            gCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACt                                                                                                          >  2:141993/1‑139 (MQ=255)                     aGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGa                                                                                                 >  2:532674/1‑139 (MQ=255)                      gATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACCGTCGGTATTAGTGCTCATACGCACTCCGTTGCGAt                                                                                                <  2:826627/139‑1 (MQ=255)                         tGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTgg                                                                                             <  2:135892/139‑1 (MQ=255)                          gTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGt                                                                                            >  2:123951/1‑139 (MQ=255)                                tgtgCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCt                                                                                                                    >  1:756349/1‑109 (MQ=255)                                tgtgCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTACGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCt                                                                                                                    <  2:756349/109‑1 (MQ=255)                                 gtgCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATg                                                                                     >  1:399286/1‑139 (MQ=255)                                       aaTCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGacaca                                                                               >  1:94103/1‑139 (MQ=255)                                         tCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACAcc                                                                             >  2:564487/1‑139 (MQ=255)                                           aaaaCTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCAt                                                                           <  2:94103/139‑1 (MQ=255)                                            aaaCTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATc                                                                          >  1:129699/1‑139 (MQ=255)                                                 aTGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGt                                                                     >  2:376895/1‑139 (MQ=255)                                                  tGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGtt                                                                                           <  1:981071/116‑1 (MQ=255)                                                  tGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGtt                                                                                           >  2:981071/1‑116 (MQ=255)                                                   gCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTg                                                                                                                      >  1:175709/1‑88 (MQ=255)                                                   gCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTg                                                                                                                      <  2:175709/88‑1 (MQ=255)                                                   gCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTa                                                                   >  2:898645/1‑139 (MQ=255)                                                    cAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTaa                                                                  <  2:226202/139‑1 (MQ=255)                                                      aCATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACg                                                                <  1:564487/139‑1 (MQ=255)                                                            tcTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTgg                                                                                             <  2:194056/104‑1 (MQ=255)                                                            tcTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTgg                                                                                             >  1:194056/1‑104 (MQ=255)                                                            tcTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGctgctg                                                          >  2:908756/1‑139 (MQ=255)                                                               gctgGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTa                                                       <  1:318901/139‑1 (MQ=255)                                                                 tgGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAAc                                                     >  2:670850/1‑139 (MQ=255)                                                                  ggCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGGTGCTGGTAAcg                                                    >  2:726500/1‑139 (MQ=255)                                                                    cgcACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATg                                                                                     <  1:437427/104‑1 (MQ=255)                                                                    cgcACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATg                                                                                     >  2:437427/1‑104 (MQ=255)                                                                          tACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAac                                            <  1:958566/139‑1 (MQ=255)                                                                                ctcactGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTc                                      >  1:890824/1‑139 (MQ=255)                                                                                  cactGTCCGGCACTGCTGCAAGCGCAGGTGCATACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCaa                                    >  1:664646/1‑139 (MQ=255)                                                                                     tGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCAc                                                                        <  1:162022/100‑1 (MQ=255)                                                                                     tGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCAc                                                                        >  2:162022/1‑100 (MQ=255)                                                                                     tGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAAc                                                     >  1:196148/1‑119 (MQ=255)                                                                                     tGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAAc                                                     <  2:196148/119‑1 (MQ=255)                                                                                     tGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCaaaaa                                 >  2:172627/1‑139 (MQ=255)                                                                                     tGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATGAGTGCTCATACACACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCAGCTGGTAACGCGGAAAACACCGTCaaaaa                                 >  1:723424/1‑139 (MQ=255)                                                                                       tCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACa                               >  1:26343/1‑139 (MQ=255)                                                                                       tCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACa                               <  1:532674/139‑1 (MQ=255)                                                                                       tCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACa                               >  1:1046822/1‑139 (MQ=255)                                                                                         cGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATc                             >  2:335788/1‑139 (MQ=255)                                                                                          ggCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTAACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCg                            <  2:664646/139‑1 (MQ=255)                                                                                           gCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACG‑GGAAAACACCGTCAAAAACATCGCa                          <  1:726500/139‑1 (MQ=255)                                                                                            cACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCa                          <  2:26343/139‑1 (MQ=255)                                                                                            cACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCa                          <  2:129699/139‑1 (MQ=255)                                                                                              ctgctgCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCAtt                        <  1:172627/139‑1 (MQ=255)                                                                                              ctgctgCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCAtt                        <  1:141993/139‑1 (MQ=255)                                                                                              ctgctgCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCAtt                        <  2:844416/139‑1 (MQ=255)                                                                                                 ctgcAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTaa                     >  1:354397/1‑139 (MQ=255)                                                                                                   gcAAGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACt                   <  2:491692/139‑1 (MQ=255)                                                                                                      aGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGacaca                                                                               <  2:382269/76‑1 (MQ=255)                                                                                                      aGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGacaca                                                                               >  1:382269/1‑76 (MQ=255)                                                                                                      aGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGacaca                                                                               <  2:920500/76‑1 (MQ=255)                                                                                                      aGCGCAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGacaca                                                                               >  1:920500/1‑76 (MQ=255)                                                                                                          cAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATtgt            <  1:123951/139‑1 (MQ=255)                                                                                                          cAGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATtgt            <  2:1046822/139‑1 (MQ=255)                                                                                                           aGGTGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATtgtg           >  1:545148/1‑139 (MQ=255)                                                                                                              tGCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATTGTGAgg        <  2:1111061/139‑1 (MQ=38)                                                                                                               gCACACGCGCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATTGTGAGGc       <  2:890824/139‑1 (MQ=37)                                                                                                                    cgcgCATACTGTCGGTATTAGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATTGTGAGGCTTGCa  <  2:1052404/139‑1 (MQ=18)                                                                                                                                       |                                                                                                                        GTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATTGTGAGGCTTGCA  >  NC_000913/1432155‑1432408

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