Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	638,614	G→T	T14K (ACA→AAA)	thiD ←	phosphomethylpyrimidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	638,614	0	G	T	100.0%	90.5 / NA	28	T14K (ACA→AAA)	thiD	phosphomethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (20/8);  total (20/8)

CATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTT  >  CP000730/638477‑638745                                                                                                                                          |                                                                                                                                    cATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCACACCTGCGCCAGCACTTttgt                                                                                                                                  <  1:26459/141‑1 (MQ=255)       ttGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAAc                                                                                                                            >  1:134920/1‑141 (MQ=255)             aTGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCaa                                                                                                                      >  1:54688/1‑141 (MQ=255)                        tCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAAct                                                                                                           <  1:20720/141‑1 (MQ=255)                         cGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGACCGTTTTCAAATCTGCTTGCATGCTTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAActt                                                                                                          >  1:42097/1‑141 (MQ=255)                         cGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAActt                                                                                                          >  1:278009/1‑141 (MQ=255)                         cGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAActt                                                                                                          >  1:292439/1‑141 (MQ=255)                                    gTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGcc                                                                                               <  2:130479/141‑1 (MQ=255)                                     tATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCa                                                                                              <  1:278805/141‑1 (MQ=255)                                                  aGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTc                                                                                 >  1:38639/1‑141 (MQ=255)                                                         aTGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACata                                                                          >  2:151409/1‑141 (MQ=255)                                                             cAGTTAAAGCGACAATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACataataa                                                                      >  1:90494/1‑141 (MQ=255)                                                                 tAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATAtt                                                                  <  1:33446/141‑1 (MQ=255)                                                                 tAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATAtt                                                                  >  2:231572/1‑141 (MQ=255)                                                                   aaGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATAtttt                                                                >  2:11871/1‑139 (MQ=255)                                                                   aaGCGACCATGCCATACGTATCTAATTCTTGGAACGGTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGt                                                                >  2:113805/1‑141 (MQ=255)                                                                       gACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCGAGCTCTTTTGTCAGAACCGGCATTTTTAAAAAGTTTCTTTAAATCCATTGAGCTTCACTCCGACATAATAATATTGTATtc                                                            >  1:173718/1‑141 (MQ=255)                                                                        aCCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATtca                                                           >  1:252956/1‑141 (MQ=255)                                                                        aCCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATtca                                                           >  1:50378/1‑141 (MQ=255)                                                                        aCCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATtca                                                           >  1:84791/1‑141 (MQ=255)                                                                             gCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcata                                                      <  1:64084/141‑1 (MQ=255)                                                                             gCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcata                                                      <  1:90069/141‑1 (MQ=255)                                                                                      aTCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTTTTTTTATTTTCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTa                                             >  2:229757/1‑141 (MQ=255)                                                                                         tAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAAcc                                          >  2:127591/1‑141 (MQ=255)                                                                                          aaTTCTTGGGACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCt                                         <  2:28466/141‑1 (MQ=255)                                                                                             tCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAAGATTGTATTCATCATATCATTTTTAACCTAAt                                      <  2:271658/141‑1 (MQ=255)                                                                                                tGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGa                                   >  2:624/1‑141 (MQ=255)                                                                                                                 tGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCa                  >  2:40980/1‑141 (MQ=255)                                                                                                                                 ccAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttgtt  >  1:284274/1‑141 (MQ=255)                                                                                                                                          |                                                                                                                                    CATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTT  >  CP000730/638477‑638745

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