Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,151,177 G→A 51.4% G284G (GGC→GGT ftsH ← protease, ATP‑dependent zinc‑metallo

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,151,1770GA51.4% ‑2.9 / 64.8 72G284G (GGC→GGTftsHprotease, ATP‑dependent zinc‑metallo
Reads supporting (aligned to +/- strand):  ref base G (30/5);  new base A (21/16);  total (51/21)
Fisher's exact test for biased strand distribution p-value = 9.38e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

TCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGC  >  minE/2151118‑2151245
                                                           |                                                                    
tCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTACCTTCGAAGccatccat                                                              >  1:3138939/1‑68 (MQ=255)
      aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                       >  1:2717140/1‑69 (MQ=255)
      aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                       >  1:9883/1‑69 (MQ=255)
      aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                       >  1:855419/1‑69 (MQ=255)
      aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                       >  1:787225/1‑69 (MQ=255)
      aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                       >  1:3320813/1‑69 (MQ=255)
      aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                       >  1:3237479/1‑69 (MQ=255)
      aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                       >  1:3188291/1‑69 (MQ=255)
      aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                       >  1:2735019/1‑69 (MQ=255)
      aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                       >  1:1052954/1‑69 (MQ=255)
      aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                       >  1:2381855/1‑69 (MQ=255)
      aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                       >  1:2246261/1‑69 (MQ=255)
      aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                       >  1:2182857/1‑69 (MQ=255)
      aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                       >  1:1966813/1‑69 (MQ=255)
      aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                       >  1:1835606/1‑69 (MQ=255)
      aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                       >  1:1751417/1‑69 (MQ=255)
      aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                       >  1:1769921/1‑69 (MQ=255)
      aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                       >  1:1750831/1‑69 (MQ=255)
      aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                       >  1:1164299/1‑69 (MQ=255)
      aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                       >  1:1502162/1‑69 (MQ=255)
      aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                       >  1:1747011/1‑69 (MQ=255)
      aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTAACTTCAAAACCATCCATCTCAAcc                                                       >  1:723580/1‑69 (MQ=255)
                       cgcgGCGATAACGATGATACCTTCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGt                                      >  1:1132630/1‑69 (MQ=255)
                       cgcgGCGATAACGATGATACCTTCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGt                                      >  1:2283953/1‑69 (MQ=255)
                          agcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGTGTCTg                                   <  1:514356/68‑1 (MQ=255)
                          agcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                                   <  1:2580542/68‑1 (MQ=255)
                          agcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                                   <  1:106058/68‑1 (MQ=255)
                          agcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                                   <  1:2666839/68‑1 (MQ=255)
                          agcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                                   <  1:1317227/68‑1 (MQ=255)
                          agcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                                   <  1:291634/68‑1 (MQ=255)
                          agcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                                   <  1:2996348/68‑1 (MQ=255)
                          agcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                                   <  1:128805/68‑1 (MQ=255)
                          agcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                                   <  1:1225692/68‑1 (MQ=255)
                          agcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                                   <  1:3319953/68‑1 (MQ=255)
                          agcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                                   <  1:3391647/68‑1 (MQ=255)
                          agcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                                   <  1:1904903/68‑1 (MQ=255)
                           gcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                                   <  1:2414224/68‑1 (MQ=255)
                           gcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                                   <  1:736785/68‑1 (MQ=255)
                           gcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                                   <  1:286619/68‑1 (MQ=255)
                           gcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                                   <  1:1977346/68‑1 (MQ=255)
                                aaCGATGATACCTTCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACg                             >  1:2078749/1‑69 (MQ=255)
                                aaCGATGATACCTTCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACg                             >  1:1865718/1‑69 (MQ=255)
                                aaCGATGATACCTTCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACg                             >  1:780464/1‑69 (MQ=255)
                                aaCGATGATACCTTCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACg                             >  1:742770/1‑69 (MQ=255)
                                aaCGATGATACCTTCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACg                             >  1:2747977/1‑69 (MQ=255)
                                aaCGATGATACCTTCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACg                             >  1:3257493/1‑69 (MQ=255)
                                aaCGATGATACCTTCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACg                             >  1:1569426/1‑69 (MQ=255)
                                            ttCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCa                                          <  1:647831/44‑1 (MQ=255)
                                            ttCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCa                                          <  1:399415/44‑1 (MQ=255)
                                            ttCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCa                                          <  1:2924486/44‑1 (MQ=255)
                                            ttCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCa                                          <  1:1879273/44‑1 (MQ=255)
                                            ttCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCa                                          <  1:1121647/44‑1 (MQ=255)
                                            ttCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGacc                 >  1:3431740/1‑69 (MQ=255)
                                            ttCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGacc                 >  1:1841751/1‑69 (MQ=255)
                                            ttCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGacc                 >  1:3290172/1‑69 (MQ=255)
                                            ttCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGacc                 >  1:1870606/1‑69 (MQ=255)
                                            ttCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGacc                 >  1:3229992/1‑69 (MQ=255)
                                            ttCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGacc                 >  1:1267628/1‑69 (MQ=255)
                                            ttCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGacc                 >  1:3040937/1‑69 (MQ=255)
                                            ttCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGacc                 >  1:290586/1‑69 (MQ=255)
                                            ttCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGac                  >  1:799119/1‑68 (MQ=255)
                                            ttCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGac                  >  1:2846306/1‑68 (MQ=255)
                                            ttCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGac                  >  1:2684997/1‑68 (MQ=255)
                                             tCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGacca                >  1:3424118/1‑69 (MQ=255)
                                             tCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGacca                >  1:2907735/1‑69 (MQ=255)
                                             tCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGacc                 >  1:2078202/1‑68 (MQ=255)
                                                       cGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGc            >  1:671812/1‑63 (MQ=255)
                                                       cGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCa         >  1:2263385/1‑66 (MQ=255)
                                                       cGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAg        >  1:2473318/1‑67 (MQ=255)
                                                       cGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAg        >  1:922113/1‑67 (MQ=255)
                                                       cGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAg        >  1:2641368/1‑67 (MQ=255)
                                                       cGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAg        >  1:2134256/1‑67 (MQ=255)
                                                           gCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGccc          >  1:2356245/1‑61 (MQ=255)
                                                           gCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAgc  >  1:1308837/1‑69 (MQ=255)
                                                           gCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAgc  >  1:294652/1‑69 (MQ=255)
                                                           gCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAgc  >  1:903311/1‑69 (MQ=255)
                                                           gCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAgc  >  1:3195190/1‑69 (MQ=255)
                                                           gCCATCAATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCa         >  1:3038741/1‑62 (MQ=255)
                                                           |                                                                    
TCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGC  >  minE/2151118‑2151245

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: