Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	542,530	A→G	T215A (ACG→GCG)	glxK →	glycerate kinase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	542,530	0	A	G	100.0%	126.7 / NA	36	T215A (ACG→GCG)	glxK	glycerate kinase II
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/19);  total (17/19)

GCAAGCATTGAGATGGATGAATGCGATCCGCGTCTGGCGAATTGCCATATTGAAGTAGCATGTGACGTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTGCGGGCGGT  >  NC_000913/542401‑542643                                                                                                                                  |                                                                                                                  gCAAGCATTGAGATGGATGAATGCGATCCGCGTCTGGCGAATTGCCATATTGAAGTAGCATGTGACGTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGgaga                                                                                                          <  1:593257/139‑1 (MQ=255)  cAAGCATTGAGATGGATGAATGCGATCCGCGTCTGGCGAATTGCCATATTGAAGTAGCATGTGACGTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGAt                                                                                                         <  2:636002/139‑1 (MQ=255)        ttGAGATGGATGAATGCGATCCGCGTCTGGCGAATTGCCATATTGAAGTAGCATGTGACGTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCga                                                                                                   <  1:641438/139‑1 (MQ=255)          gagaTGGATGAATGCGATCCGCGTCTGGCGAATTGCCATATTGAAGTAGCATGTGACGTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCgaag                                                                                                 <  1:643559/139‑1 (MQ=255)            gatggatgAATGCGATCCGCGTCTGGCGAATTGCTATATTGAAGTAGCATGTGACGTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCgaagaa                                                                                               >  2:355832/1‑139 (MQ=255)                      tGCGATCCGCGTCTGGCGAATTGCCATATTGAAGTAGCATGTGACGTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCGAAGAACTTGAACAgg                                                                                     <  2:252856/139‑1 (MQ=255)                          aTCCGCGTCTGGCGAATTGCCATATTGAAGTAGCATGTGACGTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCGAAGAACTTGAACAGGGGCt                                                                                 >  1:488447/1‑139 (MQ=255)                           tCCGCGTCTGGCGAATTGCCATATTGAAGTAGCATGTGACGTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCGAAGAACTTGAACAGGGGCTg                                                                                >  1:226179/1‑139 (MQ=255)                                         aTTGCTATATTGAAGTAGCATGTGACGTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCGAAGAACTTGAACAGGGGCTGCAAAATTACGCCCg                                                                  <  1:355832/139‑1 (MQ=255)                                         aTTGCCATATTGAAGTAGCATGTGACGTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCGAAGAACTTGAACAGGGGCTGCAAAATTACGCCCg                                                                  <  1:196005/139‑1 (MQ=255)                                               atatTGAAGTAGCATGTGACGTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCgaag                                                                                                 >  2:499632/1‑102 (MQ=255)                                               atatTGAAGTAGCATGTGACGTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCgaag                                                                                                 <  1:499632/102‑1 (MQ=255)                                                          gCATGTGACGTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCGAAGAACTTGAACAGGGGCTGCAAAATTACGCCCGTGTTTTACAACAGCaaa                                                 >  2:692373/1‑139 (MQ=255)                                                           cATGTGACGTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCGAAGAACTTGAACAGGGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAAc                                                >  2:47861/1‑139 (MQ=255)                                                           cATGTGACGTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCGAAGAACTTGAACAGGGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAAc                                                >  2:117859/1‑139 (MQ=255)                                                            aTGTGACGTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCGAAGAACTTGAACAGGGGCTGCAAAATTACGc                                                                     >  1:188991/1‑117 (MQ=255)                                                            aTGTGACGTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCGAAGAACTTGAACAGGGGCTGCAAAATTACGc                                                                     <  2:188991/117‑1 (MQ=255)                                                                   gTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGt                                                                                                      <  1:193572/77‑1 (MQ=255)                                                                   gTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGt                                                                                                      >  2:193572/1‑77 (MQ=255)                                                                    ttGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCGAAGAACTTGAACAGGGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTaa                                       <  1:206995/139‑1 (MQ=255)                                                                      gATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCGAAGAACTTGAACAGGGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATg                                     >  2:512528/1‑139 (MQ=255)                                                                            ccGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCGAAGAACTTGAACAGGGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGcc                               >  2:406314/1‑139 (MQ=255)                                                                               cTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCGAAGAACTTGAACAGGGGCTGCAAAATTACGccc                                                                   <  1:583702/100‑1 (MQ=255)                                                                               cTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCGAAGAACTTGAACAGGGGCTGCAAAATTACGccc                                                                   >  2:583702/1‑100 (MQ=255)                                                                                 ggTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCGAAGAACTTGAACAGGGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATg                          >  1:154120/1‑139 (MQ=255)                                                                                      gggCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCGAAGAACTTGAACAGGGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATggcggg                     >  2:480191/1‑139 (MQ=255)                                                                                            gcggcgcggcGGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCGAAGAACTTGAACAGGGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGggcggcgg               >  2:213485/1‑139 (MQ=255)                                                                                             cggcgcggcgGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCGAAGAACTTGAACAGGGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGt                                    <  1:77717/117‑1 (MQ=255)                                                                                             cggcgcggcgGCGGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCGAAGAACTTGAACAGGGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGt                                    >  2:77717/1‑117 (MQ=255)                                                                                                 gcggcggcgGTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCGAAGAACTTGAACAGGGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTg          <  1:512528/139‑1 (MQ=255)                                                                                                    gcggcggTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCGAAGAACTTGAACAGGGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTgcgg       <  2:226179/139‑1 (MQ=255)                                                                                                    gcggcggTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCGAAGAACTTGAACAGGGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTgcgg       <  1:47861/139‑1 (MQ=255)                                                                                                     cggcggTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCGAAGAACTTGAACAGGGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTgcggg      <  1:406314/139‑1 (MQ=255)                                                                                                         ggTGTTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCGAAGAACTTGAACAGGGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTGCGGGCGGt  <  1:178536/139‑1 (MQ=255)                                                                                                          gtgtTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCGAAGAAc                                                                                              <  2:166262/46‑1 (MQ=255)                                                                                                          gtgtTTGGCCCGCAAAAAGGGGCAGCGCCGGAGATGGTCGAAGAAc                                                                                              >  1:166262/1‑46 (MQ=255)                                                                                                                                  |                                                                                                                  GCAAGCATTGAGATGGATGAATGCGATCCGCGTCTGGCGAATTGCCATATTGAAGTAGCATGTGACGTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTGCGGGCGGT  >  NC_000913/542401‑542643

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