Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 88,344 A→T 38.2% E150D (GAA→GAT ftsQ → membrane anchored protein involved in growth of wall at septum

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE88,3440AT38.2% 49.3 / 56.3 68E150D (GAA→GATftsQmembrane anchored protein involved in growth of wall at septum
Reads supporting (aligned to +/- strand):  ref base A (16/26);  new base T (18/8);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 2.38e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGG  >  minE/88300‑88411
                                            |                                                                   
cATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTcc                                              <  1:2674681/68‑1 (MQ=255)
cATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTcc                                              <  1:1926438/68‑1 (MQ=255)
cATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTcc                                              <  1:1865816/68‑1 (MQ=255)
cATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTcc                                              <  1:1593664/68‑1 (MQ=255)
cATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTcc                                              <  1:1537226/68‑1 (MQ=255)
cATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTcc                                              <  1:1523333/68‑1 (MQ=255)
cATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTcc                                              <  1:687570/68‑1 (MQ=255)
cATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTcc                                              <  1:3046713/68‑1 (MQ=255)
    tGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCaa                                            >  1:1885906/1‑66 (MQ=255)
    tGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCaa                                            >  1:2078076/1‑66 (MQ=255)
    tGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCaa                                            >  1:1885567/1‑66 (MQ=255)
    tGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCaa                                            >  1:2285531/1‑66 (MQ=255)
    tGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCaa                                            >  1:3028822/1‑66 (MQ=255)
    tGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCaa                                            >  1:1846563/1‑66 (MQ=255)
    tGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCaa                                            >  1:1861058/1‑66 (MQ=255)
    tGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCaa                                            >  1:3358549/1‑66 (MQ=255)
    tGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGc                                         >  1:164286/1‑69 (MQ=255)
          aCGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTcc                                              <  1:2310553/58‑1 (MQ=255)
            gcgGAAGGAAACACCTTCAGCGTGCCGCCAGATCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATgg                                  >  1:1236687/1‑68 (MQ=255)
            gcgGAAGGAAACACCTTCAGCGTGCCGCCAGATCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGc                                 >  1:230162/1‑69 (MQ=255)
            gcgGAAGGAAACACCTTCAGCGTGCCGCCAGATCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGc                                 >  1:217177/1‑69 (MQ=255)
            gcgGAAGGAAACACCTTCAGCGTGCCGCCAGATCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGc                                 >  1:2225801/1‑69 (MQ=255)
            gcgGAAGGAAACACCTTCAGCGTGCCGCCAGATCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGc                                 >  1:1950438/1‑69 (MQ=255)
            gcgGAAGGAAACACCTTCAGCGTGCCGCCAGATCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGc                                 >  1:549523/1‑69 (MQ=255)
            gcgGAAGGAAACACCTTCAGCGTGCCGCCAGATCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGc                                 >  1:2313267/1‑69 (MQ=255)
            gcgGAAGGAAACACCTTCAGCGTGCCGCCAGATCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGc                                 >  1:2408389/1‑69 (MQ=255)
            gcgGAAGGAAACACCTTCAGCGTGCCGCCAGATCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGc                                 >  1:159688/1‑69 (MQ=255)
            gcgGAAGGAAACACCTTCAGCGTGCCGCCAGATCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGc                                 >  1:1293200/1‑69 (MQ=255)
            gcgGAAGGAAACACCTTCAGCGTGCCGCCAGATCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGc                                 >  1:3315326/1‑69 (MQ=255)
                         ccTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAgc                     >  1:2278623/1‑68 (MQ=255)
                             cAGCGTGCCGCCAGATCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCa                 >  1:2875644/1‑68 (MQ=255)
                             cAGCGTGCCGCCAGATCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCa                 >  1:2996687/1‑68 (MQ=255)
                             cAGCGTGCCGCCAGATCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCa                 >  1:3093026/1‑68 (MQ=255)
                             cAGCGTGCCGCCAGATCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCa                 >  1:1117705/1‑68 (MQ=255)
                             cAGCGTGCCGCCAGATCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCa                 >  1:2209703/1‑68 (MQ=255)
                             cAGCGTGCCGCCAGATCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCa                 >  1:1456143/1‑68 (MQ=255)
                             cAGCGTGCCGCCAGATCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCa                 >  1:1942256/1‑68 (MQ=255)
                                   gccgccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCATTGAAgt          <  1:751/69‑1 (MQ=255)
                                   gccgccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAgt          <  1:3270355/69‑1 (MQ=255)
                                   gccgccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAgt          <  1:1012108/69‑1 (MQ=255)
                                   gccgccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAgt          <  1:2335554/69‑1 (MQ=255)
                                   gccgccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAgt          <  1:1299671/69‑1 (MQ=255)
                                    ccgccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAgtg         >  1:45938/1‑69 (MQ=255)
                                    ccgccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAgtg         >  1:2119347/1‑69 (MQ=255)
                                    ccgccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAgtg         <  1:1882673/69‑1 (MQ=255)
                                    ccgccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAgtg         >  1:448934/1‑69 (MQ=255)
                                    ccgccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAgtg         >  1:2697342/1‑69 (MQ=255)
                                    ccgccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAgtg         >  1:2742514/1‑69 (MQ=255)
                                    ccgccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAgtg         >  1:1406839/1‑69 (MQ=255)
                                    ccgccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAgt          <  1:561656/68‑1 (MQ=255)
                                    ccgccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAgt          <  1:2163311/68‑1 (MQ=255)
                                      gccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAgtgt        <  1:1378797/68‑1 (MQ=255)
                                      gccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAgtgt        <  1:2198626/68‑1 (MQ=255)
                                      gccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAgtgt        <  1:2993689/68‑1 (MQ=255)
                                      gccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAgtgt        <  1:2901104/68‑1 (MQ=255)
                                      gccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAgtgt        <  1:2800546/68‑1 (MQ=255)
                                      gccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAgtgt        <  1:2531541/68‑1 (MQ=255)
                                      gccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAgtgt        <  1:251642/68‑1 (MQ=255)
                                      gccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAgtgt        <  1:55268/68‑1 (MQ=255)
                                      gccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAgtgt        <  1:490847/68‑1 (MQ=255)
                                         aGATCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCa    <  1:607064/69‑1 (MQ=255)
                                         aGATCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCa    <  1:891483/69‑1 (MQ=255)
                                         aGATCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCa    <  1:3284115/69‑1 (MQ=255)
                                         aGATCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCa    <  1:2282014/69‑1 (MQ=255)
                                         aGATCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCa    <  1:2035441/69‑1 (MQ=255)
                                         aGATCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCa    <  1:1906030/69‑1 (MQ=255)
                                         aGATCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCa    <  1:145240/69‑1 (MQ=255)
                                         aGATCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCa    <  1:1132102/69‑1 (MQ=255)
                                           aaCGCACCAGCAAGCTGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAgg  <  1:228187/69‑1 (MQ=255)
                                           aaCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAgg  <  1:3091856/69‑1 (MQ=255)
                                           aaCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAgg  <  1:5593/69‑1 (MQ=255)
                                           aaCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAgg  <  1:1528287/69‑1 (MQ=255)
                                           aaCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAgg  <  1:1481740/69‑1 (MQ=255)
                                           aaCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAgg  <  1:886876/69‑1 (MQ=255)
                                           aaCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCACCGCCAATGAAGTGTTGCAgg  <  1:1945778/69‑1 (MQ=255)
                                            |                                                                   
CATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGG  >  minE/88300‑88411

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: