Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,204,203 A→G 30.7% T433T (ACT→ACC yoaE ← fused predicted membrane proteins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,204,2030AG30.7% 31.4 / 15.8 39T433T (ACT→ACCyoaEfused predicted membrane proteins
Reads supporting (aligned to +/- strand):  ref base A (9/18);  new base G (6/6);  total (15/24)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTAATGATTTCCGGCGTTTCGTCAGCGTCCGGGAATTCACCCGCAATGGCTTCCAGCACA  >  minE/1204141‑1204263
                                                              |                                                            
tCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTAATg                                                        <  1:3102566/69‑1 (MQ=255)
tCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTAATg                                                        <  1:2454376/69‑1 (MQ=255)
tCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTAATg                                                        <  1:2778373/69‑1 (MQ=255)
tCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTAATg                                                        <  1:2865806/69‑1 (MQ=255)
tCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTAATg                                                        <  1:3014538/69‑1 (MQ=255)
tCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTAATg                                                        <  1:1908506/69‑1 (MQ=255)
tCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTAATg                                                        <  1:1799039/69‑1 (MQ=255)
tCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTAATg                                                        <  1:117858/69‑1 (MQ=255)
                       cAAATCTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTAATGATTTCCGGCGTTTCGTCAGCGTc                                 >  1:1245710/1‑69 (MQ=255)
                       cAAATCTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTAATGATTTCCGGCGTTTCGTCAGCGTc                                 >  1:1645457/1‑69 (MQ=255)
                       cAAATCTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTAATGATTTCCGGCGTTTCGTCAGCGTc                                 >  1:3226005/1‑69 (MQ=255)
                       cAAATCTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTAATGATTTCCGGCGTTTCGTCAGCGTc                                 >  1:1688342/1‑69 (MQ=255)
                           tCTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTAATGATTTCCGGCGTTTCGTCAGCGTCCggg                             >  1:1832699/1‑69 (MQ=255)
                            cTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTAAt                                                         >  1:1237549/1‑40 (MQ=255)
                            cTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTAATGATTTCCGGCGTTTCGTCAGCGt                                  >  1:1487531/1‑63 (MQ=255)
                            cTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTAATGATTTCCGGCGTTTCGTCAGCGt                                  >  1:1719805/1‑63 (MQ=255)
                            cTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTAATGATTTCCGGCGTTTCGTCAGCGt                                  >  1:3170049/1‑63 (MQ=255)
                                   gCCTTTTACCAGCCAGCCGTCACCATCAGTAATGATTTCCGGCGTTTCGTCAGCGTCCGGGAATTCAc                      <  1:1651824/68‑1 (MQ=255)
                                   gCCTTTTACCAGCCAGCCGTCACCATCAGTAATGATTTCCGGCGTTTCGTCAGCGTCCGGGAATTCAc                      <  1:1636519/68‑1 (MQ=255)
                                   gCCTTTTACCAGCCAGCCGTCACCATCAGTAATGATTTCCGGCGTTTCGTCAGCGTCCGGGAATTCAc                      <  1:1979833/68‑1 (MQ=255)
                                   gCCTTTTACCAGCCAGCCGTCACCATCAGTAATGATTTCCGGCGTTTCGTCAGCGTCCGGGAATTCAc                      <  1:3143396/68‑1 (MQ=255)
                                   gCCTTTTACCAGCCAGCCGTCACCATCAGTAATGATTTCCGGCGTTTCGTCAGCGTCCGGGAATTCAc                      <  1:101539/68‑1 (MQ=255)
                                   gCCTTTTACCAGCAAGCCGTCACCATCAGTAATGATTTCCGGCGTTTCGTCAGCGTCCGGGAATTCAc                      <  1:2260514/68‑1 (MQ=255)
                                      ttttACCAGCCAGCCGTCACCATCAGTAATGATTTCCGGCGTTTCGTCAGCGTCCGGGAATTCACCCGc                  <  1:175917/69‑1 (MQ=255)
                                       tttACCAGCCAGCCGTCACCATCAGTAATGATTTCCGGCGTTTCGTCAGCGTCCGGGAATTCACCCGc                  <  1:3229738/68‑1 (MQ=255)
                                          aCCAGCCAGCCGTCACCGTCGGTAATGATTTCCGGCGTTTCGTCAGCGTCCGGGAATTCACCCGCaa                <  1:499194/67‑1 (MQ=255)
                                          aCCAGCCAGCCGTCACCGTCGGTAATGATTTCCGGCGTTTCGTCAGCGTCCGGGAATTCACCCGCaa                <  1:3028778/67‑1 (MQ=255)
                                          aCCAGCCAGCCGTCACCGTCGGTAATGATTTCCGGCGTTTCGTCAGCGTCCGGGAATTCACCCGCaa                <  1:1341004/67‑1 (MQ=255)
                                          aCCAGCCAGCCGTCACCGTCGGTAATGATTTCCGGCGTTTCGTCAGCGTCCGGGAATTCACCCGCaa                <  1:1289739/67‑1 (MQ=255)
                                          aCCAGCCAGCCGTCACCGTCGGTAATGATTTCCGGCGTTTCGTCAGCGTCCGGGAATTCACCCGCaa                <  1:1149597/67‑1 (MQ=255)
                                          aCCAGCCAGCCGTCACCATCAGTAATGATTTCCGGCGTTTCGTCAGCGTCCGGGAATTCACCCGc                  <  1:16243/65‑1 (MQ=255)
                                           ccagccagCCGTCACCGTCGGTAATGATTTCCGGCGTTTCGTCAGCGTCCGGGAATTCACCCGCaa                <  1:2280017/66‑1 (MQ=255)
                                                  gccGTCACCATCAGTAATGATTTCCGGCGTTTCGGCAGCGTCCGGGAATTCACCCGCAATGGCTTCCAg      <  1:2093909/69‑1 (MQ=255)
                                                       cACCGTCGGTAATGATTTCCGGCGTTTCGTCAGCGTCCGGGAATTCACCCGCAATGGCTTc         >  1:2781794/1‑61 (MQ=255)
                                                       cACCGTCGGTAATGATTTCCGGCGTTTCGTCAGCGTCCGGGAATTCACCCGCAATGGCTTCAAGcaca  >  1:2991921/1‑68 (MQ=255)
                                                       cACCGTCGGTAATGATTTCCGGCGTTTCGTCAGCGTCCGGGAATTCACCCGCAATGGCTTCAAGcaca  >  1:2092520/1‑68 (MQ=255)
                                                       cACCGTCGGTAATGATTTCCGGCGTTTCGTCAGCGTCCGGGAATTCACCCGCAATGGCTTCAAGcaca  >  1:193428/1‑68 (MQ=255)
                                                       cACCGTCGGTAATGATTTCCGGCGTTTCGTCAGCGTCCGGGAATTCACCCGCAATGGCTTCAAGcaca  >  1:1681331/1‑68 (MQ=255)
                                                       cACCGTCGGTAATGATTTCCGGCGTTTCGTCAGCGTCCGGGAATTCACCCGCAATGGCTTCAAGcaca  >  1:1062105/1‑68 (MQ=255)
                                                              |                                                            
TCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTAATGATTTCCGGCGTTTCGTCAGCGTCCGGGAATTCACCCGCAATGGCTTCCAGCACA  >  minE/1204141‑1204263

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: