Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_013198	2,164,871	A→G	D219D (GAT→GAC)	LGG_02097	type II restriction enzyme methyltransferase

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_013198	2,164,871	0	A→G	100.0%	10.7	47	D219D (GAT→GAC)	LGG_02097	type II restriction enzyme methyltransferase
Reads supporting (aligned to +/- strand):  new base (22/3):  ref base (14/8):  total (36/11)

ACCCAGTCCGTCGTAAACAGTTGGGTGGCGGCAGGGATGTCGTTTGTCTTGATGGGACCACCATTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGGAGATATTGATTTAGCGCGTTCGCTTGCTT  .  NC_013198/2164722‑2165011                                                                                                                                                      |                                                                                                                                             aCCCAGTCCGTCGTAAACAGTTGGGTGGCGGCAGGGATGTCGTTTGTCTTGATGGGACCACCATTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACAt                                                                                                                                             >  2:892676/1‑151 aCCCAGTCCGTCGTAAACAGTTGGGTGGCGGCAGGGATGTCGTTTGTCTTGATGGGACCACCATTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACAt                                                                                                                                             >  1:393323/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                                                                                                <  1:1006582/151‑1                                              gTCTTGATGGGACCACCATTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGAccc                                                                                                >  2:53761/1‑151                                                  ataagagacagccaTTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAAc                                                                                            <  1:274902/140‑1                                                            accaTTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCtt                                                                                  >  1:462999/1‑151                                                             ccaTTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAAcctgtctcttat                                                                                 >  2:274908/1‑141                                                             ccaTTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTg                                                                                 >  1:246241/1‑151                                                             ccaTTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTg                                                                                 >  2:425100/1‑151                                                                ttGGTGTTGACTTCTTTTTCGTGCATGTCAATGTGGTAATACTGATACCGCCACCCAATAATCTCAACCTGACCGCACTTGTCTACATCAAAATCATCTGCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAg                                                                              <  2:362501/151‑1                                                                ttGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAg                                                                              <  2:624568/151‑1                                                                       tGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACaa                                                                       >  1:667079/1‑151                                                                            actTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACaaaaaga                                                                  >  2:377527/1‑146                                                                            actTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACaaaaaga                                                                  >  2:412251/1‑146                                                                               ttATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACaaaaagaaag                                                               <  1:53761/151‑9                                                                                              gtattgtAATACTGATACAGCCACCCAATAATCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGttt                                                >  2:508984/4‑151                                                                                              gtattgtAATACTGATACAGCCACCCAATAATCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGttt                                                >  1:338026/4‑151                                                                                              gtattgtAATACTGATACAGCCACCCAATAATCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGttt                                                >  1:1003256/4‑151                                                                                                 ttgtAATACTGATACAGCCACCCAATAATCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGtttttt                                             >  2:661692/1‑151                                                                                                 ttgtAATACTGATACAGCCACCCAATAATCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGtttttt                                             >  1:244456/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                       >  1:611361/1‑151                                                                                                                       cAATAATCTCAACCTGACCGCCCTTGTCCACGTCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGGAGATATTg                       >  2:25802/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                                                                                                                                                      |                                                                                                                                             ACCCAGTCCGTCGTAAACAGTTGGGTGGCGGCAGGGATGTCGTTTGTCTTGATGGGACCACCATTGATGTTGACTACTTTATCGTGCGGTTCAATGTTGTAATACTGATACAGCCACCCAATAATCTCAACCTGACCACCCTTGTCTACATCAAAATCATCTTCACGAATATCATCCACAAGATCACGAATGACCCCAACCGAATAGCTTGGAGTGAACAATAATTGCATGAAATCGTCTGTTTTTTCGAAAAGTTCAGGGAGATATTGATTTAGCGCGTTCGCTTGCTT  .  NC_013198/2164722‑2165011 Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG