Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 983,528 C→T 14.6% intergenic (+440/‑128) 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,5280CT14.6% 117.5 / 10.8 48intergenic (+440/‑128)ydfC/dicBhypothetical protein/cell division inhibition protein
Reads supporting (aligned to +/- strand):  ref base C (23/18);  new base T (2/5);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTT  >  minE/983461‑983597
                                                                   |                                                                     
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:2745282/1‑69 (MQ=255)
                   ccTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCaaaaa                                                   >  1:3708794/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:2612737/70‑1 (MQ=255)
                                            tggtgAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTAt                         <  1:2617353/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)
                                                             cgcACGCGGAACAGTTAAAACCAAAAACAGTGTTATgg                                        <  1:237180/38‑1 (MQ=255)
                                                             cgcACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAAc       <  1:622245/71‑1 (MQ=255)
                                                              gcACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGt                           >  1:1550826/1‑50 (MQ=255)
                                                               cACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGt                                      >  1:129334/1‑38 (MQ=255)
                                                                aCGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTgg    >  1:1613779/1‑71 (MQ=255)
                                                                aCGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTgg    <  1:3387694/71‑1 (MQ=255)
                                                                aCGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTgg    >  1:1055554/1‑71 (MQ=255)
                                                                   cgGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGtt  >  1:2365705/1‑70 (MQ=255)
                                                                   |                                                                     
TCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTT  >  minE/983461‑983597

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: