Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 915,579 C→T 30.8% D43D (GAC→GAT ycjF → conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE915,5790CT30.8% 66.0 / 31.4 71D43D (GAC→GATycjFconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (25/24);  new base T (4/18);  total (29/42)
Fisher's exact test for biased strand distribution p-value = 1.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCC  >  minE/915517‑915646
                                                              |                                                                   
ttCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAgcg                                                               <  1:210617/69‑1 (MQ=255)
ttCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAgcg                                                               <  1:1177808/69‑1 (MQ=255)
ttCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAgcg                                                               <  1:1228726/69‑1 (MQ=255)
ttCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAgcg                                                               <  1:3363662/69‑1 (MQ=255)
       cgcAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGga                                                          >  1:1659562/1‑67 (MQ=255)
       cgcAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGga                                                          >  1:188140/1‑67 (MQ=255)
         cagcagACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGgaagaa                                                      <  1:3224182/69‑1 (MQ=255)
         cagcagACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGgaagaa                                                      <  1:2381961/69‑1 (MQ=255)
         cagcagACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGgaagaa                                                      <  1:169710/69‑1 (MQ=255)
         cagcagACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGgaagaa                                                      <  1:2724656/69‑1 (MQ=255)
         cagcagACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGgaagaa                                                      <  1:1556509/69‑1 (MQ=255)
         cagcagACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGgaagaa                                                      <  1:422238/69‑1 (MQ=255)
         cagcagACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGgaagaa                                                      <  1:467506/69‑1 (MQ=255)
              gACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAgg                                                  >  1:3029474/1‑68 (MQ=255)
              gACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAgg                                                  >  1:3009578/1‑68 (MQ=255)
              gACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAgg                                                  >  1:2276362/1‑68 (MQ=255)
              gACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAgg                                                  >  1:2592850/1‑68 (MQ=255)
              gACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAgg                                                  >  1:2542177/1‑68 (MQ=255)
              gACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAgg                                                  >  1:1960599/1‑68 (MQ=255)
              gACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAgg                                                  >  1:2331430/1‑68 (MQ=255)
                 cTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGACACGCTCGACGAAGCGCAGGAAGAAGAGGGGCa                                              >  1:3210374/1‑69 (MQ=255)
                      aCGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGt                                           >  1:144512/1‑67 (MQ=255)
                      aCGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGt                                           >  1:45633/1‑67 (MQ=255)
                        gAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGa                                        <  1:3240454/68‑1 (MQ=255)
                        gAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGa                                        <  1:1841183/68‑1 (MQ=255)
                           aaTCAGGCGCAAAATTTTGCCCCAGCCACGCTCGATGAAGCGCAGGAAGAAGAGGGGCAAGTTGAAGCg                                    >  1:3044757/1‑69 (MQ=255)
                           aaTCAGGCGCAAAATTTTGCCCCAGCCACGCTCGATGAAGCGCAGGAAGAAGAGGGGCAAGTTGAAGCg                                    >  1:3369714/1‑69 (MQ=255)
                           aaTCAGGCGCAAAATTTTGCCCCAGCCACGCTCGATGAAGCGCAGGAAGAAGAGGGGCAAGTTGAAGCg                                    >  1:472637/1‑69 (MQ=255)
                           aaTCAGGCGCAAAATTTTGCCCCAGCCACGCTCGATGAAGCGCAGGAAGAAGAGGGGCAAGTTGAAGCg                                    >  1:713080/1‑69 (MQ=255)
                             tCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGt                                  <  1:2435398/69‑1 (MQ=255)
                             tCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGt                                  <  1:2221946/69‑1 (MQ=255)
                             tCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGt                                  <  1:977760/69‑1 (MQ=255)
                             tCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGt                                  <  1:1597176/69‑1 (MQ=255)
                             tCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGt                                  <  1:1125791/69‑1 (MQ=255)
                              cAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGaagaag                                                     <  1:1905728/49‑1 (MQ=255)
                              cAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGaagaag                                                     <  1:2495635/49‑1 (MQ=255)
                              cAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGaagaag                                                     <  1:2967937/49‑1 (MQ=255)
                               aGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAgg                                                           >  1:1111680/1‑42 (MQ=255)
                               aGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTaa                                >  1:1046424/1‑69 (MQ=255)
                               aGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTaa                                >  1:37611/1‑69 (MQ=255)
                               aGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTaa                                >  1:1444269/1‑69 (MQ=255)
                               aGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTaa                                >  1:3105065/1‑69 (MQ=255)
                               aGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTaa                                >  1:2980140/1‑69 (MQ=255)
                               aGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTaa                                >  1:649640/1‑69 (MQ=255)
                               aGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTaa                                >  1:2709471/1‑69 (MQ=255)
                               aGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTaa                                >  1:260296/1‑69 (MQ=255)
                               aGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTaa                                >  1:172796/1‑69 (MQ=255)
                               aGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTaa                                >  1:945360/1‑69 (MQ=255)
                               aGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTaa                                >  1:1790034/1‑69 (MQ=255)
                               aGGCGCAAAATTTTGCCC‑‑‑CCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATgg                             >  1:1662482/1‑69 (MQ=255)
                                 gcgcAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCACGAAGAAGAGGGGCAAGTCGAAGCGGTAATg                              <  1:1793776/69‑1 (MQ=255)
                                                    ccACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGa           <  1:2567360/69‑1 (MQ=255)
                                                    ccACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGa           <  1:2441553/69‑1 (MQ=255)
                                                       cGCTCGATGAAGCGCAGGAAGAAGAGGGGCAAGTTGAAGCGTTGATGGACGCAGCGTTACGTCCGAAAc        <  1:815063/69‑1 (MQ=255)
                                                       cGCTCGATGAAGCGCAGGAAGAAGAGGGGCAAGTTGAAGCGGTGATGGACGCAGCGTTACGTCCGAAAc        <  1:964576/69‑1 (MQ=255)
                                                       cGCTCGATGAAGCGCAGGAAGAAGAGGGGCAAGTTGAAGCGGTGATGGACGCAGCGTTACGTCCGAAAc        <  1:1328745/69‑1 (MQ=255)
                                                       cGCTCGATGAAGCGCAGGAAGAAGAGGGGCAAGTTGAAGCGGTGATGGACGCAGCGTTACGTCCGAAAc        <  1:1589825/69‑1 (MQ=255)
                                                       cGCTCGATGAAGCGCAGGAAGAAGAGGGGCAAGTTGAAGCGGTGATGGACGCAGCGTTACGTCCGAAAc        <  1:3426919/69‑1 (MQ=255)
                                                       cGCTCGATGAAGCGCAGGAAGAAGAGGGGCAAGTTGAAGCGGTGATGGACGCAGCGTTACGTCCGAAAc        <  1:1480938/69‑1 (MQ=255)
                                                       cGCTCGATGAAGCGCAGGAAGAAGAGGGGCAAGTTGAAGCGGTGATGGACGCAGCGTTACGTCCGAAAc        <  1:1504681/69‑1 (MQ=255)
                                                       cGCTCGATGAAGCGCAGGAAGAAGAGGGGCAAGTTGAAGCGGTGATGGACGCAGCGTTACGTCCGAAAc        <  1:1507664/69‑1 (MQ=255)
                                                       cGCTCGATGAAGCGCAGGAAGAAGAGGGGCAAGTTGAAGCGGTGATGGACGCAGCGTTACGTCCGAAAc        <  1:1913508/69‑1 (MQ=255)
                                                       cGCTCGATGAAGCGCAGGAAGAAGAGGGGCAAGTTGAAGCGGTGATGGACGCAGCGTTACGTCCGAAAc        <  1:2993510/69‑1 (MQ=255)
                                                       cGCTCGATGAAGCGCAGGAAGAAGAGGGGCAAGTTGAAGCGGTGATGGACGCAGCGTTACGTCCGAAAc        <  1:1684446/69‑1 (MQ=255)
                                                       cGCTCGATGAAGCGCAGGAAGAAGAGGGGCAAGTTGAAGCGGTGATGGACGCAGCGTTACGTCCGAAAc        <  1:266895/69‑1 (MQ=255)
                                                       cGCTCGATGAAGCGCAGGAAGAAGAGGGGCAAGTTGAAGCGGTGATGGACGCAGCGTTACGTCCGAAAc        <  1:1758869/69‑1 (MQ=255)
                                                       cGCTCGATGAAGCGCAGGAAGAAGAGGGGCAAGTTGAAGCGGTGATGGACGCAGCGTTACGTCCGAAAc        <  1:1011695/69‑1 (MQ=255)
                                                       cGCTCGATGAAGCGCAGGAAGAAGAGGGGCAAGTTGAAGCGGTGATGGACGCAGCGTTACGTCCGAAAc        <  1:2237889/69‑1 (MQ=255)
                                                       cGCTCGATGAAGCGCAGGAAGAAGAGGGGCAAGTTGAAGCGGTATGGACGCAGCGTTACGTCCGAAAc        <  1:1559644/68‑1 (MQ=255)
                                                       cGCTCGATGAAGCGCAGGAAGAAGAGGGGCAAGATGAAGCGGTGATGGACGCAGCGTTACGTCCGAAAc        <  1:1516707/69‑1 (MQ=255)
                                                        gCTCGATGAAGCGCAGGAAGAAGAGGGGCAAGTTGAAGCGGTGATGGACGCAGCGTTACGTCCGAAAc        <  1:1133543/68‑1 (MQ=255)
                                                              cgaAGCTCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGcc  <  1:2770059/68‑1 (MQ=255)
                                                              cgaAGCGCATGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGGTACGTCCGAAACGCAGcc  <  1:3217680/68‑1 (MQ=255)
                                                              cgaAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGcc  <  1:676010/68‑1 (MQ=255)
                                                              |                                                                   
TTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCC  >  minE/915517‑915646

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