Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,713,323	Δ1 bp	coding (232/564 nt)	tag →	3‑methyl‑adenine DNA glycosylase I, constitutive

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,713,321	0	A	.	100.0%	185.8 / NA	44	coding (230/564 nt)	tag	3‑methyl‑adenine DNA glycosylase I, constitutive
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (18/26);  total (18/26)

AACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTA  >  NC_000913/3713246‑3713396                                                                            |                                                                            aaCGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGaagact                                                                        >  2:725061/1‑77 (MQ=255)    gcgAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGt                                                                    <  1:724986/81‑1 (MQ=255)      gAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTa                                                                   >  2:725055/1‑80 (MQ=255)        aaaCTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAg                                                                >  2:725066/1‑81 (MQ=255)         aaCTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAg                                                                >  1:725045/1‑80 (MQ=255)          aCTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGa                                                              <  2:724987/81‑1 (MQ=255)           cTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGAc                                                             <  1:724989/81‑1 (MQ=255)           cTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGAc                                                             <  2:724988/81‑1 (MQ=255)            tATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACg                                                            >  2:725065/1‑81 (MQ=255)              tCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGcc                                                          <  1:724990/81‑1 (MQ=255)                gcgcCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCg                                                         <  1:724991/80‑1 (MQ=255)                gcgcCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCg                                                         <  1:724992/80‑1 (MQ=255)                 cgcCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCggg                                                       <  1:724993/81‑1 (MQ=255)                  gcCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGa                                                      >  2:725043/1‑81 (MQ=255)                    cTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGAtt                                                    <  2:724994/81‑1 (MQ=255)                       cTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATc                                                 >  1:725060/1‑81 (MQ=255)                        tttCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATcc                                                >  1:725054/1‑81 (MQ=255)                         ttCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCg                                               >  1:725051/1‑81 (MQ=255)                          tcatcaGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATcc                                                <  1:724995/79‑1 (MQ=255)                            atcaGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCa                                            <  1:724996/81‑1 (MQ=255)                              caGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATc                                          >  1:725068/1‑81 (MQ=255)                               aGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATc                                          <  1:724999/80‑1 (MQ=255)                               aGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATc                                          <  1:724998/80‑1 (MQ=255)                               aGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCg                                         <  2:724997/81‑1 (MQ=255)                                gTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATc                                          >  2:725050/1‑79 (MQ=255)                                gTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGa                                        <  1:725000/81‑1 (MQ=255)                                 ttCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAg                                       <  2:725001/81‑1 (MQ=255)                                         ggTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAAtt                               >  1:725062/1‑81 (MQ=255)                                         ggTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAAtt                               >  2:725053/1‑81 (MQ=255)                                         ggTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAAtt                               <  1:725002/81‑1 (MQ=255)                                              aGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGc                          <  2:725003/81‑1 (MQ=255)                                               ggTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCa                         <  1:725004/81‑1 (MQ=255)                                                   gcagcaATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAatta                     >  1:725070/1‑81 (MQ=255)                                                    cagcaATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAattat                    <  1:725005/81‑1 (MQ=255)                                                     agcaATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAattatt                   >  1:725069/1‑81 (MQ=255)                                                      gcaATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTg                  <  2:725006/81‑1 (MQ=255)                                                         aTGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGt                >  2:725058/1‑80 (MQ=255)                                                            cAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATg            >  2:725064/1‑81 (MQ=255)                                                               gAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATgcgc         >  1:725072/1‑81 (MQ=255)                                                                 agagGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCgg       <  2:725007/81‑1 (MQ=255)                                                                     gATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGt   <  1:725010/81‑1 (MQ=255)                                                                     gATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGt   <  1:725011/81‑1 (MQ=255)                                                                     gATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGt   <  2:725009/81‑1 (MQ=255)                                                                     gATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGt   <  2:725008/81‑1 (MQ=255)                                                                      aTGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTa  >  2:725049/1‑81 (MQ=255)                                                                            |                                                                            AACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTA  >  NC_000913/3713246‑3713396

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