Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,155,638 T→C 16.4% intergenic (+149/+37) 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
*Exported2,155,6380TC16.4% 194.7 / 25.4 79intergenic (+149/+37)dacB/obgED‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base T (48/18);  new base C (4/9);  total (52/27)
Fisher's exact test for biased strand distribution p-value = 7.99e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

CTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCG  >  Exported/2155573‑2155708
                                                                  |                                                                      
cTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAgcg                                                                    >  1:1525361/1‑70 (MQ=255)
   aTCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCAt                                                                 >  1:2995949/1‑70 (MQ=255)
   aTCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAt                                                                 >  1:2438123/1‑71 (MQ=11)
    tCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATc                                                                <  1:218615/71‑1 (MQ=11)
        gCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGc                                                            <  1:4763606/71‑1 (MQ=11)
        gCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGc                                                            <  1:2681351/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:4678403/4‑61 (MQ=2)
             cataaTCTTGAAAATTCAATATATTGCAGTTATCTGGTAGGCCAGGTAAGCGTAGCGCATc                                                                >  1:3137942/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:1090193/1‑70 (MQ=255)
              aaaaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGAttt                                                      >  1:1547663/1‑70 (MQ=255)
              aaaaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTg                                                     >  1:3810780/1‑71 (MQ=255)
              aaaaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTg                                                     >  1:3142566/1‑71 (MQ=255)
                aaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCAt                                                                 >  1:2970338/1‑57 (MQ=255)
                aaTCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGc                                                   >  1:4833084/1‑71 (MQ=255)
                 aatgTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCg                                                  <  1:4549032/68‑1 (MQ=2)
                  tCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGAtt                                                       <  1:1816328/65‑1 (MQ=255)
                       cAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCa                                                               <  1:830321/53‑1 (MQ=11)
                       cAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGttgcc                                             <  1:1749085/70‑4 (MQ=2)
                        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)
                          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)
                           ttCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCt                                                           <  1:1600126/52‑1 (MQ=255)
                            tCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAg                                       >  1:5107934/1‑71 (MQ=255)
                            tCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAg                                       >  1:5525137/1‑71 (MQ=255)
                                cATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTa                                              <  1:1990295/60‑1 (MQ=255)
                                cATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTa                                              <  1:5125964/60‑1 (MQ=255)
                                 atatTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGAtt                                  <  1:1322188/71‑1 (MQ=255)
                                   atTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCa                                                                  <  1:2561492/37‑1 (MQ=14)
                                   atTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCa                                                                  <  1:2887096/37‑1 (MQ=14)
                                    ttGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAAc                               >  1:2779820/1‑71 (MQ=255)
                                      gCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCa                                                               >  1:422365/1‑37 (MQ=255)
                                      gCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCa                                                               >  1:3699390/1‑37 (MQ=255)
                                        aaTTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGa                                                         <  1:3559613/41‑1 (MQ=255)
                                          ttCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTatcat                                          <  1:2978480/54‑1 (MQ=255)
                                           tctcTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCt                            >  1:4014597/1‑67 (MQ=255)
                                            ctctTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGt                         >  1:2199740/1‑69 (MQ=255)
                                            ctctTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTa                        >  1:368561/1‑70 (MQ=255)
                                            ctctTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTa                        >  1:1341972/1‑70 (MQ=255)
                                            ctctTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTa                        >  1:4052852/1‑70 (MQ=255)
                                               ctTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTatca                                           <  1:3920681/49‑1 (MQ=255)
                                                ttGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTaa                                >  1:1400714/1‑59 (MQ=255)
                                                ttGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTaa                       <  1:1620921/68‑1 (MQ=255)
                                                ttGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTaa                       <  1:1140017/68‑1 (MQ=255)
                                                    aGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTAAATGAAc                 >  1:653652/1‑70 (MQ=255)
                                                    aGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTAAATGAACt                >  1:5032469/1‑71 (MQ=255)
                                                     ggCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTAAATGAACTc               <  1:5274481/71‑1 (MQ=255)
                                                      gCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTAAATGAACTCa              <  1:5155150/71‑1 (MQ=255)
                                                         tGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTAAATGAACTCa              <  1:3616268/68‑1 (MQ=255)
                                                           aTAAGCGTAGCGCATCAGGCTGATTTGGCGTTTAtcatca                                        >  1:5010749/1‑40 (MQ=255)
                                                           aTAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGt                                      <  1:1475873/42‑1 (MQ=255)
                                                           aTAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGt                                      <  1:2800478/42‑1 (MQ=255)
                                                           aTAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGAtt                                  >  1:4029900/1‑46 (MQ=255)
                                                           aTAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTa                                 >  1:813473/1‑47 (MQ=255)
                                                           aTAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGcc         >  1:1989766/1‑71 (MQ=255)
                                                            tAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTAAATg                    <  1:4430770/59‑1 (MQ=255)
                                                            tAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCct        >  1:5529151/1‑71 (MQ=255)
                                                             aaGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGa                                    >  1:4503272/1‑42 (MQ=255)
                                                             aaGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGa                                    >  1:5526966/1‑42 (MQ=255)
                                                              aGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTAAATGaa                  >  1:4364609/1‑59 (MQ=255)
                                                              aGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTAAATGaa                  >  1:1988480/1‑59 (MQ=255)
                                                              aGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTAAATGa                   >  1:2557608/1‑58 (MQ=255)
                                                              aGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCcttc      >  1:990033/1‑71 (MQ=255)
                                                                cGTAGCGCATCAGGCTGATTTGGCGTTTAtcatca                                        >  1:5451966/1‑35 (MQ=255)
                                                                cGTAGCGCATCAGGCTGATTTGGCGTTTAtcatca                                        >  1:4749567/1‑35 (MQ=255)
                                                                cGTAGCGCATCAGGCTGATTTGGCGTTTAtcatca                                        >  1:1892441/1‑35 (MQ=255)
                                                                cGTAGCGCATCAGGCTGATTTGGCGTTTAtcatca                                        >  1:896941/1‑35 (MQ=255)
                                                                 gTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTaaa                      >  1:608927/1‑52 (MQ=255)
                                                                  tAGCGCATCAGGCTGATTTGGCGTTTATCATCAGt                                      <  1:2459959/35‑1 (MQ=255)
                                                                  tAGCGCATCAGGCTGATTTGGCGTTTATCATCAGt                                      <  1:1245485/35‑1 (MQ=255)
                                                                  tAGCGCATCAGGCTGATTTGGCGTTTATCATCAGt                                      <  1:113179/35‑1 (MQ=255)
                                                                  tAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGa                                    >  1:5050123/1‑37 (MQ=255)
                                                                  tAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGa                                    >  1:4166956/1‑37 (MQ=255)
                                                                  tAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGa                                    >  1:3358836/1‑37 (MQ=255)
                                                                  tAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCcttctt    >  1:5283982/1‑69 (MQ=255)
                                                                  tAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCttcttc   >  1:358710/1‑70 (MQ=255)
                                                                  tAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTtcg  <  1:2021258/71‑1 (MQ=255)
                                                                  |                                                                      
CTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTAAATGAACTCAACGCCTTCTTCG  >  Exported/2155573‑2155708

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: