Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,155,605 C→T 41.7% intergenic (+116/+70) 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,6050CT41.7% 105.3 / 77.3 72intergenic (+116/+70)dacB/obgED‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base C (31/11);  new base T (16/14);  total (47/25)
Fisher's exact test for biased strand distribution p-value = 8.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

TTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAG  >  minE/2155545‑2155671
                                                            |                                                                   
ttGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAAt                                                           >  1:1546328/1‑71 (MQ=255)
 tGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACatat                                                                 >  1:5692738/1‑64 (MQ=40)
 tGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACatat                                                                 >  1:276101/1‑64 (MQ=40)
       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)
                     tAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGc                                                              <  1:5697768/47‑1 (MQ=255)
                     tAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTg                                                   <  1:1200361/57‑1 (MQ=255)
                       aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTa                                                 >  1:3049852/1‑57 (MQ=255)
                       aaCGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCg                                    >  1:4316326/1‑71 (MQ=11)
                        aCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAACATATTGCAATTCTCTTg                                                   <  1:3097118/54‑1 (MQ=255)
                        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:5235795/1‑58 (MQ=255)
                                      cTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAg                                 >  1:5233610/1‑58 (MQ=255)
                                       tACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGcaat                    >  1:5363373/1‑68 (MQ=0)
                                        aCAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGAtt                  >  1:381283/1‑71 (MQ=255)
                                         cataaTCTTGAAAATTCAATATATTGCAGTTATCTGGTAGGCCAGGTAAGCGTAGCGCATc                           >  1:3137942/4‑61 (MQ=2)
                                         cataaTCTTGAAAATTCAATATATTGCAGTTATCTGGTAGGCCAGGTAAGCGTAGCGCATc                           >  1:4678403/4‑61 (MQ=2)
                                         cataaTCTTGAAAATTCAATATATTGCAGTTATCTGGTAGGCCAGGTAAGCGTAGCGCATc                           >  1:2272322/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)
                                             aatgTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCg             <  1:4549032/68‑1 (MQ=2)
                                             aTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGc                                              >  1:724251/1‑38 (MQ=255)
                                              tCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGAtt                  <  1:1816328/65‑1 (MQ=255)
                                                 cgCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGc                                     >  1:4031576/2‑43 (MQ=11)
                                                 cgCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCg                                    >  1:4014754/2‑44 (MQ=11)
                                                 cgCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCg                                    >  1:3318039/2‑44 (MQ=11)
                                                   cAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt                                         >  1:1325812/1‑38 (MQ=9)
                                                   cAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATaa                                       >  1:4732087/1‑40 (MQ=9)
                                                   cAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCa                          <  1:830321/53‑1 (MQ=11)
                                                   cAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGttgcc        <  1:1749085/70‑4 (MQ=2)
                                                   cAAATTAAATAAATTGCAATTCAACTTGTAGGCCTGATaa                                       >  1:4538888/1‑40 (MQ=9)
                                                    aaaTTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGtt           <  1:163315/66‑1 (MQ=2)
                                                    aaaTTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTatca      >  1:4118040/1‑71 (MQ=255)
                                                    aaaTTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTatca      >  1:3386257/1‑71 (MQ=255)
                                                     aaTTCAATAAATTGCAATTCATCTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGtt           <  1:2665475/65‑1 (MQ=11)
                                                      aTTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGc                                     <  1:495751/39‑1 (MQ=2)
                                                      aTTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGc                                     <  1:5051387/39‑1 (MQ=2)
                                                      aTTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGcaa                     >  1:3344038/1‑53 (MQ=0)
                                                      aTTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGttgc         >  1:1678013/1‑64 (MQ=2)
                                                      aTTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTatcat     >  1:4919494/1‑70 (MQ=255)
                                                      aTTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTatcat     >  1:5502185/1‑70 (MQ=255)
                                                       ttCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCg                                    >  1:3941578/1‑39 (MQ=2)
                                                       ttCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCt                      <  1:1600126/52‑1 (MQ=255)
                                                        tCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAg  >  1:5525137/1‑71 (MQ=255)
                                                        tCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAg  >  1:5107934/1‑71 (MQ=255)
                                                            cATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTa         <  1:5125964/60‑1 (MQ=255)
                                                            cATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTa         <  1:1990295/60‑1 (MQ=255)
                                                            |                                                                   
TTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAG  >  minE/2155545‑2155671

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: