Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_013198	2,164,859	0	A→G	100.0%	0.5	50	G223G (GGT→GGC)	LGG_02097	type 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