Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,577,317 C→T T214T (ACG→ACA)  mcrC ← 5‑methylcytosine‑specific restriction enzyme subunit McrC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,577,3170CT100.0% 166.0 / NA 48T214T (ACG→ACA) mcrC5‑methylcytosine‑specific restriction enzyme subunit McrC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/23);  total (25/23)

ATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTG  >  NC_000913/4577181‑4577452
                                                                                                                                        |                                                                                                                                       
aTTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGtt                                                                                                                                      <  1:665931/140‑1 (MQ=255)
 tttGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAGGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTa                                                                                                                                     >  2:383413/1‑140 (MQ=255)
   tGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAAt                                                                                                                                   >  2:706886/1‑140 (MQ=255)
   tGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAAt                                                                                                                                   >  1:362387/1‑140 (MQ=255)
       gTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCTGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCac                                                                                                                               >  1:72679/1‑140 (MQ=255)
         cAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCacga                                                                                                                             >  1:89844/1‑140 (MQ=255)
         cAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCacga                                                                                                                             >  1:71016/1‑140 (MQ=255)
         cAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCacg                                                                                                                              <  1:183344/139‑1 (MQ=255)
         cAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCacg                                                                                                                              >  2:183344/1‑139 (MQ=255)
                           cTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGa                                                                                                           <  2:534735/140‑1 (MQ=255)
                                   cGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTg                                                                                                   >  1:319725/1‑140 (MQ=255)
                                    gAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGa                                                                                                  <  2:31357/140‑1 (MQ=255)
                                     aaTGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGAt                                                                                                 <  1:466329/140‑1 (MQ=255)
                                      aTGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATa                                                                                                <  2:362387/140‑1 (MQ=255)
                                            aTGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTaa                                                                                          >  1:472952/1‑140 (MQ=255)
                                             tGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAAt                                                                                         <  2:588632/140‑1 (MQ=255)
                                                 aGTTTCCATTCGAGGTGACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACa                                                                                     >  2:479079/1‑140 (MQ=255)
                                                   tttCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTCAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACAtc                                                                                   >  1:342954/1‑140 (MQ=255)
                                                    ttCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACAtct                                                                                  >  1:498070/1‑140 (MQ=255)
                                                        aTTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCttt                                                                              >  2:629905/1‑140 (MQ=255)
                                                                  aCAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGtttctt                                                                    >  1:462126/1‑140 (MQ=255)
                                                                     aaTTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCa                                                                 <  1:706886/140‑1 (MQ=255)
                                                                         aaGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAAt                                                             <  2:174753/140‑1 (MQ=255)
                                                                          aGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATc                                                            <  1:595173/140‑1 (MQ=255)
                                                                          aGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATc                                                            >  2:595173/1‑140 (MQ=255)
                                                                          aGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATc                                                            <  2:498070/140‑1 (MQ=255)
                                                                            tGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCAt                                                          <  2:472952/140‑1 (MQ=255)
                                                                                 gATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGaaa                                                     <  1:479079/140‑1 (MQ=255)
                                                                                  aTCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAAc                                                    >  1:233915/1‑140 (MQ=255)
                                                                                        tataCTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAgt                                              <  1:383413/140‑1 (MQ=255)
                                                                                            cTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATc                                                            <  2:127345/122‑1 (MQ=255)
                                                                                            cTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATc                                                            >  1:127345/1‑122 (MQ=255)
                                                                                            cTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTcc                                          >  1:160095/1‑140 (MQ=255)
                                                                                             tCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCt                                         <  1:341910/140‑1 (MQ=255)
                                                                                             tCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCt                                         <  1:529057/140‑1 (MQ=255)
                                                                                             tCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCt                                         >  2:469577/1‑140 (MQ=255)
                                                                                               gATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGaaa                                                     <  1:428028/126‑1 (MQ=255)
                                                                                               gATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGaaa                                                     >  2:428028/1‑126 (MQ=255)
                                                                                                 tGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGt                                     >  1:186206/1‑140 (MQ=255)
                                                                                                      cccATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACCTCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGttt                                   >  1:78782/1‑137 (MQ=255)
                                                                                                      cccATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGttt                                   <  2:78782/137‑1 (MQ=255)
                                                                                                      cccATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGa                                >  1:208242/1‑140 (MQ=255)
                                                                                                       ccATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGAc                               <  2:273086/140‑1 (MQ=255)
                                                                                                        cATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGAcc                              <  2:423235/140‑1 (MQ=255)
                                                                                                        cATTTTAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGAcc                              <  2:71016/140‑1 (MQ=255)
                                                                                                           tttAAATAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTgg                           >  1:609188/1‑140 (MQ=255)
                                                                                                               aaTAAGAGCGGGTTGTGTTTGCAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATa                       <  2:89844/140‑1 (MQ=255)
                                                                                                                                    cAGATGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTg  >  2:606577/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                       
ATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTG  >  NC_000913/4577181‑4577452

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑