Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,190,510	T→G	intergenic (+185/+225)	yjaZ → / ← thiH	stationary phase growth adaptation protein/tyrosine lyase, involved in thiamine‑thiazole moiety synthesis

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,190,510	0	T	G	100.0%	179.3 / NA	47	intergenic (+185/+225)	yjaZ/thiH	stationary phase growth adaptation protein/tyrosine lyase, involved in thiamine‑thiazole moiety synthesis
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (23/24);  total (23/24)

ATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTTACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAA  >  NC_000913/4190362‑4190656                                                                                                                                                     |                                                                                                                                                   aTCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGAc                                                                                                                                                  >  1:615797/1‑151 (MQ=255)  tctcTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGAct                                                                                                                                                 >  2:434045/1‑151 (MQ=255)  tctcTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGATGAGTATTGAct                                                                                                                                                 >  1:610874/1‑151 (MQ=255)       ggTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGActctct                                                                                                                                             <  1:159679/150‑1 (MQ=255)       ggTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGActctc                                                                                                                                              <  2:341767/149‑1 (MQ=255)       ggTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGActctc                                                                                                                                              <  1:201538/149‑1 (MQ=255)       ggTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGActct                                                                                                                                               <  1:750156/148‑1 (MQ=255)            tCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACaa                                                                                                                                        <  2:448804/150‑1 (MQ=255)                  cGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACAATCAGAtc                                                                                                                                 <  1:181917/151‑1 (MQ=255)                    ccATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACAATCAGAtctc                                                                                                                               >  2:441462/1‑151 (MQ=255)                    ccATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACAATCAGAtctc                                                                                                                               <  2:676409/151‑1 (MQ=255)                    ccATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACAATCAGAtctc                                                                                                                               >  1:122897/1‑151 (MQ=255)                      aTTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACAATCAGATCTCtt                                                                                                                             <  2:365531/151‑1 (MQ=255)                         aCCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACAATCAGATCTCTTTc                                                                                                                           <  2:624751/150‑1 (MQ=255)                             gACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACAATCAGATCTCTTTCATTg                                                                                                                       >  2:591596/1‑150 (MQ=255)                                         aCCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGAt                                                                                                          <  2:433381/151‑1 (MQ=255)                                             tGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGc                                                                                                       <  2:516939/150‑1 (MQ=255)                                               cGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTc                                                                                                    >  2:787452/1‑151 (MQ=255)                                                 aTAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAg                                                                                                  <  1:265647/151‑1 (MQ=255)                                                        cAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGc                                                                                           >  1:837296/1‑151 (MQ=255)                                                          ggTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGCAt                                                                                         >  2:515698/1‑151 (MQ=255)                                                                     gTCTGAATTTCTTCAGGCTGCTGCATCCTGGAAGATGAGAACCTGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACCATCAGATCTCTTTCATTGCTCAACAGGCGCTGGCTTCAGACTTTGCATTACGGAAtttt                                                                              >  2:208327/1‑151 (MQ=255)                                                                         gAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGCATTACGGAAttttt                                                                             <  2:203405/148‑1 (MQ=255)                                                                              tcttcAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGCATTACGGAATTTTTAAGAAAg                                                                      <  2:410595/150‑1 (MQ=255)                                                                                  cAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGCATTACGGAATTTTTAAGAAAGGCAgg                                                                 >  1:374000/1‑151 (MQ=255)                                                                                  cAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGCATTACGGAATTTTTAAGAAAGGCAg                                                                  <  1:589142/150‑1 (MQ=255)                                                                                   aGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGCATTACGGAATTTTTAAGAAAGGCAggg                                                                >  2:690402/1‑151 (MQ=255)                                                                                         ctgcATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAAc                                                          >  2:707177/1‑151 (MQ=255)                                                                                           gcATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGa                                                        <  1:703200/151‑1 (MQ=255)                                                                                               ccTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGgaa                                                    <  1:301238/151‑1 (MQ=255)                                                                                                 tGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGga                                                     >  1:271056/1‑148 (MQ=255)                                                                                                  ggAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGgaagaa                                                 <  1:66497/151‑1 (MQ=255)                                                                                                         gagaACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCtttt                                           <  1:787342/150‑1 (MQ=255)                                                                                                           gaACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGc                                        <  2:328479/151‑1 (MQ=255)                                                                                                                    ttCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGaca                                >  2:672941/1‑150 (MQ=255)                                                                                                                     tCTTATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGaca                                >  2:615540/1‑149 (MQ=255)                                                                                                                       ttATTTTCGTCTCTATCATAGTTGAGTATTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACAt                            >  1:420143/1‑151 (MQ=255)                                                                                                                                 ctctATCATAGTTGAGTATTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGAc                  >  1:566268/1‑151 (MQ=255)                                                                                                                                  tctATCATAGTTGAGTATTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACg                 >  1:668668/1‑151 (MQ=255)                                                                                                                                   ctATCATAGTTGAGTATTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACgg                >  1:319611/1‑151 (MQ=255)                                                                                                                                     aTCATAGTTGAGTATTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGc               <  2:144805/150‑1 (MQ=255)                                                                                                                                       cATAGTTGAGTATTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACt            >  2:471337/1‑151 (MQ=255)                                                                                                                                           gTTGAGTATTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTAcc         >  2:171839/1‑150 (MQ=255)                                                                                                                                                gTATTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaaaa   <  1:415386/151‑1 (MQ=255)                                                                                                                                                  aTTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaaaaa  >  1:843797/1‑150 (MQ=255)                                                                                                                                                  aTTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaaaaa  <  1:780588/150‑1 (MQ=255)                                                                                                                                                  aTTGACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaaaaa  <  1:120577/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   ATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTCTATCATAGTTGAGTATTTACTCTCTTACAATCAGATCTCTTTCATTGCTCAACAGGCGATGGCTTCAGACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAA  >  NC_000913/4190362‑4190656

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