Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,823,162 C→T H255Y (CAC→TAC)  spoT → bifunctional (p)ppGpp synthase/hydrolase SpoT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,823,1620CT96.3% 88.8 / ‑6.0 27H255Y (CAC→TAC) spoTbifunctional (p)ppGpp synthase/hydrolase SpoT
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (12/14);  minor base A (0/1);  total (12/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGG‑CGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAA‑TGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACA‑GCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATG  >  NC_000913/3822915‑3823389
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                     
tttAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGG‑CGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACt                                                                                                                                                                                                                                    >  1:222160/1‑251 (MQ=255)
    ggTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGG‑CGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGAt                                                                                                                                                                                                                                <  1:147082/251‑1 (MQ=255)
       aTCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGG‑CGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCa                                                                                                                                                                                                                              <  1:141339/250‑1 (MQ=255)
       aTCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGG‑CGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCa                                                                                                                                                                                                                              <  1:67262/250‑1 (MQ=255)
       aTCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGG‑CGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCAt                                                                                                                                                                                                                             >  2:99207/1‑251 (MQ=255)
                        cGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGG‑CGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTcc                                                                                                                                                                                                            <  1:51219/251‑1 (MQ=255)
                         gAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGG‑CGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCcg                                                                                                                                                                                                           >  2:46763/1‑251 (MQ=255)
                          aaCTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGG‑CATTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCcgc                                                                                                                                                                                                          <  2:222372/251‑1 (MQ=255)
                                                    cTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGG‑CGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAA‑TGATTCTGACACCTg                                                                                                                                                                               >  1:216663/1‑251 (MQ=255)
                                                                   tatcgtatactgatgcaggtggagAAACCCACGCGCGGCACCCGTAAGGACATGATCCCGAAGATTCTTTCTGCCATCGCAGGG‑CGTTTGCAGCAAGCGGGAATACTGTGCCGCGGCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAACTGGTGCTAAAAAAGCAGCGTTTTTACTCCATCATCGACATCTACGCTTTCCGCGCGATCGTCAA‑TGATTCTGCCACCTGTTATCGCGTGCTgcg                                                                                                                                                                >  2:208278/19‑249 (MQ=255)
                                                                              cAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGG‑CGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAA‑TGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGcac                                                                                                                                                      <  1:154989/250‑1 (MQ=255)
                                                                                           aaagcaccgcgcgGAAACCGCAAAGGGATGATCCAGAAGATTCTTTCTGAAATCACAGGG‑CGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTACTCGATTTACTGCAAAATGGTGCTCACAGAGCAGCGTTTTTACTCGCTCATGGACATCTACGCTTTCCGCGAGATCG‑CAACTGATTCTGACACCTGTTGTCGCGTGCTGGGCCAGATGC‑CACGCCTGTACACGcc                                                                                                                                       >  2:306198/8‑251 (MQ=255)
                                                                                           aaaGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGG‑CGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAA‑TGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACA‑GCCTGTACAAGcc                                                                                                                                       <  1:317189/251‑1 (MQ=255)
                                                                                                      gggCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGG‑CGTTTGCAGGAATCGGGAATACCGTGCCGAGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAA‑TGATTCTGACACCTGTTATCGCGTGAGGGGCCAGATGCACA‑GCCTGTACACGCCGCGTCCGGGcc                                                                                                                            <  2:210473/250‑1 (MQ=255)
                                                                                                             cGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGG‑CGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAA‑TGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACA‑GCCTGTACAAGCCGCGTCTGGGCCGCGTGa                                                                                                                      >  2:9873/1‑250 (MQ=255)
                                                                                                                      atgatCCAGAAGATTCTTTCTGAAATCGAAGGG‑CGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAA‑TGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACA‑GCCTGTACAAGCCGCGTCCGGGCCGCGTGAACGACTatat                                                                                                            >  2:234422/1‑251 (MQ=255)
                                                                                                                        gatCCAGAAGATTCTTTCTGAAATCGAAGGG‑CGTTTGCAGGACGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAA‑TGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACA‑GCCTGTACAAGCCGCGTCCGGGCCGCGAGAAAGACTATATCg                                                                                                          >  2:278163/1‑251 (MQ=255)
                                                                                                                        gatCCAGAAGATTCTTTCTGAAATCGAAGGG‑CGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAA‑TGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACA‑GCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCg                                                                                                          <  2:251006/251‑1 (MQ=255)
                                                                                                                              gaataTTCTTTCGAAAATCGAAGGGGTGTTTG‑GGGAGGCGGGAATACGGGGCCGCGTCAGGGGTCGGGAGAAGCATCTTTATTCGATTTACTGAAAAATGGTGCTCAAAGAGCAGCGTTTTAACTCGATAATGGACCTCTACGCTTTCCGAGTGATCGTCAA‑TGATTCGGACACCTGTTATCGCGTGCTGGGACAGATGCACA‑GCCGATACAAGCCGCGTCCGGGCCGCGTGAAAGACTatatttccat                                                                                                      <  2:164126/245‑7 (MQ=255)
                                                                                                                                         ccgAAATCGAAGGG‑CGTGTGCAGGAAGCGGGAATACCGTGCCGTGTCAGTGGTCGCGAGAAGCATGTTTATTCGATTTTCTTCAAAATGGTGCTCAAAGATCAGCGTTTTTACTCTATCATGGACATCTACGCTTTCCGCGTGATCGTCAA‑TGATTCTGACACCTGTTAGCGCGTGCTGGGCCAGATTCACA‑GCTATTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAAc                                                                                         <  2:118193/249‑1 (MQ=255)
                                                                                                                                                      g‑CGTGTGCAGGAAGCGGGAATACCGTGCCGCGTCATTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAA‑TGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACA‑GCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCtt                                                                            <  1:274144/251‑1 (MQ=255)
                                                                                                                                                                                         aGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAA‑TGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACA‑GCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGcc                                          >  2:37695/1‑251 (MQ=255)
                                                                                                                                                                                          gcggTCGCGAGACGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCCTGGACATCTACGCTTTCCGCGTGATCGTCAA‑TGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACA‑GCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCCCGGCGTGCCg                                         >  2:29801/3‑251 (MQ=255)
                                                                                                                                                                                                              tATTCGAGTTACGGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCTTGATCGTCAA‑TGATTCGGACACCTGTTATCGCGTGCTGGGCCATATGCACA‑GCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTa                      <  2:13172/250‑1 (MQ=255)
                                                                                                                                                                                                                   gATTTACTGCAAAATGGGGCTCAAAGATCAGCTTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAA‑TGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACA‑GCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAg                <  2:179927/251‑1 (MQ=255)
                                                                                                                                                                                                                           gCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAA‑TGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACA‑GCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGAc        >  1:159828/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ggTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAA‑TGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACA‑GCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATg  <  1:28689/250‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                     
TTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGG‑CGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAA‑TGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACA‑GCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATG  >  NC_000913/3822915‑3823389

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