Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,882,378 A→T 17.7% R252R (CGT→CGA 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,3780AT17.7% 199.7 / 30.5 101R252R (CGT→CGArelA(p)ppGpp synthetase I/GTP pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base A (25/58);  new base T (12/6);  total (37/64)
Fisher's exact test for biased strand distribution p-value = 6.03e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

AAACAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGA  >  minE/1882312‑1882439
                                                                  |                                                             
aaaCAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGgacg                                                             >  1:1990615/1‑69 (MQ=255)
aaaCAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGgacg                                                             >  1:2301953/1‑69 (MQ=255)
aaaCAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGgacg                                                             >  1:253589/1‑69 (MQ=255)
aaaCAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGgacg                                                             >  1:2616926/1‑69 (MQ=255)
aaaCAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGgacg                                                             >  1:1273125/1‑69 (MQ=255)
aaaCAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGgacg                                                             >  1:2628525/1‑69 (MQ=255)
aaaCAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGgacg                                                             >  1:1568022/1‑69 (MQ=255)
      cTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATa                                                       <  1:2687033/69‑1 (MQ=255)
      cTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATa                                                       <  1:1185577/69‑1 (MQ=255)
           caAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATacac                                                    >  1:876905/1‑67 (MQ=255)
                  gAGGTTCTTTTTCTGCATTTTACGCCATATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                           <  1:2534461/69‑1 (MQ=255)
                  gAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                           <  1:1440192/69‑1 (MQ=255)
                  gAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                           <  1:45524/69‑1 (MQ=255)
                  gAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                           <  1:1467626/69‑1 (MQ=255)
                  gAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                           <  1:1626772/69‑1 (MQ=255)
                  gAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                           <  1:825567/69‑1 (MQ=255)
                  gAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                           <  1:2491062/69‑1 (MQ=255)
                  gAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                           <  1:1942484/69‑1 (MQ=255)
                  gAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                           <  1:1423512/69‑1 (MQ=255)
                  gAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                           <  1:1322730/69‑1 (MQ=255)
                  gAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                           <  1:2982653/69‑1 (MQ=255)
                  gAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                           <  1:3160928/69‑1 (MQ=255)
                  gAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                           <  1:1123705/69‑1 (MQ=255)
                   aGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                           <  1:2839903/68‑1 (MQ=255)
                        cTTTTTCTGCATTTTAGGCCAGATGCTGTAGATGTGTTTCGGACGACCATa                                                       >  1:1560124/1‑51 (MQ=255)
                        cTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCAt                                                        >  1:2044405/1‑50 (MQ=255)
                        cTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATa                                                       >  1:1950669/1‑51 (MQ=255)
                        cTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATa                                                       >  1:1924082/1‑51 (MQ=255)
                        cTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATa                                                       >  1:715573/1‑51 (MQ=255)
                        cTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATa                                                       >  1:1099655/1‑51 (MQ=255)
                        cTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATa                                                       >  1:111086/1‑51 (MQ=255)
                        cTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATa                                                       >  1:1037081/1‑51 (MQ=255)
                        cTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATa                                                       >  1:2977402/1‑51 (MQ=255)
                        cTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATa                                                       >  1:1297442/1‑51 (MQ=255)
                        cTTTTTCTGCATTTTACGCCAGATGCCGTAGATGTGTTTCGGACGACCATa                                                       >  1:2637611/1‑51 (MQ=255)
                           tttCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                           <  1:3219854/60‑1 (MQ=255)
                             tCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCa                                 <  1:2200608/68‑1 (MQ=255)
                             tCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCa                                 <  1:2900846/68‑1 (MQ=255)
                             tCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCa                                 <  1:643976/68‑1 (MQ=255)
                             tCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCa                                 <  1:1078158/68‑1 (MQ=255)
                             tCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCa                                 <  1:196666/68‑1 (MQ=255)
                             tCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCa                                 <  1:1710023/68‑1 (MQ=255)
                             tCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCa                                 <  1:809253/68‑1 (MQ=255)
                             tCTGAATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGAACATACACTTCCGCTTTAACGCCTTCa                                 <  1:2664282/68‑1 (MQ=255)
                              cTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAg                                <  1:2521838/68‑1 (MQ=255)
                              cTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAg                                <  1:3031318/68‑1 (MQ=255)
                              cTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAg                                <  1:688625/68‑1 (MQ=255)
                                   ttttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTc                           <  1:2172306/68‑1 (MQ=255)
                                   ttttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTc                           <  1:1082936/68‑1 (MQ=255)
                                   ttttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTc                           <  1:199599/68‑1 (MQ=255)
                                   ttttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCa                          >  1:2781318/1‑69 (MQ=255)
                                   ttttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCa                          >  1:247757/1‑69 (MQ=255)
                                   ttttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCa                          >  1:2714847/1‑69 (MQ=255)
                                   ttttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCa                          >  1:1578815/1‑69 (MQ=255)
                                   ttttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCa                          >  1:1583189/1‑69 (MQ=255)
                                     ttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATacac                                                    >  1:1136484/1‑41 (MQ=255)
                                     ttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCAt                         >  1:2252410/1‑68 (MQ=255)
                                     ttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCAt                         >  1:2415964/1‑68 (MQ=255)
                                     ttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCAt                         >  1:1484043/1‑68 (MQ=255)
                                     ttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCAt                         >  1:2563348/1‑68 (MQ=255)
                                          ccAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAgc                   >  1:3351888/1‑69 (MQ=255)
                                          ccAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAgc                   >  1:3315393/1‑69 (MQ=255)
                                             gATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCttt                                           <  1:3281604/42‑1 (MQ=255)
                                              aTGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAg                                <  1:1875347/52‑1 (MQ=255)
                                               tGCTGTAGATGTGTTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCa              <  1:2594890/69‑1 (MQ=255)
                                               tGCTGTAGATGTGTTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCa              <  1:863065/69‑1 (MQ=255)
                                               tGCTGTAGATGTGTTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCa              <  1:531738/69‑1 (MQ=255)
                                               tGCTGTAGATGTGTTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCa              <  1:2499759/69‑1 (MQ=255)
                                               tGCTGTAGATGTGTTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCa              <  1:1429997/69‑1 (MQ=255)
                                               tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa              <  1:3024413/69‑1 (MQ=255)
                                               tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa              <  1:1882531/69‑1 (MQ=255)
                                               tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa              <  1:1278626/69‑1 (MQ=255)
                                               tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa              <  1:1285844/69‑1 (MQ=255)
                                               tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa              <  1:875258/69‑1 (MQ=255)
                                               tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa              <  1:262711/69‑1 (MQ=255)
                                               tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa              <  1:1744129/69‑1 (MQ=255)
                                               tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa              <  1:2312177/69‑1 (MQ=255)
                                               tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa              <  1:175995/69‑1 (MQ=255)
                                               tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa              <  1:1878815/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:965747/59‑1 (MQ=255)
                                                         gtgtTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa              <  1:622342/59‑1 (MQ=255)
                                                         gtgtTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa              <  1:969798/59‑1 (MQ=255)
                                                         gtgtTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa              <  1:3369300/59‑1 (MQ=255)
                                                         gtgtTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa              <  1:2339363/59‑1 (MQ=255)
                                                         gtgtTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa              <  1:1625744/59‑1 (MQ=255)
                                                           gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAcaa  >  1:3423460/1‑67 (MQ=255)
                                                           gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAcaa  >  1:3396844/1‑67 (MQ=255)
                                                           gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAcaa  >  1:3334062/1‑67 (MQ=255)
                                                           gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAcaa  >  1:2025154/1‑67 (MQ=255)
                                                           gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAcaa  >  1:163754/1‑67 (MQ=255)
                                                           gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAcaa  >  1:1486866/1‑67 (MQ=255)
                                                           gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAcaa  >  1:942004/1‑67 (MQ=255)
                                                           gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAcaa  >  1:1224723/1‑67 (MQ=255)
                                                           gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAca   >  1:3308541/1‑67 (MQ=255)
                                                           gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAca   >  1:1220239/1‑67 (MQ=255)
                                                           gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCAACTCAGCGCGCAGATGGCCGAcaa  >  1:1750404/1‑67 (MQ=255)
                                                                  acgacCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCa      <  1:1371645/58‑1 (MQ=255)
                                                                  |                                                             
AAACAGCTCATCAAAGGCGAGGTTCTTTTTCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGA  >  minE/1882312‑1882439

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: