Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,104,065	G→A	*351* (TAG→TAA)	mutY →	adenine DNA glycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,104,065	0	G	A	100.0%	113.5 / NA	37	*351* (TAG→TAA)	mutY	adenine DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/16);  total (21/16)

CGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAGCGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGCATCTATAACGAG  >  NC_000913/3103938‑3104191                                                                                                                                |                                                                                                                               cgtcATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCg                                                                                                                     <  1:676009/139‑1 (MQ=255) cgtcATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGG‑TTAACGCGTGAGTCGa                                                                                                                    >  2:530540/1‑139 (MQ=255)         aCCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAagag                                                                                                             <  2:393485/139‑1 (MQ=255)                           ggCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGaa                                                                                           >  2:121608/1‑139 (MQ=255)                           ggCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGaa                                                                                           >  1:36665/1‑139 (MQ=255)                                tGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGAtt                                                                                      >  1:86968/1‑139 (MQ=255)                                       tGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAATGCGTGa                                                                                                                         >  1:627346/1‑97 (MQ=255)                                       tGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAATGCGTGa                                                                                                                         <  2:627346/97‑1 (MQ=255)                                                 tAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCa                                                                     <  1:121608/139‑1 (MQ=255)                                                 tAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCa                                                                     >  1:338157/1‑139 (MQ=255)                                                   gcgcAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAAc                                                                   >  2:90541/1‑139 (MQ=255)                                                    cgcAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACg                                                                  >  2:70401/1‑139 (MQ=255)                                                         ccgccgTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAg                                                             >  1:586070/1‑139 (MQ=255)                                                         ccgccgTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAg                                                             >  2:649176/1‑139 (MQ=255)                                                                       ccTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGAt                                               >  1:428580/1‑139 (MQ=255)                                                                           gCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTc                                           >  1:290758/1‑139 (MQ=255)                                                                                 cccGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGt                                     >  1:453384/1‑139 (MQ=255)                                                                                     tGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTa                                                                               <  2:335833/93‑1 (MQ=255)                                                                                     tGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTa                                                                               >  1:335833/1‑93 (MQ=255)                                                                                        aGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCgg                              <  2:36665/139‑1 (MQ=255)                                                                                            tttGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGTGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGa                                                         <  2:388444/108‑1 (MQ=255)                                                                                            tttGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGa                                                         >  1:388444/1‑108 (MQ=255)                                                                                             ttGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGc                         <  1:70401/139‑1 (MQ=255)                                                                                                     gcagTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGTAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGaaaa                 <  1:90541/139‑1 (MQ=255)                                                                                                     gcagTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGaaaa                 <  2:632459/139‑1 (MQ=255)                                                                                                         ttACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGCa             >  2:708601/1‑139 (MQ=255)                                                                                                          tACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGGtttt                                            >  2:85210/1‑107 (MQ=255)                                                                                                          tACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGAtttt                                            <  1:85210/107‑1 (MQ=255)                                                                                                          tACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGCAt            >  1:433000/1‑139 (MQ=255)                                                                                                                cTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTa                                                                               <  2:621620/66‑1 (MQ=255)                                                                                                                cTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTa                                                                               >  1:621620/1‑66 (MQ=255)                                                                                                                 tGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACg                                                                                         >  1:760611/1‑55 (MQ=255)                                                                                                                 tGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACg                                                                                         <  2:760611/55‑1 (MQ=255)                                                                                                                   gcgcgcCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCa                                                                     >  2:274667/1‑73 (MQ=255)                                                                                                                   gcgcgcCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCa                                                                     <  1:274667/73‑1 (MQ=255)                                                                                                                    cgcgcCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGCATCTATAACgag  <  1:708601/139‑1 (MQ=255)                                                                                                                    cgcgcCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGCATCTATAACgag  <  1:649176/139‑1 (MQ=255)                                                                                                                                |                                                                                                                               CGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAGCGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGCATCTATAACGAG  >  NC_000913/3103938‑3104191

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