Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,858,365 G→A 7.1% T213T (ACG→ACA amiB → N‑acetylmuramoyl‑l‑alanine amidase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,858,3650GA 7.1% 214.8 / 4.9 70T213T (ACG→ACAamiBN‑acetylmuramoyl‑l‑alanine amidase II
Reads supporting (aligned to +/- strand):  ref base G (30/35);  new base A (3/2);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

ATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTA  >  minE/2858295‑2858435
                                                                      |                                                                      
aTAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTAcg                                                                        <  1:4127050/71‑1 (MQ=255)
 tAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTAcgc                                                                       >  1:243696/1‑71 (MQ=255)
  aaaaTTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTAcgcg                                                                      >  1:1019880/1‑71 (MQ=255)
  aaaaTTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTAcgcg                                                                      >  1:3055129/1‑71 (MQ=255)
  aaaaTTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTAcgcg                                                                      >  1:1732525/1‑71 (MQ=255)
  aaaaTTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTAcgcg                                                                      >  1:4312116/1‑71 (MQ=255)
   aaaTTATCATCGCTATTGATGCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCgg                                                                     <  1:3623281/70‑1 (MQ=255)
    aaTTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCggg                                                                    <  1:36626/71‑1 (MQ=255)
    aaTTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCgg                                                                     <  1:484878/70‑1 (MQ=255)
     aTTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCgg                                                                     >  1:1612290/1‑69 (MQ=255)
         tcatcGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGaaa                                                               <  1:3512436/71‑1 (MQ=255)
         tcatcGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGaaa                                                               <  1:3900724/71‑1 (MQ=255)
         tcatcGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGaa                                                                <  1:3709452/70‑1 (MQ=255)
            tcGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAAt                                                            <  1:1738574/71‑1 (MQ=255)
                  ttGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGaaa                                                               <  1:922074/62‑1 (MQ=255)
                  ttGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCAc                                                       >  1:1919142/1‑70 (MQ=255)
                  ttGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCAc                                                       >  1:3269390/1‑70 (MQ=255)
                  ttGATGCCGGACACGGCGGTCAGGACCCTGGCGATATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCAc                                                       <  1:1047205/70‑1 (MQ=255)
                   tGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCAcca                                                     >  1:109300/1‑71 (MQ=255)
                   tGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCAc                                                       >  1:2556947/1‑69 (MQ=255)
                    gATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTAcgc                                                                       >  1:3471248/1‑52 (MQ=255)
                    gATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTc                                                         >  1:1329670/1‑66 (MQ=255)
                    gATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCAccat                                                    >  1:4345185/1‑71 (MQ=255)
                     aTGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCAccatc                                                   >  1:2889014/1‑71 (MQ=255)
                        ccGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCAccatcgcc                                                <  1:3144231/71‑1 (MQ=255)
                          ggACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCAccatcgccat                                              >  1:3103203/1‑71 (MQ=255)
                                  cggTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTa                                      <  1:1395357/71‑1 (MQ=255)
                                    gTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCAccatcgcc                                                <  1:3619113/59‑1 (MQ=255)
                                    gTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTaaa                                    <  1:3825678/71‑1 (MQ=255)
                                    gTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTa                                      <  1:4096391/69‑1 (MQ=255)
                                    gTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTa                                      <  1:2634266/69‑1 (MQ=255)
                                    gTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTa                                      <  1:1706670/69‑1 (MQ=255)
                                     tCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTaaa                                    <  1:4157100/70‑1 (MQ=255)
                                      cAGGACCCTGGCGCTATCGGCCCTGGTGGTACACGGGAGAAAAATGTCAccatcgccatt                                             >  1:853284/1‑59 (MQ=255)
                                      cAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAAtt                                  >  1:2536364/1‑71 (MQ=255)
                                       aGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAAtt                                  <  1:1172053/70‑1 (MQ=255)
                                       aGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAAtt                                  <  1:1950281/70‑1 (MQ=255)
                                       aGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTg                                 >  1:1560801/1‑71 (MQ=255)
                                        ggACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGc                                <  1:3356518/71‑1 (MQ=255)
                                         gACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCg                               >  1:515132/1‑71 (MQ=255)
                                          aCCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCg                               >  1:2300869/1‑70 (MQ=255)
                                            ccTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTAc                            <  1:2040667/71‑1 (MQ=255)
                                               ggCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACttt                         >  1:2894024/1‑71 (MQ=255)
                                               ggCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACttt                         >  1:1435308/1‑71 (MQ=255)
                                                gcgcTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACttt                         <  1:3903485/70‑1 (MQ=255)
                                                 cgcTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGc                       >  1:934213/1‑71 (MQ=255)
                                                   cTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTc                     <  1:4288368/71‑1 (MQ=255)
                                                   cTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTc                     <  1:2292727/71‑1 (MQ=255)
                                                    tATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTc                     >  1:4072431/1‑70 (MQ=255)
                                                     aTCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTa                             >  1:830566/1‑61 (MQ=255)
                                                     aTCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTa                             >  1:3930457/1‑61 (MQ=255)
                                                     aTCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCaa                   >  1:3197568/1‑71 (MQ=255)
                                                     aTCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCaa                   >  1:252820/1‑71 (MQ=255)
                                                      tCGGCCCTGGTGGTACACGGGAGAAAAATGTCACCATCGCCATTGCGCGTAAAttac                                <  1:1286003/57‑3 (MQ=255)
                                                      tCGGCCCTGGTGGTACACGGGAGAAAAATGTCACCATCGCCATTGCGCGTAAAttac                                <  1:2154138/57‑3 (MQ=255)
                                                      tCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCaa                   <  1:1314699/70‑1 (MQ=255)
                                                      tCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCa                    <  1:2891506/69‑1 (MQ=255)
                                                       cGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAAt                  >  1:400803/1‑70 (MQ=255)
                                                       cGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAAt                  >  1:2408133/1‑70 (MQ=255)
                                                        ggCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACtt                          >  1:4341446/1‑61 (MQ=255)
                                                        ggCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACtt                          >  1:2423720/1‑61 (MQ=255)
                                                        ggCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATg                 <  1:414418/70‑1 (MQ=255)
                                                            ctggtggtACACGGGAGAAAAATGTCACCATCGCCATTGCGCGTAAATTACGTACCTTGCTCAATAACGa             >  1:1779576/3‑70 (MQ=255)
                                                            ctggtggtACACGGGAGAAAAATGTCACCATCGCCATTGCGCGTAAATTACGTACCTTGCTCAATAACGAt            >  1:4387568/3‑71 (MQ=255)
                                                            ccGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGAt            <  1:2343633/71‑1 (MQ=255)
                                                             cGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGAt            <  1:1258117/70‑1 (MQ=255)
                                                             cGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGAt            <  1:3320091/70‑1 (MQ=255)
                                                              ggtggtACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTaaa                                    >  1:2900503/1‑45 (MQ=255)
                                                              ggtggtACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTATATTGAGTACTTTGCTCAATGACGATcc          >  1:1251470/1‑71 (MQ=255)
                                                              ggtggtACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATcc          <  1:603375/71‑1 (MQ=255)
                                                              ggtggtACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATcc          <  1:2190640/71‑1 (MQ=255)
                                                               gtggtACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATcc          >  1:2583043/1‑70 (MQ=255)
                                                               gtggtACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATcc          >  1:2290739/1‑70 (MQ=255)
                                                                tggtACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGa        <  1:2689720/71‑1 (MQ=255)
                                                                tggtACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGa        <  1:3836976/71‑1 (MQ=255)
                                                                 ggtACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGAt       >  1:3787403/1‑71 (MQ=255)
                                                                   tACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGa        <  1:490466/68‑1 (MQ=255)
                                                                     cgcgGGAGAAAAATGTCACCATCGCCATCGCGCGTAAAt                                   >  1:3351669/1‑39 (MQ=255)
                                                                     cgcgGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTc                     <  1:2071487/53‑1 (MQ=255)
                                                                     cgcgGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGAt       >  1:3385930/1‑67 (MQ=255)
                                                                     cgcgGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGtt    >  1:2096148/1‑70 (MQ=255)
                                                                     cgcgGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGtt    >  1:2528724/1‑70 (MQ=255)
                                                                      gcgGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTa  <  1:1518090/71‑1 (MQ=255)
                                                                      gcgGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTa  <  1:2125246/71‑1 (MQ=255)
                                                                      |                                                                      
ATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTA  >  minE/2858295‑2858435

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: