Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_vch_tpiA	2,817,170	C→A	E112* (GAA→TAA)	yqaB ←	fructose‑1‑P and 6‑phosphogluconate phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_vch_tpiA	2,817,170	0	C	A	100.0%	51.6 / NA	19	E112* (GAA→TAA)	yqaB	fructose‑1‑P and 6‑phosphogluconate phosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (9/10);  total (9/10)

CACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGAT  >  NC_000913_3_vch_tpiA/2817033‑2817310                                                                                                                                          |                                                                                                                                             cACGCCCATACGCTGCGCGCACAACAAAACTGTGTCTGGCGCTTGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTc                                                                                                                                   <  1:10295/149‑1 (MQ=255)     cccATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGGGATCGACAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAACTGCGCCAGCAACGCATCAGCGATGGAGCATTAACACCCCGTTCCTAc                                                                                                                               >  1:201629/1‑149 (MQ=255)               gcgcgcACAACAAAAATGTGTCGGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTggg                                                                                                                     <  2:10007/149‑1 (MQ=255)                  cgcACAACAAAAATGTGTCGGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGTCGCTTTAACTCCCCGTTCCTACAGCCATTGGGcga                                                                                                                  <  2:148889/149‑1 (MQ=255)                      caacaaAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGATGCGCCAGCAATGCATCAACGATGACGCTTTAACTCCCCCTTCCTACAGCCATTGGGCgacgac                                                                                                              >  1:53462/1‑149 (MQ=255)                          aaaaaTGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATg                                                                                                          >  1:131954/1‑149 (MQ=255)                                    cTGGGGCGGGTTTATGGTGTTTGGCGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTc                                                                                                <  2:344126/149‑1 (MQ=255)                                    cTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCAAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTc                                                                                                <  2:376858/149‑1 (MQ=255)                                                       tatGACATGATCGGCAGCGCCGACGGCGTCAAAAGAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGCCGCTTTAACTCACCGCTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGaa                                                                             <  2:291436/147‑1 (MQ=255)                                                          gACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCg                                                                          <  2:152662/149‑1 (MQ=255)                                                             gTGATCGGCAGTGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGtt                                                                       <  1:162803/149‑1 (MQ=255)                                                                                    cAAAATAATGGCGTAATACCAGGTGCGCCAGCAATGCCTAAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCACCGACCATGCCCACTTTTCACCCCATCACCAAGAGGACGTGGCTCGACGCTATCCACCAACATAcat                                                >  2:352375/1‑147 (MQ=255)                                                                                     aaaaTAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATActtc                                               >  1:31078/1‑149 (MQ=255)                                                                                      aaaTAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATActtct                                              >  1:366308/1‑149 (MQ=255)                                                                                                     cccAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCGCTCCACCCTGTCGATGCCGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGtttt                               <  2:292299/149‑1 (MQ=255)                                                                                                               ccAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACgcgcg                     >  1:180010/1‑148 (MQ=255)                                                                                                                      tGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCAt              <  1:266021/149‑1 (MQ=255)                                                                                                                          tCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGCTCGACGCTATCCAGCAGCAAACCACAAACTGCTCCTGTTTTTTCACGCGCTAACGCACGCgg          >  1:63631/1‑149 (MQ=255)                                                                                                                                  ggCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACTACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCGGCATACTTCTTGATGCTTCTGGTTTTTCACGCACTAACGCATGCGGAGCGAGAt  >  1:206608/1‑149 (MQ=255)                                                                                                                                          |                                                                                                                                             CACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGAT  >  NC_000913_3_vch_tpiA/2817033‑2817310

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑