Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,155,592 C→T 18.7% intergenic (+103/+83) 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,5920CT18.7% 177.0 / 29.7 81intergenic (+103/+83)dacB/obgED‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base C (45/20);  new base T (8/7);  total (54/27)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01

CAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGC  >  minE/2155523‑2155659
                                                                     |                                                                    
cAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCg                                                                     >  1:4382492/1‑71 (MQ=255)
   attgttGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                   >  1:474442/1‑70 (MQ=255)
     tgttGACAACGTACGCATTGTTCATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAACAGTGCa                                                                 >  1:534654/3‑70 (MQ=11)
     tgttGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                  >  1:339141/3‑69 (MQ=255)
     tgctGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCa                                                                 >  1:4908578/4‑70 (MQ=11)
      gctGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAt                                                              >  1:2276917/3‑71 (MQ=11)
       ctGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAAt                                                              >  1:4521451/2‑71 (MQ=255)
        aaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                             >  1:297689/3‑71 (MQ=11)
         gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                             >  1:1553547/1‑70 (MQ=255)
          aCAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCaa                                                          >  1:2350845/1‑71 (MQ=11)
          aCAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCaa                                                          >  1:4820455/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:4553164/68‑3 (MQ=11)
                tGCGCTTTGTTTATGCCAGATGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCaata                                                        <  1:4563233/68‑3 (MQ=11)
                    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                 <  1:970542/55‑1 (MQ=255)
                      ttGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAAt                                               >  1:1546328/1‑71 (MQ=255)
                       tGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACatat                                                     >  1:5692738/1‑64 (MQ=40)
                       tGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACatat                                                     >  1:276101/1‑64 (MQ=40)
                            aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                          <  1:2334484/54‑1 (MQ=255)
                             tGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTg                                       <  1:4561491/71‑1 (MQ=255)
                                  gATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGc                                                  >  1:444215/1‑56 (MQ=255)
                                  gATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCa                                                 >  1:3612334/1‑57 (MQ=255)
                                    tGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtct                                           <  1:891784/60‑1 (MQ=255)
                                    tGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCt                                <  1:5120457/71‑1 (MQ=255)
                                      cggcgTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCt                                >  1:253713/1‑70 (MQ=11)
                                      cggcgTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg                               <  1:2459078/71‑1 (MQ=14)
                                         cgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTg                                       <  1:695466/59‑1 (MQ=255)
                                         cgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGa                              >  1:2469852/1‑68 (MQ=255)
                                         cgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGAt                             >  1:2812655/1‑69 (MQ=255)
                                         cgTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa                            <  1:1772080/71‑1 (MQ=14)
                                           tAAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                                              <  1:614982/35‑1 (MQ=40)
                                           tAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                            <  1:3992907/37‑1 (MQ=40)
                                           tAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGc                                                  <  1:5697768/47‑1 (MQ=255)
                                           tAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTg                                       <  1:1200361/57‑1 (MQ=255)
                                            aaaCGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                            <  1:79329/36‑1 (MQ=40)
                                            aaaCGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                            <  1:3829410/36‑1 (MQ=40)
                                             aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                          >  1:1040965/1‑37 (MQ=37)
                                             aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                          >  1:909019/1‑37 (MQ=37)
                                             aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                          >  1:726113/1‑37 (MQ=37)
                                             aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                          >  1:1087824/1‑37 (MQ=37)
                                             aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                          >  1:1094565/1‑37 (MQ=37)
                                             aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                          >  1:2181420/1‑37 (MQ=37)
                                             aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTa                                     >  1:3049852/1‑57 (MQ=255)
                                             aaCGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCg                        >  1:4316326/1‑71 (MQ=11)
                                              aCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAACATATTGCAATTCTCTTg                                       <  1:3097118/54‑1 (MQ=255)
                                              aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                          >  1:3463578/1‑36 (MQ=33)
                                              aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTg                                       <  1:5603832/54‑1 (MQ=255)
                                              aCGCCTTATCCGGCCTACAAAATCCTGCAAATTCAACATATTGCAATTCTCTTg                                       <  1:3922463/54‑1 (MQ=255)
                                               cGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCgca                      <  1:775285/71‑3 (MQ=11)
                                                  cTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAgcg                   >  1:1525361/1‑70 (MQ=255)
                                                     aTCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtct                                           >  1:4672749/1‑43 (MQ=255)
                                                     aTCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCAt                >  1:2995949/1‑70 (MQ=255)
                                                     aTCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAt                >  1:2438123/1‑71 (MQ=11)
                                                      tCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtct                                           >  1:1891833/1‑42 (MQ=255)
                                                      tCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATc               <  1:218615/71‑1 (MQ=11)
                                                       ccGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGt                                      >  1:4635591/1‑46 (MQ=255)
                                                          gCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGc           <  1:2681351/71‑1 (MQ=11)
                                                          gCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGc           <  1:4763606/71‑1 (MQ=11)
                                                            cTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAg                     >  1:5233610/1‑58 (MQ=255)
                                                            cTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAg                     >  1:5235795/1‑58 (MQ=255)
                                                             tACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGcaat        >  1:5363373/1‑68 (MQ=0)
                                                              aCAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGAtt      >  1:381283/1‑71 (MQ=255)
                                                               cataaTCTTGAAAATTCAATATATTGCAGTTATCTGGTAGGCCAGGTAAGCGTAGCGCATc               >  1:2272322/4‑61 (MQ=2)
                                                               cataaTCTTGAAAATTCAATATATTGCAGTTATCTGGTAGGCCAGGTAAGCGTAGCGCATc               >  1:3137942/4‑61 (MQ=2)
                                                               cataaTCTTGAAAATTCAATATATTGCAGTTATCTGGTAGGCCAGGTAAGCGTAGCGCATc               >  1:4678403/4‑61 (MQ=2)
                                                               cAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGAttt     >  1:623795/1‑71 (MQ=255)
                                                               cAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGcaat        <  1:1641788/69‑4 (MQ=0)
                                                                aaaaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGAttt     >  1:1547663/1‑70 (MQ=255)
                                                                aaaaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGAttt     >  1:1090193/1‑70 (MQ=255)
                                                                aaaaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTg    >  1:3142566/1‑71 (MQ=255)
                                                                aaaaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTg    >  1:3810780/1‑71 (MQ=255)
                                                                  aaTCGTGCAAATTCAACATATTGCAATTCTCTTGt                                      >  1:3145127/1‑35 (MQ=255)
                                                                  aaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCAt                >  1:2970338/1‑57 (MQ=255)
                                                                  aaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGc  >  1:4833084/1‑71 (MQ=255)
                                                                   aTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGc                                  >  1:724251/1‑38 (MQ=255)
                                                                    tCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGAtt      <  1:1816328/65‑1 (MQ=255)
                                                                     |                                                                    
CAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGC  >  minE/2155523‑2155659

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: