Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,447,108 G→A 50.7% A91A (GCC→GCT rcsC ← hybrid sensory kinase in two‑component regulatory system with RcsB and YojN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,447,1080GA50.7% ‑5.6 / 62.8 69A91A (GCC→GCTrcsChybrid sensory kinase in two‑component regulatory system with RcsB and YojN
Reads supporting (aligned to +/- strand):  ref base G (24/10);  new base A (16/19);  total (40/29)
Fisher's exact test for biased strand distribution p-value = 5.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAG  >  minE/1447046‑1447167
                                                              |                                                           
aaCCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAAc                                                       >  1:1353289/1‑69 (MQ=255)
                    agaACCTCGCCAGGTGTTACTCATTGCGTAACAATCAGAGTCAGCa                                                          >  1:1435195/1‑46 (MQ=255)
                    agaACCTCGCCAGGTGTTACTCATTGCGGAACAATCAGAGTCAGCa                                                          >  1:914588/1‑46 (MQ=255)
                    agaACCTCGCCAGGTGTTACTCATTGCGGAACAATCAGAGTCAGCa                                                          >  1:2353234/1‑46 (MQ=255)
                    agaACCTCGCCAGGTGTTACTCATTGCGGAACAATCAGAGTCAGCa                                                          >  1:2462702/1‑46 (MQ=255)
                    agaACCTCGCCAGGTGTTACTCATTGCGGAACAATCAGAGTCAGCa                                                          >  1:2552663/1‑46 (MQ=255)
                    agaACCTCGCCAGGTGTTACTCATTGCGGAACAATCAGAGTCAGCa                                                          >  1:2648383/1‑46 (MQ=255)
                    agaACCTCGCCAGGTGTTACTCATTGCGGAACAATCAGAGTCAGCa                                                          >  1:3000207/1‑46 (MQ=255)
                    agaACCTCGCCAGGTGTTACTCATTGCGGAACAATCAGAGTCAGCa                                                          >  1:3099771/1‑46 (MQ=255)
                    agaACCTCGCCAGGTGTTACTCATTGCGGAACAATCAGAGTCAGCa                                                          >  1:3139961/1‑46 (MQ=255)
                    agaACCTCGCCAGGTGTTACTCATTGCGGAACAATCAGAGTCAGCa                                                          >  1:3389063/1‑46 (MQ=255)
                    agaACCTCGCCAGGTGTTACTCATTGCGGAACAATCAGAGTCAGCa                                                          >  1:1102557/1‑46 (MQ=255)
                    agaACCTCGCCAGGTGTTACTCATTGCGGAACAATCAGAGTCAGCa                                                          >  1:476507/1‑46 (MQ=255)
                    agaACCTCGCCAGGTGTTACTCATTGCGGAACAATCAGAGTCAGCa                                                          >  1:1276390/1‑46 (MQ=255)
                    agaACCTCGCCAGGTGTTACTCATTGCGGAACAATCAGAGTCAGCa                                                          >  1:1104498/1‑46 (MQ=255)
                    agaACCTCGCCAGGTGTTACTCATTGCGGAACAATCAGAGTCAGCa                                                          >  1:591522/1‑46 (MQ=255)
                    agaACCTCGCCAGGTGTTACTCATTGCGGAACAATCAGAGTCAGCa                                                          >  1:1338076/1‑46 (MQ=255)
                     gaACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGcaca                                  <  1:1135378/69‑1 (MQ=255)
                                 gtgtTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGtt                      >  1:2286443/1‑69 (MQ=255)
                                    ttACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAgg                                      >  1:3060690/1‑50 (MQ=255)
                                    ttACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAgg                                      >  1:2049580/1‑50 (MQ=255)
                                    ttACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAgg                                      >  1:2476187/1‑50 (MQ=255)
                                    ttACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAgg                                      >  1:2670086/1‑50 (MQ=255)
                                    ttACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAgg                                      >  1:3109511/1‑50 (MQ=255)
                                    ttACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAgg                                      >  1:372388/1‑50 (MQ=255)
                                    ttACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTtct                   >  1:1640377/1‑69 (MQ=255)
                                      aCTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTtctc                  >  1:1304028/1‑68 (MQ=255)
                                      aCTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCg                 >  1:2969626/1‑69 (MQ=255)
                                      aCTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCg                 >  1:105578/1‑69 (MQ=255)
                                      aCTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCg                 >  1:2558580/1‑69 (MQ=255)
                                      aCTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCg                 >  1:2662038/1‑69 (MQ=255)
                                      aCTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCg                 >  1:3377700/1‑69 (MQ=255)
                                      aCTCATTGCGGAACAATCAGAGTCAGCAAACAGCGGTTCAAACGCAGGCACATCCGTCTGCGTTtctc                  <  1:2571120/68‑1 (MQ=255)
                                      aCTCATTGCGGAACAATCAGAGTCAGCAAACAGCGGTTCAAACGCAGGCACATCCGTCTGCGTTtctc                  <  1:1043567/68‑1 (MQ=255)
                                      aCTCATTGCGGAACAATCAGAGTCAGCAAACAGCGGTTCAAACGCAGGCACATCCGTCTGCGTTtctc                  <  1:2209797/68‑1 (MQ=255)
                                      aCTCATTGCGGAACAATCAGAGTCAGCAAACAGCGGTTCAAACGCAGGCACATCCGTCTGCGTTtctc                  <  1:1287694/68‑1 (MQ=255)
                                      aCTCATTGCGGAACAATCAGAGTCAGCAAACAGCGGTTCAAACGCAGGCACATCCGTCTGCGTTtctc                  <  1:1495072/68‑1 (MQ=255)
                                      aCTCATTGCGGAACAATCAGAGTCAGCAAACAGCGGTTCAAACGCAGGCACATCCGTCTGCGTTtctc                  <  1:1730008/68‑1 (MQ=255)
                                      aCTCATTGCGGAACAATCAGAGTCAGCAAACAGCGGTTCAAACGCAGGCACATCCGTCTGCGTTtctc                  <  1:1778702/68‑1 (MQ=255)
                                      aCTCATTGCGGAACAATCAGAGTCAGCAAACAGCGGTTCAAACGCAGGCACATCCGTCTGCGTTtctc                  <  1:1937929/68‑1 (MQ=255)
                                      aCTCATTGCGGAACAATCAGAGTCAGCAAACAGCGGTTCAAACGCAGGCACATCCGTCTGCGTTtctc                  <  1:2007526/68‑1 (MQ=255)
                                      aCTCATTGCGGAACAATCAGAGTCAGCAAACAGCGGTTCAAACGCAGGCACATCCGTCTGCGTTtctc                  <  1:2087898/68‑1 (MQ=255)
                                      aCTCATTGCGGAACAATCAGAGACAGCAAACAGCGGTTCAAACGCAGGCACATCCGTCTGCGTTtctc                  <  1:2966387/68‑1 (MQ=255)
                                        tCATTGCGGAACAATCAGAGTCAGCAAACAGCGGTTCAAACGCAGGAACATCCGTCTGCGTTtctc                  <  1:1394483/66‑1 (MQ=255)
                                             gCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCAcgc          >  1:936754/1‑69 (MQ=255)
                                             gCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCAcgc          >  1:3422171/1‑69 (MQ=255)
                                             gCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCAcgc          >  1:1714166/1‑69 (MQ=255)
                                             gCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCAcgc          >  1:2865695/1‑69 (MQ=255)
                                             gCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCAcgc          >  1:510145/1‑69 (MQ=255)
                                                gAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCAc            >  1:2254146/1‑64 (MQ=255)
                                                gAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCggg       <  1:172697/69‑1 (MQ=255)
                                                gAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCggg       <  1:3259934/69‑1 (MQ=255)
                                                gAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCgg        >  1:1821227/1‑68 (MQ=255)
                                                gAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCgg        >  1:424191/1‑68 (MQ=255)
                                                gAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGAggg       <  1:1028301/69‑1 (MQ=255)
                                                   cAATCGGAGTCGGCAAACAGCGGTTCAAACGCAgg                                      >  1:54277/1‑35 (MQ=255)
                                                     aTCGGAGTCGTCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAg  <  1:1894969/69‑1 (MQ=255)
                                                     aTCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAg  <  1:1407785/69‑1 (MQ=255)
                                                     aTCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAg  <  1:1390964/69‑1 (MQ=255)
                                                     aTCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAg  <  1:674552/69‑1 (MQ=255)
                                                     aTCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAg  <  1:373632/69‑1 (MQ=255)
                                                     aTCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAg  <  1:876906/69‑1 (MQ=255)
                                                     aTCAGAGTCAGCAAACAGCTGTTCAAACGCAGGCACATCCGTCTGCGTTTCTCGTCCACGCGGGGAAAg  <  1:2648891/69‑1 (MQ=255)
                                                     aTCAGAGTCAGCAAACAGCGGTTCAAACGCAGGCACATCCGTCTGCGTTTCTCGTCCACGCGGGGAAAg  <  1:1114141/69‑1 (MQ=255)
                                                     aTCAGAGTCAGCAAACAGCGGTTCAAACGCAGGCACATCCGTCTGCGTTTCTCGTCCACGCGGGGAAAg  <  1:3183933/69‑1 (MQ=255)
                                                     aTCAGAGTCAGCAAACAGCGGTTCAAACGCAGGCACATCCGTCTGCGTTTCTCGTCCACGCGGGGAAAg  <  1:1389268/69‑1 (MQ=255)
                                                     aTCAGAGTCAGCAAACAGCGGTTCAAACGCAGGCACATCCGTCTGCGTTTCTCGTCCACGCGGGGAAAg  <  1:2533046/69‑1 (MQ=255)
                                                     aTCAGAGTCAGCAAACAGCGGTTCAAACGCAGGCACATCCGTCTGCGTTTCTCGTCCACGCGGGGAAAg  <  1:1681163/69‑1 (MQ=255)
                                                     aTCAGAGTCAGCAAACAGCGGTTCAAACGCAGGCACATCCGTCTGCGTTTCTCGTCCACGCGGGGAAAg  <  1:2637172/69‑1 (MQ=255)
                                                              ggCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAg  <  1:2874284/60‑1 (MQ=255)
                                                              |                                                           
AACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAG  >  minE/1447046‑1447167

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: