Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 983,519 G→A 14.9% intergenic (+431/‑137) ydfC → / → dicB hypothetical protein/cell division inhibition protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE983,5190GA14.9% 115.3 / 11.2 47intergenic (+431/‑137)ydfC/dicBhypothetical protein/cell division inhibition protein
Reads supporting (aligned to +/- strand):  ref base G (22/18);  new base A (2/5);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

CAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGT  >  minE/983452‑983588
                                                                   |                                                                     
cAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAgc                                                                    >  1:2105614/1‑71 (MQ=255)
    tttACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCAc                                                                >  1:2003373/1‑71 (MQ=255)
     ttACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCAcg                                                               >  1:1446336/1‑71 (MQ=255)
     ttACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCAcg                                                               >  1:3508298/1‑71 (MQ=255)
     ttACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCAc                                                                <  1:1242156/70‑1 (MQ=255)
      tACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCAcg                                                               <  1:1626949/70‑1 (MQ=255)
         tCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCgg                                                            >  1:703380/1‑70 (MQ=255)
           cGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGaa                                                          >  1:1744422/1‑70 (MQ=255)
             tGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGATTGCGGCTGTGCGCACGCGGAACa                                                        <  1:3273197/70‑1 (MQ=255)
                cTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGGACGCGGAACAGTTa                                                    <  1:1872369/71‑1 (MQ=255)
                   ctctGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTaaaa                                                 >  1:2045059/1‑71 (MQ=255)
                   ctctGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTaa                                                   <  1:1253196/69‑1 (MQ=255)
                   ctctGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTaa                                                   <  1:3112833/69‑1 (MQ=255)
                       gccgccCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCa                                              >  1:2303951/1‑70 (MQ=255)
                           cccTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCAcgc                                                              <  1:3592838/50‑1 (MQ=255)
                            ccTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCaaaaa                                          >  1:3708794/1‑69 (MQ=255)
                            ccTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCaaaaa                                          >  1:2745282/1‑69 (MQ=255)
                             cTTATTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTaa                                                   <  1:2741076/59‑1 (MQ=255)
                              tttttAAAGTGAATTTTGTGATGCGGTGAATGCGGCTAAGCGCACGTGGCACAGTTaaaa                                                 <  1:2502227/60‑1 (MQ=255)
                              tttttAAAGTGAATTTTGTCATGCGGTGAATGCGGCTAAGCGCACGTGGCACAGTTaaaa                                                 <  1:1748936/60‑1 (MQ=255)
                               ttttAAAGTGAATTTTGTGATGCGGTGAATGCGGCTAAGCGCACGTGGCACAGTTaaaagt                                               <  1:447940/61‑3 (MQ=255)
                                tttAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGAGCACGCgg                                                            >  1:3257778/1‑47 (MQ=255)
                                tttAAAGTGAATTTTGTGATGCGGTGAATGCGGCTAAGCGCACGTggc                                                           >  1:800831/1‑47 (MQ=37)
                                  tAAAGTGAATTTTGTGATGCGGTGAATGCGGCTAAGCGCACGTGGCACAGt                                                      <  1:113572/51‑1 (MQ=37)
                                    aaGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTAt                                 <  1:368105/70‑1 (MQ=255)
                                    aaGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTAt                                 <  1:608316/70‑1 (MQ=255)
                                    aaGTGAATTTTGTGATGCGGTGAATGCGGCTAAGCGCACGTGGCACAGTTa                                                    >  1:1951489/1‑51 (MQ=37)
                                       tGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCgg                                                            <  1:991038/40‑1 (MQ=255)
                                       tGAATTTTGTGATGCGGTGAATGCGGCTAAGCGCACGTGGCACAg                                                       <  1:2499119/45‑1 (MQ=25)
                                         aaTTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATgg                               <  1:77106/67‑1 (MQ=255)
                                         aaTTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTg                            <  1:2978353/70‑1 (MQ=255)
                                           ttttGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGAt                         >  1:3282629/1‑71 (MQ=255)
                                             ttGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTa                                  >  1:459098/1‑60 (MQ=255)
                                               gtgATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATtctc                     >  1:1177374/1‑71 (MQ=255)
                                                tgATGTGGTGAATGCGGCTGAGCGCACGCGGAACa                                                        >  1:2955667/1‑35 (MQ=255)
                                                 gATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTg                            >  1:1738694/1‑62 (MQ=255)
                                                   tgtgGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTa                 <  1:1642272/71‑1 (MQ=255)
                                                     tggtgAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTAt                <  1:2617353/70‑1 (MQ=255)
                                                     tggtgAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTAt                <  1:2612737/70‑1 (MQ=255)
                                                       gtgAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTgg                           >  1:1505269/1‑57 (MQ=255)
                                                       gtgAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTg                   >  1:201841/1‑65 (MQ=255)
                                                       gtgAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCg             >  1:545256/1‑71 (MQ=255)
                                                        tgAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCgg            <  1:1623228/71‑1 (MQ=255)
                                                         gAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAAc                                         >  1:3534131/1‑41 (MQ=255)
                                                         gAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGt                             >  1:1533084/1‑53 (MQ=255)
                                                          aaTGCGGCTGAGCGCACGCGGAACAGTTAAAACAAAAAACAGTGTTATGGGTGGAGTctct                    >  1:483760/1‑61 (MQ=255)
                                                                 cTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGt                             <  1:2666997/45‑1 (MQ=255)
                                                                   gAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGt  >  1:3049490/1‑70 (MQ=255)
                                                                   |                                                                     
CAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGT  >  minE/983452‑983588

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: