Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,155,616 T→G 6.9% intergenic (+127/+59) 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,6160TG 6.9% 143.0 / 7.1 72intergenic (+127/+59)dacB/obgED‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base T (40/27);  new base G (3/2);  total (43/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCT  >  Exported/2155552‑2155682
                                                                |                                                                   
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: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)
                                      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:3318039/2‑44 (MQ=11)
                                          cgCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCg                                               >  1:4014754/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:5051387/39‑1 (MQ=2)
                                               aTTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGc                                                <  1:495751/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: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:2887096/37‑1 (MQ=14)
                                                        atTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCa                                        <  1:2561492/37‑1 (MQ=14)
                                                         ttGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAAc     >  1:2779820/1‑71 (MQ=255)
                                                           gCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCa                                     >  1:3699390/1‑37 (MQ=255)
                                                           gCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCa                                     >  1:422365/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)
                                                                |                                                                   
TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATTAACGCT  >  Exported/2155552‑2155682

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: