Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_013198	2,158,710	0	G→A	100.0%	9.5	41	V25V (GTG→GTA)	LGG_02095	type II restriction enzyme methylase subunit
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (8/13):  ref base (11/9):  total (19/22)

TAACGCACCGGCTAGTATTCTAGCCTCTAGTATAACTAACAGTTAAAGCCGTGGCTATAATCCTACTTTCTTTTTGTTCACTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGACGTGGACAAGGGCGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACAATACCGAGCCCAAAGACAGTGCCTTTAAGAAAAAGAAGTATTTGACTTCTGATATTCCCCCAGTGACGCAACTTTTTACGCCAGACTGGATCGTCAAGTAC  .  NC_013198/2158563‑2158860                                                                                                                                                    |                                                                                                                                                       tAACGCACCGGCTAGTATTCTAGCCTCTAGTATAACTAACAGTTAAAGCCGTGGCTATAATCCTACTTTCTTTTTGTTCACTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGATGTAGAc                                                                                                                                                     <  2:178522/151‑1  aaCGCACCGGCTAGTATTCTAGCCTCTAGTATAACTAACAGTTAAAGCCGTGGCTATAATCCTACTTTCTTTTTGTTCACTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGACGTGGACa                                                                                                                                                    <  1:1172711/151‑1  aaCGCACCGGCTAGTATTCTAGCCTCTAGTATAACTAACAGTTAAAGCCGTGGCTATAATCCTACTTTCTTTTTGTTCACTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGACGTGGACa                                                                                                                                                    <  2:490443/151‑1      cACCGGCTATTATTCTAGCCTCTAGTATAACTAACAGTTAAAGCCGTGGCTATAATCCTACTTTCTTTTTGTTCACTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGATGTAGACAAggg                                                                                                                                                <  2:481671/151‑1      cACCGGCTAGTATTCTAGCCTCTAGTATAACTAACAGTTAAAGCCGTGGCTATAATCCTACTTTCTTTTTGTTCACTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGATGTAGACAAggg                                                                                                                                                <  1:601662/151‑1      cACCGGCTAGTATTCTAGCCTCTAGTATAACTAACAGTTAAAGCCGTGGCTATAATCCTACTTTCTTTTTGTTCACTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGATGTAGACAAggg                                                                                                                                                <  1:559590/151‑1      acagggCTAGTATTCTAGCCTCTAGTATAACTAACAGTTAAAGCCGTGGCTATAATCCTACTTTCTTTTTGTTCACTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGATGTAGACAAggg                                                                                                                                                <  2:1245371/147‑1          ggCTAGTATTCTAGCCTCTAGTATAACTAACAGTTAAAGCCGTGGCTATAATCCTACTTTCTTTTTGTTCACTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGATGTAGACAAGGGctgt                                                                                                                                            >  1:1245737/1‑148               gTATTCTAGCCTCTAGTATAACTAACAGTTAAAGCCGTGGCTATAATCCTACTTTCTTTTTGTTCACTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGATGTAGACAAGGGTGGTCAGGt                                                                                                                                       >  2:961701/1‑151                    cTAGCCTCTAGTATAACTAACAGTTAAAGCCGTGGCTATAATCCTACTTTCTTTTTGTTCACTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGATGTAGACAAGGGTGGTCAGGTTgaga                                                                                                                                  <  1:416055/151‑1                     tAGCCTCTAGTATAACTAACAGTTAAAGCCGTGGCTATAATCCTACTTTCTTTTTGTTCACTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGATGTAGACAAGGGTGGTCAGGTTGAGat                                                                                                                                 <  1:73167/151‑1                     tAGCCTCTAGTATAACTAACAGTTAAAGCCGTGGCTATAATCCTACTTTCTTTTTGTTCACTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGATGTAGACAAGGGTGGTCAGGTTGAGat                                                                                                                                 <  1:499124/151‑1                     tAGCCTCTAGTATAACTAACAGTTAAAGCCGTGGCTATAATCCTACTTTCTTTTTGTTCACTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGATGTAGACAAGGGTGGTCAGGTTGAGat                                                                                                                                 <  1:295091/151‑1                      aGCCTCTAGTATAACTAACAGTTAAAGCCGTGGCTATAATCCTACTTTCTTTTTGTTCACTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGCTATTCGTGAAGATGATTTTGATGTAGACAAGGGTGGTGAGGTTGAGatt                                                                                                                                >  2:983659/1‑151                                      aacAGTTAAAGCCGTGGCTATAATCCTACTTTCTTTTTGTTCACTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGATGTAGACAAGGGTGGTCAGGTTGAGATTATTGGGTGGCTGTATc                                                                                                                <  2:1067084/151‑1                                          gTTAAAGCCGTGGCTATAATCCTACTTTCTTTTTGTTCACTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGATGTAGACAAGGGTGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTa                                                                                                            <  1:76193/151‑1                                                gacagggCTATAATCCTACTTTCTTTTTGTTCACTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGATGTAGACAAGGGTGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACaa                                                                                                      <  1:770522/146‑1                                                gCCGTGGCTATAATCCTACTTTCTTTTTGTTCACTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGACGTGGACAAGGGCGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACaa                                                                                                      <  1:440130/151‑1                                                 ccGTGGCTATAATCCTACTTTCTTTTTGTTCACTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGATGTAGACAAGGGTGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACaac                                                                                                     <  2:680808/151‑2                                                     ggCTATAATCCTACTTTCTTTTTGTTCACTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGATGTAGACAAGGGTGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACaactgtc                                                                                                 >  2:770340/1‑146                                                     ggCTATAATCCTACTTTCTTTTTGTTCACTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGATGTAGACAAGGGTGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACaacattg                                                                                                 >  1:844752/1‑146                                                       cTATAATCCTACTTTCTTTTTGTTCACTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGATGTAGACAAGGGTGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACaacattgac                                                                                               >  2:407952/1‑144                                                                catgcaattattGTTCACTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGACGTGGACAAGGGCGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACAATACCGAGCCCAAAGAc                                                                                      <  1:179375/141‑1                                                                catgcaattattGTTCACTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGACGTGGACAAGGGCGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACAATACCGAGCCCAAAGAc                                                                                      <  2:140401/141‑1                                                                                aCTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGACGTGGACACGGGCGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACAATACCGAGCCCAAAGACAGTGCCTTTAAGaaaa                                                                      >  2:354193/1‑151                                                                                 cTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGACGTGGACAAGGGCGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACAATACCGAGCCCAAAGACAGTGCCTTTAAGaaaaa                                                                     >  1:1040128/1‑151                                                                                 cTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGACGTGGACAAGGGCGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACAATACCGAGCCCAAAGACAGTGCCTTTAAGaaaaa                                                                     >  2:1135/1‑151                                                                                 cTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGACGTGGACAAGGGCGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACAATACCGAGCCCAAAGACAGTGCCTTTAAGaaaaa                                                                     >  1:1074734/1‑151                                                                                     aaGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGACGTGGACAAGGGCGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACAATACCGAGCCCAAAGACAGTGCCTTTAAGAAAaagaag                                                                 >  1:332617/1‑151                                                                                             cGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGACGTGGACAAGGGCGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACAATACCGAGCCCAAAGACAGTGCCTTTAAGAAAAAGAAGTATTTGAc                                                         <  1:248949/151‑1                                                                                             cGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGACGTGGACAAGGGCGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACAATACCGAGCCCAAAGACAGTGCCTTTAAGAAAAAGAAGTATTTGAc                                                         <  2:478959/151‑1                                                                                                ttGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGACGTGGACAAGGGCGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACAATACCGAGCCCAAAGACAGTGCCTTTAAGAAAAAGAAGTATTTGACTTc                                                      <  2:636665/151‑1                                                                                                    ggTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGACGTGGACAAGGGCGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACAATACCGAGCCCAAAGACAGTGCCTTTAAGAAAAAGAAGTATTTGACTTCTGat                                                  >  1:608968/1‑151                                                                                                    ggTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGACGTGGACAAGGGCGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACAATACCGAGCCCAAAGACAGTGCCTTTAAGAAAAAGAAGTATTTGACTTCTGat                                                  >  2:1241415/1‑151                                                                                                         ttCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGACGTGGACAAGGGCGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACAATACCGAGCCCAAAGACAGTGCCTTTAAGAAAAAGAAGTATTTGACTTCTGATATTcc                                             >  2:326421/1‑151                                                                                                         ttCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGACGTGGACAAGGGCGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACAATACCGAGCCCAAAGACAGTGCCTTTAAGAAAAAGAAGTATTTGACTTCTGATATTcc                                             >  2:980957/1‑151                                                                                                          tCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGACGTGGACAAGGGCGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACAATACCGAGCCCAAAGACAGTGCCTTTAAGAAAAAGAAGTATTTGACTTCTGATATTccc                                            <  1:293209/151‑1                                                                                                                       gatgatATTCGTGAAGATGATTTTGACGTGGACAAGGGCGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACAATACCGAGCCCAAAGACAGTGCCTTTAAGAAAAAGAAGTATTTGACTTCTGATATTCCCCCAGTGACGCAAc                               >  1:963700/1‑151                                                                                                                         tgatATTCGTGAAGATGATTTTGATGTAGACAAGGGTGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACAATACCGAGCCCAAAGACAGTGCCTTTAAGGAAAAGAAGTATTTGACTTCTGATATTCCCCCAGTGACGCAACtt                             >  2:1010145/1‑151                                                                                                                                          atTTTGATGTAGACAAGGGCGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACAATACCGAGCCCAAAGACAGTGCCTTTAAGAAAAAGAAGTATTTGACTTCTGATATTCCCCCAGTGACGCAACTTTTTACGCCAGACTGGAt            >  2:54364/1‑151                                                                                                                                          atTTTGACGTGGACAAGGGCGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACAATACCGAGCCCAAAGACAGTGCCTTTAAGAAAAAGAAGTATTTGACTTCTGATATTCCCCCAGTGACGCAACTTTTTACGCCAGACTGGAt            >  2:792839/1‑151                                                                                                                                                    ggACAAGGGCGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACAATACCGAGCCCAAAGACAGTGCCTTTAAGAAAAAGAAGTATTTGACTTCTGATATTCCCCCAGTGACGCAACTTTTTACGCCAGACTGGATCGTCAAGTAc  >  1:684897/1‑151                                                                                                                                                    ggACAAGGGCGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACAATACCGAGCCCAAAGACAGTGCCTTTAAGAAAAAGAAGTATTTGACTTCTGATATTCCCCCAGTGACGCAACTTTTTACGCCAGACTGGATCGTCAAGTAc  <  1:11029/151‑1                                                                                                                                                    ggACAAGGGCGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACAATACCGAGCCCAAAGACAGTGCCTTTAAGAAAAAGAAGTATTTGACTTCTGATATTCCCCCAGTGACGCAACTTTTTACGCCAGACTGGATCGTCAAGTAc  >  2:944518/1‑151                                                                                                                                                    |                                                                                                                                                       TAACGCACCGGCTAGTATTCTAGCCTCTAGTATAACTAACAGTTAAAGCCGTGGCTATAATCCTACTTTCTTTTTGTTCACTCCAAGCTATTCGGTTGGGGTCATTCGTGATCTTGTGGATGATATTCGTGAAGATGATTTTGACGTGGACAAGGGCGGTCAGGTTGAGATTATTGGGTGGCTGTATCAGTATTACAATACCGAGCCCAAAGACAGTGCCTTTAAGAAAAAGAAGTATTTGACTTCTGATATTCCCCCAGTGACGCAACTTTTTACGCCAGACTGGATCGTCAAGTAC  .  NC_013198/2158563‑2158860 Base quality scores:  ATCG  < 3 ≤  ATCG  < 23 ≤  ATCG  < 30 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG