Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,155,567 A→T 16.3% intergenic (+78/+108) dacB → / ← obgE D‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,155,5670AT16.3% 94.9 / 13.8 43intergenic (+78/+108)dacB/obgED‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base A (22/14);  new base T (5/2);  total (27/16)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

CCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAA  >  minE/2155502‑2155634
                                                                 |                                                                    
cccGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGc                                                                 >  1:1121781/1‑71 (MQ=255)
        aTGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCg                                                         >  1:1369230/1‑71 (MQ=255)
           gCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                      >  1:1195792/1‑71 (MQ=255)
           gCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                       <  1:2048745/70‑1 (MQ=255)
              ggggCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                    >  1:43025/1‑70 (MQ=255)
                 gCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaa                                                <  1:2269870/71‑1 (MQ=255)
                    acAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATc                                             >  1:1908706/2‑71 (MQ=255)
                     cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATc                                             <  1:2314111/70‑1 (MQ=255)
                     cAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCg                                            >  1:556950/1‑71 (MQ=255)
                     cAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCg                                            >  1:320434/1‑71 (MQ=255)
                       gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                          >  1:986703/1‑71 (MQ=255)
                       gATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                       <  1:2841771/58‑1 (MQ=255)
                         ttgttGACAACGTACGCATTGTTCATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAACAGTGCa                                        <  1:3369653/68‑1 (MQ=1)
                         ttattGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                     >  1:531414/1‑58 (MQ=255)
                           attGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTaaa                                                                    >  1:1615768/1‑41 (MQ=255)
                            ttGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                     >  1:653638/1‑71 (MQ=255)
                            ttGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCaaa                                      <  1:2719167/70‑1 (MQ=255)
                             tGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAAtt                                    >  1:178360/1‑71 (MQ=255)
                               aCAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCaa                                 >  1:1098698/1‑71 (MQ=11)
                               aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                    <  1:1707583/69‑1 (MQ=39)
                                 acaGTGCGCTTTGTTTATGCCAGATGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCaat                                <  1:3620551/69‑2 (MQ=11)
                                 aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                 >  1:2247584/1‑70 (MQ=25)
                                 aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                 >  1:2252239/1‑70 (MQ=25)
                                 aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                 >  1:1649628/1‑70 (MQ=35)
                                 aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                 >  1:1840022/1‑70 (MQ=35)
                                    gTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACata                             <  1:1502857/71‑1 (MQ=39)
                                     tGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACatat                            >  1:3163074/1‑71 (MQ=39)
                                       cgcTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTg                          >  1:592977/1‑71 (MQ=39)
                                              tttATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtc                    <  1:112347/70‑1 (MQ=255)
                                                tATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTctct                >  1:160016/1‑71 (MQ=255)
                                                    ccGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                     >  1:1427629/1‑47 (MQ=38)
                                                    ccGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTa            <  1:757168/71‑1 (MQ=255)
                                                     cGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                   <  1:3889732/48‑1 (MQ=255)
                                                      ggATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAg           <  1:2820287/70‑1 (MQ=255)
                                                       gATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg          >  1:378758/1‑71 (MQ=11)
                                                         tGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCt       >  1:2547970/1‑71 (MQ=255)
                                                          gcggcgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCt       >  1:1473641/1‑70 (MQ=255)
                                                            ggcgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                    >  1:924232/1‑40 (MQ=40)
                                                            ggcgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGa     >  1:1435700/1‑70 (MQ=255)
                                                            ggcgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGAt    <  1:2525497/71‑1 (MQ=255)
                                                            ggcgTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGa     >  1:3172492/1‑71 (MQ=14)
                                                              cgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTg                          <  1:204757/48‑1 (MQ=255)
                                                              cgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATaa  <  1:1881545/71‑1 (MQ=255)
                                                               gTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                   >  1:1065880/1‑38 (MQ=40)
                                                                 |                                                                    
CCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAA  >  minE/2155502‑2155634

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