Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,577,958	C→T	M1M (GTG→ATG) † *460* (TAG→TAA)	mcrC ← mcrB ←	5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC 5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,577,958	0	C	T	100.0%	220.6 / NA	69	M1M (GTG→ATG) † *460* (TAG→TAA)	mcrC mcrB	5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC 5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (39/30);  total (39/30)

ACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCA  >  NC_000913/4577829‑4578090                                                                                                                                  |                                                                                                                                     acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTcc                                                                                                                             >  2:612022/1‑139 (MQ=255) acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTcc                                                                                                                             >  2:530901/1‑139 (MQ=255)      aGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCtaa                                                                                                                        >  2:1195589/1‑139 (MQ=255)      aGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCtaa                                                                                                                        >  2:1326475/1‑139 (MQ=255)           atCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAAtt                                                                                                                   >  2:1009949/1‑139 (MQ=255)                  agaTTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTccc                                                                                                                            <  2:1253015/123‑1 (MQ=255)                  agaTTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTccc                                                                                                                            >  1:1253015/1‑123 (MQ=255)                  agaTTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGt                                                                                                            >  2:246453/1‑139 (MQ=255)                   gaTTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTc                                                                                                           <  1:246453/139‑1 (MQ=255)                      tGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAg                                                                                                                                >  1:1292945/1‑115 (MQ=255)                      tGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAg                                                                                                                                <  2:1292945/115‑1 (MQ=255)                        ttACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATttgt                                                                                                                <  1:489691/129‑1 (MQ=255)                        ttACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATttgt                                                                                                                >  2:489691/1‑129 (MQ=255)                             gggTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt                                                                                                 <  2:437203/139‑1 (MQ=255)                                       tCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAggg                                                                                       <  2:390534/139‑1 (MQ=255)                                        cAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAgggg                                                                                      >  1:462584/1‑139 (MQ=255)                                         aaGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGt                                                                                     >  1:141417/1‑139 (MQ=255)                                          aGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGtc                                                                                    <  2:462584/139‑1 (MQ=255)                                              ttGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGTGGtcatca                                                                                >  2:885137/1‑139 (MQ=255)                                               tGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCt                                                                                                                          <  2:1018211/96‑1 (MQ=255)                                               tGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCt                                                                                                                          >  1:1018211/1‑96 (MQ=255)                                                 ccTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAg                                                                             <  1:530901/139‑1 (MQ=255)                                                    gCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGaaa                                                                          <  1:1326475/139‑1 (MQ=255)                                                                         cccATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCg                                                     >  1:528541/1‑139 (MQ=255)                                                                           cATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGa                                                                                                                                 >  1:968609/1‑61 (MQ=255)                                                                           cATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGa                                                                                                                                 <  2:968609/61‑1 (MQ=255)                                                                           cATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGat                                                   >  2:1144593/1‑139 (MQ=255)                                                                             tGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGatat                                                 <  1:1009949/139‑1 (MQ=255)                                                                             tGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGatat                                                 >  1:1003205/1‑139 (MQ=255)                                                                                  aGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTAt                                                                                                                                   >  2:541142/1‑52 (MQ=255)                                                                                  aGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTAt                                                                                                                                   <  1:541142/52‑1 (MQ=255)                                                                                  aGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTc                                                                                                                              >  2:1056083/1‑57 (MQ=255)                                                                                  aGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTc                                                                                                                              <  1:1056083/57‑1 (MQ=255)                                                                                      aaGTATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGtcatca                                                                                <  2:1041677/99‑1 (MQ=255)                                                                                      aaGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGtcatca                                                                                >  1:1041677/1‑99 (MQ=255)                                                                                       aGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATcac                                       >  2:1162176/1‑139 (MQ=255)                                                                                        gCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATcaca                                      <  2:1003205/139‑1 (MQ=255)                                                                                          aTGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAAt                                    <  1:612022/139‑1 (MQ=255)                                                                                            gTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCt                                                                                                    >  1:366497/1‑73 (MQ=255)                                                                                            gTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCt                                                                                                    <  2:366497/73‑1 (MQ=255)                                                                                                   atatTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAg                           >  1:1032927/1‑139 (MQ=255)                                                                                                    tatTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTGGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGc                          >  1:257906/1‑139 (MQ=255)                                                                                                      ttACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg                                                                                                   <  2:652996/64‑1 (MQ=255)                                                                                                      ttACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg                                                                                                   >  1:652996/1‑64 (MQ=255)                                                                                                      ttACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAggg                                                                                       <  2:1031292/76‑1 (MQ=255)                                                                                                      ttACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAggg                                                                                       >  1:1031292/1‑76 (MQ=255)                                                                                                       tACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATttgt                                                                                                                <  1:321803/50‑1 (MQ=255)                                                                                                       tACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATttgt                                                                                                                >  2:321803/1‑50 (MQ=255)                                                                                                       tACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCAt                       <  2:1032927/139‑1 (MQ=255)                                                                                                       tACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCAt                       >  2:1191858/1‑139 (MQ=255)                                                                                                        aCGGACAGGTATCACGGGCTGTTCCATTATGATTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCAtt                      >  1:310010/1‑139 (MQ=255)                                                                                                            aCAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGt                  <  1:1162176/139‑1 (MQ=255)                                                                                                              aGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAgggg                                                       <  1:175883/100‑1 (MQ=255)                                                                                                              aGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAgggg                                                       >  2:175883/1‑100 (MQ=255)                                                                                                                gTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt                                                                                                 >  1:138858/1‑56 (MQ=255)                                                                                                                gTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt                                                                                                 <  2:138858/56‑1 (MQ=255)                                                                                                                gTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATcc              >  1:457437/1‑139 (MQ=255)                                                                                                                   tCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGga           <  2:310010/139‑1 (MQ=255)                                                                                                                      cGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGagag        >  2:451993/1‑139 (MQ=255)                                                                                                                      cGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGagag        >  1:734308/1‑139 (MQ=255)                                                                                                                        ggCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGt      >  2:1159614/1‑139 (MQ=255)                                                                                                                          cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGtcatca                                                                                <  2:631921/63‑1 (MQ=255)                                                                                                                          cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGtcatca                                                                                >  1:631921/1‑63 (MQ=255)                                                                                                                          cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGtcatca                                        <  2:333454/103‑1 (MQ=255)                                                                                                                          cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGtcatca                                        >  1:333454/1‑103 (MQ=255)                                                                                                                           tGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGcc   >  1:159016/1‑139 (MQ=255)                                                                                                                            gTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCa  <  1:1191858/139‑1 (MQ=255)                                                                                                                             ttCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGa                                                                            <  1:742601/64‑1 (MQ=255)                                                                                                                             ttCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGa                                                                            >  2:742601/1‑64 (MQ=255)                                                                                                                                  |                                                                                                                                     ACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCA  >  NC_000913/4577829‑4578090

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