Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,740,926	G→A	E408K (GAG→AAG)	avtA →	valine‑‑pyruvate aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,740,926	0	G	A	100.0%	67.9 / NA	20	E408K (GAG→AAG)	avtA	valine‑‑pyruvate aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (7/13);  total (7/13)

TACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTC  >  NC_000913/3740850‑3741000                                                                             |                                                                           tACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAAAGAt                                                                        <  1:669628/81‑1 (MQ=255)       tCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAAAGATAGAAag                                                                  <  1:669629/81‑1 (MQ=255)              aTGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAAAGATAGAAAGAGCCTgg                                                           >  1:669656/1‑81 (MQ=255)               tGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAAAGATAGAAAGAGCCTggg                                                          <  1:669630/81‑1 (MQ=255)                       aaCTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAAAGATAGAAAGAGCCTGGGCTGAAAGt                                                  <  1:669631/81‑1 (MQ=255)                            cGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAAAGATAGAAAGAGCCTGGGCTGAAAGTCAc                                               <  1:669632/79‑1 (MQ=255)                                  agagCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAAAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAAtt                                       <  2:669633/81‑1 (MQ=255)                                        ggAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAAAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGa                                 >  2:669653/1‑81 (MQ=255)                                        ggAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAAAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGa                                 >  1:669654/1‑81 (MQ=255)                                           gaAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAAAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATgc                              <  1:669634/81‑1 (MQ=255)                                            aaaaTTGAGGCGGGGGTGAAGATTCTGGCGGAAAAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATgcg                             >  1:669658/1‑81 (MQ=255)                                                ttGAGGCGGGGGTGAAGATTCTGGCGGAAAAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTAc                         <  2:669635/81‑1 (MQ=255)                                                  gAGGCGGGGGTGAAGATTCTGGCGGAAAAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGc                       >  2:669661/1‑81 (MQ=255)                                                      cGGGGGTGAAGATTCTGGCGGAAAAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTAt                   >  2:669657/1‑81 (MQ=255)                                                       gggggTGAAGATTCTGGCGGAAAAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATc                  <  2:669636/81‑1 (MQ=255)                                                        ggggTGAAGATTCTGGCGGAAAAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCa                 <  2:669637/81‑1 (MQ=255)                                                              aaGATTCTGGCGGAAAAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCt            <  1:669638/80‑1 (MQ=255)                                                                  ttCTGGCGGAAAAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACAt        >  2:669655/1‑80 (MQ=255)                                                                  ttCTGGCGGAAAAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATg       <  2:669639/81‑1 (MQ=255)                                                                       gCGGAAAAGATAGAAAGAGCCTGTGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTc  <  1:669640/81‑1 (MQ=255)                                                                             |                                                                           TACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTC  >  NC_000913/3740850‑3741000

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