Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,103,116	T→C	L35P (CTC→CCC)	mutY →	adenine DNA glycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,103,116	0	T	C	100.0%	173.1 / NA	50	L35P (CTC→CCC)	mutY	adenine DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/25);  total (25/25)

TCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTGTGGACCGGGC  >  NC_000913/3102981‑3103250                                                                                                                                        |                                                                                                                                       tCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCa                                                                                                                                    <  2:153632/140‑1 (MQ=255)   ttGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGa                                                                                                                                  >  1:664867/1‑140 (MQ=255)     gCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAg                                                                                                                                >  1:772665/1‑140 (MQ=255)          caacaaCAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTAAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATg                                                                                                                           <  2:772665/140‑1 (MQ=255)            acaacaGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGtt                                                                                                                         <  1:28385/140‑1 (MQ=255)            acaacaGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGtt                                                                                                                         <  1:435465/140‑1 (MQ=255)                     aaTTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAAc                                                                                                                >  1:154463/1‑140 (MQ=255)                      aTTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACt                                                                                                               <  1:26420/140‑1 (MQ=255)                       ttCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTc                                                                                                              >  1:279082/1‑140 (MQ=255)                        tCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCATACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCa                                                                                                             >  1:416336/1‑140 (MQ=255)                                cATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGa                                                                                                     <  1:696168/140‑1 (MQ=255)                                  tGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGAcc                                                                                                   >  2:636625/1‑140 (MQ=255)                                   gCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCg                                                                                                  <  2:279082/140‑1 (MQ=255)                                        cGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATc                                                                                             >  1:166913/1‑140 (MQ=255)                                         gTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATcc                                                                                            >  2:444833/1‑140 (MQ=255)                                          tCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATccc                                                                                           >  1:365828/1‑140 (MQ=255)                                            gCAATCTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCt                                                                                         <  1:236561/140‑1 (MQ=255)                                                     aGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGa                                                                                  <  1:77726/138‑1 (MQ=255)                                                     aGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGa                                                                                  >  2:77726/1‑138 (MQ=255)                                                            gTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTAt                                                                         >  2:615909/1‑140 (MQ=255)                                                                  gACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGgcgcg                                                                   <  1:529773/140‑1 (MQ=255)                                                                       gTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGgcgcg                                                                   <  1:318215/135‑1 (MQ=255)                                                                       gTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGgcgcg                                                                   >  2:318215/1‑135 (MQ=255)                                                                       gTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTcc                                                              >  1:612091/1‑140 (MQ=255)                                                                           gATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCATACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGAc                                                          <  2:416336/140‑1 (MQ=255)                                                                                 tACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGAc                                                    >  1:195252/1‑140 (MQ=255)                                                                                 tACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGAc                                                    >  2:77695/1‑140 (MQ=255)                                                                                  aCGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGAcc                                                   >  1:670510/1‑140 (MQ=255)                                                                                     ggCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGAt                                                <  2:366244/140‑1 (MQ=255)                                                                                      gCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGAtc                                               >  2:440684/1‑140 (MQ=255)                                                                                      gCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGAtc                                               >  1:399306/1‑140 (MQ=255)                                                                                      gCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGAtc                                               >  2:3197/1‑140 (MQ=255)                                                                                       cGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGAtct                                              <  1:678765/140‑1 (MQ=255)                                                                                       cGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGAtct                                              >  1:720338/1‑140 (MQ=255)                                                                                            aaCTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCa                                         <  2:601428/140‑1 (MQ=255)                                                                                                ctGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATgc                                     <  2:195252/140‑1 (MQ=255)                                                                                                 tGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATgcg                                    <  1:3197/140‑1 (MQ=255)                                                                                                        gCAAATTGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTAACCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTcg                             >  2:575211/1‑140 (MQ=255)                                                                                                             ttGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGaa                        >  2:506874/1‑140 (MQ=255)                                                                                                              tGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAg                       <  1:440684/140‑1 (MQ=255)                                                                                                              tGACAAGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAg                       <  2:685157/140‑1 (MQ=255)                                                                                                                   aGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTtctc                  <  2:664867/140‑1 (MQ=255)                                                                                                                   aGACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTtctc                  <  2:154463/140‑1 (MQ=255)                                                                                                                    gACGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTcc                 <  2:612091/140‑1 (MQ=255)                                                                                                                     aCGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCa                <  2:720338/140‑1 (MQ=255)                                                                                                                      cGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTtctc                  >  2:688641/1‑137 (MQ=255)                                                                                                                      cGCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTtctc                  <  1:688641/137‑1 (MQ=255)                                                                                                                       gCCCTACAAAGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACt              >  1:443878/1‑140 (MQ=255)                                                                                                                              aaaGTATGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTGTGGAc       <  1:615909/140‑1 (MQ=255)                                                                                                                                   aTGGCCCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTGTGGACCGGGc  >  2:602280/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                       TCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTGTGGACCGGGC  >  NC_000913/3102981‑3103250

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: ‑