Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,882,393 C→A 19.3% A247A (GCG→GCT relA ← (p)ppGpp synthetase I/GTP pyrophosphokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,882,3930CA19.3% 188.4 / 28.5 93A247A (GCG→GCTrelA(p)ppGpp synthetase I/GTP pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base C (19/56);  new base A (12/6);  total (31/62)
Fisher's exact test for biased strand distribution p-value = 1.65e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.36e-01

GAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCG  >  minE/1882330‑1882459
                                                               |                                                                  
gAGGTTCTTTTTCTGCATTTTACGCCATATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                                               <  1:2534461/69‑1 (MQ=255)
gAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                                               <  1:2491062/69‑1 (MQ=255)
gAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                                               <  1:154682/69‑1 (MQ=255)
gAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                                               <  1:2982653/69‑1 (MQ=255)
gAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                                               <  1:1626772/69‑1 (MQ=255)
gAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                                               <  1:1467626/69‑1 (MQ=255)
gAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                                               <  1:1440192/69‑1 (MQ=255)
gAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                                               <  1:3160928/69‑1 (MQ=255)
gAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                                               <  1:1423512/69‑1 (MQ=255)
gAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                                               <  1:1322730/69‑1 (MQ=255)
gAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                                               <  1:1942484/69‑1 (MQ=255)
gAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                                               <  1:45524/69‑1 (MQ=255)
gAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                                               <  1:268112/69‑1 (MQ=255)
gAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                                               <  1:825567/69‑1 (MQ=255)
gAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                                               <  1:1123705/69‑1 (MQ=255)
 aGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                                               <  1:635906/68‑1 (MQ=255)
 aGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                                               <  1:2839903/68‑1 (MQ=255)
         tttCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                                               <  1:3219854/60‑1 (MQ=255)
           tCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCa                                                     <  1:1710023/68‑1 (MQ=255)
           tCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCa                                                     <  1:2200608/68‑1 (MQ=255)
           tCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCa                                                     <  1:1078158/68‑1 (MQ=255)
           tCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCa                                                     <  1:2900846/68‑1 (MQ=255)
           tCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCa                                                     <  1:643976/68‑1 (MQ=255)
           tCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCa                                                     <  1:809253/68‑1 (MQ=255)
           tCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCa                                                     <  1:196666/68‑1 (MQ=255)
           tCTGAATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGAACATACACTTCCGCTTTAACGCCTTCa                                                     <  1:2664282/68‑1 (MQ=255)
            cTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAg                                                    <  1:3031318/68‑1 (MQ=255)
            cTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAg                                                    <  1:688625/68‑1 (MQ=255)
            cTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAg                                                    <  1:2521838/68‑1 (MQ=255)
                 ttttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTc                                               <  1:2172306/68‑1 (MQ=255)
                 ttttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTc                                               <  1:199599/68‑1 (MQ=255)
                 ttttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTc                                               <  1:3267771/68‑1 (MQ=255)
                 ttttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTc                                               <  1:1082936/68‑1 (MQ=255)
                 ttttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCa                                              >  1:2714847/1‑69 (MQ=255)
                 ttttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCa                                              >  1:2781318/1‑69 (MQ=255)
                 ttttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCa                                              >  1:247757/1‑69 (MQ=255)
                 ttttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCa                                              >  1:1583189/1‑69 (MQ=255)
                 ttttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCa                                              >  1:1578815/1‑69 (MQ=255)
                   ttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCAt                                             >  1:2563348/1‑68 (MQ=255)
                   ttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCAt                                             >  1:1484043/1‑68 (MQ=255)
                   ttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCAt                                             >  1:2252410/1‑68 (MQ=255)
                   ttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCAt                                             >  1:2415964/1‑68 (MQ=255)
                        ccAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAgc                                       >  1:3315393/1‑69 (MQ=255)
                        ccAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAgc                                       >  1:1276040/1‑69 (MQ=255)
                        ccAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAgc                                       >  1:3351888/1‑69 (MQ=255)
                           gATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                                               <  1:3281604/42‑1 (MQ=255)
                            aTGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAg                                                    <  1:1875347/52‑1 (MQ=255)
                             tGCTGTAGATGTGTTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCa                                  <  1:2499759/69‑1 (MQ=255)
                             tGCTGTAGATGTGTTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCa                                  <  1:863065/69‑1 (MQ=255)
                             tGCTGTAGATGTGTTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCa                                  <  1:531738/69‑1 (MQ=255)
                             tGCTGTAGATGTGTTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCa                                  <  1:1429997/69‑1 (MQ=255)
                             tGCTGTAGATGTGTTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCa                                  <  1:2594890/69‑1 (MQ=255)
                             tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                  <  1:1285844/69‑1 (MQ=255)
                             tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                  <  1:1278626/69‑1 (MQ=255)
                             tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                  <  1:3024413/69‑1 (MQ=255)
                             tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                  <  1:1744129/69‑1 (MQ=255)
                             tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                  <  1:262711/69‑1 (MQ=255)
                             tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                  <  1:175995/69‑1 (MQ=255)
                             tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                  <  1:1798290/69‑1 (MQ=255)
                             tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                  <  1:2312177/69‑1 (MQ=255)
                             tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                  <  1:1878815/69‑1 (MQ=255)
                             tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                  <  1:1882531/69‑1 (MQ=255)
                             tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                  <  1:875258/69‑1 (MQ=255)
                             tGCTGTAGATGTGTTTCGGACGACCATACACTTACGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                  <  1:2426769/69‑1 (MQ=255)
                              gCTGTAGATGTGTTTCGGTCGACCATACACTTCAGCTTTAACACCTACAGCTTTCATCTCAGCGCGCa                                  <  1:1157505/68‑1 (MQ=255)
                              gCTGTAGATGTGTTTCGGACGACCATAGACTTCCGCTTTAACGCCTTCa                                                     <  1:2582384/49‑1 (MQ=255)
                                       gtgtTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                  <  1:622342/59‑1 (MQ=255)
                                       gtgtTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                  <  1:3369300/59‑1 (MQ=255)
                                       gtgtTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                  <  1:965747/59‑1 (MQ=255)
                                       gtgtTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                  <  1:969798/59‑1 (MQ=255)
                                       gtgtTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                  <  1:2339363/59‑1 (MQ=255)
                                       gtgtTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                  <  1:2515459/59‑1 (MQ=255)
                                       gtgtTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                  <  1:1625744/59‑1 (MQ=255)
                                         gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAcaa                      >  1:2025154/1‑67 (MQ=255)
                                         gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAcaa                      >  1:2155422/1‑67 (MQ=255)
                                         gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAcaa                      >  1:942004/1‑67 (MQ=255)
                                         gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAcaa                      >  1:1224723/1‑67 (MQ=255)
                                         gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAcaa                      >  1:163754/1‑67 (MQ=255)
                                         gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAcaa                      >  1:3423460/1‑67 (MQ=255)
                                         gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAcaa                      >  1:3396844/1‑67 (MQ=255)
                                         gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAcaa                      >  1:1486866/1‑67 (MQ=255)
                                         gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAcaa                      >  1:3334062/1‑67 (MQ=255)
                                         gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAca                       >  1:1220239/1‑67 (MQ=255)
                                         gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAca                       >  1:3308541/1‑67 (MQ=255)
                                         gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCAACTCAGCGCGCAGATGGCCGAcaa                      >  1:1750404/1‑67 (MQ=255)
                                                acgacCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCa                          <  1:1371645/58‑1 (MQ=255)
                                                         cacTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAgtg      >  1:1543971/1‑69 (MQ=255)
                                                         cacTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAgtg      >  1:923619/1‑69 (MQ=255)
                                                             tCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTcg  >  1:2114627/1‑69 (MQ=255)
                                                             tCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTcg  >  1:2560801/1‑69 (MQ=255)
                                                             tCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTcg  >  1:1205101/1‑69 (MQ=255)
                                                             tCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTcg  >  1:3053387/1‑69 (MQ=255)
                                                             tCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTcg  >  1:3196501/1‑69 (MQ=255)
                                                               |                                                                  
GAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCG  >  minE/1882330‑1882459

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