Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,270,398	A→G	N54S (AAC→AGC)	yjbQ →	thiamine phosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,270,398	0	A	G	100.0%	174.7 / NA	47	N54S (AAC→AGC)	yjbQ	thiamine phosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/16);  total (31/16)

CGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGCTGGGAACATCGCTTGTAT  >  NC_000913/4270254‑4270541                                                                                                                                                 |                                                                                                                                                cGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCg                                                                                                                                            >  1:91166/1‑150 (MQ=255) cGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGa                                                                                                                                           <  1:168781/151‑1 (MQ=255)  gCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCCTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGAt                                                                                                                                          >  1:68504/1‑151 (MQ=255)       cGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCAc                                                                                                                                     <  2:6753/151‑1 (MQ=255)         cTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCg                                                                                                                                   >  2:553910/1‑151 (MQ=255)               gcCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGc                                                                                                                              >  2:681615/1‑150 (MQ=255)                    aTCTCGTGGGGTTCATCTGGTAACGGCTGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTCACATTGGCTTACTGCATCTGCTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACCCCACGAc                                                                                                                        >  2:44129/1‑151 (MQ=255)                            ggTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCg                                                                                                                <  1:452894/151‑1 (MQ=255)                                 cATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGttttt                                                                                                           <  2:718570/151‑1 (MQ=255)                                   tCTGGTAACGGGTGAAATTCTGAATCAGCGGGCGGAGATGCCGCGCGGTAACATCGGCGTACTGCTGCTGTTGCTGCAACATACCTCCGCCGCGCTGACACTTAATGAGAGGTGCGATCCCACCGTGCGCCACGACATGGAGCGtgtttt                                                                                                          >  2:291601/1‑150 (MQ=255)                                             gATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCAcc                                                                                                >  1:513766/1‑150 (MQ=255)                                              aTGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGt                                                                                              >  2:344727/1‑151 (MQ=255)                                                gAAATTCTGAATCAGGTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCATCATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGtt                                                                                             <  1:465319/150‑1 (MQ=255)                                                          aTCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACgg                                                                                  >  2:732513/1‑151 (MQ=255)                                                            cAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGa                                                                                 <  1:121781/150‑1 (MQ=255)                                                            cAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCACCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGa                                                                                 >  1:164845/1‑150 (MQ=255)                                                             aGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGaa                                                                                >  1:155991/1‑150 (MQ=255)                                                                tggctgATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAAtt                                                                             >  2:364277/1‑150 (MQ=255)                                                                tggctgATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAAtt                                                                             >  1:675502/1‑150 (MQ=255)                                                                 ggctgATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTAt                                                                           >  1:602715/1‑151 (MQ=255)                                                                 ggctgATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTAt                                                                           >  1:381644/1‑151 (MQ=255)                                                                     gATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAg                                                                        <  2:615068/150‑1 (MQ=255)                                                                     gATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGc                                                                       >  2:700249/1‑151 (MQ=255)                                                                      atatGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGc                                                                       >  1:432614/1‑150 (MQ=255)                                                                        atGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCAt                                                                     <  2:249093/150‑1 (MQ=255)                                                                        atGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATg                                                                    >  1:512179/1‑151 (MQ=255)                                                                          gCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGAc                                                                  <  1:5993/151‑1 (MQ=255)                                                                           ccGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACt                                                                 <  2:256189/151‑1 (MQ=255)                                                                           ccGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACt                                                                 >  2:310388/1‑151 (MQ=255)                                                                             gcgcgTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTAt                                                               >  1:588037/1‑151 (MQ=255)                                                                               gcgTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATg                                                              >  1:502825/1‑150 (MQ=255)                                                                                          ggCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACg                                                  <  1:740812/151‑1 (MQ=255)                                                                                                 tGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGcc                                           >  2:66558/1‑151 (MQ=255)                                                                                                   cATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCtt                                         >  1:20965/1‑151 (MQ=255)                                                                                                        gTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCat                                    >  2:709029/1‑151 (MQ=255)                                                                                                        gTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCat                                    >  2:605584/1‑151 (MQ=255)                                                                                                            cttcAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATc                                 <  2:427038/147‑1 (MQ=255)                                                                                                               cAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCaa                               <  1:153071/149‑1 (MQ=255)                                                                                                               cAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAAt                             <  1:690654/151‑1 (MQ=255)                                                                                                                     aCCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCa                        <  2:393981/150‑1 (MQ=255)                                                                                                                        tCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGc                    >  1:299834/1‑151 (MQ=255)                                                                                                                          cGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGCTg                  >  1:648442/1‑151 (MQ=255)                                                                                                                          cGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGCTg                  >  1:315058/1‑151 (MQ=255)                                                                                                                                   gACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGCTGGGAACATc          <  2:685842/150‑1 (MQ=255)                                                                                                                                      acTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGCTGGGAACATCGCtt      >  1:407309/1‑151 (MQ=255)                                                                                                                                       cTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGCTGGGAACATCGCTTg     >  1:596715/1‑151 (MQ=255)                                                                                                                                          aaTGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGCTGGGAACATCGCTTGTAt  >  2:204613/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                CGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGCTGGGAACATCGCTTGTAT  >  NC_000913/4270254‑4270541

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