Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,918,220	T→A	N16Y (AAC→TAC)	atpG ←	ATP synthase F1 complex subunit gamma

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,918,220	0	T	A	100.0%	147.0 / NA	42	N16Y (AAC→TAC)	atpG	ATP synthase F1 complex subunit gamma
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (23/19);  total (23/19)

TGGCCGCCATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCC  >  NC_000913/3918142‑3918291                                                                               |                                                                        tGGCCGCCATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACt                                                                       >  2:705765/1‑81 (MQ=255) tGGCCGCCATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACt                                                                       >  1:705755/1‑81 (MQ=255)  ggCCGCCATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACTg                                                                      >  1:705766/1‑81 (MQ=255)   gccgccATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACt                                                                       <  1:705697/79‑1 (MQ=255)   gccgccATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACt                                                                       <  2:705698/79‑1 (MQ=255)   gccgccATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACt                                                                       <  2:705696/79‑1 (MQ=255)     cgccATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACTGGAc                                                                   <  1:705699/81‑1 (MQ=255)     cgccATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACTGGAc                                                                   <  1:705700/81‑1 (MQ=255)      gccATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCTATCGCTTTAGTGATCTTTTGCGTGTACTGGACg                                                                  >  1:705734/1‑81 (MQ=255)      gccATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACTGGACg                                                                  >  2:705763/1‑81 (MQ=255)      gccATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACTGGACg                                                                  >  1:705764/1‑81 (MQ=255)         aTGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACTGGACGCtt                                                               <  2:705701/81‑1 (MQ=255)         aTGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACTGGACGCtt                                                               >  2:705744/1‑81 (MQ=255)                ccTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACTGGACGCTTGCGATCt                                                        >  1:705743/1‑81 (MQ=255)                 cTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACTGGACGCTTGCGATCtt                                                       <  2:705702/81‑1 (MQ=255)                       tttACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACTGGACGCTTGCGATCTTACTACg                                                 >  1:705736/1‑81 (MQ=255)                       tttACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACTGGACGCTTGCGATCTTACTACg                                                 >  2:705758/1‑81 (MQ=255)                         tACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACTGGACGCTTGCGATCTTACTACGt                                                <  2:705703/80‑1 (MQ=255)                            gCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACTGGACGCTTGCGATCTTACTACGTAtct                                            <  1:705704/81‑1 (MQ=255)                                 ttGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACTGGACGCTTGCGATCTTACTACGTATCTCtttt                                       >  1:705770/1‑81 (MQ=255)                                 ttGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACTGGACGCTTGCGATCTTACTACGTATCTCttt                                        <  2:705705/80‑1 (MQ=255)                                  tGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACTGGACGCTTGCGATCTTACTACGTATCTCtttt                                       <  1:705706/80‑1 (MQ=255)                                   ggAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACTGGACGCTTGCGATCTTACTACGTATCTCTTTTg                                      <  2:705707/80‑1 (MQ=255)                                     aaGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACTGGACGCTTGCGATCTTACTACGTATCTCTTTTgcgc                                   >  1:705773/1‑81 (MQ=255)                                     aaGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACTGGACGCTTGCGATCTTACTACGTATCTCTTTTgcgc                                   >  2:705772/1‑81 (MQ=255)                                      aGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGcc                                  <  2:705708/81‑1 (MQ=255)                                         ggcgACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGc                               >  1:705752/1‑81 (MQ=255)                                          gcgACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGcc                              <  1:705709/81‑1 (MQ=255)                                          gcgACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGc                               >  2:705746/1‑80 (MQ=255)                                            gACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTACTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCAt                            >  1:705725/1‑81 (MQ=255)                                              ccATCTCCATCGCTTTAGTGATCTTTTGCGTGTACTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCAt                            >  2:705753/1‑79 (MQ=255)                                              ccATCTCCATCGCTTTAGTGATCTTTTGCGTGTACTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGa                          <  1:705710/81‑1 (MQ=255)                                                     cATCGCTTTAGTGATCTTTTGCGTGTACTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTtctc                   >  1:705757/1‑81 (MQ=255)                                                       tCGCTTTAGTGATCTTTTGCGTGTACTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTctcct                 <  2:705711/81‑1 (MQ=255)                                                        cGCTTTAGTGATCTTTTGCGTGTACTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTctcct                 >  2:705759/1‑80 (MQ=255)                                                        cGCTTTAGTGATCTTTTGCGTGTACTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTctcct                 >  2:705782/1‑80 (MQ=255)                                                             tAGTGATCTTTTGCGTGTACTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAAt             <  2:705712/79‑1 (MQ=255)                                                              aGTGATCTTTTGCGTGTACTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGcc          >  2:705729/1‑81 (MQ=255)                                                               gTGATCTTTTGCGTGTACTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGcc          <  2:705713/80‑1 (MQ=255)                                                                    cTTTTGCGTGTACTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCgg    >  2:705771/1‑81 (MQ=255)                                                                      tttGCGTGTACTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGcc  >  2:705738/1‑81 (MQ=255)                                                                       ttGCGTGTACTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGc   <  1:705714/79‑1 (MQ=255)                                                                               |                                                                        TGGCCGCCATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCC  >  NC_000913/3918142‑3918291

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