Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,155,567 A→T 19.3% intergenic (+78/+108) 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,5670AT19.3% 98.8 / 24.9 57intergenic (+78/+108)dacB/obgED‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base A (26/20);  new base T (5/6);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 5.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATA  >  Exported/2155497‑2155633
                                                                      |                                                                   
gAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTa                                                                     <  1:3823536/71‑1 (MQ=255)
        ggCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCt                                                              >  1:4772423/1‑70 (MQ=255)
           acaTGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATc                                                          >  1:3364822/1‑71 (MQ=255)
           acaTGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATc                                                          >  1:4800180/1‑71 (MQ=255)
            caTGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTa                                                                     <  1:1301565/59‑1 (MQ=255)
              tGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTa                                                                     <  1:1927413/57‑1 (MQ=255)
              tGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTa                                                                     <  1:2977235/57‑1 (MQ=255)
              tGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTa                                                                     <  1:2921788/57‑1 (MQ=255)
                gCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                      <  1:4255506/70‑1 (MQ=255)
                    cggCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                  >  1:3673691/2‑70 (MQ=11)
                     ggCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaa                                                <  1:1553601/71‑1 (MQ=11)
                     ggCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                 >  1:3050819/1‑70 (MQ=11)
                     ggCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                 >  1:5520019/1‑70 (MQ=255)
                       cTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                              <  1:3245084/71‑1 (MQ=11)
                       cTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                              <  1:3841057/71‑1 (MQ=11)
                       cTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                              <  1:3285390/71‑1 (MQ=255)
                        tacAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                              <  1:3603031/52‑1 (MQ=11)
                        tacAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaaa                                              <  1:4236095/68‑1 (MQ=11)
                        tacAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaaa                                              <  1:1356754/68‑1 (MQ=11)
                          cAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCg                                           >  1:4382492/1‑71 (MQ=255)
                           aGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTaaa                                                                   >  1:1573383/1‑46 (MQ=255)
                             attgttGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                         >  1:474442/1‑70 (MQ=255)
                               tgttGACAACGTACGCATTGTTCATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAACAGTGCa                                       >  1:534654/3‑70 (MQ=11)
                               tgttGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                            >  1:4728969/3‑49 (MQ=11)
                               tgttGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                        >  1:339141/3‑69 (MQ=255)
                               tgctGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCa                                       >  1:4908578/4‑70 (MQ=11)
                                gctGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAt                                    >  1:2276917/3‑71 (MQ=11)
                                 ttGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                  <  1:2908284/57‑1 (MQ=40)
                                 ctGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAAt                                    >  1:4521451/2‑71 (MQ=255)
                                  tGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                         >  1:4614147/1‑48 (MQ=11)
                                  tGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTa                                                            <  1:3205179/46‑1 (MQ=40)
                                  aaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                   >  1:297689/3‑71 (MQ=11)
                                   gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                   >  1:1553547/1‑70 (MQ=255)
                                    aCAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCaa                                >  1:4820455/1‑71 (MQ=11)
                                    aCAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCaa                                >  1:2350845/1‑71 (MQ=11)
                                    aCAAATTGCGCTTTATGTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCaa                                >  1:4762957/1‑71 (MQ=11)
                                    aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                  >  1:1072952/1‑70 (MQ=39)
                                     aacaGTGCGCTTTGTTTATGCCAGATGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCaa                                >  1:4747307/4‑71 (MQ=11)
                                      aaaGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCaat                               <  1:422614/71‑2 (MQ=255)
                                        aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaatac                             <  1:4900152/71‑4 (MQ=9)
                                          tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                >  1:2710053/1‑66 (MQ=255)
                                          tGCGCTTTGTTTATGCCAGATGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCaata                              <  1:4553164/68‑3 (MQ=11)
                                          tGCGCTTTGTTTATGCCAGATGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCaata                              <  1:4563233/68‑3 (MQ=11)
                                              cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                       <  1:970542/55‑1 (MQ=255)
                                                ttGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAAt                     >  1:1546328/1‑71 (MQ=255)
                                                 tGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACatat                           >  1:5692738/1‑64 (MQ=40)
                                                 tGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACatat                           >  1:276101/1‑64 (MQ=40)
                                                      aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                <  1:2334484/54‑1 (MQ=255)
                                                       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)
                                                                     tAAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                    <  1:614982/35‑1 (MQ=40)
                                                                     tAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                  <  1:3992907/37‑1 (MQ=40)
                                                                     tAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGc                        <  1:5697768/47‑1 (MQ=255)
                                                                     tAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTg             <  1:1200361/57‑1 (MQ=255)
                                                                      aaaCGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                  <  1:3829410/36‑1 (MQ=40)
                                                                      aaaCGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                  <  1:79329/36‑1 (MQ=40)
                                                                      |                                                                   
GAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATA  >  Exported/2155497‑2155633

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: