Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,823,664	T→C	V422A (GTG→GCG)	spoT →	bifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,823,664	0	T	C	100.0%	112.5 / NA	31	V422A (GTG→GCG)	spoT	bifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/15);  total (16/15)

GAATTTATCGAGAGCGTTAAATCCGATCTCTTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCG  >  NC_000913/3823519‑3823811                                                                                                                                                  |                                                                                                                                                    gAATTTATCGAGAGCGTTAAATCCGATCTCTTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGgc                                                                                                                                                 <  1:619141/150‑1 (MQ=255)   aTTTATCGAGAGCGTTAAATCCGATCTCTTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCa                                                                                                                                              <  2:455340/151‑1 (MQ=255)              gCGTTAAATCCGATCTCTTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAACTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCg                                                                                                                                   <  2:464363/151‑1 (MQ=255)               cGTTAAATCCGATCTCTTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGc                                                                                                                                  >  1:53282/1‑151 (MQ=255)                            ctctTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGc                                                                                                                     >  1:614407/1‑151 (MQ=255)                              ctTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTg                                                                                                                   >  1:562584/1‑151 (MQ=255)                                tCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCt                                                                                                                    <  2:408189/148‑1 (MQ=255)                                tCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGt                                                                                                                  <  1:44834/150‑1 (MQ=255)                                   cGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCa                                                                                                              <  1:64664/151‑1 (MQ=255)                                    ggATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGTAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAg                                                                                                             <  2:153786/151‑1 (MQ=255)                                                    ttCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTc                                                                                             >  2:421993/1‑151 (MQ=255)                                                       acacCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCaa                                                                                           >  2:549204/1‑150 (MQ=255)                                                          ccGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCg                                                                                       >  2:338518/1‑151 (MQ=255)                                                                   cgcATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCa                                                                              <  1:248723/151‑1 (MQ=255)                                                                             aGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTcc                                                                    <  2:206852/151‑1 (MQ=255)                                                                               ctgcctgcCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCgg                                                                  >  1:323999/1‑151 (MQ=255)                                                                                   ctgccGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGgcgc                                                              >  2:627394/1‑151 (MQ=255)                                                                                         gTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGccc                                                        <  1:164778/151‑1 (MQ=255)                                                                                               cGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGc                                                  <  1:202400/151‑1 (MQ=255)                                                                                                 cccGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCg                                                <  2:75088/151‑1 (MQ=255)                                                                                                      cGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTgg                                           >  2:196896/1‑151 (MQ=255)                                                                                                      cGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTg                                            <  1:351164/150‑1 (MQ=255)                                                                                                         cTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTg                                        >  2:338074/1‑151 (MQ=255)                                                                                                         cTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTg                                        >  1:593561/1‑151 (MQ=255)                                                                                                          tacGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGAACAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTg                                        >  2:155757/3‑150 (MQ=255)                                                                                                                tATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTg                                 >  2:273991/1‑151 (MQ=255)                                                                                                                    cAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGt                             >  2:218076/1‑151 (MQ=255)                                                                                                                      gTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTa                           <  1:589284/151‑1 (MQ=255)                                                                                                                      gTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTa                           <  1:445166/151‑1 (MQ=255)                                                                                                                                        cATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGcc          >  1:232618/1‑150 (MQ=255)                                                                                                                                         aTGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCAGCTGTCGCAGCCGCTTACCAGCGGTCAAAACGTTGAAATCATTACCGCTCCGGGCGCTAGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCa         >  2:339545/1‑150 (MQ=255)                                                                                                                                               gcgcGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTc   <  1:325054/150‑1 (MQ=255)                                                                                                                                               gcgcGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCg  >  1:430832/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                    GAATTTATCGAGAGCGTTAAATCCGATCTCTTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCG  >  NC_000913/3823519‑3823811

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