Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,184,452	G→A	A22V (GCC→GTC)	thiD ←	hydroxy‑methylpyrimidine kinase and hydroxy‑phosphomethylpyrimidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,184,452	0	G	A	100.0%	161.8 / NA	47	A22V (GCC→GTC)	thiD	hydroxy‑methylpyrimidine kinase and hydroxy‑phosphomethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/25);  total (22/25)

CCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGGCCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATGCCTGCACCTCCTGCGTCAATTGCCAGAGCGCATCA  >  NC_000913/2184318‑2184552                                                                                                                                       |                                                                                                     ccGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTg                                                                                                  <  2:1095440/139‑1 (MQ=255) ccGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTg                                                                                                  <  1:805245/139‑1 (MQ=255)  cGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGa                                                                                                 >  1:96127/1‑139 (MQ=255)    aaTCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAAt                                                                                               <  2:175147/139‑1 (MQ=255)    aaTCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAAt                                                                                               >  1:29291/1‑139 (MQ=255)              ggcggcGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGcacca                                                                                     <  2:373472/139‑1 (MQ=255)                      cAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCgg                                                                             <  2:212872/139‑1 (MQ=255)                        aaaTCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGAt                                                                           <  1:398255/139‑1 (MQ=255)                              ggCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTAc                                                                     >  2:93989/1‑139 (MQ=255)                                          tACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTATGATCGACCTTAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCa                                                         <  2:681771/139‑1 (MQ=255)                                          tACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCa                                                         <  1:1008696/139‑1 (MQ=255)                                                 aCTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGtt                                                  >  2:229097/1‑139 (MQ=255)                                                     tACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAAtt                                              <  2:29291/139‑1 (MQ=255)                                                          cACGGGTATTTTTCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTTCCTGGGAGGTTTTAAGATCCACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGttt                                         <  2:531370/139‑1 (MQ=255)                                                            cGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAAtt                                              >  2:881038/1‑132 (MQ=255)                                                            cGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAAtt                                              <  1:881038/132‑1 (MQ=255)                                                               gTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGAt                                                                           <  1:1008800/100‑1 (MQ=255)                                                               gTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGAt                                                                           >  2:1008800/1‑100 (MQ=255)                                                               gTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATGc                                    >  1:454729/1‑139 (MQ=255)                                                                    ttGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATGCCTGCa                               >  1:656572/1‑139 (MQ=255)                                                                                       ataaCTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATGCCTGCACCTCCTGCGTCAATTGCCa            >  1:54635/1‑139 (MQ=255)                                                                                       ataaCTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGAACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATGCCTGCACCTCCTGCg                      >  1:312277/1‑129 (MQ=255)                                                                                       ataaCTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGAACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATGCCTGCACCTCCTGCg                      <  2:312277/129‑1 (MQ=255)                                                                                          aCTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGc                                                                 >  1:72473/1‑83 (MQ=255)                                                                                          aCTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGc                                                                 <  2:72473/83‑1 (MQ=255)                                                                                          aCTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATGCCTGCACCTCCTGCGTCAATTGCCagag         >  2:158150/1‑139 (MQ=255)                                                                                          aCTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATGCCTGCACCTCCTGCGTCAATTGCCagag         <  1:158150/139‑1 (MQ=255)                                                                                          aCTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATGCCTGCACCTCCTGCGTCAATTGCCagag         <  1:229597/139‑1 (MQ=255)                                                                                          aCTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATGCCTGCACCTCCTGCGTCAATTGCCagag         >  1:903362/1‑139 (MQ=255)                                                                                               gCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAATATCGACCTGAATCCCCg                                                                                          <  2:633951/53‑1 (MQ=255)                                                                                               gCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCg                                                                                          >  1:633951/1‑53 (MQ=255)                                                                                               gCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATGCCTGCACCTCCTGCg                      >  2:843428/1‑121 (MQ=255)                                                                                               gCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATGCCTGCACCTCCTGCg                      <  1:843428/121‑1 (MQ=255)                                                                                               gCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATGCCTGCACCTCCTGCGTCAATTGCCAGAGCGCAt    >  2:408630/1‑139 (MQ=255)                                                                                                 aaCCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACt                                                                                <  1:601899/61‑1 (MQ=255)                                                                                                 aaCCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACt                                                                                >  2:601899/1‑61 (MQ=255)                                                                                                 aaCCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATGCCTGCACCTCCTGCGTCAATTGCCAGAGCGCATCa  <  1:93989/139‑1 (MQ=255)                                                                                                   ccATAAGCGCCAAGTGCCGAGAAGGTTTTAAGCTCGACCTGAATCCCCGCACCACCACTCGGATCa                                                                         >  1:867848/1‑66 (MQ=255)                                                                                                   ccATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCTGATCa                                                                         <  2:867848/66‑1 (MQ=255)                                                                                                   ccATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATGCCTGCACCTCCTGCGTCAATTGCCagag         >  1:938158/1‑130 (MQ=255)                                                                                                   ccATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATGCCTGCACCTCCTGCGTCAATTGCCagag         <  2:938158/130‑1 (MQ=255)                                                                                                             cAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGAt                                                                           >  1:206844/1‑54 (MQ=255)                                                                                                             cAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGAt                                                                           <  2:206844/54‑1 (MQ=255)                                                                                                             cAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTAc                                                                     >  1:516473/1‑60 (MQ=255)                                                                                                             cAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTAc                                                                     <  2:516473/60‑1 (MQ=255)                                                                                                                   ccGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACTGTCAATCGTCAGAGCGTTAATTCGTTTCATGCCTGCAcc                             >  2:177514/1‑94 (MQ=255)                                                                                                                   ccGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATGCCTGCAcc                             <  1:177514/94‑1 (MQ=255)                                                                                                                                       |                                                                                                     CCGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGGCCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATGCCTGCACCTCCTGCGTCAATTGCCAGAGCGCATCA  >  NC_000913/2184318‑2184552

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