Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,155,617:1 +A 7.1% intergenic (+128/+58) 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,6171.A 7.1% 270.3 / 17.2 70intergenic (+128/+58)dacB/obgED‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base . (43/22);  new base A (3/2);  total (46/24)
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

TTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTA  >  minE/2155548‑2155686
                                                                      |                                                                     
tttATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtc                                                                        <  1:112347/70‑1 (MQ=255)
  tATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTctct                                                                    >  1:160016/1‑71 (MQ=255)
      ccGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTa                                                                <  1:757168/71‑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)
              ggcgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGa                                                         >  1:1435700/1‑70 (MQ=255)
              ggcgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGAt                                                        <  1:2525497/71‑1 (MQ=255)
              ggcgTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGa                                                         >  1:3172492/1‑71 (MQ=14)
                cgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATaa                                                      <  1:1881545/71‑1 (MQ=255)
                    aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGc                                                             >  1:3349055/1‑60 (MQ=255)
                    aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCg                                                   >  1:1649692/1‑70 (MQ=255)
                    aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCg                                                   >  1:1818822/1‑70 (MQ=255)
                    aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCg                                                   >  1:2966096/1‑70 (MQ=255)
                     aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCCACATATTGCAATTCTCTTGTAGGCCTGATAAGCGtc                                                 >  1:900512/1‑70 (MQ=255)
                     aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCg                                                   >  1:505943/1‑69 (MQ=255)
                        ccTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAgcg                                              >  1:1178601/1‑71 (MQ=255)
                         cTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCgca                                                 <  1:1723047/68‑3 (MQ=11)
                           tATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCAt                                           >  1:3104888/1‑71 (MQ=255)
                            aTCCGGCCTACAAAATCGTGCAAATTCAACATATTGCTATTCTCTTGTAGGCCTGATAAGCGTAGCGCAt                                           >  1:828670/1‑70 (MQ=255)
                               cGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAgg                                       <  1:669121/71‑1 (MQ=255)
                                 gCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCt                                     >  1:1363920/1‑71 (MQ=255)
                                  ccTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGc                                      >  1:627270/1‑69 (MQ=255)
                                  ccTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCt                                     >  1:1387815/1‑70 (MQ=255)
                                  ccTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCt                                     >  1:1067304/1‑70 (MQ=255)
                                   cTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTg                                                                  <  1:2097864/40‑1 (MQ=255)
                                   cTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGcaa                                    >  1:1006971/1‑69 (MQ=2)
                                      cAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGt                                                                 >  1:2639210/1‑38 (MQ=255)
                                      cAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGAt                                                        >  1:15816/1‑47 (MQ=255)
                                      cAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGAttt                                >  1:380661/1‑71 (MQ=255)
                                      cAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGcaat                                   <  1:2393310/69‑4 (MQ=0)
                                       aaaaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATa                                                       <  1:3633402/47‑1 (MQ=255)
                                       aaaaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCg                                                   <  1:3844684/51‑1 (MQ=255)
                                       aaaaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCa                                         <  1:3614853/61‑1 (MQ=255)
                                       aaaaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCa                                         <  1:1851692/61‑1 (MQ=255)
                                       aaaaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCa                                         <  1:962781/61‑1 (MQ=255)
                                         aaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGc                             >  1:2928902/1‑71 (MQ=255)
                                         aaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGc                             >  1:983841/1‑71 (MQ=255)
                                         aaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGc                             >  1:227844/1‑71 (MQ=255)
                                           tCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGt                           >  1:725990/1‑71 (MQ=255)
                                           tCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGt                           >  1:2420039/1‑71 (MQ=255)
                                           tCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGt                           <  1:3388900/71‑1 (MQ=255)
                                            cGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGtt                          >  1:1361377/1‑71 (MQ=255)
                                            cGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGtt                          >  1:2834110/1‑71 (MQ=255)
                                               gCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTat                       >  1:36949/1‑71 (MQ=255)
                                                 aaaTTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTa                                                 <  1:1476321/43‑1 (MQ=255)
                                                 aaaTTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGttt                         <  1:265077/67‑1 (MQ=255)
                                                  aaTTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGAtt                                 >  1:2471532/1‑58 (MQ=255)
                                                  aaTTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGttt                         <  1:2116251/66‑1 (MQ=255)
                                                  aaTTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTatcat                    >  1:354385/1‑71 (MQ=255)
                                                   aTTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCg                                                   <  1:3824181/39‑1 (MQ=255)
                                                   aTTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTatcat                    >  1:2639318/1‑70 (MQ=255)
                                                   aTTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTatcat                    >  1:643573/1‑70 (MQ=255)
                                                   aTTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTatcat                    >  1:641304/1‑70 (MQ=255)
                                                    ttCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGTCGTTTAtcatca                  <  1:3049201/71‑1 (MQ=255)
                                                     tCAACATATTGCAATTCTCTTGTAGGCCTGATAAGc                                                    >  1:2589245/1‑36 (MQ=255)
                                                     tCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTa                                                 >  1:1665264/1‑39 (MQ=255)
                                                       aaCATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGt                <  1:3788766/70‑1 (MQ=255)
                                                        aCATATTGCAATTCTCTAGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCACTg               >  1:2637782/1‑70 (MQ=255)
                                                         cATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTatcat                    <  1:1192717/64‑1 (MQ=255)
                                                         cATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGAt             >  1:1179388/1‑71 (MQ=255)
                                                         cATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGAt             >  1:3320652/1‑71 (MQ=255)
                                                          atatTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGAtt            <  1:2267692/71‑1 (MQ=255)
                                                          atatTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGAtt            <  1:522826/71‑1 (MQ=255)
                                                               gCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACg        >  1:2054445/1‑70 (MQ=255)
                                                               gCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGc       >  1:2981596/1‑71 (MQ=255)
                                                                cAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCt      >  1:2716851/1‑71 (MQ=255)
                                                                 aaTTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGc       >  1:1095751/1‑69 (MQ=255)
                                                                  aTTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTATATCATCAGTGATTAACGCTTg    >  1:1558789/1‑71 (MQ=255)
                                                                   ttCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGc                             >  1:2536611/1‑45 (MQ=255)
                                                                   ttCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGt                <  1:2213412/58‑1 (MQ=255)
                                                                    tctcTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTa           <  1:1864194/62‑1 (MQ=255)
                                                                     ctctTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTa  <  1:1529365/70‑1 (MQ=255)
                                                                     ctctTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTa  >  1:1478130/1‑70 (MQ=255)
                                                                     ctctTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTa  >  1:3797727/1‑70 (MQ=255)
                                                                      |                                                                     
TTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTA  >  minE/2155548‑2155686

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: