Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,103,479	G→A	G156E (GGG→GAG)	mutY →	adenine DNA glycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,103,479	0	G	A	100.0%	177.0 / NA	50	G156E (GGG→GAG)	mutY	adenine DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/26);  total (24/26)

GGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTC  >  NC_000913/3103357‑3103615                                                                                                                           |                                                                                                                                         gggcgtcgggcgtTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCgaga                                                                                                                         <  1:459034/140‑1 (MQ=255)  ggcgtcgggcgtTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGaa                                                                                                                        >  1:593018/1‑140 (MQ=255)    cgtcgggcgtTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATa                                                                                                                      >  1:417280/1‑140 (MQ=255)    cgtcgggcgtTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATa                                                                                                                      <  2:787438/140‑1 (MQ=255)             ttCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGa                                                                                                             >  2:1129497/1‑140 (MQ=255)             ttCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGa                                                                                                             >  1:446017/1‑140 (MQ=255)                    gCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGCGAATAAATTATGGAGTTTgag                                                                                                      <  2:829094/140‑1 (MQ=255)                       ggCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTgagcga                                                                                                   >  1:607156/1‑140 (MQ=255)                       ggCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTgag                                                                                                      <  1:581093/137‑1 (MQ=255)                       ggCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTgag                                                                                                      >  2:581093/1‑137 (MQ=255)                        gcgcGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTgagcgag                                                                                                  >  2:848813/1‑140 (MQ=255)                                     gCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCcg                                                                                     <  2:869045/140‑1 (MQ=255)                                        ttCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCgg                                                                                  <  2:740435/140‑1 (MQ=255)                                        ttCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCgg                                                                                  <  1:185747/140‑1 (MQ=255)                                         tctcTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGt                                                                                 <  2:584808/140‑1 (MQ=255)                                                gTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTg                                                                          >  2:833178/1‑140 (MQ=255)                                                  aaGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGa                                                                        <  2:446017/140‑1 (MQ=255)                                                   aGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGaa                                                                       >  1:1114447/1‑140 (MQ=255)                                                    gCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAAc                                                                      >  1:228454/1‑140 (MQ=255)                                                     cACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACg                                                                     >  1:883567/1‑140 (MQ=255)                                                       cTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGt                                                                   <  1:848813/140‑1 (MQ=255)                                                         ttCCGATTCTCGACGGAAACGTCAAGCGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGATTAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCACGGTTGGCGTGGAACGGttt                                                                 >  2:492303/1‑140 (MQ=255)                                                          tCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGGGCGAGCAGGTGACGCCCGCGGTTGGTGTGGAACGGTTTa                                                                >  2:350758/1‑140 (MQ=255)                                                          tCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTa                                                                >  1:809910/1‑140 (MQ=255)                                                          tCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTa                                                                >  2:64126/1‑140 (MQ=255)                                                               ttCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAg                                                           >  1:549961/1‑140 (MQ=255)                                                               ttCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAg                                                           <  1:64126/140‑1 (MQ=255)                                                                   cGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCga                                                       <  2:557055/140‑1 (MQ=255)                                                                   cGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCga                                                       <  1:1129497/140‑1 (MQ=255)                                                                       ggTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCgatgat                                                   <  2:883567/140‑1 (MQ=255)                                                                         tAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATgg                                                 >  1:252395/1‑140 (MQ=255)                                                                           aCGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGAt                                               <  1:390178/140‑1 (MQ=255)                                                                           aCGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGAt                                               <  2:124950/140‑1 (MQ=255)                                                                            cGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGAtt                                              >  2:763882/1‑140 (MQ=255)                                                                                aaaCGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGATTTggg                                          >  1:214422/1‑140 (MQ=255)                                                                                 aaCGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGGGCGAGCAGGTGACGCCCGCGGTTGGTGTGGAACGGTTTAATCAGGCGATGATGGATTTGGGt                                         <  1:350758/140‑1 (MQ=255)                                                                                  aCGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGATTTGGGTg                                        >  2:947296/1‑140 (MQ=255)                                                                                    gcgTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGATTTGGGTGCg                                      <  1:833178/140‑1 (MQ=255)                                                                                     cgTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGATTTGGGTGCga                                     >  2:109660/1‑140 (MQ=255)                                                                                              gcgcTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGATTTGGGt                                         <  1:644121/127‑1 (MQ=255)                                                                                              gcgcTGCTATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGATTTGGGt                                         >  2:644121/1‑127 (MQ=255)                                                                                                    ctATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGATTTGGGTGCGATGATTTGTACGCGCt                      <  1:947296/140‑1 (MQ=255)                                                                                                    ctATGCTGTAAGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGATTTGGGTGCGATGATTTGTACGCGCt                      <  1:109660/140‑1 (MQ=255)                                                                                                             aaGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTAGCGTGGAACGGTTTAATCAGGCGATGATGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGa             <  2:549961/140‑1 (MQ=255)                                                                                                              aGCGGCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGaa            >  2:601262/1‑140 (MQ=255)                                                                                                                  gCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGt        >  1:1124385/1‑140 (MQ=255)                                                                                                                  gCTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGt        <  1:763882/140‑1 (MQ=255)                                                                                                                   cTGGCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGtt       >  1:720966/1‑140 (MQ=255)                                                                                                                     ggCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCg     <  1:674982/140‑1 (MQ=255)                                                                                                                     ggCCTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCg     <  2:1114447/140‑1 (MQ=255)                                                                                                                        cTGAGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGctc  <  2:252395/140‑1 (MQ=255)                                                                                                                           |                                                                                                                                         GGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTC  >  NC_000913/3103357‑3103615

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