Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,728,012	A→G	Y59C (TAC→TGC)	gloA →	glyoxalase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,728,012	0	A	G	100.0%	106.2 / NA	30	Y59C (TAC→TGC)	gloA	glyoxalase I
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/19);  total (11/19)

GATTTTTATACCAAAGTGCTGGGCATGAAACTGCTGCGTACCAGCGAAAACCCGGAATACAAATACTCACTGGCGTTTGTTGGCTACGGCCCGGAAACCGAAGAAGCGGTGATTGAACTGACCTACAACTGGGGCGTGGATAAATACGAACTCGGCACTGCTTATGGTCACATCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTG  >  NC_000913/1727888‑1728124                                                                                                                             |                                                                                                                 gATTTTTATACCAAAGTGCTGGGCATGAAACTGCTGCGTACCAGCGAAAACCCGGAATACAAATACTCACTGGCGTTTGTTGGCTACGGCCCGGAAACCGAAGAAGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGa                                                                                                   <  1:725038/140‑1 (MQ=255)         taCCAAAGTGCTGGGCATGAAACTGCTGCGTACCAGCGAAAACCCGGAATACAAATACTCACTGGCGTTTGTTGGCTACGGCCCGGAAACCGAAGAAGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGATAAATACg                                                                                           >  2:700099/1‑140 (MQ=255)             aaaGTGCTGGGCATGAAACTGCTGCGTACCAGCGAAAACCCGGAATACAAATACTCACTGGCGTTTGTTGGCTACGGCCCGGAAACCGAAGAAGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGATAAATACGAACt                                                                                       <  2:39699/140‑1 (MQ=255)               aGTGCTGGGCATGAAACTGCTGCGTACCAGCGAAAACCCGGAATACAAATACTCACTGGCGTTTGTTGGCTACGGCCCGGAAACCGAAGAAGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGATAAATACGAACTCg                                                                                     >  1:320839/1‑140 (MQ=255)                  gCTGGGCATGAAACTGCTGCGTACCAGCGAAAACCCGGAATACAAATACTCACTGGCGTTTGTTGGCTACGGCCCGGAAACCGAAGAAGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGATAAATACGAACTCTGCa                                                                                  <  1:736138/140‑1 (MQ=255)                  gCTGGGCATGAAACTGCTGCGTACCAGCGAAAACCCGGAATACAAATACTCACTGGCGTTTGTTGGCTACGGCCCGGAAACCGAAGAAGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGATAAATACGAACTCGGCa                                                                                  <  2:2485/140‑1 (MQ=255)                   cTGGGCATGAAACTGCTGCGTACCAGCGAAAACCCGGAATACAAATACTCACTGGCGTTTGTTGGCTACGGCCCGGAAACCGAAGAAGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGATAAATACGAACTCGGCAc                                                                                 <  1:171976/140‑1 (MQ=255)                           gAAACTGCTGCGTACCAGCGAAAACCCGGAATACAAATACTCACTGGCGTTTGTTGGCTACGGCCCGGAAACCGAAGAAGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGATAAATACGAACTCGGCACTGCTTATg                                                                         >  1:277536/1‑140 (MQ=255)                           gAAACTGCTGCGTACCAGCGAAAACCCGGAATACAAATACTCACTGGCGTTTGTTGGCTACGGCCCGGAAACCGAAGAAGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGATAAATACGAACTCGGCACTGCTTATg                                                                         <  2:127525/140‑1 (MQ=255)                                 gctgcGTACCAGCGAAAACCCGGAATACAAATACTCACTGGCGTTTGTTGGCTACGGCCCGGAAACCGAAGAAGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGATAAATACGAACTCGGc                                                                                   >  1:383547/1‑124 (MQ=255)                                 gctgcGTACCAGCGAAAACCCGGAATACAAATACTCACTGGCGTTTGTTGGCTACGGCCCGGAAACCGAAGAAGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGATAAATACGAACTCGGc                                                                                   <  2:383547/124‑1 (MQ=255)                                  ctgcGTACCAGCGAAAACCCGGAATACAAATACTCACTGGCGTTTGTTGGCTACGGCCCGGAAACCGAAGAAGCGGTGATTGAACTGACCTGGAACTGGGGCGTGGATAAATACGAACTCGGCACTGCTTATGGTCACAt                                                                  <  1:134334/140‑1 (MQ=255)                                   tgcGTACCAGCGAAAACCCGGAATACAAATACTCACTGGCGTTTGTTGGCTACGGCCCGGAAACCGAAGAAGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGATAAATACGAACTCGGCACTGCTTATGGTCACATc                                                                 <  2:412201/140‑1 (MQ=255)                                       tACCAGCGAAAACCCGGAATACAAATACTCACTGGCGTTTGTTGGCTACGGCCCGGAAACCGAAGAAGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGATAAATACGAACTCGGCACTGCTTATGGTCACATCgcgc                                                             <  2:922201/140‑1 (MQ=255)                                       tACCAGCGAAAACCCGGAATACAAATACTCACTGGCGTTTGTTGGCTACGGCCCGGAAACCGAAGAAGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGATAAATACGAACTCGGCACTGCTTATGGTCACATCgcgc                                                             >  2:710482/1‑140 (MQ=255)                                                    ccGGAATACAAATACTCACTGGCGTTTGTTGGCTACGGCCCGGAAACCGAAGAAGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGATAAATACGAACTCGGCACTGCTTATGGTCACATCGCGCTTAGCGTAGATaa                                                <  1:57151/140‑1 (MQ=255)                                                    ccGGAATACAAATACTCACTGGCGTTTGTTGGCTACGGCCCGGAAACCGAAGAAGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGATAAATACGAACTCGGCACTGCTTATGGTCACATCGCGCTTAGCGTAGATaa                                                <  2:84478/140‑1 (MQ=255)                                                      ggAATACAAATACTCACTGGCGTTTGTTGGCTACGGCCCGGAAACCGAAGAAGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGATAAATACGAACTCGGCACTGCTTATGGTCACATCGCGCTTAGCGTAGATAAcg                                              >  1:875825/1‑140 (MQ=255)                                                                  cTCACTGGCGTTTGTTGGCTACGGCCCGGAAACCGAAGAAGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGATAAATACGAACTCGGCACTGCTTATGGTCATATCGCGCTTAGCGTAGATaa                                                <  2:663071/126‑1 (MQ=255)                                                                  cTCACTGGCGTTTGTTGGCTACGGCCCGGAAACCGAAGAAGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGATAAATACGAACTCGGCACTGCTTATGGTCACATCGCGCTTAGCGTAGATaa                                                >  1:663071/1‑126 (MQ=255)                                                                            tttGTTGGCTACGGCCCGGAAACCGAAGAAGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGATAAATACGAACTCGGCACTGCTTATGGTCACATCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAAt                        >  2:485769/1‑140 (MQ=255)                                                                            tttGTTGGCTACGGCCCGGAAACCGAAGAAGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGATAAATACGAACTCGGCACTGCTTATGGTCACATCGCGCTTAGCATAGATAACGCCGCTGAAGCGTGCGAAAAAAt                        >  1:758288/1‑140 (MQ=255)                                                                                ttgGCTACGGCCCGGAAACCGAAGAAGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGATAAATACGAACTCGGCACTGCTTATGGTCACATCGCGCTTAGCATAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGt                    <  2:758288/140‑1 (MQ=255)                                                                                   gCTACGGCCCGGAAACCGAAGAAGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGATAAATACGAACTCGGCACTGCTTATGGTCACATCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCaa                 <  1:485769/140‑1 (MQ=255)                                                                                     tACGGCCCGGAAACCGAAGAAGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGATAAATACGAACTCGGCACTGCTTATGGTCACATCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCaaaa               <  2:277536/140‑1 (MQ=255)                                                                                       cGGCCCGGAAACAGAAGAAGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGATAAATACGAACTCGGCACTGCTTATGGTCACATCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACg             >  2:224427/1‑140 (MQ=255)                                                                                          cccGGAAACCGAAGAAGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGATAAATACGAACTCGGCACTGCTTATGGTCACATCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACgggg          <  2:320839/140‑1 (MQ=255)                                                                                          cccGGAAACCGAAGAAGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGATAAATACGAACTCGGCACTGCTTATGGTCACATCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACgggg          <  1:710482/140‑1 (MQ=255)                                                                                          cccGGAAACCGAAGAAGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGATAAATACGAACTCGGCACTGCTTATGGTCACATCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACgggg          >  2:403741/1‑140 (MQ=255)                                                                                                  cagaagaaGCGGTGATTGAACTGACCTGCAACTGGGGCGTGGATAAATACGAACTCGGCACTGCTTATGGTCACATCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTg  <  1:224427/138‑1 (MQ=255)                                                                                                                             |                                                                                                                 GATTTTTATACCAAAGTGCTGGGCATGAAACTGCTGCGTACCAGCGAAAACCCGGAATACAAATACTCACTGGCGTTTGTTGGCTACGGCCCGGAAACCGAAGAAGCGGTGATTGAACTGACCTACAACTGGGGCGTGGATAAATACGAACTCGGCACTGCTTATGGTCACATCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTG  >  NC_000913/1727888‑1728124

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