Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,912,885 A→G I256T (ATC→ACC)  relA ← (p)ppGpp synthetase I/GTP pyrophosphokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,912,8850AG97.2% 124.2 / ‑3.8 36I256T (ATC→ACC) relA(p)ppGpp synthetase I/GTP pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (16/19);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TGGCGATAGTGAGTGTGCACTATCCCCAGTGCGGCATAGCAATCCTGTAAACGCTCGGCGACAATACGTACCGCACGCACATCAAACAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGG  >  NC_000913/2912747‑2913034
                                                                                                                                          |                                                                                                                                                     
tGGCGATAGTGAGTGTGCACTATCCCCAGTGCGGCATAGCAATCCTGTAAACGCTCGGCGACAATACGTACCGCACGCACATCAAACAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGGTGTGTTTCGgac                                                                                                                                           >  2:93333/1‑151 (MQ=255)
      tAGTGAGTGTGCACTATCCCCAGTGCGGCATAGCAATCCTGTAAACGCTCGGCGACAATACGTACCGCACGCACATCAAACAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGGTGTGTTTCGGACGACCAt                                                                                                                                     >  1:317449/1‑151 (MQ=255)
        gtgagtgTGCACTATCCCCAGTGCGGCATAGCAATCCTGTAAACGCTCGGCGACAATACGTACCGCACGCACATCAAACAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGGTGTGTTTCGGACGACCATac                                                                                                                                   >  2:453769/1‑151 (MQ=255)
            gtgtGCACTATCCCCAGTGCGGCATAGCAATCCTGTAAACGCTCGGCGACAATACGTACCGCACGCACATCAAACAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGGTGTGTTTCGGACGACCATACACtt                                                                                                                               <  2:249134/151‑1 (MQ=255)
                  aCTATCCCCAGTGCGGCATAGCAATCCTGTAAACGCTCGGCGACAATACGTACCGCACGCACATCAAACAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCtt                                                                                                                         <  2:147262/151‑1 (MQ=255)
                       ccccAGTGCGGCATAGCAATCCTGTAAACGCTCGGCGACAATACGTACCGCACGCACATCAAACAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAAc                                                                                                                     <  2:131868/150‑1 (MQ=255)
                       ccccAGTGCGGCATAGCAATCCTGTAAACGCTCGGCGACAATACGTACCGCACGCACATCAAACAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAAc                                                                                                                     <  1:93333/150‑1 (MQ=255)
                              gCGGCATAGCAATCCTGTAAACGCTCGGCGACAATACGTACCGCACGCACATCAAACAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAg                                                                                                             <  1:328700/151‑1 (MQ=255)
                                  cATAGCAATCCTGTAAACGCTCGGCGACAATACGTACCGCACGCACATCAAACAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCttt                                                                                                         <  1:583496/151‑1 (MQ=255)
                                               tAAACGCTCGGCGACAATACGTACCGCACGCACATCAAACAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAgcgcgc                                                                                            >  1:300582/1‑151 (MQ=255)
                                                   cGCTCGGCGACAATACGTACCGCACGCACATCAAACAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAgcgcgc                                                                                            >  1:108268/1‑147 (MQ=255)
                                                   cGCTCGGCGACAATACGTACCGCACGCACATCAAACAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGAt                                                                                        <  1:427899/151‑1 (MQ=255)
                                                        ggCGACAATACGTACCGCACGCACATCAAACAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCa                                                                                   <  2:185852/151‑1 (MQ=255)
                                                           gACAATACGTACCGCACGCACATCAAACAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACg                                                                                <  2:190978/151‑1 (MQ=255)
                                                              aaTACGTACCGCACGCACATCAAACAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAAc                                                                             >  2:192120/1‑151 (MQ=255)
                                                                   gTACCGCACGCACATCAAACAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTc                                                                        <  2:135268/151‑1 (MQ=255)
                                                                     aCCGCACGCACATCAAACAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCg                                                                       <  2:574510/150‑1 (MQ=255)
                                                                           cgcacATCAAACAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAgt                                                                >  2:53866/1‑151 (MQ=255)
                                                                                   aaaCAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCgg                                                        >  1:439417/1‑151 (MQ=255)
                                                                                   aaaCAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCgg                                                        >  2:404014/1‑151 (MQ=255)
                                                                                       aGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCg                                                     >  1:483035/1‑150 (MQ=255)
                                                                                        gCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCgag                                                   <  1:148714/151‑1 (MQ=255)
                                                                                         cTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCg                                                     >  2:77253/1‑148 (MQ=255)
                                                                                               aaaGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGt                                            >  2:614469/1‑151 (MQ=255)
                                                                                                        gTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATgcagc                                   <  1:374826/151‑1 (MQ=255)
                                                                                                          tCTTTTTCTGCATTTTACGCCAGATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATgcagca                                  <  1:301267/150‑1 (MQ=255)
                                                                                                                tCTGCATTTTACGCCAGATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTgg                            <  1:162899/150‑1 (MQ=255)
                                                                                                                       tttACGCCAGATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGt                    >  1:341643/1‑151 (MQ=255)
                                                                                                                        ttACGCCAGATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGt                    <  1:611154/150‑1 (MQ=255)
                                                                                                                         tACGCCAGATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGgtt                  >  2:301031/1‑151 (MQ=255)
                                                                                                                             ccAGATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTAt              <  2:134482/151‑1 (MQ=255)
                                                                                                                               aGATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTc            <  2:47204/151‑1 (MQ=255)
                                                                                                                                gATGCTGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCg           >  1:19933/1‑151 (MQ=255)
                                                                                                                                     tGTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTtgg      >  1:463004/1‑151 (MQ=255)
                                                                                                                                      gTAGGTGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTtg       <  2:16609/149‑1 (MQ=255)
                                                                                                                                         gATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTtggatgg  >  1:347396/1‑151 (MQ=255)
                                                                                                                                          |                                                                                                                                                     
TGGCGATAGTGAGTGTGCACTATCCCCAGTGCGGCATAGCAATCCTGTAAACGCTCGGCGACAATACGTACCGCACGCACATCAAACAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGG  >  NC_000913/2912747‑2913034

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