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%	166.8 / NA	52	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 (30/22);  total (30/22)

AATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATC  >  NC_000913/4577822‑4578092                                                                                                                                         |                                                                                                                                       aaTACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTa                                                                                                                                      >  2:694109/1‑139 (MQ=255)      atacCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGt                                                                                                                                 >  1:191021/1‑139 (MQ=255)        acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTcc                                                                                                                               >  2:453257/1‑139 (MQ=255)        acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTcc                                                                                                                               >  2:179780/1‑139 (MQ=255)                    cAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATg                                                                                                                                    >  2:574013/1‑122 (MQ=255)                    cAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATg                                                                                                                                    <  1:574013/122‑1 (MQ=255)                    cAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATttg                                                                                                                   >  2:623502/1‑139 (MQ=255)                    cAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATttg                                                                                                                   >  1:323429/1‑139 (MQ=255)                    cAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATttg                                                                                                                   >  2:765162/1‑139 (MQ=255)                     aagaagATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATttgt                                                                                                                  <  1:623502/139‑1 (MQ=255)                          gaTTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTc                                                                                                             <  2:191021/139‑1 (MQ=255)                           aTTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTcc                                                                                                            >  1:992892/1‑139 (MQ=255)                               ttACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCAttt                                                                                                        <  1:694109/139‑1 (MQ=255)                                 aCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCt                                                                                                      <  2:87376/139‑1 (MQ=255)                                  ccGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg                                                                                                     <  1:1036417/139‑1 (MQ=255)                                    gggTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt                                                                                                   <  1:179780/139‑1 (MQ=255)                                     ggTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg                                                                                                     <  1:626994/136‑1 (MQ=255)                                     ggTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg                                                                                                     >  2:626994/1‑136 (MQ=255)                                         tGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTcc                                                                                                                               <  1:614360/106‑1 (MQ=255)                                         tGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTcc                                                                                                                               >  2:614360/1‑106 (MQ=255)                                                aaGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGa                                                                                                                                   <  2:155889/95‑1 (MQ=255)                                                aaGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGa                                                                                                                                   >  1:155889/1‑95 (MQ=255)                                                           gCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGaaa                                                                            <  2:323429/139‑1 (MQ=255)                                                            cTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAat                                                                           >  2:1144573/1‑139 (MQ=255)                                                                           aTAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAg                                                            <  2:963860/139‑1 (MQ=255)                                                                                  cATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTc                                                                                                             <  1:971812/83‑1 (MQ=255)                                                                                  cATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTc                                                                                                             >  2:971812/1‑83 (MQ=255)                                                                                           gTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGtcat                                            >  1:91904/1‑139 (MQ=255)                                                                                                cATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCAttt                                                                                                        <  1:521617/74‑1 (MQ=255)                                                                                                cATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCAttt                                                                                                        >  2:521617/1‑74 (MQ=255)                                                                                                cATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACaa                                       >  2:94315/1‑139 (MQ=255)                                                                                                 aTGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAAt                                      >  1:451294/1‑139 (MQ=255)                                                                                                             ttACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTtata                                                                                            >  1:545369/1‑73 (MQ=255)                                                                                                             ttACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTtata                                                                                            <  2:545369/73‑1 (MQ=255)                                                                                                              tACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCAt                         >  2:1231441/1‑139 (MQ=255)                                                                                                                 ggACAGGTATGACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGc                      <  2:992892/139‑1 (MQ=255)                                                                                                                 ggACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGc                      <  1:1144573/139‑1 (MQ=255)                                                                                                                   aCAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGt                    >  2:459877/1‑139 (MQ=255)                                                                                                                       gTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt                                                                                                   >  2:63910/1‑56 (MQ=255)                                                                                                                       gTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt                                                                                                   <  1:63910/56‑1 (MQ=255)                                                                                                                       gTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATcc                >  2:517612/1‑139 (MQ=255)                                                                                                                              gggCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAgg         <  1:453257/139‑1 (MQ=255)                                                                                                                               ggCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTa                                                               <  2:433835/84‑1 (MQ=255)                                                                                                                               ggCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTa                                                               >  1:433835/1‑84 (MQ=255)                                                                                                                                 cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGc                      <  2:279417/123‑1 (MQ=255)                                                                                                                                 cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGc                      >  1:279417/1‑123 (MQ=255)                                                                                                                                 cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGc      >  2:1299934/1‑139 (MQ=255)                                                                                                                                 cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGc      >  1:677801/1‑139 (MQ=255)                                                                                                                                    ttCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTg                       >  2:1084303/1‑119 (MQ=255)                                                                                                                                    ttCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTg                       <  1:1084303/119‑1 (MQ=255)                                                                                                                                    ttCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCAt   >  1:856649/1‑139 (MQ=255)                                                                                                                                     tCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATc  >  1:375395/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                       AATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATC  >  NC_000913/4577822‑4578092

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