Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,950,945	G→A	G38G (GGG→GGA)	yecD →	isochorismatase family protein YecD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,950,945	0	G	A	100.0%	101.8 / NA	29	G38G (GGG→GGA)	yecD	isochorismatase family protein YecD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/13);  total (16/13)

GCCCGAACTCATAAAAATTAACAAAATTTGTCGTTCCGCCATCGGCTAATCGCATTAAGGTGTGAGGCACGATTTTGTTTTTTCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATC  >  NC_000913/1950738‑1951165                                                                                                                                                                                                                |                                                                                                                                                                                                                             gcacGAACTCATAAAAATTAACCAAATTTGTCGTTCCGCCATCGGCTAATCGCATTAACGTGTGAAACACGATTTTGTTTTTTCAGCAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACACGAAGGCATCTTACCTTTTGCCGGTAGTCCCCATACTGCCGATGAGGTAGTTAATCGCGCCGGAAAGCTGGCGGCGAAATTTCGCTCCAGCGGACAGCCCGTGTTTc                                                                                                                                                                                   >  2:46823/4‑251 (MQ=255) gCCCGAACTCATAAAAATTAACAAAATTTGTCGTTCCGCCATCGGCTAATCGCATTAAGGTGTGAGGCACGATTTTGTTTTTTCAGGAGTCATTATTCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGAAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGttt                                                                                                                                                                                    <  1:304447/250‑1 (MQ=255)            tAAAAATTAACAAAATTTGTCGTTCCGCCATCGGCTAATCGCATTAAGGTGTGAGGCACGATTTTGTTTTTTCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGAAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGtt                                                                                                                                                                        >  1:74701/1‑251 (MQ=255)             aaaaaTTAACAAAATTTGTCGTTCCGCCATCGGCTAATCGCATTAAGGTGTGAGGCACGATTTTGTTTTTTCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGAAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCNTTg                                                                                                                                                                       >  2:375719/1‑251 (MQ=255)              aaaaTTAACAAAATTTGTCGTTCCGCCATCGGCTAATCGCATTAAGGTGTGAGGCACGATTTTGTTTTTTCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGAAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTgg                                                                                                                                                                      >  1:123562/1‑251 (MQ=255)               aaaTTAACAAAATTTGTCGTTCCGCCATCGGCTAATCGCATTAAGGTGTGAGGCACGATTTTGTTTTTTCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGAAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGc                                                                                                                                                                     >  2:155022/1‑251 (MQ=255)                    aaCAAAATTTGTCGTTCCGCCATCGGCTAATCGCATTAAGGTGTGAGGCACGATTTTGTTTTTTCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGAAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGt                                                                                                                                                                 <  1:352845/250‑1 (MQ=255)                                    ccGCCATCGGCTAATCGCATTAAGGTGTGAGGCACGATTTTGTTTTTTCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGAAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGaa                                                                                                                                                <  1:317674/251‑1 (MQ=255)                                         aTCGGCTAATCGCATTAAGGTGTGAGGCACGATTTTGTTTTTTCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGAAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCAtt                                                                                                                                           >  1:235435/1‑251 (MQ=255)                                                        tAAGGTGTGAGGCACGATTTTGTTTTTTCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGAAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCg                                                                                                                             >  1:218032/1‑250 (MQ=255)                                                          aGGTGTGAGGCACGATTTTGTTTTTTCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTTGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGAAAGCTGGCGGCTAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTTCGCGTTGGCTGGTCTGCCTTTTACGCTGAAGCATTAAAACAGCCGGTCGATg                                                                                                                          <  1:384363/251‑1 (MQ=255)                                                          aGGTGTGAGGCACGATTTTGTTTTTTCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGAAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg                                                                                                                          >  1:157403/1‑251 (MQ=255)                                                              gtgAGGCACGATTTTGTTTTTTCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGAAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGcc                                                                                                                      <  1:375719/251‑1 (MQ=255)                                                                   gCACGATTTTGTTTTTTCAGGAGTCATTATGCTTTAACTTAATTCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGAAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTTGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCAc                                                                                                                 <  1:209675/251‑1 (MQ=255)                                                                   gCACGATTTTGTTTTTTCAGGAGTAATTATGCTTGAACTCAATGCTAAAACCACTGCACTGGTGGTGATCGATTTACAAGAAGGCCTCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGAAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTAAGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCAc                                                                                                                 >  2:412179/1‑251 (MQ=255)                                                                                                    tGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGAAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAAc                                                                                >  2:88429/1‑251 (MQ=255)                                                                                                     gAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGTCAGCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGAAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAAc                                                                                <  1:235929/250‑1 (MQ=255)                                                                                                      aaCTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAACGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGAAAGCTAGCGGCGCAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACAt                                                                              >  1:19668/1‑251 (MQ=255)                                                                                                           aaTGCTAAAACCACTGCGCTGGAGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGAAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGACCGAAAATTGGTGGCAACATCctgc                                                                         >  2:276916/1‑251 (MQ=255)                                                                                                              gCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGAAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCgg                                                                                                                                <  1:366647/193‑1 (MQ=255)                                                                                                              gCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGAAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCgg                                                                                                                                >  2:366647/1‑193 (MQ=255)                                                                                                              gCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGAAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCtgctgc                                                                      >  1:44170/1‑251 (MQ=255)                                                                                                                     cctcTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGGGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGAAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTGACCTGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGt                                                               <  2:406136/248‑1 (MQ=255)                                                                                                                                      aTCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGAAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGatat                                              >  1:88086/1‑251 (MQ=255)                                                                                                                                                    agGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGAAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAAc                                <  1:155022/251‑1 (MQ=255)                                                                                                                                                                       ggtggtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGAAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCg                                                  <  2:116342/214‑1 (MQ=255)                                                                                                                                                                       ggtggtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGAAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCg                                                  >  1:116342/1‑214 (MQ=255)                                                                                                                                                                             ccACATTCTGCCGATGAGGTAGTTAATCGCGCCGGAAAGCTGGCGTCGAAATTTCGCGCCAGCAGTCAGCCCGTGTTTCTTGTGCGCGTTGGCTGGTCTGCCCATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGTCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTtacggtac       <  1:384125/251‑1 (MQ=255)                                                                                                                                                                                   aCTGCCGAGGAGGTAGTTAATCGAGCCGGAAAGCTGGCGGCGAAATTGCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTTCGGATc  <  2:398249/250‑1 (MQ=255)                                                                                                                                                                                                                |                                                                                                                                                                                                                             GCCCGAACTCATAAAAATTAACAAAATTTGTCGTTCCGCCATCGGCTAATCGCATTAAGGTGTGAGGCACGATTTTGTTTTTTCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATC  >  NC_000913/1950738‑1951165

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