Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,394,163	G→A	L111L (CTG→TTG)	ogt ←	methylated‑DNA‑‑[protein]‑cysteine S‑methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,394,163	0	G	A	100.0%	295.9 / NA	74	L111L (CTG→TTG)	ogt	methylated‑DNA‑‑[protein]‑cysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/46);  total (28/46)

GCAAGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGC  >  NZ_CP009273/1394077‑1394247                                                                                       |                                                                                     gCAAGGTACGGCGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTg                                                                                   <  1:2428562/90‑1 (MQ=255)   aaGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGc                                                                                 <  1:2192221/90‑1 (MQ=255)    aGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGcc                                                                                >  1:2614520/1‑90 (MQ=255)    aGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGcc                                                                                >  1:1757790/1‑90 (MQ=255)    aGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGcc                                                                                >  1:694113/1‑90 (MQ=255)    aGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGcc                                                                                >  2:1465715/1‑90 (MQ=255)       tACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTa                                                                             <  2:1864163/90‑1 (MQ=255)        acgacgATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTaa                                                                            >  2:954915/1‑90 (MQ=255)        acgacgATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTaa                                                                            >  1:1241255/1‑90 (MQ=255)           acgaTGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGc                                                                         <  1:1465715/90‑1 (MQ=255)           acgaTGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGc                                                                         >  2:1470623/1‑90 (MQ=255)           acgaTGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGc                                                                         <  1:1183983/90‑1 (MQ=255)              aTGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCAtt                                                                      >  1:411654/1‑90 (MQ=255)              aTGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCAtt                                                                      >  2:2082008/1‑90 (MQ=255)              aTGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAACTGCTCAGCCAATTGGCCGTAATGCAtt                                                                      >  2:1826410/1‑90 (MQ=255)               tGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTa                                                                     >  2:1929554/1‑90 (MQ=255)               tGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTa                                                                     >  1:249836/1‑90 (MQ=255)                gCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTAc                                                                    >  1:326705/1‑90 (MQ=255)                gCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTAc                                                                    >  1:111968/1‑90 (MQ=255)                gCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTAc                                                                    >  1:2130829/1‑90 (MQ=255)                gCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTAc                                                                    >  1:1476679/1‑90 (MQ=255)                   gATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTg                                                                 <  2:2612577/90‑1 (MQ=255)                   gATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTg                                                                 <  1:817259/90‑1 (MQ=255)                   gATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTg                                                                 <  1:1593023/90‑1 (MQ=255)                     tGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGcc                                                               <  2:862591/90‑1 (MQ=255)                     tGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGcc                                                               <  1:2854026/90‑1 (MQ=255)                     tGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGcc                                                               <  1:1552000/90‑1 (MQ=255)                        gATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGc                                                            >  1:1205630/1‑90 (MQ=255)                        gATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGc                                                            >  1:909546/1‑90 (MQ=255)                         aTTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCa                                                           <  1:842507/90‑1 (MQ=255)                          ttCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAg                                                          <  1:1826410/90‑1 (MQ=255)                          ttCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAg                                                          <  1:1630/90‑1 (MQ=255)                             gATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAgggg                                                       >  2:2889755/1‑90 (MQ=255)                              aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGa                                                      <  1:109688/90‑1 (MQ=255)                              aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGa                                                      <  1:2970504/90‑1 (MQ=255)                                ccGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATa                                                    <  2:326705/90‑1 (MQ=255)                                ccGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATa                                                    <  2:2694248/90‑1 (MQ=255)                                ccGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATa                                                    <  2:2563413/90‑1 (MQ=255)                                 cGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAg                                                   <  2:111968/90‑1 (MQ=255)                                 cGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAg                                                   <  2:1205630/90‑1 (MQ=255)                                 cGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAg                                                   <  2:2403725/90‑1 (MQ=255)                                 cGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAg                                                   <  2:2130829/90‑1 (MQ=255)                                 cGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAg                                                   <  2:1757790/90‑1 (MQ=255)                                   tttGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTg                                                 <  1:3188/90‑1 (MQ=255)                                       ccGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTa                                             <  1:476244/90‑1 (MQ=255)                                          cACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAgtg                                          <  2:2096138/90‑1 (MQ=255)                                          cACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAgtg                                          <  2:3022770/90‑1 (MQ=255)                                                      gcgccgcgccAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACtt                              >  2:2254076/1‑90 (MQ=255)                                                         ccgcgccAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTcgc                           <  2:694113/90‑1 (MQ=255)                                                          cgcgccAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTcgcg                          >  1:2054240/1‑90 (MQ=255)                                                           gcgccAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCgcgc                         >  2:837617/1‑90 (MQ=255)                                                           gcgccAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCgcgc                         <  2:821390/90‑1 (MQ=255)                                                           gcgccAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCgcgc                         <  2:1476679/90‑1 (MQ=255)                                                                   aCGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATgg                 >  1:120532/1‑90 (MQ=255)                                                                      gCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAg                                                          >  1:309871/1‑46 (MQ=255)                                                                      gCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAg                                                          <  2:309871/46‑1 (MQ=255)                                                                      gCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATa                                                    <  2:616209/52‑1 (MQ=255)                                                                      gCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATa                                                    >  1:616209/1‑52 (MQ=255)                                                                        ccAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTcc            <  2:2880403/90‑1 (MQ=255)                                                                        ccAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTcc            <  2:8735/90‑1 (MQ=255)                                                                        ccAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTcc            <  2:411654/90‑1 (MQ=255)                                                                        ccAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTcc            <  2:2919048/90‑1 (MQ=255)                                                                        ccAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTcc            <  2:1409855/90‑1 (MQ=255)                                                                              gCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCg      >  1:1434012/1‑90 (MQ=255)                                                                               cTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGt     <  1:2889755/90‑1 (MQ=255)                                                                               cTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGt     <  1:2082008/90‑1 (MQ=255)                                                                                tCAGCCAATTGGCCGTAATGCATTACCTGCCCGCa                                                           <  1:953729/35‑1 (MQ=255)                                                                                tCAGCCAATTGGCCGTAATGCATTACCTGCCCGCa                                                           >  2:953729/1‑35 (MQ=38)                                                                                tCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTcccc          <  2:1212377/84‑1 (MQ=255)                                                                                tCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTcccc          >  1:1212377/1‑84 (MQ=255)                                                                                  aGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGc  <  1:2547397/90‑1 (MQ=255)                                                                                  aGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGc  <  1:2254076/90‑1 (MQ=255)                                                                                  aGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGc  <  2:816585/90‑1 (MQ=255)                                                                                  aGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGc  <  1:2917009/90‑1 (MQ=255)                                                                                       |                                                                                     GCAAGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGC  >  NZ_CP009273/1394077‑1394247

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑