Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,103,953	G→A	C314Y (TGC→TAC)	mutY →	adenine DNA glycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,103,953	0	G	A	100.0%	106.6 / NA	31	C314Y (TGC→TAC)	mutY	adenine DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/17);  total (14/17)

CCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGG  >  NC_000913/3103875‑3104024                                                                               |                                                                        ccGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTACa                                                                       >  2:539863/1‑81 (MQ=255)  cgcgTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTACAt                                                                      <  1:539812/81‑1 (MQ=255)      tttCGGCATACCTTCAGCTATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTACATGGAt                                                                  >  2:539867/1‑81 (MQ=255)       ttCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTACATGGAt                                                                  >  1:539864/1‑80 (MQ=255)               aCCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTACATGGATGAAGGCAAt                                                         >  1:539870/1‑81 (MQ=255)                ccTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTACATGGATGAAGGCAAt                                                         <  1:539813/80‑1 (MQ=255)                ccTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTACATGGATGAAGGCAAt                                                         <  2:539814/80‑1 (MQ=255)                ccTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTACATGGATGAAGGCAATg                                                        <  2:539815/81‑1 (MQ=255)                 cTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTACATGGATGAAGGCAATgc                                                       >  1:539874/1‑81 (MQ=255)                           ttCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTACATGGATGAAGGCAATGCGCTCTGGta                                              >  2:539866/1‑80 (MQ=255)                           ttCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTACATGGATGAAGGCAATGCGCTCTGGta                                              >  2:539876/1‑80 (MQ=255)                            tCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTACATGGATGAAGGCAATGCGCTCTGGtat                                             <  2:539816/80‑1 (MQ=255)                                    gATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTACATGGTTGAAGGCAATGCGCTCTGGTATAACTTAgcg                                    >  2:539868/1‑81 (MQ=255)                                    gATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTACATGGATGAAGGCAATGCGCTCTGGTATAACTTAgc                                     <  1:539818/80‑1 (MQ=255)                                    gATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTACATGAATGAAGGCAATGCGCTCTGGTATAACTTGgcg                                    <  2:539817/81‑1 (MQ=255)                                          gtgCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTACATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAAccg                              >  1:539877/1‑81 (MQ=255)                                          gtgCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTACATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAAcc                               <  1:539819/80‑1 (MQ=255)                                           tgCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTACATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAAccgc                             <  1:539820/81‑1 (MQ=255)                                           tgCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTACATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAAccg                              <  1:539821/80‑1 (MQ=255)                                              cTATGTGGCTTCCCGTGTCGTCATTCACCGGCTACATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGt                          >  1:539878/1‑81 (MQ=255)                                                 tgtgGCTTCCCGTGTCGTCATTCACCGGCTACATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAg                       <  2:539822/81‑1 (MQ=255)                                                 tgtgGCTTCCCGTGTCGTCATTCACCGGCTACATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAg                       <  1:539823/81‑1 (MQ=255)                                                   tgGCTTCCCGTGTCGTCATTCACCGGCTACATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGtt                     <  2:539824/81‑1 (MQ=255)                                                     gCTTCCCGTGTCGTCATTCACCGGCTACATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTgg                   >  2:539881/1‑81 (MQ=255)                                                           cGTGTCGTCATTCACCGGCTACATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGc             <  2:539825/81‑1 (MQ=255)                                                            gtgtCGTCATTCACCGGCTACATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCg            >  1:539880/1‑81 (MQ=255)                                                              gtcgtcATTCACCGGCTACATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGc          >  2:539892/1‑81 (MQ=255)                                                                 gtcATTCACCGGCTACATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTcc       >  2:539884/1‑81 (MQ=255)                                                                  tcATTCACCGGCTACATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTccc      <  2:539826/81‑1 (MQ=255)                                                                     ttCACCGGCTACATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTg   <  1:539827/81‑1 (MQ=255)                                                                       cACCGGCTACATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTgg  <  1:539828/80‑1 (MQ=255)                                                                               |                                                                        CCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGG  >  NC_000913/3103875‑3104024

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