Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,823,683 C→T R428R (CGC→CGT 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_0009133,823,6830CT94.3% 100.8 / ‑0.7 35R428R (CGC→CGTspoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (22/11);  total (23/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.68e-01

CTTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGT  >  NC_000913/3823548‑3823812
                                                                                                                                       |                                                                                                                                 
ctTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAg                                                                                                                                >  2:107787/1‑139 (MQ=255)
 ttCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGc                                                                                                                               >  2:764805/1‑139 (MQ=255)
       gATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTAcc                                                                                                                         >  2:37648/1‑139 (MQ=255)
           agaTTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCt                                                                                                                     >  1:1085901/1‑139 (MQ=255)
                tACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGc                                                                                                                <  2:982892/139‑1 (MQ=255)
                 aCGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCa                                                                                                               >  1:954135/1‑139 (MQ=255)
                      ttCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGc                                                                                                          >  1:424552/1‑139 (MQ=255)
                                          tGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGtt                                                                                      <  2:935425/139‑1 (MQ=255)
                                                            tGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCg                                                                    >  2:417785/1‑139 (MQ=255)
                                                              cAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCggg                                                                  >  1:358668/1‑139 (MQ=255)
                                                                aCGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGgcg                                                                <  2:424552/139‑1 (MQ=255)
                                                                  gCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGc                                                                                                                >  1:1142926/1‑89 (MQ=255)
                                                                  gCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGc                                                                                                                <  2:1142926/89‑1 (MQ=255)
                                                                          aCTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGaa                                                      <  1:107787/139‑1 (MQ=255)
                                                                               gCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACTAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCg                                                 <  2:1085901/139‑1 (MQ=255)
                                                                                 ttATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCt                                               <  1:417785/139‑1 (MQ=255)
                                                                                           cATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACt                                     >  1:1215755/1‑139 (MQ=255)
                                                                                              aCCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTg                                  >  1:433382/1‑139 (MQ=255)
                                                                                                  atatCGGTCATGCCTGCGTGGACGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCTGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGt                              <  2:721654/139‑1 (MQ=255)
                                                                                                   tatCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGtt                             >  2:694768/1‑139 (MQ=255)
                                                                                                   tatCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGtt                             >  2:776649/1‑139 (MQ=255)
                                                                                                            tGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCt                                                              <  2:465097/97‑1 (MQ=255)
                                                                                                            tGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCt                                                              >  1:465097/1‑97 (MQ=255)
                                                                                                            tGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGaaa                    >  1:919316/1‑139 (MQ=255)
                                                                                                            tGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGaaa                    >  2:1074316/1‑139 (MQ=255)
                                                                                                            tGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGaaa                    >  1:1212719/1‑139 (MQ=255)
                                                                                                             gCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCt                                          >  2:261360/1‑116 (MQ=255)
                                                                                                             gCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCt                                          <  1:261360/116‑1 (MQ=255)
                                                                                                             gCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAg                   >  1:1052510/1‑139 (MQ=255)
                                                                                                                  cGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAgcgcgc              >  2:1044879/1‑139 (MQ=255)
                                                                                                                  cGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTAAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAgcgcgc              >  1:771450/1‑139 (MQ=255)
                                                                                                                  cGTGGGCGCACGCGTTGAACGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGACGCTTGGCTGAACTTTGTCGTTAGCTCGAAAgcgcgc              >  1:1011002/1‑139 (MQ=255)
                                                                                                                      ggCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCa          >  1:1071811/1‑139 (MQ=255)
                                                                                                                      ggCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCa          <  1:1021165/139‑1 (MQ=255)
                                                                                                                              cgTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGt  <  1:764805/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                 
CTTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGT  >  NC_000913/3823548‑3823812

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: