Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,155,544 T→G 14.3% intergenic (+55/+131) 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,5440TG14.3% 99.9 / 17.2 58intergenic (+55/+131)dacB/obgED‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base T (24/25);  new base G (4/4);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.38e-01

TATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACAT  >  minE/2155475‑2155607
                                                                     |                                                               
tATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCt                                                                 <  1:2648301/70‑1 (MQ=255)
                 aaaaaGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                               >  1:2760005/1‑71 (MQ=255)
                  aaaaGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATgcgg                                              >  1:146625/1‑71 (MQ=255)
                   aaaGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                             >  1:2048622/1‑71 (MQ=255)
                      gAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                           >  1:837092/1‑70 (MQ=255)
                      gAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTa                                          <  1:3823536/71‑1 (MQ=255)
                              ggCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCt                                   >  1:4772423/1‑70 (MQ=255)
                                 acaTGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGca                                                      >  1:4035108/1‑47 (MQ=255)
                                 acaTGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATc                               >  1:4800180/1‑71 (MQ=255)
                                 acaTGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATc                               >  1:3364822/1‑71 (MQ=255)
                                  caTGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGtt                                                            <  1:3065056/41‑1 (MQ=255)
                                  caTGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGtt                                                            <  1:3530875/41‑1 (MQ=255)
                                  caTGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTa                                          <  1:1301565/59‑1 (MQ=255)
                                    tGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTa                                          <  1:1927413/57‑1 (MQ=255)
                                    tGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTa                                          <  1:2977235/57‑1 (MQ=255)
                                    tGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTa                                          <  1:2921788/57‑1 (MQ=255)
                                      gCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGc                           <  1:4255506/70‑1 (MQ=255)
                                          cggCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                       >  1:3673691/2‑70 (MQ=11)
                                           ggCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaa                     <  1:1553601/71‑1 (MQ=11)
                                           ggCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                      >  1:3050819/1‑70 (MQ=11)
                                           ggCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                      >  1:5520019/1‑70 (MQ=255)
                                             cTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                   <  1:3245084/71‑1 (MQ=11)
                                             cTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                   <  1:3841057/71‑1 (MQ=11)
                                             cTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                   <  1:3285390/71‑1 (MQ=255)
                                              tacAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgc                                                <  1:255958/39‑1 (MQ=14)
                                              tacAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgc                                                <  1:5139261/39‑1 (MQ=14)
                                              tacAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgc                                                <  1:4175532/39‑1 (MQ=14)
                                              tacAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                   <  1:3603031/52‑1 (MQ=11)
                                              tacAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaaa                   <  1:1356754/68‑1 (MQ=11)
                                              tacAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaaa                   <  1:4236095/68‑1 (MQ=11)
                                               acAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgc                                                <  1:5145935/39‑1 (MQ=14)
                                                cAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCg                >  1:4382492/1‑71 (MQ=255)
                                                 aGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTaaa                                        >  1:1573383/1‑46 (MQ=255)
                                                   attgttGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGtt                                          <  1:864957/42‑2 (MQ=255)
                                                   attgttGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGtt                                          <  1:5410303/42‑2 (MQ=255)
                                                   attgttGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg              >  1:474442/1‑70 (MQ=255)
                                                   attgttGACAAAGTGCGCTTTGTTTATGCCGGATGCAGCGtt                                          <  1:1713143/42‑2 (MQ=11)
                                                     tgttGACAACGTACGCATTGTTCATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAACAGTGCa            >  1:534654/3‑70 (MQ=11)
                                                     tgttGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                 >  1:4728969/3‑49 (MQ=11)
                                                     tgttGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc             >  1:339141/3‑69 (MQ=255)
                                                     tgctGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCa            >  1:4908578/4‑70 (MQ=11)
                                                      gctGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAt         >  1:2276917/3‑71 (MQ=11)
                                                       ttGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                       <  1:2908284/57‑1 (MQ=40)
                                                       ctGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAAt         >  1:4521451/2‑71 (MQ=255)
                                                        tGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc                              >  1:4614147/1‑48 (MQ=11)
                                                        tGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTa                                 <  1:3205179/46‑1 (MQ=40)
                                                        aaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt        >  1:297689/3‑71 (MQ=11)
                                                         gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt        >  1:1553547/1‑70 (MQ=255)
                                                          aCAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCaa     >  1:4820455/1‑71 (MQ=11)
                                                          aCAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCaa     >  1:2350845/1‑71 (MQ=11)
                                                          aCAAATTGCGCTTTATGTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCaa     >  1:4762957/1‑71 (MQ=11)
                                                          aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc       >  1:1072952/1‑70 (MQ=39)
                                                           aacaGTGCGCTTTGTTTATGCCAGATGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCaa     >  1:4747307/4‑71 (MQ=11)
                                                            aaaGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCaat    <  1:422614/71‑2 (MQ=255)
                                                              aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaatac  <  1:4900152/71‑4 (MQ=9)
                                                                tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa     >  1:2710053/1‑66 (MQ=255)
                                                                tGCGCTTTGTTTATGCCAGATGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCaata   <  1:4563233/68‑3 (MQ=11)
                                                                tGCGCTTTGTTTATGCCAGATGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCaata   <  1:4553164/68‑3 (MQ=11)
                                                                    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa            <  1:970542/55‑1 (MQ=255)
                                                                     |                                                               
TATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACAT  >  minE/2155475‑2155607

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: