Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	582,081	A→G	N183S (AAC→AGC)	tfaD →	DLP12 prophage; putative tail fiber assembly protein TfaD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	582,081	0	A	G	100.0%	202.6 / NA	58	N183S (AAC→AGC)	tfaD	DLP12 prophage; putative tail fiber assembly protein TfaD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/24);  total (34/24)

GCGCCACTTCAGGATGCTGTAGATCTGGAAATCGCAACGGAGGAAGAAACCTCATTGCTGGAAGCCTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCT  >  NC_000913/581945‑582213                                                                                                                                         |                                                                                                                                     gcgcCACTTCAGGATGCTGTAGATCTGGAAATCGCAACGGAGGAAGAAACCTCATTGCTGGAAGCCTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGcc                                                                                                                                    >  1:421812/1‑139 (MQ=255)   gcCACTTCAGGATGCTGTAGATCTGGAAATCGCAACGGAGGAAGAAACCTCATTGCTGGAAGCCTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTg                                                                                                                                 >  1:223015/1‑140 (MQ=255)    ccACTTCAGGATGCTGTAGATCTGGAAATCGCAACGGAGGAAGAAACCTCATTGCTGGAAGCCTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGt                                                                                                                                >  1:91215/1‑140 (MQ=255)     cACTTCAGGATGCTGTAGATCTGGAAATCGCAACGGAGGAAGAAACCTCATTGCTGGAAGCCTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTc                                                                                                                               <  1:635635/140‑1 (MQ=255)     cACTTCAGGATGCTGTAGATCTGGAAATCG‑AACGGAGGAAGAAACCTCATTGCTGGAAGCCTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCa                                                                                                                              <  1:282701/140‑1 (MQ=255)       cTTCAGGATGCTGTAGATCTGGAAATCGCAACGGAGGAAGAAACCTCATTGCTGGAAGCCTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAg                                                                                                                             >  2:435441/1‑140 (MQ=255)       cTTCAGGATGCTGTAGATCTGGAAATCGCAACGGAGGAAGAAACCTCATTGCTGGAAGCCTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAg                                                                                                                             <  2:91215/140‑1 (MQ=255)          cAGGATGCTGTAGATCTGGAAATCGCAACGGAGGAAGAAACCTCATTGCTGGAAGCCTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGcc                                                                                                                                    >  1:308194/1‑130 (MQ=255)            ggATGCTGTAGATCTGGAAATCGCAACGGAGGAAGAAACCTCATTGCTGGAAGCCTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGt                                                                                                                        <  1:435441/140‑1 (MQ=255)             gATGCTGTAGATCTGGAAATCGCAACGGAGGAAGAAACCTCATTGCTGGAAGCCTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTa                                                                                                                       <  1:317753/140‑1 (MQ=255)                 cTGTAGATCTGGAAATCGCAACGGAGGAAGAAACCTCATTGCTGGAAGCCTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAg                                                                                                                         <  1:101334/134‑1 (MQ=255)                 cTGTAGATCTGGAAATCGCAACGGAGGAAGAAACCTCATTGCTGGAAGCCTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAg                                                                                                                         >  2:101334/1‑134 (MQ=255)                  tGTAGATCTGGAAATCGCAACGGAGGAAGAAACCTCATTGCTGGAAGCCTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCAt                                                                                                                  >  2:439380/1‑140 (MQ=255)                     aGATCTGGAAATCGCAACGGAGGAAGAAACCTCATTGCTGGAAGCCTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTgg                                                                                                               <  2:38645/140‑1 (MQ=255)                      gATCTGGAAATCGCAACGGAGGAAGAAACCTCATTGCTGGAAGCCTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTggg                                                                                                              >  1:319077/1‑140 (MQ=255)                            gAAATCGCAACGGAGGAAGAAACCTCATTGCTGGAAGCCTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATg                                                                                                        >  1:405433/1‑140 (MQ=255)                              aaTCGCAACGGAGGAAGAAACCTCATTGCTGGAAGCCTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGcc                                                                                                      >  2:174033/1‑140 (MQ=255)                                 cGCAACGGAGGAAGAAACCTCATTGCTGGAAGCCTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCgca                                                                                                   >  2:410643/1‑140 (MQ=255)                                 cGCAACGGAGGAAGAAACCTCATTGCTGGAAGCCTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACAAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCgca                                                                                                   >  1:91930/1‑140 (MQ=255)                                   cAACGGAGGAAGAAACCTCATTGCTGGAAGCCTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCgcagc                                                                                                 <  2:319077/140‑1 (MQ=255)                                          ggaAGAAACCTCATTGCTGGAAGCCTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCtt                                                                                          <  1:410643/140‑1 (MQ=255)                                                  ccTCATTGCTGGAAGCCTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGa                                                                                  <  1:446367/140‑1 (MQ=255)                                                          cTGGAAGCCTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACa                                                                               >  2:266326/1‑135 (MQ=255)                                                          cTGGAAGCCTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACa                                                                               <  1:266326/135‑1 (MQ=255)                                                             gAAGCCTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACAtgctg                                                                          >  1:212101/1‑137 (MQ=255)                                                             gAAGCCTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACAtgctg                                                                          <  2:212101/137‑1 (MQ=255)                                                               aGCCTGGAAAAAATATCGGGAGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGtt                                                                     >  2:137655/1‑140 (MQ=255)                                                                                gggTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTcc                                                    >  1:18387/1‑140 (MQ=255)                                                                                  gtgtTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTg                                                  <  1:595777/140‑1 (MQ=255)                                                                                      tGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGc                                              >  2:235488/1‑140 (MQ=255)                                                                                       gCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCg                                             <  1:174033/140‑1 (MQ=255)                                                                                          gAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGaa                                          <  1:218518/140‑1 (MQ=255)                                                                                            aCCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGaaaa                                        >  2:650856/1‑140 (MQ=255)                                                                                            aCCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGaaaa                                        >  1:86373/1‑140 (MQ=255)                                                                                            aCCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGaaaa                                        <  1:235488/140‑1 (MQ=255)                                                                                               gtgtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACAt                                     <  1:425973/140‑1 (MQ=255)                                                                                               gtgtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACAt                                     <  2:86373/140‑1 (MQ=255)                                                                                                tgtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACAtt                                    >  1:417961/1‑140 (MQ=255)                                                                                                tgtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACAtt                                    >  2:667840/1‑140 (MQ=255)                                                                                                tgtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACAtt                                    >  1:469467/1‑140 (MQ=255)                                                                                                tgtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACAtt                                    >  1:42394/1‑140 (MQ=255)                                                                                                 gtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATta                                   <  2:208686/140‑1 (MQ=255)                                                                                                     aTACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACAt                                     >  2:373693/1‑134 (MQ=255)                                                                                                     aTACATCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACAt                                     <  1:373693/134‑1 (MQ=255)                                                                                                          tCAACTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTAc                          >  2:13267/1‑140 (MQ=255)                                                                                                            aaCTGCACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCt                        >  1:317549/1‑140 (MQ=255)                                                                                                                 cACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAAc                   >  2:127247/1‑140 (MQ=255)                                                                                                                 cACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAAc                   >  1:655393/1‑140 (MQ=255)                                                                                                                 cACCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAAc                   >  2:615798/1‑140 (MQ=255)                                                                                                                  aCCTGATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAAcc                  <  1:127247/140‑1 (MQ=255)                                                                                                                      gATATTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTgg              >  2:686130/1‑140 (MQ=255)                                                                                                                       atatTGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGt             <  2:405433/140‑1 (MQ=255)                                                                                                                           tGAGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTa         >  2:93319/1‑140 (MQ=255)                                                                                                                             aGTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAAc       >  1:487702/1‑140 (MQ=255)                                                                                                                              gTGGCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACa      <  2:120672/140‑1 (MQ=255)                                                                                                                                 gCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATtc   <  1:650856/140‑1 (MQ=255)                                                                                                                                 gCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATtc   >  1:236613/1‑140 (MQ=255)                                                                                                                                 gCCTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATtc   >  1:18254/1‑140 (MQ=255)                                                                                                                                  ccTACGAGCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATtct  >  1:535935/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                     GCGCCACTTCAGGATGCTGTAGATCTGGAAATCGCAACGGAGGAAGAAACCTCATTGCTGGAAGCCTGGAAAAAATATCGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCT  >  NC_000913/581945‑582213

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