Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,913,035	A→G	L206P (CTC→CCC)	relA ←	(p)ppGpp synthetase I/GTP pyrophosphokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,913,035	0	A	G	100.0%	172.7 / NA	46	L206P (CTC→CCC)	relA	(p)ppGpp synthetase I/GTP pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/20);  total (26/20)

TTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCA  >  NC_000913/2912891‑2913181                                                                                                                                                 |                                                                                                                                                   ttCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCTTGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTaa                                                                                                                                               <  1:23054/150‑1 (MQ=255)          ccATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGt                                                                                                                                     <  1:230599/151‑1 (MQ=255)            aTACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTTGTGTTCGCGGTCGCGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAg                                                                                                                                         <  2:668126/145‑1 (MQ=255)              acacTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATc                                                                                                                                 >  1:487319/1‑151 (MQ=255)                acTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCtt                                                                                                                               <  1:177018/151‑1 (MQ=255)                   tCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTcc                                                                                                                             <  1:37163/150‑1 (MQ=255)                   tCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTc                                                                                                                              <  1:636712/149‑1 (MQ=255)                    ccGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTcc                                                                                                                             >  2:379052/1‑149 (MQ=255)                     cGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGt                                                                                                                          >  1:164038/1‑151 (MQ=255)                           aaCGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCAt                                                                                                                    <  1:540957/151‑1 (MQ=255)                                 ttCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGt                                                                                                              >  2:108219/1‑151 (MQ=255)                                           aTCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTc                                                                                                    >  1:163481/1‑151 (MQ=255)                                              tcAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCg                                                                                                  >  1:606617/1‑150 (MQ=255)                                               cAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTc                                                                                                    <  2:623514/147‑1 (MQ=255)                                                                aCGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTg                                                                               >  2:231808/1‑151 (MQ=255)                                                                 cGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGc                                                                              >  1:313758/1‑151 (MQ=255)                                                                 cGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGc                                                                              <  2:58068/151‑1 (MQ=255)                                                                  gAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCg                                                                             <  1:323971/151‑1 (MQ=255)                                                                    aCTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGt                                                                            >  2:656060/1‑150 (MQ=255)                                                                      tctTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGt                                                                            >  2:116157/1‑148 (MQ=255)                                                                               tAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTac                                                                >  2:141134/1‑151 (MQ=255)                                                                                 gtgtTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTacac                                                              >  1:169872/1‑151 (MQ=255)                                                                                      cgcgGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACActctg                                                          <  2:84828/150‑2 (MQ=255)                                                                                      cgcgGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCttt                                                         >  1:95748/1‑151 (MQ=255)                                                                                       gcgGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCtttt                                                        >  1:96254/1‑151 (MQ=255)                                                                                       gcgGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCtttt                                                        >  1:581792/1‑151 (MQ=255)                                                                                         ggTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTgc                                                      >  2:428458/1‑151 (MQ=255)                                                                                          gTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTgcc                                                     >  1:657735/1‑151 (MQ=255)                                                                                            cGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTgccgc                                                   <  2:642732/151‑1 (MQ=255)                                                                                             gagaCGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTgccgcc                                                  <  1:277732/151‑1 (MQ=255)                                                                                              agaCGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCa                                                 <  2:516946/151‑1 (MQ=255)                                                                                                       tCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGtt                                        >  2:458108/1‑151 (MQ=255)                                                                                                        cATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTc                                       <  2:204466/151‑1 (MQ=255)                                                                                                           gcagcaGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATc                                    >  2:367327/1‑151 (MQ=255)                                                                                                             agcagTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCtt                                  <  2:628325/151‑1 (MQ=255)                                                                                                                  tttGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGgc                             >  2:106895/1‑151 (MQ=255)                                                                                                                    tGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGgcgc                           >  2:655113/1‑151 (MQ=255)                                                                                                                     ggCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCa                          <  1:461257/151‑1 (MQ=255)                                                                                                                      gCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCAt                         <  1:102919/151‑1 (MQ=255)                                                                                                                       cAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATc                        >  2:523323/1‑151 (MQ=255)                                                                                                                        aaTTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCt                       <  2:174546/151‑1 (MQ=255)                                                                                                                                tGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTc               >  2:144985/1‑151 (MQ=255)                                                                                                                                 gTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTcg              <  2:260767/151‑1 (MQ=255)                                                                                                                                  tATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTcgc             >  1:133757/1‑151 (MQ=255)                                                                                                                                   aTTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTcgcg            >  1:547788/1‑151 (MQ=255)                                                                                                                                    ttCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCgcgc           >  1:282883/1‑151 (MQ=255)                                                                                                                                          tggatgggGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATg      <  1:50206/150‑1 (MQ=255)                                                                                                                                             atgggGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCa  <  1:372687/151‑1 (MQ=255)                                                                                                                                             atgggGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCa  <  2:648451/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                   TTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCA  >  NC_000913/2912891‑2913181

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