Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 634,576 C→T 46.5% L518L (CTG→TTG)  mukB → fused chromosome partitioning proteins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE634,5760CT46.5% 5.8 / 46.5 58L518L (CTG→TTG) mukBfused chromosome partitioning proteins
Reads supporting (aligned to +/- strand):  ref base C (22/9);  new base T (4/23);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 1.94e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

TCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTCTGCTGGCA  >  minE/634518‑634644
                                                          |                                                                    
tcgcGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGc                                                             <  1:1997735/68‑1 (MQ=255)
 cgcgAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCg                                                            >  1:2471483/1‑68 (MQ=255)
 cgcgAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCg                                                            >  1:2855887/1‑68 (MQ=255)
          ttGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGaa                                                  >  1:2501500/1‑69 (MQ=255)
          ttGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGaa                                                  >  1:2299641/1‑69 (MQ=255)
          ttGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGaa                                                  >  1:1119803/1‑69 (MQ=255)
          ttGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGaa                                                  >  1:3146750/1‑69 (MQ=255)
          ttGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGaa                                                  >  1:1747045/1‑69 (MQ=255)
          ttGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGaa                                                  >  1:1479720/1‑69 (MQ=255)
          ttGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGaa                                                  >  1:45615/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:1894456/68‑1 (MQ=255)
               cgAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGAACTgg                                              <  1:111115/68‑1 (MQ=255)
               cgAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGAACTgg                                              <  1:2266450/68‑1 (MQ=255)
               cgAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCCATTAAGCGAACTgg                                              <  1:3376759/68‑1 (MQ=255)
                    gggTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGAACTGGAACAGc                                        >  1:2816144/1‑69 (MQ=255)
                      gTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:2976428/69‑1 (MQ=255)
                      gTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:981558/69‑1 (MQ=255)
                      gTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:2735088/69‑1 (MQ=255)
                      gTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:633206/69‑1 (MQ=255)
                      gTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:2918566/69‑1 (MQ=255)
                      gTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:621023/69‑1 (MQ=255)
                      gTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:623343/69‑1 (MQ=255)
                      gTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:1015414/69‑1 (MQ=255)
                      gTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:3403309/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:2991857/69‑1 (MQ=255)
                                                  ttCAGCCGCTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTc          <  1:2519368/69‑1 (MQ=255)
                                                  ttCAGCCGCTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTc          <  1:1059388/69‑1 (MQ=255)
                                                   tCAGCCGCTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTct         >  1:1750044/1‑69 (MQ=255)
                                                   tCAGCCGCTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTct         >  1:3418083/1‑69 (MQ=255)
                                                   tCAGCCGCTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTct         >  1:914845/1‑69 (MQ=255)
                                                    cAGCCGTTGCGGATGCGCTTAAGCGAACTGGAACAGCGt                                      <  1:2746563/39‑1 (MQ=37)
                                                          cTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTCTGCTGGCa  <  1:2357181/69‑1 (MQ=255)
                                                          cTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTCTGCTGGCa  <  1:875047/69‑1 (MQ=255)
                                                          cTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTCTGCTGGCa  <  1:2982265/69‑1 (MQ=255)
                                                          cTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTCTGCTGGCa  <  1:753567/69‑1 (MQ=255)
                                                          cTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTCTGCTGGCa  <  1:1386724/69‑1 (MQ=255)
                                                          cTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTCTGCTGGCa  <  1:2592138/69‑1 (MQ=255)
                                                          cTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTCTGCTGGCa  <  1:3224818/69‑1 (MQ=255)
                                                          |                                                                    
TCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGCGAACTGGAACAGCGTCTACGCGAGCAGCAAGAAGCTGAGCGTCTGCTGGCA  >  minE/634518‑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: