Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0131982,164,8590A→G100.0% 0.5 50G223G (GGT→GGCLGG_02097type II restriction enzyme methyltransferase
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (21/4):  ref base (16/9):  total (37/13)

ATGTACCGGACAACCCAGTCCGTCGTAAACAGTTGGGTGGCGGCAGGGATGTCGTTTGTCTTGATGGGACCACCATTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGGAGATATTGATTTAGCGCGTTCGC  .  NC_013198/2164710‑2165005
                                                                                                                                                     |                                                                                                                                                  
aTGTACCGGACAACCCAGTCCGTCGTAAACAGTTGGGTGGCGGCAGGGATGTCGTTTGTCTTGATGGGACCACCATTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGaccac                                                                                                                                                   >  2:979482/1‑151
       ggACAACCCAGTCCGTCGTAAACAGTTGGGTGGCGGCAGGGATGTCGTTTGTCTTGATGGGACCACCATTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTc                                                                                                                                            >  2:844540/1‑151
        ggcAACCCAGTCCGTCGTAAACAGTTGGGTGGCGGCAGGGATGTCGTTTGTCTTGATGGGACCACCATTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCt                                                                                                                                           <  2:991284/149‑1
        gACAACCCAGTCCGTCGTAAACAGTTGGGTGGCGGCAGGGATGTCGTTTGTCTTTATAGGACCACCATTGATGTTAACTACTTAATCGTGCGGTTCAATGTTGTAATACTGATTCAGCCACCCAATAATCTCAACCTGACCGCCCTTGTcc                                                                                                                                           <  2:1037420/151‑2
            aCCCAGTCCGTCGTAAACAGTTGGGTGGCGGCAGGGATGTCGTTTGTCTTGATGGGACCACCATTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACAt                                                                                                                                       >  2:892676/1‑151
            aCCCAGTCCGTCGTAAACAGTTGGGTGGCGGCAGGGATGTCGTTTGTCTTGATGGGACCACCATTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACAt                                                                                                                                       >  1:393323/1‑151
             cccAGTCCGTCGTAAACAGTTGGGTGGCGGCAGGGATGTCGTTTGTCTTGATGGGACCACCATTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTccacgtc                                                                                                                                      <  2:590465/151‑7
                     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:25802/1‑151
                                                                                                                                  cAATAATCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGGAGATATTg                 >  2:180367/1‑151
                                                                                                                                  cAATAATCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGGAGATATTg                 >  1:611361/1‑151
                                                                                                                                      aatCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGGAGATATTGAttt             <  1:509124/151‑1
                                                                                                                                       atCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGGAGATATTGATTTa            >  1:923809/1‑151
                                                                                                                                       atCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGGAGATATTGATTTa            >  2:205756/1‑151
                                                                                                                                       atCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGGAGATATTGATTTa            >  2:226328/1‑151
                                                                                                                                           cAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGGAGATATTGATTTAgcgc        >  1:838342/1‑151
                                                                                                                                              ccTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGATCAGGGAGATATTGATTTAgcacgct     >  2:80666/1‑146
                                                                                                                                                tGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGGAGATATTGATTTATCGCGTTCg   >  2:932217/1‑151
                                                                                                                                                 gACCACCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGGAGATATTGATTTAGCGCGTTCgc  >  1:303389/1‑151
                                                                                                                                                     |                                                                                                                                                  
ATGTACCGGACAACCCAGTCCGTCGTAAACAGTTGGGTGGCGGCAGGGATGTCGTTTGTCTTGATGGGACCACCATTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGGAGATATTGATTTAGCGCGTTCGC  .  NC_013198/2164710‑2165005

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 27 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG