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_0009134,577,9580CT100.0% 157.1 / NA 49M1M (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 (25/24);  total (25/24)

ACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGC  >  NC_000913/4577829‑4578088
                                                                                                                                 |                                                                                                                                  
acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAg                                                                                                                              >  2:1423807/1‑136 (MQ=255)
acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAg                                                                                                                              <  1:1423807/136‑1 (MQ=255)
  cccAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTcccc                                                                                                                         >  2:1328622/1‑139 (MQ=255)
            cAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAg                                                                                                                              >  2:1436300/1‑124 (MQ=255)
            cAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAg                                                                                                                              <  1:1436300/124‑1 (MQ=255)
                         aCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCt                                                                                                  <  2:1712017/139‑1 (MQ=255)
                               tATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTc                                                                                                                            >  1:1404039/1‑107 (MQ=255)
                               tATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTc                                                                                                                            <  2:1404039/107‑1 (MQ=255)
                                 tGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTta                                                                                          >  1:1306215/1‑139 (MQ=255)
                                      tCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAggg                                                                                     <  1:1328622/139‑1 (MQ=255)
                                       cAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAgggg                                                                                    >  2:1455251/1‑139 (MQ=255)
                                        aaGGTTTGCCTGCTTAATTTCCTGTAAATAACCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGtc                                                                                  >  2:235235/1‑139 (MQ=255)
                                            tttGCCTGCTTAATTTCCTGTAAATAACCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGtcatca                                                                              <  1:235235/139‑1 (MQ=255)
                                                    cTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAat                                                                       >  2:516506/1‑139 (MQ=255)
                                                      tAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAatat                                                                     <  1:516506/139‑1 (MQ=255)
                                                      tAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGATGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAatat                                                                     >  2:1040536/1‑139 (MQ=255)
                                                                 aaaTAACCCCATGCATAGGTAAGCATGTTATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTaa                                                          <  1:1040536/139‑1 (MQ=255)
                                                                          cATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt                                                                                               <  1:253923/93‑1 (MQ=255)
                                                                          cATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt                                                                                               >  2:253923/1‑93 (MQ=255)
                                                                                             aaTAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTaaa                                                         >  2:231181/1‑112 (MQ=255)
                                                                                             aaTAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTaaa                                                         <  1:231181/112‑1 (MQ=255)
                                                                                                  atatTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATAttcttc                                                                <  1:664800/100‑1 (MQ=255)
                                                                                                  atatTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATAttcttc                                                                >  2:664800/1‑100 (MQ=255)
                                                                                                      tACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAgggg                                                     <  2:1131043/107‑1 (MQ=255)
                                                                                                      tACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAgggg                                                     >  1:1131043/1‑107 (MQ=255)
                                                                                                        cGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGtcatc                                       <  1:442244/119‑1 (MQ=255)
                                                                                                        cGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGtcatc                                       >  2:442244/1‑119 (MQ=255)
                                                                                                         ggACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGc                  <  1:937928/139‑1 (MQ=255)
                                                                                                            cAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTa               >  2:1024403/1‑139 (MQ=255)
                                                                                                             aGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTc                                                                                                         >  2:1692675/1‑48 (MQ=255)
                                                                                                             aGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTc                                                                                                         <  1:1692675/48‑1 (MQ=255)
                                                                                                             aGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTta                                                                                          <  2:1020322/63‑1 (MQ=255)
                                                                                                             aGGTACCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTta                                                                                          >  1:1020322/1‑63 (MQ=255)
                                                                                                               gTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTaaa                                                         <  1:523153/94‑1 (MQ=255)
                                                                                                               gTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTaaa                                                         >  2:523153/1‑94 (MQ=255)
                                                                                                               gTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGcc                       <  2:1012042/128‑1 (MQ=255)
                                                                                                               gTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGcc                       >  1:1012042/1‑128 (MQ=255)
                                                                                                                        gCTGTTCTATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTg   <  1:1455251/139‑1 (MQ=255)
                                                                                                                         cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt                                                                                               <  2:1303634/46‑1 (MQ=255)
                                                                                                                         cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt                                                                                               >  1:1303634/1‑46 (MQ=255)
                                                                                                                         cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg                                                                                                 <  2:110586/44‑1 (MQ=255)
                                                                                                                         cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg                                                                                                 >  1:110586/1‑44 (MQ=255)
                                                                                                                         cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAggg                                                                                     <  2:1280491/56‑1 (MQ=255)
                                                                                                                         cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAggg                                                                                     >  1:1280491/1‑56 (MQ=255)
                                                                                                                         cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGc  >  1:586088/1‑139 (MQ=255)
                                                                                                                            ttCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTta                                                                                          <  1:360246/48‑1 (MQ=255)
                                                                                                                            ttCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTta                                                                                          >  2:360246/1‑48 (MQ=255)
                                                                                                                            ttCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAAt                                  <  1:861735/104‑1 (MQ=255)
                                                                                                                            ttCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAAt                                  >  2:861735/1‑104 (MQ=255)
                                                                                                                                 |                                                                                                                                  
ACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGC  >  NC_000913/4577829‑4578088

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