Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 634,587 A→C 42.9% R521R (CGA→CGC mukB → fused chromosome partitioning proteins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE634,5870AC42.9% 10.1 / 67.2 63R521R (CGA→CGCmukBfused chromosome partitioning proteins
Reads supporting (aligned to +/- strand):  ref base A (20/16);  new base C (4/23);  total (24/39)
Fisher's exact test for biased strand distribution p-value = 1.45e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTCTGCTGGCA  >  minE/634528‑634644
                                                           |                                                         
ttGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGaa                                                  >  1:2501500/1‑69 (MQ=255)
ttGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGaa                                                  >  1:3146750/1‑69 (MQ=255)
ttGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGaa                                                  >  1:1747045/1‑69 (MQ=255)
ttGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGaa                                                  >  1:45615/1‑69 (MQ=255)
ttGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGaa                                                  >  1:1479720/1‑69 (MQ=255)
ttGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGaa                                                  >  1:2299641/1‑69 (MQ=255)
ttGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGaa                                                  >  1:1119803/1‑69 (MQ=255)
ttGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGa                                                   >  1:3230675/1‑68 (MQ=255)
ttGCGAGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGaa                                                  >  1:3039307/1‑69 (MQ=255)
 tGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCg                                                    >  1:789670/1‑66 (MQ=255)
     cgAAGGGGTCGATCAGCGTCACCTGGCATAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGAACTgg                                              <  1:1279710/68‑1 (MQ=255)
     cgAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGAACTgg                                              <  1:2266450/68‑1 (MQ=255)
     cgAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGAACTgg                                              <  1:111115/68‑1 (MQ=255)
     cgAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGAACTgg                                              <  1:1894456/68‑1 (MQ=255)
     cgAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCCATTAAGCGAACTgg                                              <  1:3376759/68‑1 (MQ=255)
          gggTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGAACTGGAACAGc                                        >  1:2816144/1‑69 (MQ=255)
            gTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:2735088/69‑1 (MQ=255)
            gTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:3403309/69‑1 (MQ=255)
            gTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:2918566/69‑1 (MQ=255)
            gTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:2976428/69‑1 (MQ=255)
            gTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:1015414/69‑1 (MQ=255)
            gTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:621023/69‑1 (MQ=255)
            gTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:623343/69‑1 (MQ=255)
            gTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:633206/69‑1 (MQ=255)
            gTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:981558/69‑1 (MQ=255)
            gTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:1867107/69‑1 (MQ=255)
            gTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:1003077/69‑1 (MQ=255)
            gTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:1362912/69‑1 (MQ=255)
            gTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:1716397/69‑1 (MQ=255)
            gTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:1794336/69‑1 (MQ=255)
            gTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:1982910/69‑1 (MQ=255)
            gTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:2105047/69‑1 (MQ=255)
            gTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:2230329/69‑1 (MQ=255)
            gTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:235974/69‑1 (MQ=255)
             tCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:461175/68‑1 (MQ=255)
             tCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:2434888/68‑1 (MQ=255)
                     cGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:264855/60‑1 (MQ=255)
                       tCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGagc                            >  1:1034511/1‑68 (MQ=255)
                       tCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGagc                            >  1:742866/1‑68 (MQ=255)
                       tCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGagc                            >  1:1900350/1‑68 (MQ=255)
                         aCCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGTACTGGAACAGCGTCTGCGTGagca                           >  1:2218468/1‑67 (MQ=255)
                         aCCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGTCTGCGTGagca                           >  1:961619/1‑67 (MQ=255)
                         aCCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGTCTGCGTGagca                           >  1:3400591/1‑67 (MQ=255)
                         aCCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGTCTGCGTGagc                            >  1:881425/1‑66 (MQ=255)
                          ccTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:2502463/55‑1 (MQ=255)
                                    cAGGTTCAGCCGCTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGGTGAg              >  1:2987077/1‑69 (MQ=255)
                                    cAGGTTCAGCCGCTGAGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAg              >  1:3397770/1‑69 (MQ=255)
                                      ggTTCAGCCGCTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCg            >  1:2933701/1‑69 (MQ=255)
                                        ttCAGCCGCTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTc          <  1:1059388/69‑1 (MQ=255)
                                        ttCAGCCGCTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTc          <  1:2519368/69‑1 (MQ=255)
                                        ttCAGCCGCTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTc          <  1:2991857/69‑1 (MQ=255)
                                         tCAGCCGCTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTct         >  1:3418083/1‑69 (MQ=255)
                                         tCAGCCGCTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTct         >  1:1750044/1‑69 (MQ=255)
                                         tCAGCCGCTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTct         >  1:914845/1‑69 (MQ=255)
                                          cAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:2746563/39‑1 (MQ=37)
                                                cTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTCTGCTGGCa  <  1:875047/69‑1 (MQ=255)
                                                cTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTCTGCTGGCa  <  1:2357181/69‑1 (MQ=255)
                                                cTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTCTGCTGGCa  <  1:753567/69‑1 (MQ=255)
                                                cTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTCTGCTGGCa  <  1:1386724/69‑1 (MQ=255)
                                                cTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTCTGCTGGCa  <  1:3224818/69‑1 (MQ=255)
                                                cTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTCTGCTGGCa  <  1:2592138/69‑1 (MQ=255)
                                                cTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTCTGCTGGCa  <  1:2982265/69‑1 (MQ=255)
                                                  gCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTCTGCTGGCa  <  1:2496993/67‑1 (MQ=255)
                                                           |                                                         
TTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTCTGCTGGCA  >  minE/634528‑634644

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: