Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,155,607 T→A 36.0% intergenic (+118/+68) 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,6070TA36.0% 122.7 / 63.6 72intergenic (+118/+68)dacB/obgED‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base T (32/13);  new base A (10/14);  total (45/27)
Fisher's exact test for biased strand distribution p-value = 2.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATT  >  Exported/2155545‑2155676
                                                              |                                                                      
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:2272322/4‑61 (MQ=2)
                                         cataaTCTTGAAAATTCAATATATTGCAGTTATCTGGTAGGCCAGGTAAGCGTAGCGCATc                                >  1:4678403/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: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: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:5125964/60‑1 (MQ=255)
                                                            cATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTa              <  1:1990295/60‑1 (MQ=255)
                                                             atatTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGAtt  <  1:1322188/71‑1 (MQ=255)
                                                              |                                                                      
TTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATT  >  Exported/2155545‑2155676

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: