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%	188.1 / NA	48	L206P (CTC→CCC)	relA	(p)ppGpp synthetase I/GTP pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/23);  total (25/23)

GTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGA  >  NC_000913/2912887‑2913175                                                                                                                                                     |                                                                                                                                             gtgtTTCGGTCGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGgggg                                                                                                                                            <  1:52436/151‑1 (MQ=255)       cGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACgg                                                                                                                                      >  2:595742/1‑151 (MQ=255)            gacCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCa                                                                                                                                 <  2:75208/151‑1 (MQ=255)                acacacTTCCGCTCTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGGCGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTa                                                                                                                              <  2:696020/148‑1 (MQ=255)                  acacTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAAt                                                                                                                            <  1:461101/150‑1 (MQ=255)                  acacTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATc                                                                                                                           >  2:752915/1‑151 (MQ=255)                      ttCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTcc                                                                                                                       >  2:150934/1‑151 (MQ=255)                        ccGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCa                                                                                                                      <  2:339729/150‑1 (MQ=255)                               aaCGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCa                                                                                                               <  2:659972/150‑1 (MQ=255)                               aaCGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCAt                                                                                                              >  2:714464/1‑151 (MQ=255)                                      tCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGt                                                                                                        <  1:685557/150‑1 (MQ=255)                                            ttCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCg                                                                                                  <  1:545427/150‑1 (MQ=255)                                              cATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGAtt                                                                                               <  1:165380/151‑1 (MQ=255)                                              cATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGAtt                                                                                               <  1:700063/151‑1 (MQ=255)                                              cATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGAtt                                                                                               >  2:226974/1‑151 (MQ=255)                                                tctcAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTcc                                                                                             >  1:567046/1‑151 (MQ=255)                                                   cAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTc                                                                                              <  2:696032/147‑1 (MQ=255)                                                       gcgcAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACgg                                                                                      >  1:205252/1‑151 (MQ=255)                                                               gACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCa                                                                               <  2:913600/150‑1 (MQ=255)                                                                    aCGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTg                                                                         >  1:43543/1‑151 (MQ=255)                                                                      gAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCg                                                                       >  2:325399/1‑151 (MQ=255)                                                                          tctTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGt                                                                      >  1:174400/1‑148 (MQ=255)                                                                             tCGATGTAGTGTTCGCGGTCGGGGCGACGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAGGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATg                                                                 <  1:622906/150‑1 (MQ=255)                                                                             tCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATg                                                                 <  1:33476/150‑1 (MQ=255)                                                                                aTGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTgg                                                             >  1:49625/1‑151 (MQ=255)                                                                                     gtgtTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTacac                                                        >  2:31813/1‑151 (MQ=255)                                                                                            cgGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGGGCGTAGATGTTGGTACACTCTTTTg                                                 >  2:283339/1‑151 (MQ=255)                                                                                             ggTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTgc                                                >  2:118798/1‑151 (MQ=255)                                                                                                 gagaCGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTgccgcc                                            >  1:296082/1‑151 (MQ=255)                                                                                                 gagaCGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTgccg                                              <  1:773866/149‑1 (MQ=255)                                                                                                        cGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTacac                                     >  1:562328/1‑151 (MQ=255)                                                                                                        cGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTacac                                     >  1:354122/1‑151 (MQ=255)                                                                                                               gcagcaGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCAt                               <  1:382294/150‑1 (MQ=255)                                                                                                               gcagcaGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATc                              <  1:280256/151‑1 (MQ=255)                                                                                                                  gcagTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTc                           <  1:246564/151‑1 (MQ=255)                                                                                                                       ttGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGgc                       >  1:159271/1‑150 (MQ=255)                                                                                                                        tGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGgcg                      <  2:844310/150‑1 (MQ=255)                                                                                                                         ggCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGgcgc                     <  1:390481/150‑1 (MQ=255)                                                                                                                              ttCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCttt               >  2:558044/1‑151 (MQ=255)                                                                                                                                cGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCtttt              >  1:832056/1‑150 (MQ=255)                                                                                                                                 gTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTAc            >  2:589765/1‑151 (MQ=255)                                                                                                                                   ttGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACtt          >  1:68015/1‑151 (MQ=255)                                                                                                                                    tGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTc         >  1:845374/1‑151 (MQ=255)                                                                                                                                     gTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTcg        <  2:758507/151‑1 (MQ=255)                                                                                                                                     gTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTc         <  1:283373/150‑1 (MQ=255)                                                                                                                                         tCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCgcgc     >  2:553254/1‑150 (MQ=255)                                                                                                                                         tCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCgcgc     <  1:660968/150‑1 (MQ=255)                                                                                                                                          cGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCgcgc     >  2:710584/1‑149 (MQ=255)                                                                                                                                           ggTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGa  <  1:865183/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                             GTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGA  >  NC_000913/2912887‑2913175

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: ‑