Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,032,423	G→A	V219V (GTC→GTT)	yedJ ←	putative HD superfamily phosphohydrolase YedJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,032,423	0	G	A	100.0%	104.0 / NA	31	V219V (GTC→GTT)	yedJ	putative HD superfamily phosphohydrolase YedJ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/15);  total (16/15)

ATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCT‑CCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGG‑CCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACA  >  NC_000913/2032204‑2032661                                                                                                                                                                                                                              |                                                                                                                                                                                                                                               aTTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCT‑CCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGG‑CCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACCCCTTCATTCTCCCCCGCCAGTTCGGCAc                                                                                                                                                                                                                 <  1:202881/251‑1 (MQ=255)      gCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCT‑CCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGG‑CCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGc                                                                                                                                                                                                            >  1:257481/1‑251 (MQ=255)       cTGCGCAACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCT‑CCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGG‑CCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAGAACGCATCTATCACCTTGTCATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCt                                                                                                                                                                                                           >  1:331998/1‑251 (MQ=255)        tGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCT‑CCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCAGCGG‑CCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCtt                                                                                                                                                                                                          >  2:104172/1‑251 (MQ=255)         gcgcCACCAACGTTTTAACATCAAAAATTTGCAGCAGTTTCT‑CCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGG‑CCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCGCCCCCGCCAGTTCGGCACTGAGCt                                                                                                                                                                                                           <  2:30826/249‑1 (MQ=255)          cgcCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCT‑CCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGG‑CCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTcccc                                                                                                                                                                                                                               >  2:256553/1‑228 (MQ=255)          cgcCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCT‑CCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGG‑CCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTcccc                                                                                                                                                                                                                               <  1:256553/228‑1 (MQ=255)            caaccaACGTTTTAACATCATAAATTTGCAGAAGTTTCT‑CCAGCATTGCCTGTGCGGAATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGG‑CCATCTGAGAGGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCCCCTTGTGATCAACACCTGCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGcc                                                                                                                                                                                                      <  2:257481/249‑1 (MQ=255)                  acggtttAACATCATAAATTTGCAGCAGTTTCT‑CCAGCATTGCCTGTGCGGCATTCCCGGTGATGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGG‑CCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGAGGAAAACGCCTCTATCACCTTGTGATCACCACCTTCATTCTCCCCCTCCAGTTCGCCACTGAGCTTCGCCATatcc                                                                                                                                                                                                >  2:58/5‑248 (MQ=255)                      tttAACATCATAAATTTGCAGC‑GTTTTTGCCAGCATTGCATGTGCGGCATTCCCGGGTCTGTAGGAATCGGTTACTGATATATTTTCCTGCATTTATTTC‑CCGGTCCAGCTGGTATGTCCGACAGTTTACCATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCTTCTATCACCTTGTTTTCAACACCTTCATGCTCCCCCTCCAGTTCGGCTCTGAGCTTTGCCAGAAACTCg                                                                                                                                                                                             <  2:404551/250‑1 (MQ=255)                           catcatAAATTTGCAGCAGTTTCT‑CCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATGTTTTCCTGCATAGATTTCACCGG‑CCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTaa                                                                                                                                                                                        >  2:212135/1‑250 (MQ=255)                                 aaaTTTGCAGCAGTTTCT‑CCAGCATTGCCTGTGCGGCATTCCCGGGGCTGTATGAATCGGTTACTGATATATTTTCCTGCTTAGATTTCCCCGG‑CCATCTGAGATGGCCTACAGTTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACCCATCTATCACCTTGTGATCAACACCTTCATCCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATgc                                                                                                                                                                                 <  1:229874/251‑1 (MQ=255)                                                                      cATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGG‑CCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCGCCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAgctgc                                                                                                                                                            <  1:402396/236‑1 (MQ=255)                                                                      cATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGG‑CCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAgctgc                                                                                                                                                            >  2:402396/1‑236 (MQ=255)                                                                                   gTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGG‑CCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCg                                                                                                                                >  1:317603/1‑251 (MQ=255)                                                                                          tCGGTTACTGATATATTTTCCTGCATAGATTTCACCGG‑CCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCG‑TAGCTTc                                                                                                                        >  1:196934/1‑251 (MQ=255)                                                                                                             ccTGCATAGATTTCCCCGG‑CCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCCCCTTGTGATCAACCCCTTCATTCCCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGaaa                                                                                                      <  2:388307/251‑1 (MQ=255)                                                                                                             ccTGCATAGATTTCACCGG‑CCATCTGAGATGGCCGACAGGTTAACATAATTCCAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGaaa                                                                                                      <  1:223921/251‑1 (MQ=255)                                                                                                                  aTAGATTTCACCGG‑CCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACCCCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTg                                                                                                                                           <  1:7334/209‑1 (MQ=255)                                                                                                                  aTAGATTTCACCGG‑CCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTg                                                                                                                                           >  2:7334/1‑209 (MQ=255)                                                                                                                         tCACCGG‑CCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCCCGCGCTCAGCCAGCGGATGACAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGcgcg                                                                                          >  2:202771/1‑251 (MQ=255)                                                                                                                                  cATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGttt                                                                                  >  2:124109/1‑251 (MQ=255)                                                                                                                                        agaTGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGAGCTCAGCCAGCGGAAGAAAACGCATCTATAACCTTGTGATCAACACCTTCATTATCCCCCGCCAGTTCGGCACTGAGATTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTatcatc                                                                            >  1:43141/1‑251 (MQ=255)                                                                                                                                            ggCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAg                                                                         <  1:278891/250‑1 (MQ=255)                                                                                                                                                   aGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCAAGCGGATGAAAACGCATCTATCCCCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCAACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGAt                                                                 <  2:280530/251‑1 (MQ=255)                                                                                                                                                    ggTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGAGGAAAACGCATCTATCCCCTTGTGATCAACCCCTTCATTCCCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATg                                                                <  1:237764/251‑1 (MQ=255)                                                                                                                                                                     ccATTTATGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTAATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGAGCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGAGAACGGATCtt                                               >  1:109863/1‑251 (MQ=255)                                                                                                                                                                        tttAGGGCCACCCGCTCAGCCAGCGGATGAAAACGCATCTATCCCCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCG‑TAGCTTTAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCAc                                            <  2:196934/250‑1 (MQ=255)                                                                                                                                                                                                          gCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAAc          >  1:327276/1‑251 (MQ=255)                                                                                                                                                                                                              cTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGc       >  2:155396/1‑250 (MQ=255)                                                                                                                                                                                                                  cACCTTGTGATCAACCCCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAaca  <  1:155396/251‑1 (MQ=255)                                                                                                                                                                                                                              |                                                                                                                                                                                                                                               ATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCT‑CCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGG‑CCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACA  >  NC_000913/2032204‑2032661

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