Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 780,059 G→T 36.2% L55L (CTG→CTT ycgL → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE780,0590GT36.2% 77.5 / 56.9 80L55L (CTG→CTTycgLconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (16/35);  new base T (13/16);  total (29/51)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGCCGATATTGAAAAGGTTA  >  minE/780007‑780108
                                                    |                                                 
ttCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGc                                      >  1:320078/1‑66 (MQ=255)
ttCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGc                                      >  1:3192518/1‑66 (MQ=255)
ttCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGc                                      >  1:1785771/1‑66 (MQ=255)
ttCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGc                                      >  1:2688190/1‑66 (MQ=255)
ttCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGc                                      >  1:360571/1‑66 (MQ=255)
    tGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGgaaga                                <  1:393838/68‑1 (MQ=255)
    tGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGgaaga                                <  1:1533219/68‑1 (MQ=255)
    tGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGgaaga                                <  1:2920274/68‑1 (MQ=255)
     gAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGgaaga                                <  1:963900/67‑1 (MQ=255)
     gAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGgaaga                                <  1:1704694/67‑1 (MQ=255)
        gAACTGATGAAAGGTTTTGGTCAGCCTCGGTTAGCGATGATTCTTCCGCTGGATGGGCGGAAGAAAc                             <  1:2978832/67‑1 (MQ=255)
        gAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTTCCGCTGGATGGGCGGAAGAAAc                             <  1:2786820/67‑1 (MQ=255)
        gAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTTCCGCTGGATGGGCGGAAGAAAc                             <  1:1176247/67‑1 (MQ=255)
        gAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTTCCGCTGGATGGGCGGAAGAAAc                             <  1:404821/67‑1 (MQ=255)
        gAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTTCCGCTGGATGGGCGGAAGAAAc                             <  1:2261387/67‑1 (MQ=255)
        gAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTTCCGCTGGATGGGCGGAAGAAAc                             <  1:2164138/67‑1 (MQ=255)
        gAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTTCCGCTGGATGGGCGGAAGAAAc                             <  1:772752/67‑1 (MQ=255)
        gAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTTCCGCTGGATGGGCGGAAGAAAc                             <  1:831094/67‑1 (MQ=255)
        gAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTTCCGCTGGATGGGCGGAAGAAAc                             <  1:2000046/67‑1 (MQ=255)
        gAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTTCCGCTGGATGGGCGGAAGAAAc                             <  1:1722095/67‑1 (MQ=255)
        gAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTTCCGCTGGATGGGCGGAAGAAAc                             <  1:948223/67‑1 (MQ=255)
        gAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTTCCGCTGGATGGGCGGAAGAAAc                             <  1:1662009/67‑1 (MQ=255)
        gAACTGATCAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTTCCGCTGGATGGGCGGAAGAAAc                             <  1:1858466/67‑1 (MQ=255)
         aaCTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTTCCGCTGGATGGGCGGAAGAAAc                             <  1:2787151/66‑1 (MQ=255)
           cTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTc                        <  1:2245492/69‑1 (MQ=255)
           cTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTc                        <  1:250878/69‑1 (MQ=255)
           cTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTc                        <  1:2458780/69‑1 (MQ=255)
           cTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTc                        <  1:2387096/69‑1 (MQ=255)
           cTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTc                        <  1:2998128/69‑1 (MQ=255)
           cTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTc                        <  1:333829/69‑1 (MQ=255)
           cTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTc                        <  1:2054059/69‑1 (MQ=255)
            tgatgaAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTc                        <  1:1649650/68‑1 (MQ=255)
            tgatgaAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTc                        <  1:67156/68‑1 (MQ=255)
            tgatgaAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTc                        <  1:629721/68‑1 (MQ=255)
            tgatgaAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTc                        <  1:569006/68‑1 (MQ=255)
                       tttGGTCAGCCTCAGTTAGCGATGATTCTTCCGCTGGATGGGc                                      >  1:1264929/1‑43 (MQ=255)
                       tttGGTCAGCCTCAGTTAGCGATGATTCTTCCGCTGGATGGGc                                      >  1:2528919/1‑43 (MQ=255)
                       tttGGTCAGCCTCAGTTAGCGATGATTCTTCCGCTGGATGGGc                                      >  1:1224496/1‑43 (MQ=255)
                       tttGGTCAGCCTCAGTTAGCGATGATTCTTCCGCTGGATGGGc                                      >  1:1206930/1‑43 (MQ=255)
                       tttGGTCAGCCTCAGTTAGCGATGATTCTTCCGCTGGATGGGc                                      >  1:1677863/1‑43 (MQ=255)
                       tttGGTCAGCCTCAGTTAGCGATGATTCTTCCGCTGGATGGGc                                      >  1:2104937/1‑43 (MQ=255)
                       tttGGTCAGCCTCAGTTAGCGATGATTCTTCCGCTGGATGGGc                                      >  1:860141/1‑43 (MQ=255)
                       tttGGTCAGCCTCAGTTAGCGATGATTCTTCCGCTGGATGGGc                                      >  1:2314430/1‑43 (MQ=255)
                       tttGGTCAGCCTCAGTTAGCGATGATTCTTCCGCTGGATGGGc                                      >  1:1743260/1‑43 (MQ=255)
                       tttGGTCAGCCTCAGTTAGCGATGATTCTTCCGCTGGATGGGc                                      >  1:2244999/1‑43 (MQ=255)
                       tttGGTCAGCCTCAGTTAGCGATGATTCTTCCGCTGGATGGGc                                      >  1:412362/1‑43 (MQ=255)
                       tttGGTCAGCCTCAGTTAGCGATGATTCTTCCGCTGGATGGGc                                      >  1:684095/1‑43 (MQ=255)
                       tttGGTCAGCCTCAGTTAGCGATGATTCTTCCGCTGGATGGGc                                      >  1:175763/1‑43 (MQ=255)
                        ttGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAAt                     <  1:1037554/59‑1 (MQ=255)
                        ttGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAAt                     <  1:3417994/59‑1 (MQ=255)
                        ttGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAAt                     <  1:2163857/59‑1 (MQ=255)
                        ttGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAAt                     <  1:742392/59‑1 (MQ=255)
                          ggTCAGCCTCAGTTAGCGATGATTCTTCCGCTGGATGGGCGGAAGAAAc                             <  1:2651508/49‑1 (MQ=255)
                           gTCAGCCTCAGTTAGCGATGATTCTTCCGCTGGATGGGCGGAAGAAAc                             <  1:3079647/48‑1 (MQ=255)
                              aGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGCCGATATTGaaaa       >  1:3422973/1‑67 (MQ=255)
                              aGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGCCGATATTGAAAAg      >  1:530306/1‑68 (MQ=255)
                              aGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGCCGATATTGAAAAg      >  1:994184/1‑68 (MQ=255)
                              aGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGCCGATATTGAAAAg      >  1:1867905/1‑68 (MQ=255)
                              aGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGCCGATATTGAAAAg      >  1:237192/1‑68 (MQ=255)
                                 cTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGCCGATATTGAAAAGGTTa  <  1:2454915/69‑1 (MQ=255)
                                 cTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGCCGATATTGAAAAGGTTa  <  1:1370015/69‑1 (MQ=255)
                                 cTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGCCGATATTGAAAAGGTTa  <  1:1590877/69‑1 (MQ=255)
                                 cTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGCCGATATTGAAAAGGTTa  <  1:869429/69‑1 (MQ=255)
                                 cTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGCCGATATTGAAAAGGTTa  <  1:2112776/69‑1 (MQ=255)
                                 cTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGCCGATATTGAAAAGGTTa  <  1:2119105/69‑1 (MQ=255)
                                 cTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGCCGATATTGAAAAGGTTa  <  1:2129169/69‑1 (MQ=255)
                                 cTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGCCGATATTGAAAAGGTTa  <  1:2413215/69‑1 (MQ=255)
                                 cTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGCCGATATTGAAAAGGTTa  <  1:2622906/69‑1 (MQ=255)
                                 cTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGCCGATATTGAAAAGGTTa  <  1:2855000/69‑1 (MQ=255)
                                 cTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGCCGATATTGAAAAGGTTa  <  1:3382594/69‑1 (MQ=255)
                                        aGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGcc                  >  1:3316118/1‑46 (MQ=255)
                                        aGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGcc                  >  1:2631613/1‑46 (MQ=255)
                                        aGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGcc                  >  1:2686270/1‑46 (MQ=255)
                                        aGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGcc                  <  1:222026/46‑1 (MQ=255)
                                        aGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGcc                  >  1:2697058/1‑46 (MQ=255)
                                        aGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGcc                  <  1:2752342/46‑1 (MQ=255)
                                        aGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGcc                  <  1:3426945/46‑1 (MQ=255)
                                        aGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGcc                  >  1:1914635/1‑46 (MQ=255)
                                        aGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGcc                  <  1:2826623/46‑1 (MQ=255)
                                        aGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGcc                  >  1:310622/1‑46 (MQ=255)
                                                    |                                                 
TTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGCCGATATTGAAAAGGTTA  >  minE/780007‑780108

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: