Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,151,174 G→A 48.6% F285F (TTC→TTT ftsH ← protease, ATP‑dependent zinc‑metallo

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,151,1740GA48.6% 2.9 / 66.1 76F285F (TTC→TTTftsHprotease, ATP‑dependent zinc‑metallo
Reads supporting (aligned to +/- strand):  ref base G (34/5);  new base A (21/16);  total (55/21)
Fisher's exact test for biased strand distribution p-value = 4.40e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.00e-01

AGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAG  >  minE/2151109‑2151239
                                                                 |                                                                 
agGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTACCTTCGaa                                                                 >  1:88011/1‑68 (MQ=255)
agGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTACCTTCGaa                                                                 >  1:219704/1‑68 (MQ=255)
agGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTACCTTCGaa                                                                 >  1:2841183/1‑68 (MQ=255)
agGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTACCTTCGaa                                                                 >  1:1945970/1‑68 (MQ=255)
         tCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTACCTTCGAAGccatccat                                                        >  1:3138939/1‑68 (MQ=255)
               aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                 >  1:1052954/1‑69 (MQ=255)
               aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                 >  1:2717140/1‑69 (MQ=255)
               aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                 >  1:2381855/1‑69 (MQ=255)
               aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                 >  1:2735019/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:3188291/1‑69 (MQ=255)
               aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                 >  1:3237479/1‑69 (MQ=255)
               aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                 >  1:3320813/1‑69 (MQ=255)
               aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                 >  1:1769921/1‑69 (MQ=255)
               aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                 >  1:1164299/1‑69 (MQ=255)
               aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                 >  1:1751417/1‑69 (MQ=255)
               aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                 >  1:1750831/1‑69 (MQ=255)
               aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                 >  1:1747011/1‑69 (MQ=255)
               aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                 >  1:855419/1‑69 (MQ=255)
               aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                 >  1:1502162/1‑69 (MQ=255)
               aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                 >  1:9883/1‑69 (MQ=255)
               aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                 >  1:1835606/1‑69 (MQ=255)
               aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                                                 >  1:787225/1‑69 (MQ=255)
               aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTAACTTCAAAACCATCCATCTCAAcc                                                 >  1:723580/1‑69 (MQ=255)
                                cgcgGCGATAACGATGATACCTTCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGt                                >  1:2283953/1‑69 (MQ=255)
                                cgcgGCGATAACGATGATACCTTCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGt                                >  1:1132630/1‑69 (MQ=255)
                                   agcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGTGTCTg                             <  1:514356/68‑1 (MQ=255)
                                   agcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                             <  1:106058/68‑1 (MQ=255)
                                   agcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                             <  1:2666839/68‑1 (MQ=255)
                                   agcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                             <  1:291634/68‑1 (MQ=255)
                                   agcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                             <  1:2996348/68‑1 (MQ=255)
                                   agcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                             <  1:3319953/68‑1 (MQ=255)
                                   agcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                             <  1:3391647/68‑1 (MQ=255)
                                   agcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                             <  1:2580542/68‑1 (MQ=255)
                                   agcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                             <  1:1225692/68‑1 (MQ=255)
                                   agcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                             <  1:128805/68‑1 (MQ=255)
                                   agcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                             <  1:1904903/68‑1 (MQ=255)
                                   agcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                             <  1:1317227/68‑1 (MQ=255)
                                    gcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                             <  1:286619/68‑1 (MQ=255)
                                    gcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                             <  1:1977346/68‑1 (MQ=255)
                                    gcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                             <  1:736785/68‑1 (MQ=255)
                                    gcgATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAACCAGCATCTGGTTCAGAGTCTg                             <  1:2414224/68‑1 (MQ=255)
                                         aaCGATGATACCTTCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACg                       >  1:742770/1‑69 (MQ=255)
                                         aaCGATGATACCTTCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACg                       >  1:2747977/1‑69 (MQ=255)
                                         aaCGATGATACCTTCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACg                       >  1:2078749/1‑69 (MQ=255)
                                         aaCGATGATACCTTCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACg                       >  1:1569426/1‑69 (MQ=255)
                                         aaCGATGATACCTTCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACg                       >  1:1865718/1‑69 (MQ=255)
                                         aaCGATGATACCTTCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACg                       >  1:780464/1‑69 (MQ=255)
                                         aaCGATGATACCTTCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACg                       >  1:3257493/1‑69 (MQ=255)
                                                     ttCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCa                                    <  1:1121647/44‑1 (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)
                                                     ttCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGacc           >  1:3229992/1‑69 (MQ=255)
                                                     ttCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGacc           >  1:3290172/1‑69 (MQ=255)
                                                     ttCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGacc           >  1:1870606/1‑69 (MQ=255)
                                                     ttCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGacc           >  1:3040937/1‑69 (MQ=255)
                                                     ttCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGacc           >  1:1841751/1‑69 (MQ=255)
                                                     ttCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGacc           >  1:3431740/1‑69 (MQ=255)
                                                     ttCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGacc           >  1:1267628/1‑69 (MQ=255)
                                                     ttCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGacc           >  1:290586/1‑69 (MQ=255)
                                                     ttCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGac            >  1:2684997/1‑68 (MQ=255)
                                                     ttCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGac            >  1:799119/1‑68 (MQ=255)
                                                     ttCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGac            >  1:2846306/1‑68 (MQ=255)
                                                      tCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGacca          >  1:2907735/1‑69 (MQ=255)
                                                      tCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGacca          >  1:3424118/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:2134256/1‑67 (MQ=255)
                                                                cGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAg  >  1:2473318/1‑67 (MQ=255)
                                                                cGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAg  >  1:922113/1‑67 (MQ=255)
                                                                cGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAg  >  1:2641368/1‑67 (MQ=255)
                                                                 |                                                                 
AGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAG  >  minE/2151109‑2151239

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: