Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,155,618 T→A 38.3% intergenic (+129/+57) 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,6180TA38.3% 117.8 / 78.0 73intergenic (+129/+57)dacB/obgED‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base T (30/15);  new base A (13/15);  total (43/30)
Fisher's exact test for biased strand distribution p-value = 1.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTA  >  Exported/2155552‑2155686
                                                                   |                                                                    
tGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTg                                                                  <  1:4561491/71‑1 (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)
              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:4763606/71‑1 (MQ=11)
                             gCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGc                                      <  1:2681351/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:4678403/4‑61 (MQ=2)
                                  cataaTCTTGAAAATTCAATATATTGCAGTTATCTGGTAGGCCAGGTAAGCGTAGCGCATc                                          >  1:2272322/4‑61 (MQ=2)
                                  cataaTCTTGAAAATTCAATATATTGCAGTTATCTGGTAGGCCAGGTAAGCGTAGCGCATc                                          >  1:3137942/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: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:3386257/1‑71 (MQ=255)
                                             aaaTTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTatca                     >  1:4118040/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:5502185/1‑70 (MQ=255)
                                               aTTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTatcat                    >  1:4919494/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)
                                                      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:1341972/1‑70 (MQ=255)
                                                                 ctctTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTa  >  1:4052852/1‑70 (MQ=255)
                                                                 ctctTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTa  >  1:368561/1‑70 (MQ=255)
                                                                   |                                                                    
TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCTTGTA  >  Exported/2155552‑2155686

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: