Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,447,102 G→A 50.0% S93S (TCC→TCT 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,1020GA50.0% ‑3.7 / 63.1 70S93S (TCC→TCTrcsChybrid sensory kinase in two‑component regulatory system with RcsB and YojN
Reads supporting (aligned to +/- strand):  ref base G (25/10);  new base A (16/19);  total (41/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.61e-01

GCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAG  >  minE/1447039‑1447167
                                                               |                                                                 
gcgcATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTc                                                              >  1:1546252/1‑69 (MQ=255)
       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:1338076/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)
                            gaACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGcaca                                  <  1:1135378/69‑1 (MQ=255)
                                        gtgtTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGtt                      >  1:2286443/1‑69 (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:372388/1‑50 (MQ=255)
                                           ttACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAgg                                      >  1:3060690/1‑50 (MQ=255)
                                           ttACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAgg                                      >  1:3109511/1‑50 (MQ=255)
                                           ttACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTtct                   >  1:1640377/1‑69 (MQ=255)
                                             aCTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTtctc                  >  1:1304028/1‑68 (MQ=255)
                                             aCTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCg                 >  1:105578/1‑69 (MQ=255)
                                             aCTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCg                 >  1:2969626/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:2007526/68‑1 (MQ=255)
                                             aCTCATTGCGGAACAATCAGAGTCAGCAAACAGCGGTTCAAACGCAGGCACATCCGTCTGCGTTtctc                  <  1:2087898/68‑1 (MQ=255)
                                             aCTCATTGCGGAACAATCAGAGTCAGCAAACAGCGGTTCAAACGCAGGCACATCCGTCTGCGTTtctc                  <  1:1937929/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:510145/1‑69 (MQ=255)
                                                    gCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCAcgc          >  1:3422171/1‑69 (MQ=255)
                                                    gCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCAcgc          >  1:2865695/1‑69 (MQ=255)
                                                    gCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCAcgc          >  1:1714166/1‑69 (MQ=255)
                                                       gAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCAc            >  1:2254146/1‑64 (MQ=255)
                                                       gAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCggg       <  1:3259934/69‑1 (MQ=255)
                                                       gAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCggg       <  1:172697/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:876906/69‑1 (MQ=255)
                                                            aTCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAg  <  1:674552/69‑1 (MQ=255)
                                                            aTCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAg  <  1:373632/69‑1 (MQ=255)
                                                            aTCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAg  <  1:1390964/69‑1 (MQ=255)
                                                            aTCAGAGTCAGCAAACAGCTGTTCAAACGCAGGCACATCCGTCTGCGTTTCTCGTCCACGCGGGGAAAg  <  1:2648891/69‑1 (MQ=255)
                                                            aTCAGAGTCAGCAAACAGCGGTTCAAACGCAGGCACATCCGTCTGCGTTTCTCGTCCACGCGGGGAAAg  <  1:1389268/69‑1 (MQ=255)
                                                            aTCAGAGTCAGCAAACAGCGGTTCAAACGCAGGCACATCCGTCTGCGTTTCTCGTCCACGCGGGGAAAg  <  1:1114141/69‑1 (MQ=255)
                                                            aTCAGAGTCAGCAAACAGCGGTTCAAACGCAGGCACATCCGTCTGCGTTTCTCGTCCACGCGGGGAAAg  <  1:1681163/69‑1 (MQ=255)
                                                            aTCAGAGTCAGCAAACAGCGGTTCAAACGCAGGCACATCCGTCTGCGTTTCTCGTCCACGCGGGGAAAg  <  1:3183933/69‑1 (MQ=255)
                                                            aTCAGAGTCAGCAAACAGCGGTTCAAACGCAGGCACATCCGTCTGCGTTTCTCGTCCACGCGGGGAAAg  <  1:2533046/69‑1 (MQ=255)
                                                            aTCAGAGTCAGCAAACAGCGGTTCAAACGCAGGCACATCCGTCTGCGTTTCTCGTCCACGCGGGGAAAg  <  1:2637172/69‑1 (MQ=255)
                                                               |                                                                 
GCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAG  >  minE/1447039‑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: