Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,501,821 A→G 42.8% D269D (GAT→GAC lrhA ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,501,8210AG42.8% 21.3 / 65.1 84D269D (GAT→GAClrhADNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (20/28);  new base G (13/23);  total (33/51)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

GCCACGGATTATGGTAGCTCTCCATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACA  >  minE/1501754‑1501884
                                                                   |                                                               
gCCACGGATTATGGTAGCTCTCCATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGa                                                                 >  1:3364785/1‑68 (MQ=255)
gCCACGGATTATGGTAGCTCTCCATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGAt                                                                >  1:108943/1‑69 (MQ=255)
gCCACGGATTATGGTAGCTCTCCATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGAt                                                                >  1:1173902/1‑69 (MQ=255)
gCCACGGATTATGGTAGCTCTCCATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGAt                                                                >  1:2532044/1‑69 (MQ=255)
gCCACGGATTATGGTAGCTCTCCATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGAt                                                                >  1:2473215/1‑69 (MQ=255)
gCCACGGATTATGGTAGCTCTCCATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGAt                                                                >  1:2941318/1‑69 (MQ=255)
gCCACGGATTATGGTAGCTCTCCATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGAt                                                                >  1:1601262/1‑69 (MQ=255)
   aCGGATTATGGTAGCTCTCCATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCa                                                              <  1:3167774/68‑1 (MQ=255)
   aCGGATTATGGTAGCTCTCCATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCa                                                              <  1:750877/68‑1 (MQ=255)
   aCGGATTATGGTAGCTCTCCATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCa                                                              <  1:1155920/68‑1 (MQ=255)
   aCGGATTATGGTAGCTCTCCATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCa                                                              <  1:659328/68‑1 (MQ=255)
   aCGGATTATGGTAGCTCTCCATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCa                                                              <  1:2742074/68‑1 (MQ=255)
   aCGGATTATGGTAGCTCTCCATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCa                                                              <  1:3211390/68‑1 (MQ=255)
   aCGGATTATGGTAGCTCTCCATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCa                                                              <  1:2959263/68‑1 (MQ=255)
    cGGATTATGGTAGCTCTCCATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCa                                                              <  1:1916778/67‑1 (MQ=255)
             gTAGCTCTCCATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCaa                                                   <  1:52533/69‑1 (MQ=255)
                        tccccTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGGTCATAGCACAGCAAATATTCGGt                                          <  1:944526/64‑1 (MQ=255)
                        tCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGGTCATAGCACAGCAAATATTCGGt                                          <  1:2998632/67‑1 (MQ=255)
                        tCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGGTCATAGCACAGCAAATATTCGGt                                          <  1:2556292/67‑1 (MQ=255)
                        tCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGGTCATAGCACAGCAAATATTCGGt                                          <  1:3064543/67‑1 (MQ=255)
                        tCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGGTCATAGCACAGCAAATATTCGGt                                          <  1:2250581/67‑1 (MQ=255)
                        tCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGGTCATAGCACAGCAAATATTCGGt                                          <  1:2059549/67‑1 (MQ=255)
                        tCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGGTCATAGCACAGCAAATATTCGGt                                          <  1:2789200/67‑1 (MQ=255)
                        tCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGGTCATAGCACAGCAAATATTCGGt                                          <  1:2558467/67‑1 (MQ=255)
                        tCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGGTCATAGCACAGCAAATATTCGGt                                          <  1:3180164/67‑1 (MQ=255)
                        tCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGGTCATAGCACAGCAAATATTCGGt                                          <  1:1856510/67‑1 (MQ=255)
                        tCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGGTCATAGCACAGCAAATATTCGGt                                          <  1:3215207/67‑1 (MQ=255)
                        tCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGGTCATAGCACAGCAAATATTCGGt                                          <  1:1703122/67‑1 (MQ=255)
                        tCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGGTCATAGCACAGCAAATATTCGGt                                          <  1:1698666/67‑1 (MQ=255)
                        tCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGGTCATAGCACAGCAAATATTCGGt                                          <  1:1575933/67‑1 (MQ=255)
                        tCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGGTCATAGCACAGCAAATATTCGGt                                          <  1:41255/67‑1 (MQ=255)
                        tCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGGTCATAGCACAGCAAATATTCGGt                                          <  1:604448/67‑1 (MQ=255)
                        tCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGGTCATAGCACAGCAAATATTCGGt                                          <  1:14586/67‑1 (MQ=255)
                        tCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGGTCATAGCACAGCAAATATTCGGt                                          <  1:1431077/67‑1 (MQ=255)
                        tCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGGTCATAGCACAGCAAATATTCGGt                                          <  1:767444/67‑1 (MQ=255)
                        tCGCCTGATAGATCACCTGGGCCACCTCATTATTACTGGACGGGTCATAGCACAGCAAATATTCGGt                                          <  1:3428560/67‑1 (MQ=255)
                             tGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGt                                   <  1:2912638/69‑1 (MQ=255)
                             tGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGt                                   <  1:622933/69‑1 (MQ=255)
                             tGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGt                                   <  1:1018272/69‑1 (MQ=255)
                             tGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGt                                   <  1:2541288/69‑1 (MQ=255)
                             tGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGt                                   <  1:2443683/69‑1 (MQ=255)
                             tGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGt                                   <  1:1044813/69‑1 (MQ=255)
                              gATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGt                                   <  1:1496412/68‑1 (MQ=255)
                                      aCCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGt                                   <  1:52171/60‑1 (MQ=255)
                                         tGGGCCAGCTCATTATTACTGGACGGGTCATAGCACAGCAAATATTCGGt                                          <  1:3197562/50‑1 (MQ=255)
                                         tGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAg                       <  1:1469481/69‑1 (MQ=255)
                                         tGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAg                       <  1:1153099/69‑1 (MQ=255)
                                         tGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAg                       <  1:3400043/69‑1 (MQ=255)
                                         tGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAg                       <  1:760586/69‑1 (MQ=255)
                                         tGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAg                       <  1:2345736/69‑1 (MQ=255)
                                         tGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAg                       <  1:2876208/69‑1 (MQ=255)
                                         tGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAg                       <  1:2909868/69‑1 (MQ=255)
                                           ggCCAGCTCATTATTACTGGACGGGTCATAGCACAGCAAATATTCGGt                                          <  1:2021387/48‑1 (MQ=255)
                                              cAGCTCATTATTACTGGACGGGTCATAGCACAGCAAATATTCGGt                                          <  1:1630418/45‑1 (MQ=255)
                                              cAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATc                  <  1:197540/69‑1 (MQ=255)
                                              cAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATc                  <  1:1856314/69‑1 (MQ=255)
                                              cAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATc                  <  1:2363978/69‑1 (MQ=255)
                                                 cTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTAc               >  1:589796/1‑69 (MQ=255)
                                                  tCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAg                       <  1:1009487/60‑1 (MQ=255)
                                                   cATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGTCGGTAAGCCATCTACTc             >  1:1347767/1‑69 (MQ=255)
                                                   cATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTc             >  1:2343343/1‑69 (MQ=255)
                                                   cATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTc             >  1:594984/1‑69 (MQ=255)
                                                   cATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTc             >  1:2527420/1‑69 (MQ=255)
                                                   cATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTc             >  1:363125/1‑69 (MQ=255)
                                                   cATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTc             >  1:339637/1‑69 (MQ=255)
                                                   cATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTc             >  1:3271636/1‑69 (MQ=255)
                                                   cATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTc             >  1:3254110/1‑69 (MQ=255)
                                                   cATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTc             >  1:1844917/1‑69 (MQ=255)
                                                   cATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTc             >  1:3205434/1‑69 (MQ=255)
                                                   cATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTc             >  1:94226/1‑69 (MQ=255)
                                                       attaCTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCa           >  1:3435504/1‑67 (MQ=255)
                                                       attaCTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCa           >  1:467996/1‑67 (MQ=255)
                                                                cGGGTCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCATcaca  >  1:1282923/1‑67 (MQ=255)
                                                                cGGGTCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGtaca  >  1:592293/1‑67 (MQ=255)
                                                                cGGGTCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGcaca  >  1:751635/1‑67 (MQ=255)
                                                                cGGGTCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGcaca  >  1:965869/1‑67 (MQ=255)
                                                                cGGGTCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGcaca  >  1:1471243/1‑67 (MQ=255)
                                                                cGGGTCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGcaca  >  1:419204/1‑67 (MQ=255)
                                                                cGGGTCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGcaca  >  1:1512434/1‑67 (MQ=255)
                                                                cGGGTCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGcaca  >  1:3235023/1‑67 (MQ=255)
                                                                cGGGTCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGcaca  >  1:2311275/1‑67 (MQ=255)
                                                                cGGGTCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGcaca  >  1:2312253/1‑67 (MQ=255)
                                                                cGGGTCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGcaca  >  1:2792506/1‑67 (MQ=255)
                                                                cGGGTCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGcaca  >  1:2643981/1‑67 (MQ=255)
                                                                cGGGTCATAGCACAGCAAATATTAGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGcaca  >  1:512031/1‑67 (MQ=255)
                                                                   |                                                               
GCCACGGATTATGGTAGCTCTCCATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACA  >  minE/1501754‑1501884

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: