Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,738,392 C→T 33.2% K473K (AAG→AAA lepA ← GTP binding membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,738,3920CT33.2% 55.2 / 34.9 69K473K (AAG→AAAlepAGTP binding membrane protein
Reads supporting (aligned to +/- strand):  ref base C (38/8);  new base T (3/20);  total (41/28)
Fisher's exact test for biased strand distribution p-value = 2.84e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

AGCGCATCAACACGTTCACCGTTGATTAATACGTCTACACGTACCATGTCGGACGCCTGGAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCGAAGAAGTCGAGCAC  >  minE/1738327‑1738460
                                                                 |                                                                    
aGCGCATCAACACGTTCACCGTTGATTAATACGTCTACACGTACCATGTCGGACGCCTGGAAGCGCTTg                                                                   <  1:3148109/69‑1 (MQ=255)
aGCGCATCAACACGTTCACCGTTGATTAATACGTCTACACGTACCATGTCGGACGCCTGGAAGCGCTTg                                                                   <  1:3409814/69‑1 (MQ=255)
aGCGCATCAACACGTTCACCGTTGATTAATACGTCTACACGTACCATGTCGGACGCCTGGAAGCGCTTg                                                                   <  1:746266/69‑1 (MQ=255)
aGCGCATCAACACGTTCACCGTTGATTAATACGTCTACACGTACCATGTCGGACGCCTGGAAGCGCTTg                                                                   <  1:3018140/69‑1 (MQ=255)
    cATCAACACGTTCACCGTTGATTAATACGTCTACACGTACCATGTCGGACGCCTGGAAGCGCTTGAAGt                                                               >  1:2319424/1‑69 (MQ=255)
    cATCAACACGTTCACCGTTGATTAATACGTCTACACGTACCATGTCGGACGCCTGGAAGCGCTTGAAGt                                                               >  1:2374842/1‑69 (MQ=255)
    cATCAACACGTTCACCGTTGATTAATACGTCTACACGTACCATGTCGGACGCCTGGAAGCGCTTGAAGt                                                               >  1:3143414/1‑69 (MQ=255)
    cATCAACACGTTCACCGTTGATTAATACGTCTACACGTACCATGTCGGACGCCTGGAAGCGCTTGAAGt                                                               >  1:3229680/1‑69 (MQ=255)
    cATCAACACGTTCACCGTTGATTAATACGTCTACACGTACCATGTCGGACGCCTGGAAGCGCTTGAAGt                                                               >  1:2803100/1‑69 (MQ=255)
    cATCAACACGTTCACCGTTGATTAATACGTCTACACGTACCATGTCGGACGCCTGGAAGCGCTTGAAGt                                                               >  1:600140/1‑69 (MQ=255)
    cATCAACACGTTCACCGTTGATTAATACGTCTACACGTACCATGTCGGACGCCTGGAAGCGCTAGAAGt                                                               >  1:3120888/1‑69 (MQ=255)
                cACCGTTGATTAATACGTCTACACGTACCATGTCGGACGCCTGGAAGCGCTTGGAGTTGTAATCCaga                                                    >  1:2916932/1‑68 (MQ=255)
                cACCGTTGATTAATACGTCTACACGTACCATGTCGGACGCCTGGAAGCGCTTGAAGTTGTAATCCagag                                                   >  1:1813169/1‑69 (MQ=255)
                cACCGTTGATTAATACGTCTACACGTACCATGTCGGACGCCTGGAAGCGCTTGAAGTTGTAATCCagag                                                   >  1:2199581/1‑69 (MQ=255)
                 aCCGTTGATTAATACGTCTACACGTACCATGTCGGACGCCTGGAAGCGTTTGAAGTTGTAATCCAgaga                                                  >  1:2509221/1‑69 (MQ=255)
                 aCCGTTGATTAATACGTCTACACGTACCATGTCGGACGCCTGGAAGCGTTTGAAGTTGTAATCCAgaga                                                  >  1:3255747/1‑69 (MQ=255)
                 aCCGTTGATTAATACGTCTACACGTACCACGTCGGACGCCTGGAAGCGTTTGAAGTTGTAATCCAgaga                                                  >  1:1228804/1‑69 (MQ=255)
                               cGTCTACACGTACCATGTCGGCCGCCTGGAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCAcgcg                                    >  1:3402960/1‑69 (MQ=255)
                                     cacGTACCATGTCGGACGCCTGGAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAg                              >  1:2007342/1‑69 (MQ=255)
                                     cacGTACCATGTCGGACGCCTGGAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAg                              >  1:1675402/1‑69 (MQ=255)
                                     cacGTACCATGTCGGACGCCTGGAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAg                              >  1:3073136/1‑69 (MQ=255)
                                          aCCATGTCGGACGCCTGGAATCGTTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTc                         <  1:2629210/69‑1 (MQ=255)
                                          aCCATGTCGGACGCCTGGAATCGTTTGAAGTTGTAATCCAGAGACGCAAAACCACGCGAGGTAGATTTc                         <  1:585329/69‑1 (MQ=255)
                                          aCCATGTCGGACGCCTGGAAGCGTTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTc                         <  1:2063199/69‑1 (MQ=255)
                                          aCCATGTCGGACGCCTGGAAGCGTTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTc                         <  1:2516946/69‑1 (MQ=255)
                                          aCCATGTCGGACGCCTGGAAGCGTTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTc                         <  1:2776598/69‑1 (MQ=255)
                                          aCCATGTCGGACGCCTGGAAGCGTTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTc                         <  1:2919761/69‑1 (MQ=255)
                                          aCCATGTCGGACGCCTGGAAGCGTTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTc                         <  1:2944213/69‑1 (MQ=255)
                                          aCCATGTCGGACGCCTGGAAGCGTTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTc                         <  1:154736/69‑1 (MQ=255)
                                          aCCATGTCGGACGCCTGGAAGCGTTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTc                         <  1:1055105/69‑1 (MQ=255)
                                          aCCATGTCGGACGCCTGGAAGCGTTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTc                         <  1:881738/69‑1 (MQ=255)
                                          aCCATGTCGGACGCCTGGAAGCGTTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTc                         <  1:830954/69‑1 (MQ=255)
                                          aCCATGTCGGACGCCTGGAAGCGTTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTc                         <  1:1333298/69‑1 (MQ=255)
                                          aCCATGTCGGACGCCTGGAAGCGTTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTc                         <  1:800160/69‑1 (MQ=255)
                                          aCCATGTCGGACGCCTGGAAGCGTTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTc                         <  1:731941/69‑1 (MQ=255)
                                          aCCATGTCGGACGCCTGGAAGCGTTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTc                         <  1:1389420/69‑1 (MQ=255)
                                          aCCATGTCGGACGCCTGGAAGCGTTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTc                         <  1:1658920/69‑1 (MQ=255)
                                          aCCATGTCGGACGCCTGGAAGCGTTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTc                         <  1:529117/69‑1 (MQ=255)
                                          aCCATGTCGGACGCCTGGAAGCGTTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTc                         <  1:526102/69‑1 (MQ=255)
                                          aCCATGTCGGACGCCTGGAAGCGTTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTc                         <  1:1784528/69‑1 (MQ=255)
                                          aCCATGTCGGACGCCTGGAAGCGTTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTc                         <  1:178648/69‑1 (MQ=255)
                                                tCGGACGCCTGGAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGa                   >  1:180993/1‑69 (MQ=255)
                                                tCGGACGCCTGGAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGa                   >  1:3287570/1‑69 (MQ=255)
                                                    aCGCCTGGAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGAt                  >  1:366241/1‑66 (MQ=255)
                                                    aCGCCTGGAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCg                <  1:3197176/68‑1 (MQ=255)
                                                    aCGCCTGGAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCg                <  1:641384/68‑1 (MQ=255)
                                                    aCGCCTGGAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCg                <  1:1421934/68‑1 (MQ=255)
                                                    aCGCCTGGAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCg                <  1:2442223/68‑1 (MQ=255)
                                                     cGCCTGGAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCga               >  1:3212981/1‑68 (MQ=255)
                                                     cGCCTGGAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCga               >  1:928303/1‑68 (MQ=255)
                                                     cGCCTGGAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCga               >  1:3173825/1‑68 (MQ=255)
                                                     cGCCTGGAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCga               >  1:1009286/1‑68 (MQ=255)
                                                     cGCCTGGAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCga               >  1:1086184/1‑68 (MQ=255)
                                                     cGCCTGGAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCga               >  1:1647337/1‑68 (MQ=255)
                                                     cGCCTGGAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCga               >  1:1729862/1‑68 (MQ=255)
                                                     cGCCTGGAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCga               >  1:1497403/1‑68 (MQ=255)
                                                     cGCCTGGAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCga               >  1:1959828/1‑68 (MQ=255)
                                                     cGCCTGGAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCga               >  1:2134673/1‑68 (MQ=255)
                                                     cGCCTGGAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCga               >  1:2488569/1‑68 (MQ=255)
                                                     cGCCTGGAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCga               >  1:1365279/1‑68 (MQ=255)
                                                     cGCCTGGAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCga               >  1:1092595/1‑68 (MQ=255)
                                                     cGCCTGGAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCg                >  1:2194455/1‑67 (MQ=255)
                                                           gAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCGaagaag          >  1:789336/1‑67 (MQ=255)
                                                           gAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCGaagaag          >  1:2165633/1‑67 (MQ=255)
                                                           gAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCGaagaag          >  1:2336272/1‑67 (MQ=255)
                                                           gAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCGaagaag          >  1:1372589/1‑67 (MQ=255)
                                                           gAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCGaagaag          >  1:448362/1‑67 (MQ=255)
                                                           gAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCGaagaag          >  1:1637568/1‑67 (MQ=255)
                                                             aGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCGAAGAAGTCg       >  1:464608/1‑68 (MQ=255)
                                                                 cTTGATGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCGAAGAAGTCGAGcac  <  1:138501/69‑1 (MQ=255)
                                                                 cTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCGAAGAAGTCGAGcac  <  1:1780116/69‑1 (MQ=255)
                                                                 cTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCGAAGAAGTCGAGcac  <  1:2630512/69‑1 (MQ=255)
                                                                 cTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCGAAGAAGTCGAGcac  <  1:3435604/69‑1 (MQ=255)
                                                                 cTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCGAAGAAGTCGAGcac  <  1:3345830/69‑1 (MQ=255)
                                                                 cTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCGAAGAAGTCGAGcac  <  1:2513860/69‑1 (MQ=255)
                                                                 cTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCGAAGAAGTCGAGcac  <  1:3207575/69‑1 (MQ=255)
                                                                 cTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCGAAGAAGTCGAGcac  <  1:1088795/69‑1 (MQ=255)
                                                                 cTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCGAAGAAGTCGAGcac  <  1:926070/69‑1 (MQ=255)
                                                                 cTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCGAAGAAGTCGAGcac  <  1:2536204/69‑1 (MQ=255)
                                                                 |                                                                    
AGCGCATCAACACGTTCACCGTTGATTAATACGTCTACACGTACCATGTCGGACGCCTGGAAGCGCTTGAAGTTGTAATCCAGAGACGCATAACCACGCGAGGTAGATTTCAGGCGATCGAAGAAGTCGAGCAC  >  minE/1738327‑1738460

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: