Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,554,719	C→T	Q317Q (CAG→CAA)	maeA ←	malate dehydrogenase (oxaloacetate‑decarboxylating)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,554,719	0	C	T	100.0%	43.1 / NA	16	Q317Q (CAG→CAA)	maeA	malate dehydrogenase (oxaloacetate‑decarboxylating)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (10/6);  total (10/6)

CGTAAACAGCCCGGTCTGTCCTGAGACGCCAATCAGAATATCTGGTTTTACATTGCGCACCACATCCAGCAGTGACAGCACATCGCTGTCGGTATCCCAGTCACTGAGGTTTTCGCGCTTCTGCACCAGTTTGGTCTGGAAAGGCAGCAGGTTCGGCATCTTGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAGTGTGCCGACTGTTACCGCCGCAGTGCCCTGAATGTCATCGTTAAAAGAACAAATTTCATTGCGATAGCGGTTAAGTAACGGCATCGCATTTTTTTGAGCAAAGTCTTCAAACTGCAACAGCACGTCTGG  >  NC_000913/1554479‑1554955                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                             cGTAAACAGCCCGGTCTGTCCTGAGACGCCAATCAGAATATCTGGTTTTACATTGCGCACCACATCCAGCAGTGACAGCACATCGCTGTCCGTCTCCCCGTAACTGAGGTTTTCGCGCTTCTGCACCCGTTTGGTCTGGAAAGGCAGCAGGTTCGGCATCTTGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCACTGGGTTTGGGAgatg                                                                                                                                                                                                                                     >  1:442805/1‑250 (MQ=255)                                                                                                       aCTGAGGTTTTCGCGCTTCTGCACCAGTTTGGTCTGGAAAGGCAGCAGGGTCGGCATCTTGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTTCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTTTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTTAACCTGCGCCAATGAAGACGATTTTTGTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAGtgtgt                                                                                                                              <  1:409369/251‑2 (MQ=255)                                                                                                                            naCCAGTTTGGTCTGGAAAGGCAGCAGGTTCGGCATCTTGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTTTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAGTGTGCCGACTGTTACCGCCGCAGTGc                                                                                                         >  1:405790/2‑251 (MQ=255)                                                                                                                                                                       gtcagCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTTTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGcccacc                                                                                                                                                           >  1:151843/1‑158 (MQ=255)                                                                                                                                                                               gCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCATTCGCGCTGGGTTTGGGGGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTCCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAGTGTGCCGAATGTTACCGCCGCAGTGCCCTGAAt                                                                                                  >  1:221205/1‑207 (MQ=255)                                                                                                                                                                                     gCGATCGACCATAAAGACTTTCTGNCGCGCCGCTNNCTNGCTTANTNNNTNNNGCNGGGTTTGGGAGATGATCATNTCNNCAATGCNGCATCCCGCTGAACCTGCGCCAAGGa                                                                                                                                                                                          <  2:118178/113‑1 (MQ=255)                                                                                                                                                                                     gCGATCGAACATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTTTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGa                                                                                                                                                                                          >  1:118178/1‑113 (MQ=255)                                                                                                                                                                                           gACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTTTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGANCACNNNNngc                                                                                                                                                     >  2:311379/1‑144 (MQ=255)                                                                                                                                                                                                    gACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTTTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTNNGCCAAgnn                                                                                                                                                                                          >  2:8660/1‑96 (MQ=255)                                                                                                                                                                                                    gACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTTTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAGTGTGcc                                                                                                                             <  1:8660/159‑1 (MQ=255)                                                                                                                                                                                                              cgcgccgcTTCCTCGCTTAATCCTTCGCGCTGGGTTTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGc                                                                                                                                                        <  1:444807/122‑1 (MQ=255)                                                                                                                                                                                                              cgcgccgcTTCCTCGCTTAATCCTTCGCGCTGGGTTTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAGTGTGCCGACTGTTACCGCCGCAGTGCCCTGAATGTCATCGTTAAAAGAACAAATTTCATTGCGATAGCGGTTAAGTAACGGCATCGCATTTTTTTGAGCAAAGTCTTc                       >  1:390185/1‑251 (MQ=255)                                                                                                                                                                                                                  ccgcTTCCTCGCTTAATCCTTCGCGCTGGGTTTGGGAGATGATCATTTCGGCAAGGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGAGTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAGTGTGCCGACTGTTACCGCCGCAGTGCCCTGAATGTCATCGTTAAAAGAACAAATTTCATTGCGATAGCGGTTAAGTAACGGCATCGCATTTTTTTGAGCAAAGTCTTCaaan                   <  1:311379/251‑2 (MQ=255)                                                                                                                                                                                                                        ccTCGCTTAATCCTTCGCGCTGGGTTTGGGAGATGATCATTTCGGCAAGGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAGTGTGCCGACTGTTGCCGCCGCAGTGCCCTGAATGTCATCGTTAAAAGAACAAATTTCATTGCGATAGCTGTTAAGTAACGGCATCGCATTTTTTTGAGCAAAGTCTTCAAACTGCAAc             <  1:229786/251‑1 (MQ=255)                                                                                                                                                                                                                                   ccTTCGCGCTGGGTTTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAGTGTGCCGACTGTTACCGCCGCAGTGCCCTGAATGTCATCGTTAAAAGAACAAATTTCATTGCGATAGCGGTTAAGTAACGGCATCGCATTTTTTTGAGCAAAGTCTTCAAACTGCAACAGCACGTCTgg  >  1:227961/1‑251 (MQ=255)                                                                                                                                                                                                                                     ttCGCGCTGGGTTTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAGTGTGCCGACTGTTAcc                                                                                                                   >  1:169728/1‑136 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                             CGTAAACAGCCCGGTCTGTCCTGAGACGCCAATCAGAATATCTGGTTTTACATTGCGCACCACATCCAGCAGTGACAGCACATCGCTGTCGGTATCCCAGTCACTGAGGTTTTCGCGCTTCTGCACCAGTTTGGTCTGGAAAGGCAGCAGGTTCGGCATCTTGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAGTGTGCCGACTGTTACCGCCGCAGTGCCCTGAATGTCATCGTTAAAAGAACAAATTTCATTGCGATAGCGGTTAAGTAACGGCATCGCATTTTTTTGAGCAAAGTCTTCAAACTGCAACAGCACGTCTGG  >  NC_000913/1554479‑1554955

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