Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,910,539 C→A 81.2% R590L (CGC→CTC)  relA ← (p)ppGpp synthetase I/GTP pyrophosphokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,910,5390CA81.2% 85.5 / 27.4 48R590L (CGC→CTC) relA(p)ppGpp synthetase I/GTP pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base C (3/6);  new base A (7/32);  total (10/38)
Fisher's exact test for biased strand distribution p-value = 3.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGG  >  W3110S.gb/2910492‑2910602
                                               |                                                               
tCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGAGGTTTTGCGGCGTGTAGCtttt                                           >  1:1549147/1‑70 (MQ=255)
      tGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGAGGTTTTGCGGCGTGTAGCttttt                                          >  1:1689572/1‑65 (MQ=255)
      tGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGAGGTTTTGCGGCGTGTAGCttttt                                          >  1:325817/1‑65 (MQ=255)
      tGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTTCTg                                     >  1:830298/1‑70 (MQ=255)
      tGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTg                                     >  1:1643466/1‑70 (MQ=255)
      tGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTg                                     >  1:1946063/1‑70 (MQ=255)
      tGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTg                                     >  1:278124/1‑70 (MQ=255)
             acCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCt                             >  1:1094080/1‑71 (MQ=255)
             acCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCt                             >  1:1133181/1‑71 (MQ=255)
             acCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCt                             >  1:27068/1‑71 (MQ=255)
                           cgcgACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCgcggcggcg               <  1:687347/71‑1 (MQ=255)
                           cgcgACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCgcggcggcg               <  1:559509/71‑1 (MQ=255)
                           cgcgACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCgcggcggcg               <  1:1410514/71‑1 (MQ=255)
                           cgcgACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCgcggcggcg               <  1:300127/71‑1 (MQ=255)
                            gcgACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCgcggcggcg               <  1:1312127/70‑1 (MQ=255)
                            gcgACCGTTATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGTCGGCGt              <  1:302713/71‑1 (MQ=255)
                            gcgACCGTTATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGTCGGCGt              <  1:3003/71‑1 (MQ=255)
                            gcgACCGTTATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGt              <  1:1915545/71‑1 (MQ=255)
                            gcgACCGTTATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGt              <  1:32134/71‑1 (MQ=255)
                            gcgACCGTTATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGt              <  1:681046/71‑1 (MQ=255)
                            gcgACCGTTATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGt              <  1:1937369/71‑1 (MQ=255)
                            gcgACCGTTATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGt              <  1:946191/71‑1 (MQ=255)
                            gcgACCGTTATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGt              <  1:1012041/71‑1 (MQ=255)
                            gcgACCGTTATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGt              <  1:1878176/71‑1 (MQ=255)
                            gcgACCGTTATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGt              <  1:1647183/71‑1 (MQ=255)
                            gcgACCGTTATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGt              <  1:1611395/71‑1 (MQ=255)
                            gcgACCGTTATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGt              <  1:1446981/71‑1 (MQ=255)
                            gcgACCGTTATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGt              <  1:1296136/71‑1 (MQ=255)
                            gcgACCGTTATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGt              <  1:1293874/71‑1 (MQ=255)
                            gcgACCGTTATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGt              <  1:1136480/71‑1 (MQ=255)
                             cgACCGTTATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGt              <  1:1100024/70‑1 (MQ=255)
                                  gTTATCTTTACTGAGGTTTTTCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGt              <  1:952479/65‑1 (MQ=255)
                                  gTTATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGt              <  1:18577/65‑1 (MQ=255)
                                  gTTATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGt              <  1:967346/65‑1 (MQ=255)
                                  gTTATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGt              <  1:936646/65‑1 (MQ=255)
                                  gTTATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGt              <  1:925010/65‑1 (MQ=255)
                                  gTTATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGt              <  1:1341582/65‑1 (MQ=255)
                                  gTTATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGt              <  1:1521309/65‑1 (MQ=255)
                                  gTTATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGt              <  1:1573676/65‑1 (MQ=255)
                                  gTTATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGAGt              <  1:1435759/65‑1 (MQ=255)
                                   ttATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGt              <  1:838123/64‑1 (MQ=255)
                                   ttATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGctc       <  1:1279138/71‑1 (MQ=255)
                                   ttATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGctc       <  1:678114/71‑1 (MQ=255)
                                   ttATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGctc       <  1:519833/71‑1 (MQ=255)
                                   ttATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGctc       <  1:109908/71‑1 (MQ=255)
                                   ttATCTTTACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGAGGCGTCCTGctc       <  1:123169/71‑1 (MQ=255)
                                          tACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCgg  <  1:322649/69‑1 (MQ=255)
                                          tACTGAGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGt              <  1:609807/57‑1 (MQ=255)
                                               |                                                               
TCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGG  >  W3110S.gb/2910492‑2910602

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