Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,155,567 A→T 7.5% 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,5670AT 7.5% 144.0 / 4.4 53intergenic (+78/+108)dacB/obgED‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base A (24/25);  new base T (2/2);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGC  >  minE/2155504‑2155636
                                                               |                                                                      
cGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTa                                                                        <  1:4033543/64‑1 (MQ=255)
 ggCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCtt                                                                <  1:710142/71‑1 (MQ=255)
   cacaTGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTaaa                                                                      >  1:673511/1‑63 (MQ=255)
   cacaTGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTAt                                                              <  1:1407639/71‑1 (MQ=255)
     caTGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATcc                                                            <  1:2461240/71‑1 (MQ=255)
      aTGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCg                                                           <  1:169438/71‑1 (MQ=255)
       tGGCCGGGGCTTAAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCgg                                                          >  1:1481361/1‑71 (MQ=255)
            ggggCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                      >  1:1901834/1‑70 (MQ=255)
                    aGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGt                                             >  1:3904086/1‑71 (MQ=255)
                    aGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGt                                             <  1:1843593/71‑1 (MQ=255)
                      attgttGACAAGGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                            >  1:3074179/1‑70 (MQ=11)
                        tattGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTaaa                                                                      >  1:1460514/1‑42 (MQ=255)
                        tattGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTaaa                                                                      >  1:138373/1‑42 (MQ=255)
                          ttGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                       <  1:3978830/71‑1 (MQ=255)
                            gACAAATTGCGCTTTATGTTGCCGGATGCGGCGTa                                                                        >  1:1374726/1‑35 (MQ=0)
                            gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                      >  1:1963225/1‑70 (MQ=40)
                             aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                     <  1:2485088/70‑1 (MQ=255)
                             aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                     <  1:411988/70‑1 (MQ=39)
                             aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                     <  1:300113/70‑1 (MQ=39)
                              cAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCaa                                   >  1:242103/1‑70 (MQ=11)
                              cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                   <  1:625065/71‑1 (MQ=25)
                              cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                   >  1:2413976/1‑71 (MQ=35)
                               acaGTGCGCTTTGTTTATGCCAGATGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCaa                                   >  1:3598848/3‑70 (MQ=11)
                                aaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACa                                 >  1:85039/1‑71 (MQ=38)
                                 aGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACa                                 >  1:1464419/1‑70 (MQ=39)
                                      gcTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGc                           <  1:1571595/71‑1 (MQ=39)
                                       cTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGc                           >  1:187516/1‑70 (MQ=40)
                                         ttGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATAtt                             >  1:3317833/1‑66 (MQ=40)
                                          tGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAAtt                       >  1:294861/1‑71 (MQ=255)
                                           gTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtc                      >  1:2543323/1‑71 (MQ=255)
                                           gTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtc                      >  1:580429/1‑71 (MQ=255)
                                            tttATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCaata                                 <  1:3599264/59‑3 (MQ=255)
                                              tATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACatat                              >  1:1554861/1‑60 (MQ=255)
                                              tATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtctc                   <  1:2417866/70‑1 (MQ=255)
                                               aTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTctct                  <  1:109204/70‑1 (MQ=255)
                                                 gCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTg                >  1:2310589/1‑70 (MQ=255)
                                                  ccGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTg                >  1:3414943/1‑69 (MQ=255)
                                                   cGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtc                      <  1:1244988/63‑1 (MQ=255)
                                                    ggATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAgg            >  1:1199472/1‑71 (MQ=255)
                                                     gATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACata                               >  1:1385463/1‑52 (MQ=255)
                                                     gATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTg                <  1:2005899/66‑1 (MQ=255)
                                                      aTGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCtt                 <  1:2525949/64‑1 (MQ=255)
                                                      aTGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCtt                 <  1:2520956/64‑1 (MQ=255)
                                                      aTGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCtt                 <  1:2042141/64‑1 (MQ=255)
                                                      aTGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTg                <  1:3420709/65‑1 (MQ=255)
                                                       tGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                     <  1:2928581/44‑1 (MQ=38)
                                                       tGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCt         >  1:286787/1‑71 (MQ=255)
                                                         cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCaaa                                        >  1:3315927/1‑39 (MQ=255)
                                                         cggcgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGa       >  1:53033/1‑71 (MQ=255)
                                                         cggcgTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCt         <  1:1238079/70‑1 (MQ=11)
                                                          ggcgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACa                                 <  1:858955/45‑1 (MQ=255)
                                                          ggcgTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGa       >  1:2890438/1‑71 (MQ=14)
                                                           gcgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGAt      >  1:924625/1‑70 (MQ=255)
                                                            cgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGAt      >  1:4094102/1‑69 (MQ=255)
                                                              tAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACa                                 <  1:1949693/41‑1 (MQ=255)
                                                              tAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACa                                 <  1:3514097/41‑1 (MQ=255)
                                                              tAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACa                                 <  1:3162022/41‑1 (MQ=255)
                                                              tAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGc  <  1:2910537/71‑1 (MQ=255)
                                                               |                                                                      
CGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGC  >  minE/2155504‑2155636

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: