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%	151.9 / NA	36	coding (230/564 nt)	tag	3‑methyl‑adenine DNA glycosylase I, constitutive
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (16/20);  total (16/20)

AAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTAC  >  NC_000913/3713244‑3713397                                                                              |                                                                             aaaaCGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGaaga                                                                           >  2:658139/1‑79 (MQ=255)  aaaCGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGaagac                                                                          >  2:658149/1‑78 (MQ=255)  aaaCGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGaagac                                                                          >  2:658110/1‑78 (MQ=255)  aaaCGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGaagac                                                                          >  1:658134/1‑78 (MQ=255)   aaCGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGaagact                                                                         >  1:658147/1‑77 (MQ=255)          aaaCTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAg                                                                 <  2:658090/81‑1 (MQ=255)           aaCTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAgg                                                                <  2:658091/81‑1 (MQ=255)             cTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGAc                                                              >  1:658146/1‑81 (MQ=255)               aTCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGc                                                            >  2:658148/1‑81 (MQ=255)                 cgcgCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCg                                                          >  1:658169/1‑81 (MQ=255)                   cgcCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCgg                                                         >  2:658163/1‑80 (MQ=255)                     ccTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGAt                                                      >  1:658137/1‑81 (MQ=255)                           ttCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCg                                                >  1:658155/1‑81 (MQ=255)                               tcaGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCAt                                            <  1:658093/81‑1 (MQ=255)                                 aGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATc                                           <  2:658094/80‑1 (MQ=255)                                 aGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATc                                           <  2:658093/80‑1 (MQ=255)                                    tCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAgg                                       >  2:658159/1‑81 (MQ=255)                                      gATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAggg                                      <  2:658096/80‑1 (MQ=255)                                      gATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGa                                     <  1:658096/81‑1 (MQ=255)                                         ccGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGaaaa                                  >  1:658151/1‑81 (MQ=255)                                          cGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAAt                                 <  1:658098/81‑1 (MQ=255)                                          cGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAAt                                 <  1:658099/81‑1 (MQ=255)                                           ggTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAAtt                                >  1:658153/1‑81 (MQ=255)                                           ggTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAAtt                                <  2:658099/81‑1 (MQ=255)                                           ggTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAAtt                                <  2:658100/81‑1 (MQ=255)                                           ggTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAAtt                                >  2:658143/1‑81 (MQ=255)                                            gTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTc                               <  2:658101/81‑1 (MQ=255)                                              gAAGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAg                             >  2:658165/1‑81 (MQ=255)                                               aaGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAgg                            >  2:658140/1‑81 (MQ=255)                                                aGGTCGCAGCAATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGc                           <  1:658103/81‑1 (MQ=255)                                                     gcagcaATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAatta                      >  2:658162/1‑81 (MQ=255)                                                     gcagcaATGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAatt                       <  2:658103/80‑1 (MQ=255)                                                          aaTGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTgg                  >  1:658165/1‑80 (MQ=255)                                                            tGCAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTaa               <  1:658105/81‑1 (MQ=255)                                                              cAGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAAt              <  2:658105/80‑1 (MQ=255)                                                               aGGAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATgc            <  1:658107/81‑1 (MQ=255)                                                                 gAAGAGGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATgcgc          <  2:658107/81‑1 (MQ=255)                                                                     agGATGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGc      <  2:658108/81‑1 (MQ=255)                                                                        aTGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAATTTCAGGCAATTATTGGTAATGCGCGGGCGTa   >  2:658160/1‑81 (MQ=255)                                                                         tGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTAc  <  1:658110/81‑1 (MQ=255)                                                                         tGTCG‑AAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTAc  >  2:658155/1‑81 (MQ=255)                                                                              |                                                                             AAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTAC  >  NC_000913/3713244‑3713397

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