Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,706,510	A→G	E20G (GAG→GGG)	tag →	DNA‑3‑methyladenine glycosylase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,706,510	0	A	G	100.0%	164.0 / NA	44	E20G (GAG→GGG)	tag	DNA‑3‑methyladenine glycosylase I
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/18);  total (26/18)

TTAAATATTGACTTTTTCACCGATGCGTCAAGAAAAGCGGCTGAAATTTTTACGATCGGGTACATAGCGAGGGAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGAGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGT  >  NZ_CP009273/3706374‑3706638                                                                                                                                         |                                                                                                                                 ttAAATATTGACTTTTTCACCGATGCGTCAAGAAAAGCGGCTGAAATTTTTACGATCGGGTACATAGCGAGGGAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGt                                                                                                                                >  2:420521/1‑139 (MQ=255) ttAAATATTGACTTTTTCACCGATGCGTCAAGAAAAGCGGCTGAAATTTTTACGATCGGGTACATAGCGAGGGAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGt                                                                                                                                >  1:918072/1‑139 (MQ=255)    aaTATTGACTTTTTCACCGATGCGTCAAGAAAAGCGGCTGAAATTTTTACGATCGGGTACATAGCGAGGGAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAAtgggtggg                                                                                                                             >  2:349450/1‑139 (MQ=255)    aaTATTGACTTTTTCACCGATGCGTCAAGAAAAGCGGCTGAAATTTTTACGATCGGGTACATAGCGAGGGAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAAtgggtggg                                                                                                                             >  2:333545/1‑139 (MQ=255)         tGACTTTTTCACCGATGCGTCAAGAAAAGCGGCTGAAATTTTTACGATCGGGTCCATAGCGAGGGAAAGTATGGAGCGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGTGGGGCGTg                                                                                                                        >  2:116109/1‑139 (MQ=255)         tGACTTTTTCACCGATGCGTCAAGAAAAGCGGCTGAAATTTTTACGATCGGGTACATAGCGAGGGAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGTGGGGCGTg                                                                                                                        >  2:991627/1‑139 (MQ=255)         tGACTTTTTCACCGATGCGTCAAGAAAAGCGGCTGAAATTTTTACGATCGGGTACATAGCGAGGGAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGTGGGGCGTg                                                                                                                        >  1:712590/1‑139 (MQ=255)           aCTTTTTCACCGATGCGTCAAGAAAAGCGGCTGAAATTTTTACGATCGGGTACATAGCGGGGGAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGTGGGGAGTGcc                                                                                                                      >  2:594976/1‑139 (MQ=255)                    ccGATGCGTCAAGAAAAGCGGCTGAAATTTTTACGATCGGGTACATAGCGAGGGAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGa                                                                                                             <  1:743496/139‑1 (MQ=255)                                  aaaGCGGCTGAAATTTTTACGATCGGGTACATAGCGAGGGAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGt                                                                                               <  1:1230146/139‑1 (MQ=255)                                        gCTGAAATTTTTACGATCGGGTACATAGCGAGGGAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGaaa                                                                                         >  1:1097405/1‑139 (MQ=255)                                            aaaTTTTTACGATCGGGTACATAGCGAGGGAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGAt                                                                                     <  1:535620/139‑1 (MQ=255)                                            aaaTTTTTACGATCGGGTACATAGCGAGGGAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGAt                                                                                     <  2:1069125/139‑1 (MQ=255)                                                           ggTACATAGCGAGGGAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGgca                                                                      <  2:436264/139‑1 (MQ=255)                                                                aTAGCGAGGGAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAgg                                                                 <  2:751868/139‑1 (MQ=255)                                                                 tAGCGAGGGAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGc                                                                <  1:349450/139‑1 (MQ=255)                                                                  aGCGAGGGAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCt                                                               >  1:69768/1‑139 (MQ=255)                                                                    cGAGGGAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTg                                                              >  2:828055/1‑138 (MQ=255)                                                                    cGAGGGAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTg                                                              >  1:360437/1‑138 (MQ=255)                                                                     gAGGGAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGa                                                            <  2:1053827/139‑1 (MQ=255)                                                                     gAGGGAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGa                                                            <  1:490482/139‑1 (MQ=255)                                                                     gAGGGAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGa                                                            <  1:700701/139‑1 (MQ=255)                                                                      aGGGAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGAt                                                                                     <  1:839543/113‑1 (MQ=255)                                                                      aGGGAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGAt                                                                                     >  2:839543/1‑113 (MQ=255)                                                                      aGGGAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGa                                                            <  1:890827/138‑1 (MQ=255)                                                                      aGGGAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGa                                                            >  2:890827/1‑138 (MQ=255)                                                                      aGGGAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGAt                                                           >  2:561201/1‑139 (MQ=255)                                                                         gAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGATTAt                                                        >  1:590625/1‑139 (MQ=255)                                                                                        tGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGa                                                                                                             <  2:25732/71‑1 (MQ=255)                                                                                        tGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGa                                                                                                             >  1:25732/1‑71 (MQ=255)                                                                                        tGCGGCTGGGTGAGTCAGGACACGCTTTATATTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGATTATCGTGGAACACCGTcc                                         >  1:1281421/1‑139 (MQ=255)                                                                                                  tGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGATTATCGTGGATCACCGTCCTCAAAAAAcg                               >  2:363098/1‑139 (MQ=255)                                                                                                      tCAGGACCCGCTTTATATTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGa                            >  2:604375/1‑138 (MQ=255)                                                                                                                cgttATATAGCCAATCATGATAATGGGGGGGGCGTGCCTGTAAATGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGTCTGGGTTATCGTGGATAACCGTCCTCAAAAAACGCGAAAACTATcgcg                 >  2:83015/3‑139 (MQ=255)                                                                                                                cTTTATATTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTg                                                                                  <  2:717437/74‑1 (MQ=255)                                                                                                                cTTTATATTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTg                                                                                  >  1:717437/1‑74 (MQ=255)                                                                                                                   tataTTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCgcgc                <  1:991627/137‑1 (MQ=255)                                                                                                                   tataTTGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCgcgc                <  1:453879/137‑1 (MQ=255)                                                                                                                       ttGCCTACCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCtt          >  1:301348/1‑139 (MQ=255)                                                                                                                           cTACCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTtcat      <  1:469919/139‑1 (MQ=255)                                                                                                                             aCCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTtcat      >  1:211768/1‑137 (MQ=255)                                                                                                                             aCCATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTtcat      >  1:790446/1‑137 (MQ=255)                                                                                                                               cATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGt  >  2:14557/1‑139 (MQ=255)                                                                                                                               cATGATAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGt  >  2:376594/1‑139 (MQ=255)                                                                                                                                   aTAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAgg                                                                 <  2:1018164/72‑1 (MQ=255)                                                                                                                                   aTAATGGGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAgg                                                                 >  1:1018164/1‑72 (MQ=255)                                                                                                                                         |                                                                                                                                 TTAAATATTGACTTTTTCACCGATGCGTCAAGAAAAGCGGCTGAAATTTTTACGATCGGGTACATAGCGAGGGAAAGTATGGAACGTTGCGGCTGGGTGAGTCAGGACCCGCTTTATATTGCCTACCATGATAATGAGTGGGGCGTGCCTGAAACTGACAGTAAAAAACTGTTCGAAATGATCTGCCTTGAAGGGCAGCAGGCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGT  >  NZ_CP009273/3706374‑3706638

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