Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,518,841 A→T 7.2% T78T (ACT→ACA spoT ← bifunctional (p)ppGpp synthetase II and guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,518,8410AT 7.2% 238.1 / 9.7 83T78T (ACT→ACAspoTbifunctional (p)ppGpp synthetase II and guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base A (53/24);  new base T (4/2);  total (57/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CGACACCCCCTCTACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAG  >  minE/2518772‑2518909
                                                                     |                                                                    
cGACACCCCCTCTACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAg                                                                     >  1:3815817/1‑71 (MQ=255)
cGACACCCCCTCTACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAg                                                                     >  1:3425813/1‑71 (MQ=255)
 gACACCCCCTCTACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGt                                                                    <  1:3278142/71‑1 (MQ=255)
 gACACCCCCTCTACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGt                                                                    <  1:4274837/71‑1 (MQ=255)
        ccTCTACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTc                                                              >  1:3362916/1‑70 (MQ=255)
        ccTCTACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTc                                                              >  1:2871096/1‑70 (MQ=255)
          tctACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCa                                                             <  1:1580948/69‑1 (MQ=255)
          tctACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAAt                                                           >  1:3201122/1‑71 (MQ=255)
            tACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAAt                                                           <  1:934720/69‑1 (MQ=255)
             aCCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCAc                                                        >  1:3060684/1‑71 (MQ=255)
             aCCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCAc                                                        >  1:842526/1‑71 (MQ=255)
              ccAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATct                                                         <  1:2622718/69‑2 (MQ=255)
              ccAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCa                                                         <  1:3798712/69‑1 (MQ=255)
               cAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGt                                                      >  1:2977191/1‑71 (MQ=255)
                aGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCaa                                                            >  1:1465233/1‑64 (MQ=255)
                aGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCaa                                                            >  1:3358767/1‑64 (MQ=255)
                aGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGt                                                      >  1:2964796/1‑70 (MQ=255)
                 gCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCa                                                    <  1:1435286/71‑1 (MQ=255)
                 gCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCa                                                    <  1:2189596/71‑1 (MQ=255)
                   tCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATg                                                  >  1:3332155/1‑71 (MQ=255)
                    cGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATg                                                  >  1:2715722/1‑70 (MQ=255)
                     ggCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATgca                                                <  1:605792/71‑1 (MQ=255)
                       cGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATgcag                                               >  1:2694040/1‑70 (MQ=255)
                          cGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATgca                                                >  1:3796510/1‑66 (MQ=255)
                          cGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCagcagc                                           >  1:52649/1‑71 (MQ=255)
                          cGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCagcagc                                           >  1:828202/1‑71 (MQ=255)
                           gCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAgcg                                          >  1:940297/1‑71 (MQ=255)
                            cTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAAt                                                           >  1:3726286/1‑53 (MQ=255)
                            cTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAAt                                                           >  1:955951/1‑53 (MQ=255)
                            cTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAAt                                                           >  1:109650/1‑53 (MQ=255)
                            cTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAAt                                                           >  1:3283891/1‑53 (MQ=255)
                              tttACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCa                                                             <  1:192008/49‑1 (MQ=255)
                              tttACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCa                                                             <  1:1225846/49‑1 (MQ=255)
                                tACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGa                                                                      >  1:3918876/1‑38 (MQ=255)
                                  ccAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCAt                                   >  1:3928624/1‑71 (MQ=255)
                                  ccAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCAt                                   >  1:4001146/1‑71 (MQ=255)
                                  ccAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCAt                                   >  1:1009818/1‑71 (MQ=255)
                                          cTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTc                                                              >  1:271569/1‑36 (MQ=255)
                                          cTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTc                                                              >  1:1389926/1‑36 (MQ=255)
                                              tCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTc                          >  1:1291767/1‑68 (MQ=255)
                                               ccATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAg                      >  1:355024/1‑71 (MQ=255)
                                                cATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACg                                                       >  1:4266451/1‑37 (MQ=255)
                                                 atatCCTGGTATGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTc                    <  1:2124787/71‑1 (MQ=255)
                                                 atatCCTGGTAGGTGGCGGGAGTATCTTCAATCACGt                                                      >  1:4058352/1‑37 (MQ=255)
                                                 atatCCTGGTAGGTGGCGGGAGTATCTTCAATCACGt                                                      >  1:436797/1‑37 (MQ=255)
                                                 atatCCTGGTAGGTGGCGGGAGTATCTTCAATCACGt                                                      >  1:1637602/1‑37 (MQ=255)
                                                 atatCCTGGTAGGTGGCGGGAGTATCTTCAATCACGt                                                      >  1:3118645/1‑37 (MQ=255)
                                                 atatCCTGGTAGGTGGCGGGAGTATCTTCAATCACGt                                                      >  1:2769226/1‑37 (MQ=255)
                                                  tatCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTc                    >  1:3261600/1‑70 (MQ=255)
                                                  tatCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTc                    >  1:3861253/1‑70 (MQ=255)
                                                    tCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATgcagca                                             >  1:728838/1‑43 (MQ=255)
                                                    tCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATgca                                                >  1:2172637/1‑40 (MQ=255)
                                                    tCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGcagcag                                            >  1:2469247/1‑44 (MQ=255)
                                                    tCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGcagcag                                            >  1:1308013/1‑44 (MQ=255)
                                                     ccTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGt                >  1:492457/1‑71 (MQ=255)
                                                     ccTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGt                >  1:471984/1‑71 (MQ=255)
                                                     ccTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGt                >  1:633174/1‑71 (MQ=255)
                                                     ccTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGt                >  1:3383147/1‑71 (MQ=255)
                                                     ccTGGTAGGTGGAGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGt                >  1:2738902/1‑71 (MQ=255)
                                                      cTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGtt               <  1:3573452/71‑1 (MQ=255)
                                                      cTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGtt               <  1:945360/71‑1 (MQ=255)
                                                       tGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCg                              <  1:3272235/55‑1 (MQ=255)
                                                       tGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCg                              <  1:2022786/55‑1 (MQ=255)
                                                       tGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCg                              <  1:1956127/55‑1 (MQ=255)
                                                       tGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGtt               >  1:3466994/1‑70 (MQ=255)
                                                         gTAGGTGGCGGGTGTATCTTCAATCACATCATGCAGCAGCGCCGCCATCAGCg                              <  1:3813527/53‑1 (MQ=39)
                                                         gTAGGTGGCGGGTGTATCTTCAATCACATCATGCAGCAGCGCCGCCATCAGCg                              <  1:2014347/53‑1 (MQ=39)
                                                         gTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGttt              <  1:3314495/69‑1 (MQ=255)
                                                         gTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGttt              <  1:4276813/69‑1 (MQ=255)
                                                         gTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGttt              <  1:2595588/69‑1 (MQ=255)
                                                         gTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTc             <  1:2848010/70‑1 (MQ=255)
                                                           aGGTGGCGGGTGTATCTTCAATCACATCATGCAGCAGCGCCGCCATCAGCGtt                            >  1:3643209/1‑53 (MQ=39)
                                                           aGGTGGCGGGTGTATCTTCAATCACATCATGCAGCAGCGCCGCCATCAGCGtt                            >  1:3331001/1‑53 (MQ=39)
                                                           aGGTGGCGGGTGTATCTTCAATCACATCATGCAGCAGCGCCGCCATCAGCGtt                            >  1:3182362/1‑53 (MQ=39)
                                                           aGGTGGCGGGTGTATCTTCAATCACATCATGCAGCAGCGCCGCCATCAGCGtt                            >  1:985662/1‑53 (MQ=39)
                                                           aGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCAtc          >  1:3493867/1‑71 (MQ=255)
                                                            ggTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCAtct         >  1:2268266/1‑71 (MQ=255)
                                                              tGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGttt              <  1:3811393/64‑1 (MQ=255)
                                                              tGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGttt              <  1:904977/64‑1 (MQ=255)
                                                               ggCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCgg      >  1:3122182/1‑71 (MQ=255)
                                                                gCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCa            >  1:716077/1‑64 (MQ=255)
                                                                gCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCa            >  1:2512466/1‑64 (MQ=255)
                                                                gCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTAGAGTTTCa            >  1:1914530/1‑64 (MQ=255)
                                                                   ggAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAg  <  1:2001533/71‑1 (MQ=255)
                                                                     |                                                                    
CGACACCCCCTCTACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAG  >  minE/2518772‑2518909

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: