Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,210,758	C→T	A138V (GCG→GTG)	mcrA →	e14 prophage; 5‑methylcytosine‑specific restriction enzyme McrA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,210,758	0	C	T	100.0%	214.5 / NA	62	A138V (GCG→GTG)	mcrA	e14 prophage; 5‑methylcytosine‑specific restriction enzyme McrA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (34/28);  total (34/28)

CCCGCGATATACGCTTGAAGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAG  >  NC_000913/1210623‑1210892                                                                                                                                        |                                                                                                                                       cccGCGATATACGCTTGAAGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCt                                                                                                                                    >  2:533686/1‑140 (MQ=255)   cgcgATATACGCTTGAAGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCtat                                                                                                                                  >  2:560222/1‑140 (MQ=255)   cgcgATATACGCTTGAAGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCtat                                                                                                                                  >  2:639497/1‑140 (MQ=255)   cgcgATATACGCTTGAAGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCtat                                                                                                                                  >  2:74455/1‑140 (MQ=255)   cgcgATATACGCTTGAAGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCtat                                                                                                                                  >  2:670590/1‑140 (MQ=255)           aCGCTTGAAGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTg                                                                                                                          <  1:225901/140‑1 (MQ=255)            cGCTTGAAGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTgg                                                                                                                         <  2:789359/140‑1 (MQ=255)               ttGAAGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGgaga                                                                                                                      >  2:661375/1‑140 (MQ=255)                  aaGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACt                                                                                                                   >  1:1019504/1‑140 (MQ=255)                  aaGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACt                                                                                                                   >  2:654925/1‑140 (MQ=255)                    gATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTAt                                                                                                                 >  1:241528/1‑140 (MQ=255)                     aTGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATc                                                                                                                >  1:30984/1‑140 (MQ=255)                               aTCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTca                                                                                                      <  2:1019504/140‑1 (MQ=255)                               aTCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTca                                                                                                      <  1:74455/140‑1 (MQ=255)                               aTCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTca                                                                                                      <  2:1011723/140‑1 (MQ=255)                               aTCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTca                                                                                                      >  1:817481/1‑140 (MQ=255)                               aTCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTca                                                                                                      <  1:560222/140‑1 (MQ=255)                               aTCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTca                                                                                                      <  2:38897/140‑1 (MQ=255)                               aTCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTca                                                                                                      <  1:286536/140‑1 (MQ=255)                               aTCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTca                                                                                                      <  2:825909/140‑1 (MQ=255)                               aTCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTca                                                                                                      <  1:106023/140‑1 (MQ=255)                                  aGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAg                                                                                                   <  2:1120171/140‑1 (MQ=255)                                  aGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAg                                                                                                   >  2:188467/1‑140 (MQ=255)                                  aGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAg                                                                                                   >  1:53202/1‑140 (MQ=255)                                   ggCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGc                                                                                                  <  2:95936/140‑1 (MQ=255)                                   ggCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGc                                                                                                  <  1:1190449/140‑1 (MQ=255)                                    gCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGGTTTCACAGcc                                                                                                 >  2:1060548/1‑140 (MQ=255)                                    gCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGcc                                                                                                 <  1:533686/140‑1 (MQ=255)                                      ttATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTa                                                                                               <  1:901546/140‑1 (MQ=255)                                       tATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTAc                                                                                              >  1:872393/1‑140 (MQ=255)                                              gTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAgatgat                                                                                       >  2:586330/1‑140 (MQ=255)                                              gTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAgatgat                                                                                       >  1:778168/1‑140 (MQ=255)                                                  tACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAAt                                                                                   <  2:256092/140‑1 (MQ=255)                                                  tACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAAt                                                                                   <  1:670590/140‑1 (MQ=255)                                                   aCGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAAt                                                                                   <  2:902622/139‑1 (MQ=255)                                                   aCGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAAt                                                                                   >  1:902622/1‑139 (MQ=255)                                                              ggAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAata                                                                       <  2:241528/140‑1 (MQ=255)                                                                 aaTTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGa                                                                    >  2:302492/1‑140 (MQ=255)                                                                                            aaTATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGt                                         >  1:1011693/1‑140 (MQ=255)                                                                                                 tGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGt                                    <  2:92064/140‑1 (MQ=255)                                                                                                          ttCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAg                           <  2:550024/140‑1 (MQ=255)                                                                                                          ttCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAg                           <  1:1207113/140‑1 (MQ=255)                                                                                                          ttCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAg                           <  2:504589/140‑1 (MQ=255)                                                                                                           tCCAGGTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAgg                          >  1:558677/1‑140 (MQ=255)                                                                                                                gTATTTATAGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTcc                     >  2:204461/1‑140 (MQ=255)                                                                                                                     tataGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGaa                >  2:422914/1‑140 (MQ=255)                                                                                                                      ataGTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGaaa               <  2:30984/140‑1 (MQ=255)                                                                                                                         gTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAg                           >  2:26010/1‑125 (MQ=255)                                                                                                                         gTGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAg                           <  1:26010/125‑1 (MQ=255)                                                                                                                          tGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAg           >  2:456282/1‑140 (MQ=255)                                                                                                                          tGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAg           <  1:302492/140‑1 (MQ=255)                                                                                                                          tGACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAg           >  1:1191616/1‑140 (MQ=255)                                                                                                                          tGACAGCTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAg           >  2:530482/1‑140 (MQ=255)                                                                                                                           gACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGt          >  1:782505/1‑140 (MQ=255)                                                                                                                           gACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGt          >  1:162074/1‑140 (MQ=255)                                                                                                                           gACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGt          >  2:596369/1‑140 (MQ=255)                                                                                                                           gACAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGCAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGt          >  1:443216/1‑140 (MQ=255)                                                                                                                            aCAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGtt         >  2:939232/1‑140 (MQ=255)                                                                                                                             cAGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTg        <  1:1060548/140‑1 (MQ=255)                                                                                                                              aGTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGa       >  2:755999/1‑140 (MQ=255)                                                                                                                               gTTTTTGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAg      >  2:308960/1‑140 (MQ=255)                                                                                                                                   ttGGGTGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGCAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAg  <  2:443216/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                       CCCGCGATATACGCTTGAAGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAG  >  NC_000913/1210623‑1210892

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