Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_013198 2,164,868 T→C V220V (GTA→GTG LGG_02097 type II restriction enzyme methyltransferase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0131982,164,8680T→C100.0% 13.9 47V220V (GTA→GTGLGG_02097type II restriction enzyme methyltransferase
Reads supporting (aligned to +/- strand):  new base (22/3):  ref base (14/8):  total (36/11)

GACAACCCAGTCCGTCGTAAACAGTTGGGTGGCGGCAGGGATGTCGTTTGTCTTGATGGGACCACCATTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGGAGATATTGATTTAGCGCGTTCGCTTGCTT  .  NC_013198/2164718‑2165011
                                                                                                                                                      |                                                                                                                                               
ggcAACCCAGTCCGTCGTAAACAGTTGGGTGGCGGCAGGGATGTCGTTTGTCTTGATGGGACCACCATTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCt                                                                                                                                                 <  2:991284/149‑1
    aCCCAGTCCGTCGTAAACAGTTGGGTGGCGGCAGGGATGTCGTTTGTCTTGATGGGACCACCATTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACAt                                                                                                                                             >  1:393323/1‑151
    aCCCAGTCCGTCGTAAACAGTTGGGTGGCGGCAGGGATGTCGTTTGTCTTGATGGGACCACCATTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACAt                                                                                                                                             >  2:892676/1‑151
             gtcgtAAACAGTTGGGTGGCGGCAGGGATGTCGTTTGTCTTGATGGGACCACCATTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAatcat                                                                                                                                    >  1:35285/1‑151
                     cAGTTGGGTGGCGGCAGGGATGTCGTTTGTCTTGATGGGACCACCATTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGa                                                                                                                            >  2:406147/1‑151
                         tGGGTGGCGGCAGGGATGTCGTTTGTCTTGATGGGACCACCATTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAatat                                                                                                                        <  1:101820/151‑1
                                cggcAGGGATGTCGTTTGTCTTGATGGGACCACCATTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCcac                                                                                                                 >  1:893274/1‑151
                                        aTGTCGTTTGTCTTGATGGGACCACCATTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCAc                                                                                                         <  1:828095/151‑1
                                           tCGTTTGTCTTGATGGGACCACCATTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGaa                                                                                                      <  1:463638/151‑1
                                            cGTTTGTCTTGATGGGACCACCATTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAAt                                                                                                     >  1:1046790/1‑151
                                                 gTCTTGATGGGACCACCATTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGAccc                                                                                                >  2:53761/1‑151
                                                 gTCTTGATGGGACCACCATTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGAccc                                                                                                <  1:1006582/151‑1
                                                     ataagagacagccaTTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAAc                                                                                            <  1:274902/140‑1
                                                               accaTTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCtt                                                                                  >  1:462999/1‑151
                                                                ccaTTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAAcctgtctcttat                                                                                 >  2:274908/1‑141
                                                                ccaTTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTg                                                                                 >  2:425100/1‑151
                                                                ccaTTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTg                                                                                 >  1:246241/1‑151
                                                                   ttGGTGTTGACTTCTTTTTCGTGCATGTCAATGTGGTAATACTGATACCGCCACCCAATAATCTCAACCTGACCGCACTTGTCTACATCAAAATCATCTGCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAg                                                                              <  2:362501/151‑1
                                                                   ttGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAg                                                                              <  2:624568/151‑1
                                                                          tGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACaa                                                                       >  1:667079/1‑151
                                                                               actTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACaaaaaga                                                                  >  2:412251/1‑146
                                                                               actTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACaaaaaga                                                                  >  2:377527/1‑146
                                                                                  ttATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACaaaaagaaag                                                               <  1:53761/151‑9
                                                                                                 gtattgtAATACTGATACAGCCACCCAATAATCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGttt                                                >  2:508984/4‑151
                                                                                                 gtattgtAATACTGATACAGCCACCCAATAATCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGttt                                                >  1:1003256/4‑151
                                                                                                 gtattgtAATACTGATACAGCCACCCAATAATCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGttt                                                >  1:338026/4‑151
                                                                                                    ttgtAATACTGATACAGCCACCCAATAATCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGtttttt                                             >  1:244456/1‑151
                                                                                                    ttgtAATACTGATACAGCCACCCAATAATCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGtttttt                                             >  2:661692/1‑151
                                                                                                      gtAATANNNNNNCNNCNNNNNNNNAATCTCAACCTGACCGCNNNNNNNNNNNNNNNNNNNNNNNTANNNNATATCATCCNNNANANNACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGattttttg                                           >  2:177061/1‑149
                                                                                                      gtAATACTGATACAGCCACCCAATAATCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCg                                           >  1:371438/1‑151
                                                                                                          tACTGATACAGCCACCCAATAATCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGaaaa                                       >  2:767580/1‑151
                                                                                                           aCTGATACAGCCACCCAATAATCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTGGaaaat                                      >  2:897353/1‑150
                                                                                                            cTGATACAGCCACCCAATAATCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGt                                     >  1:851456/1‑151
                                                                                                            cTGATACAGCCACCCAATAATCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGt                                     >  1:34872/1‑151
                                                                                                                  cAGCCACCCAATAATCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAggg                               >  2:5053/1‑151
                                                                                                                   aGCCACCCAATAATCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGga                              <  2:244468/151‑1
                                                                                                                          cAATAATCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGGAGATATTg                       >  2:180367/1‑151
                                                                                                                          cAATAATCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGGAGATATTg                       >  2:25802/1‑151
                                                                                                                          cAATAATCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGGAGATATTg                       >  1:611361/1‑151
                                                                                                                              aatCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGGAGATATTGAttt                   <  1:509124/151‑1
                                                                                                                               atCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGGAGATATTGATTTa                  >  2:226328/1‑151
                                                                                                                               atCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGGAGATATTGATTTa                  >  2:205756/1‑151
                                                                                                                               atCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGGAGATATTGATTTa                  >  1:923809/1‑151
                                                                                                                                   cAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGGAGATATTGATTTAgcgc              >  1:838342/1‑151
                                                                                                                                      ccTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGATCAGGGAGATATTGATTTAgcacgct           >  2:80666/1‑146
                                                                                                                                        tGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGGAGATATTGATTTATCGCGTTCg         >  2:932217/1‑151
                                                                                                                                         gACCACCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGGAGATATTGATTTAGCGCGTTCgc        >  1:303389/1‑151
                                                                                                                                              cccTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGGAGATATTGATTTAGCGCGTTCgcttgct   >  2:242244/1‑151
                                                                                                                                               ccTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGGAGATATTGATTTAGCGCGTTCgcttgctt  <  1:730797/151‑1
                                                                                                                                                      |                                                                                                                                               
GACAACCCAGTCCGTCGTAAACAGTTGGGTGGCGGCAGGGATGTCGTTTGTCTTGATGGGACCACCATTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGGAGATATTGATTTAGCGCGTTCGCTTGCTT  .  NC_013198/2164718‑2165011

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG