Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,155,529 A→G 35.3% intergenic (+40/+146) 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,5290AG35.3% 46.4 / 27.9 34intergenic (+40/+146)dacB/obgED‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base A (12/10);  new base G (2/10);  total (14/20)
Fisher's exact test for biased strand distribution p-value = 6.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.98e-01

TTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTG  >  minE/2155459‑2155595
                                                                      |                                                                  
tttGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGatta                                                                    <  1:2278123/71‑1 (MQ=255)
       aaaGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACa                                                              >  1:2133692/1‑70 (MQ=255)
                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:3530875/41‑1 (MQ=255)
                                                  caTGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGtt                                                <  1:3065056/41‑1 (MQ=255)
                                                  caTGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTa                              <  1:1301565/59‑1 (MQ=255)
                                                    tGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTa                              <  1:2977235/57‑1 (MQ=255)
                                                    tGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTa                              <  1:2921788/57‑1 (MQ=255)
                                                    tGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTa                              <  1:1927413/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:4175532/39‑1 (MQ=14)
                                                              tacAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgc                                    <  1:5139261/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:5410303/42‑2 (MQ=255)
                                                                   attgttGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGtt                              <  1:864957/42‑2 (MQ=255)
                                                                   attgttGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg  >  1:474442/1‑70 (MQ=255)
                                                                   attgttGACAAAGTGCGCTTTGTTTATGCCGGATGCAGCGtt                              <  1:1713143/42‑2 (MQ=11)
                                                                      |                                                                  
TTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTG  >  minE/2155459‑2155595

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: