Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,518,142 G→A 41.0% I311I (ATC→ATT 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,1420GA41.0% 31.6 / 53.6 78I311I (ATC→ATTspoTbifunctional (p)ppGpp synthetase II and guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (25/21);  new base A (15/17);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.40e-01

TCCGCCATCTGGTCCATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTGGAATGGCGATA  >  minE/2518080‑2518193
                                                              |                                                   
tcCGCCATCTGGTCCATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGatca                                                 >  1:2724280/1‑67 (MQ=255)
tcCGCCATCTGGTCCATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGatca                                                 >  1:291809/1‑67 (MQ=255)
tcCGCCATCTGGTCCATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGatca                                                 >  1:2886505/1‑67 (MQ=255)
tcCGCCATCTGGTCCATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGatca                                                 >  1:153101/1‑67 (MQ=255)
tcCGCCATCTGGTCCATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGatca                                                 >  1:3018411/1‑67 (MQ=255)
tcCGCCATCTGGTCCATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGatca                                                 >  1:1455995/1‑67 (MQ=255)
tcCGCCATCTGGTCCATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGatca                                                 >  1:1895664/1‑67 (MQ=255)
tcCGCCATCTGGTCCATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGatca                                                 >  1:820750/1‑67 (MQ=255)
tcCGCCATCTGGTCCATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGatca                                                 >  1:2012655/1‑67 (MQ=255)
tcCGCCATCTGGTCCATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGatca                                                 >  1:2714629/1‑67 (MQ=255)
tcCGCCATCTGGTCCATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGatca                                                 >  1:2351983/1‑67 (MQ=255)
tcCGCCATCTGGTCCATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGatca                                                 >  1:2475577/1‑67 (MQ=255)
         tGGTCCATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGgtgt                                        >  1:866150/1‑67 (MQ=255)
         tGGTCCATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGtgtg                                       >  1:2942070/1‑68 (MQ=255)
         tGGTCCATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGtgtg                                       >  1:2832846/1‑68 (MQ=255)
         tGGTCCATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGtgtg                                       >  1:1757230/1‑68 (MQ=255)
         tGGTCCATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGtgtg                                       >  1:2476816/1‑68 (MQ=255)
         tGGTCCATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGtgtg                                       >  1:3114660/1‑68 (MQ=255)
         tGGTCCATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGtgtg                                       >  1:94013/1‑68 (MQ=255)
         tGGTCCATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGtgtg                                       >  1:946086/1‑68 (MQ=255)
             ccATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAg                                  <  1:1481377/69‑1 (MQ=255)
             ccATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAg                                  <  1:1089028/69‑1 (MQ=255)
             ccATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAg                                  <  1:936324/69‑1 (MQ=255)
                 atCTTCGGTACGGATCTGGACCTCAACCGGCACGCCGTGCGGGCCAATCATCGAGGTGTGCAAAGACt                               >  1:1032556/1‑68 (MQ=255)
                 atCTTCGGTACGGATCTGGACCTCAACCGGCACGCCGTGCGGGCCAATCATCGAGGTGTGCAAAGACt                               >  1:2959038/1‑68 (MQ=255)
                 atCTTCGGTACGGATCTGGACCTCAACCGGCACGCCGTGCGGGCCAATCATCGAGGTGTGCAAAGACt                               >  1:321963/1‑68 (MQ=255)
                 atCTTCGGTACGGATCTGGACCTCAACCGGCACGCCGTGCGGGCCAATCATCGAGGTGTGCAAAGACt                               >  1:2788601/1‑68 (MQ=255)
                 atCTTCGGTACGGATCTGGACCTCAACCGGCACGCCGTGCGGGCCAATCATCGAGGTGTGCAAAGACt                               >  1:2778625/1‑68 (MQ=255)
                 atCTTCGGTACGGATCTGGACCTCAACCGGCACGCCGTGCGGGCCAATCATCGAGGTGTGCAAAGACt                               >  1:329637/1‑68 (MQ=255)
                 atCTTCGGTACGGATCTGGACCTCAACCGGCACGCCGTGCGGGCCAATCATCGAGGTGTGCAAAGACt                               >  1:923661/1‑68 (MQ=255)
                 atCTTCGGTACGGATCTGGACCTCAACCGGCACGCCGTGCGGGCCAATCATCGAGGTGTGCAAAGACt                               >  1:254291/1‑68 (MQ=255)
                 atCTTCGGTACGGATCTGGACCTCAACCGGCACGCCGTGCGGGCCAATCATCGAGGTGTGCAAAGACt                               >  1:248084/1‑68 (MQ=255)
                 atCTTCGGTACGGATCTGGACCTCAACCGGCACGCCGTGCGGGCCAATCATCGAGGTGTGCAAAGACt                               >  1:1067210/1‑68 (MQ=255)
                 atCTTCGGTACGGATCTGGACCTCAACCGGCACGCCGTGCGGGCCAATCATCGAGGTGTGCAAAGACt                               >  1:1318826/1‑68 (MQ=255)
                 atCTTCGGTACGGATCTGGACCTCAACCGGCACGCCGTGCGGGCCAATCATCGAGGTGTGCAAAGACt                               >  1:1331348/1‑68 (MQ=255)
                 atCTTCGGTACGGATCTGGACCTCAACCGGCACGCCGTGCGGGCCAATCATCGAGGTGTGCAAAGACt                               >  1:1608496/1‑68 (MQ=255)
                 atCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTg                              <  1:1590431/69‑1 (MQ=255)
                 atCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTg                              <  1:1010292/69‑1 (MQ=255)
                 atCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTg                              <  1:536075/69‑1 (MQ=255)
                 atCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTg                              <  1:529272/69‑1 (MQ=255)
                 atCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTg                              <  1:1492805/69‑1 (MQ=255)
                 atCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTg                              <  1:1742815/69‑1 (MQ=255)
                 atCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTg                              <  1:2694475/69‑1 (MQ=255)
                 atCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTg                              <  1:2627879/69‑1 (MQ=255)
                 atCTTCGGTACGGATCTGGAACTCAACCGGCACGCCGTGCGGGCCAATCATCGAGGTGTGCAAAGACt                               >  1:2909962/1‑68 (MQ=255)
                 atCTTCGGTACGGAACTGGACCTCAACCGGCACGCCGTGCGGGCCAATCATCGAGGTGTGCAAAGACt                               >  1:2230039/1‑68 (MQ=255)
                  tCTTCGGTACGGATCTGGACCTCAACCGGCCCACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACt                               >  1:1113423/1‑67 (MQ=255)
                  tCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACt                               >  1:3381079/1‑67 (MQ=255)
                  tCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACt                               >  1:3219826/1‑67 (MQ=255)
                    ttCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTGATa                           >  1:2664769/1‑69 (MQ=255)
                    ttCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTGATa                           >  1:2460407/1‑69 (MQ=255)
                                     ccTCAACCGGCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTgg            <  1:1065618/67‑1 (MQ=255)
                                             ggCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTGGAATGGCGata  <  1:1968130/69‑1 (MQ=255)
                                             ggCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTGGAATGGCGata  <  1:428163/69‑1 (MQ=255)
                                             ggCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTGGAATGGCGata  <  1:108964/69‑1 (MQ=255)
                                             ggCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTGGAATGGCGata  <  1:750850/69‑1 (MQ=255)
                                             ggCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTGGAATGGCGata  <  1:2607002/69‑1 (MQ=255)
                                             ggCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTGGAATGGCGata  <  1:1256765/69‑1 (MQ=255)
                                             ggCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTGGAATGGCGata  <  1:3147930/69‑1 (MQ=255)
                                             ggCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTGGAATGGCGata  <  1:2574119/69‑1 (MQ=255)
                                              gCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTGGAATGGCGata  <  1:2520769/68‑1 (MQ=255)
                                                       gCGGGCCAATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTgg            <  1:113693/49‑1 (MQ=255)
                                                       gCGGGCCAATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTgg            <  1:888910/49‑1 (MQ=255)
                                                       gCGGGCCAATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTgg            <  1:2244140/49‑1 (MQ=255)
                                                       gCGGGCCAATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTgg            <  1:939550/49‑1 (MQ=255)
                                                       gCGGGCCAATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTgg            <  1:641064/49‑1 (MQ=255)
                                                       gCGGGCCAATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTgg            <  1:2694114/49‑1 (MQ=255)
                                                       gCGGGCCAATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTgg            <  1:2434844/49‑1 (MQ=255)
                                                       gCGGGCCAATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTgg            <  1:428945/49‑1 (MQ=255)
                                                       gCGGGCCAATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTgg            <  1:3329744/49‑1 (MQ=255)
                                                       gCGGGCCAATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTgg            <  1:2606735/49‑1 (MQ=255)
                                                       gCGGGCCAATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTgg            <  1:173843/49‑1 (MQ=255)
                                                       gCGGGCCAATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTgg            <  1:189584/49‑1 (MQ=255)
                                                       gCGGGCCAATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTgg            <  1:2813194/49‑1 (MQ=255)
                                                       gCGGGCCAATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTgg            <  1:1924933/49‑1 (MQ=255)
                                                       gCGGGCCAATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTgg            <  1:2043273/49‑1 (MQ=255)
                                                       gCGGGCCAATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTgg            <  1:2671073/49‑1 (MQ=255)
                                                        cGGGCCAATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTgg            <  1:1904345/48‑1 (MQ=255)
                                                              |                                                   
TCCGCCATCTGGTCCATATCTTCGGTACGGATCTGGACCTCAACCGGCACACCGTGCGGGCCGATCATCGAGGTGTGCAAAGACTGATAGCCGTTCGCTTTTGGAATGGCGATA  >  minE/2518080‑2518193

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