Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,705,446 C→T 6.5% A88A (GCC→GCT)  dam → DNA adenine methylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,705,4460CT 6.5% 1130.0 / 28.0 502A88A (GCC→GCT) damDNA adenine methylase
Reads supporting (aligned to +/- strand):  ref base C (251/217);  new base T (22/12);  total (273/229)
Fisher's exact test for biased strand distribution p-value = 2.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

Only 100 of 517 total aligned reads displayed.

AAGATGCGTACTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGCCAGGATCCGTTCCGTCGGGCGGTACTGT  >  minE/2705381‑2705513
                                                                 |                                                                   
aaGATGCGTACTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAg                                                                  >  1:1921063/1‑69 (MQ=255)
aaGATGCGTACTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAg                                                                  <  1:2076597/69‑1 (MQ=255)
  gATGCGTACTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAgg                                                                 <  1:2407572/68‑1 (MQ=255)
  gATGCGTACTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAgg                                                                 <  1:400776/68‑1 (MQ=255)
  gATGCGTACTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAgg                                                                 <  1:650137/68‑1 (MQ=255)
   aTGCGTACTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGtt                                                               <  1:1720633/69‑1 (MQ=255)
   aTGCGTACTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGtt                                                               <  1:300732/69‑1 (MQ=255)
   aTGCGTACTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGtt                                                               <  1:287963/69‑1 (MQ=255)
   aTGCGTACTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGtt                                                               <  1:1367519/69‑1 (MQ=255)
   aTGCGTACTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGtt                                                               <  1:2849704/69‑1 (MQ=255)
   aTGCGTACTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGtt                                                               <  1:2312899/69‑1 (MQ=255)
     gCGTACTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGttt                                                              >  1:1793644/1‑68 (MQ=255)
     gCGTACTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGtt                                                               >  1:58688/1‑67 (MQ=255)
     gCGTACTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGtt                                                               >  1:1946825/1‑67 (MQ=255)
     gCGTACTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGtt                                                               >  1:3271437/1‑67 (MQ=255)
     gCGTACTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGtt                                                               >  1:3121266/1‑67 (MQ=255)
      cGTACTGATGAGTACGTAGAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAg                                                                  >  1:566711/1‑63 (MQ=255)
      cGTACTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAg                                                                  >  1:156848/1‑63 (MQ=255)
      cGTACTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTAc                                                            >  1:2741587/1‑69 (MQ=255)
      cGTACTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTAc                                                            >  1:258153/1‑69 (MQ=255)
      cGTACTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTAc                                                            >  1:921719/1‑69 (MQ=255)
        tACTGATTAGTACGTACAGGCAGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCTGAGGTTTACt                                                           <  1:773317/68‑1 (MQ=255)
        tACTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCTGAGGTTTACt                                                           <  1:2526926/68‑1 (MQ=255)
        tACTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCTGAGGTTTACt                                                           <  1:3037173/68‑1 (MQ=255)
        tACTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCTGAGGTTTACt                                                           <  1:1449083/68‑1 (MQ=255)
        tACTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACt                                                           >  1:2926824/1‑68 (MQ=255)
        tACTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACt                                                           >  1:2111300/1‑68 (MQ=255)
        tACTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTa                                                          <  1:1880503/69‑1 (MQ=255)
          cTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAg                                                                  >  1:736155/1‑59 (MQ=255)
          cTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATc                                                        >  1:2619039/1‑69 (MQ=255)
          cTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATc                                                        >  1:1171396/1‑69 (MQ=255)
          cTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATc                                                        >  1:2007105/1‑69 (MQ=255)
          cTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATc                                                        <  1:2229927/69‑1 (MQ=255)
          cTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATc                                                        >  1:472012/1‑69 (MQ=255)
          cTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATc                                                        >  1:1341806/1‑69 (MQ=255)
             atgaGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGt                                                     >  1:2440512/1‑69 (MQ=255)
             atgaGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGt                                                     >  1:2786085/1‑69 (MQ=255)
             atgaGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGt                                                     >  1:2167397/1‑69 (MQ=255)
               gaGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGtt                                                    >  1:3403186/1‑68 (MQ=255)
                aGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTcc                                                  >  1:2396846/1‑69 (MQ=255)
                  tacgtacAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCcg                                                 >  1:1678006/1‑68 (MQ=255)
                  tacgtacAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCcg                                                 >  1:1217418/1‑68 (MQ=255)
                  tacgtacAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCcg                                                 >  1:2021376/1‑68 (MQ=255)
                  tacgtacAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCcg                                                 >  1:2630335/1‑68 (MQ=255)
                  tacgtacAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCcg                                                 >  1:1833550/1‑68 (MQ=255)
                  tacgtacAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCcg                                                 >  1:2964885/1‑68 (MQ=255)
                  tacgtacAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCcg                                                 >  1:1000873/1‑68 (MQ=255)
                   acgtacAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCcgcg                                               <  1:611667/69‑1 (MQ=255)
                    cgtacAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCcg                                                 >  1:987207/1‑66 (MQ=255)
                    cgtacAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCcg                                                 >  1:2355283/1‑66 (MQ=255)
                    cgtacAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCcg                                                 >  1:2475469/1‑66 (MQ=255)
                      tacAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGaa                                             <  1:896759/68‑1 (MQ=255)
                      tacAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGaa                                             <  1:3303671/68‑1 (MQ=255)
                      tacAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGaa                                             <  1:3139943/68‑1 (MQ=255)
                       acAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAaga                                           >  1:1114320/1‑69 (MQ=255)
                       acAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAaga                                           <  1:1410261/69‑1 (MQ=255)
                       acAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAaga                                           <  1:3321619/69‑1 (MQ=255)
                       acAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAaga                                           <  1:2955668/69‑1 (MQ=255)
                        cAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGt                                                     <  1:1270571/58‑1 (MQ=255)
                        cAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAagag                                          <  1:1538415/69‑1 (MQ=255)
                        cAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAagag                                          <  1:3430858/69‑1 (MQ=255)
                        cAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAagag                                          >  1:2149611/1‑69 (MQ=255)
                        cAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAagag                                          >  1:2269280/1‑69 (MQ=255)
                          ggCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGt                                         <  1:2896155/68‑1 (MQ=255)
                           gCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTc                                       >  1:147671/1‑69 (MQ=255)
                           gCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTc                                       >  1:2209824/1‑69 (MQ=255)
                           gCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTc                                       >  1:338499/1‑69 (MQ=255)
                           gCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTc                                       <  1:2689328/69‑1 (MQ=255)
                                   cgAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCcgcg                                               >  1:2048746/1‑53 (MQ=255)
                                   cgAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACaaaa                                >  1:540668/1‑68 (MQ=255)
                                    gAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGc                              <  1:1051802/69‑1 (MQ=255)
                                    gAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGc                              <  1:1506527/69‑1 (MQ=255)
                                     aGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGcc                             <  1:1638706/69‑1 (MQ=255)
                                      gCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGcc                             <  1:2865127/68‑1 (MQ=255)
                                         gtttgttCCCGAAACAAATTGCGCTGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGCCATGa                         >  1:2325335/1‑69 (MQ=255)
                                         gtttgttCCCGAAACAAATTGCGCTGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGCCAGGa                         >  1:1021932/1‑69 (MQ=255)
                                         gtttgttCCCGAAACAAATTGCGCTGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGCCAGGa                         >  1:3236221/1‑69 (MQ=255)
                                          tttgttCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGCCAGGAt                        >  1:2187623/1‑69 (MQ=255)
                                            tgttCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACaaa                                 <  1:1667351/58‑1 (MQ=255)
                                            tgttCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGcc                             <  1:963249/62‑1 (MQ=255)
                                            tgttCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGCCAGGAtcc                      >  1:1090179/1‑69 (MQ=255)
                                            tgttCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGCCAGGAtc                       <  1:1608947/68‑1 (MQ=255)
                                               tCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGCCAGGAtccgt                    <  1:3173990/68‑1 (MQ=255)
                                                  cGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGCCAGGAtccgttcc                 <  1:276725/68‑1 (MQ=255)
                                                  cGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGCCAGGAtccgttcc                 <  1:521254/68‑1 (MQ=255)
                                                  cGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGCCAGGAtccgttcc                 <  1:3339460/68‑1 (MQ=255)
                                                  cGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGCCAGGAtccgttcc                 <  1:310378/68‑1 (MQ=255)
                                                     aaCAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGCCAGGATCCGTTCCGTCg             <  1:43699/69‑1 (MQ=255)
                                                      aCAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGCCAGGATCCGTTCCGTCgg            >  1:1300062/1‑69 (MQ=255)
                                                      aCAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGCCAGGATCCGTTCCGTCg             >  1:384827/1‑68 (MQ=255)
                                                      aCAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGCCAGGATCCGTTCCGTCg             >  1:3216707/1‑68 (MQ=255)
                                                          aTTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGCCAGGATCCGTTCCGTCGGGCgg        <  1:1135830/69‑1 (MQ=255)
                                                          aTTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGCCAGGATCCGTTCCGTCGGGCgg        <  1:879463/69‑1 (MQ=255)
                                                          aTTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGCCAGGATCCGTTCCGTCGGGCgg        <  1:3056448/69‑1 (MQ=255)
                                                           ttGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGCCAGGATCCGTTCCGTCGGGCGGt       <  1:855864/69‑1 (MQ=255)
                                                           ttGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGCCAGGATCCGTTCCGTCGGGCGGt       <  1:159263/69‑1 (MQ=255)
                                                                ccGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGCCAGGATCCGTTCCGTCGGGCGGTACTg   <  1:142533/68‑1 (MQ=255)
                                                                ccGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGCCAGGATCCGTTCCGTCGGGCGGTACTg   <  1:1752819/68‑1 (MQ=255)
                                                                ccGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGCCAGGATCCGTTCCGTCGGGCGGTACTGt  <  1:719402/69‑1 (MQ=255)
                                                                ccGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGCCAGGATCCGTTCCGTCGGGCGGTACTGt  <  1:2821590/69‑1 (MQ=255)
                                                                 |                                                                   
AAGATGCGTACTGATGAGTACGTACAGGCCGCACGCGAGCTGTTTGTTCCCGAAACAAATTGCGCCGAGGTTTACTATCAGTTCCGCGAAGAGTTCAACAAAAGCCAGGATCCGTTCCGTCGGGCGGTACTGT  >  minE/2705381‑2705513

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: ‑