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%	130.0 / NA	38	A22V (GCC→GTC)	thiD	hydroxy‑methylpyrimidine kinase and hydroxy‑phosphomethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/17);  total (21/17)

CGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGGCCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATGCCTGCACCTCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACAAAACTGCCTGGCCCCT  >  NC_000913/2184319‑2184583                                                                                                                                      |                                                                                                                                    cGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGa                                                                                                                                >  1:1197975/1‑139 (MQ=255) cGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGa                                                                                                                                >  2:881910/1‑139 (MQ=255) cGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGa                                                                                                                                >  2:626846/1‑139 (MQ=255)             ggcggcGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGcacca                                                                                                                    <  1:881910/139‑1 (MQ=255)             ggcggcGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGcacca                                                                                                                    <  2:1211383/139‑1 (MQ=255)                  cgACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACt                                                                                                               <  2:180528/139‑1 (MQ=255)                  cgACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACt                                                                                                               <  1:698333/139‑1 (MQ=255)                    aCAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCg                                                                                                             >  1:613115/1‑139 (MQ=255)                      aaaaTCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGa                                                                                                           <  2:1197975/139‑1 (MQ=255)                      aaaaTCAGGCTCAATGCGATACACCGACTGTACGCCACGGGCATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCgaa                                                                                                           <  2:742764/139‑3 (MQ=255)                             ggCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTAc                                                                                                    <  1:368015/139‑1 (MQ=255)                                                aCTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGtt                                                                                 >  2:318763/1‑139 (MQ=255)                                                aCTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGtt                                                                                 >  2:1132111/1‑139 (MQ=255)                                                   gTACGCCACGGGTATTTTGCGCCACCAGTGTAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAAt                                                                              >  1:484516/1‑139 (MQ=255)                                                        ccACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGtt                                                                         >  1:847356/1‑139 (MQ=255)                                                            gggTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCAt                                                                     >  1:33743/1‑139 (MQ=255)                                                                     gcgcCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATGCCTGCAcc                                                            <  1:626846/139‑1 (MQ=255)                                                                      cgcCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGcac                                                                                                                      >  1:279865/1‑80 (MQ=255)                                                                      cgcCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGcac                                                                                                                      <  2:279865/80‑1 (MQ=255)                                                                        ccaccaGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCaa                                                                                              <  1:979216/102‑1 (MQ=255)                                                                        ccaccaGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCaa                                                                                              >  2:979216/1‑102 (MQ=255)                                                                                         aCTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATGCCTGCACCTCCTGCGTCAATTGCCagag                                        >  2:213240/1‑139 (MQ=255)                                                                                         aCTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATGCCTGCACCTCCTGCGTCAATTGCCagag                                        <  2:421321/139‑1 (MQ=255)                                                                                         aCTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATGCCTGCACCTCCTGCGTCAATTGCCagag                                        <  2:847356/139‑1 (MQ=255)                                                                                             aGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAAt                                                                                                                              <  2:980971/49‑1 (MQ=255)                                                                                             aGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAAt                                                                                                                              >  1:980971/1‑49 (MQ=255)                                                                                              gCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGcacca                                                                                                                    <  2:493589/58‑1 (MQ=255)                                                                                              gCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGcacca                                                                                                                    >  1:493589/1‑58 (MQ=255)                                                                                              gCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATg                                                                    >  2:264061/1‑106 (MQ=255)                                                                                              gCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATg                                                                    <  1:264061/106‑1 (MQ=255)                                                                                              gCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATGCCTGCACCTCCTGCGTCAATTGCCAGAGCGCAt                                   >  1:776053/1‑139 (MQ=255)                                                                                               cAACCATAAGCGCCCAGTGCCGAGAAGGGTTTAAGATCGACCTGAATCCCCGCACCACCAGGCGGGGCAGTACCGGCAATCGTCGGAGATTTAATTCGTTTCCTGCCTGCACCTCCTGCGTCAATTTCCAGATGGCATc                                  >  2:919808/1‑139 (MQ=255)                                                                                                            cAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACt                                                                                                               >  2:1055735/1‑49 (MQ=255)                                                                                                            cAAGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACt                                                                                                               <  1:1055735/49‑1 (MQ=255)                                                                                                              aGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTc                                                                                                              >  1:974748/1‑48 (MQ=255)                                                                                                              aGTGCCGAGAAGGTTTTAAGATCGACCTGAATCCCCGCACCACCACTc                                                                                                              <  2:974748/48‑1 (MQ=255)                                                                                                                            tttAAGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATGCCTGCACCTCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACAAAACTGCCTGGcc     >  1:24929/1‑139 (MQ=255)                                                                                                                               aaGATCGACCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATGCCTGCACCTCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACAAAACTGCCTGGCCCCt  >  1:1066252/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                    CGAATCGAGCTGGGCGGCGACAAAATCAGGCTCAATGCGATACACCGACTGTACGCCACGGGTATTTTGCGCCACCAGTGCAGTAATAACTGAGCAACCATAAGCGCCAAGTGCCGAGAAGGTTTTAAGATCGGCCTGAATCCCCGCACCACCACTCGGATCAGTACCGGCAATCGTCAGAGCGTTAATTCGTTTCATGCCTGCACCTCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACAAAACTGCCTGGCCCCT  >  NC_000913/2184319‑2184583

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