Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 581,876 G→A 27.2% P638P (CCG→CCA clpA → ATPase and specificity subunit of ClpA‑ClpP ATP‑dependent serine protease, chaperone activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE581,8760GA27.2% 139.4 / 57.5 81P638P (CCG→CCAclpAATPase and specificity subunit of ClpA‑ClpP ATP‑dependent serine protease, chaperone activity
Reads supporting (aligned to +/- strand):  ref base G (31/28);  new base A (10/12);  total (41/40)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

TCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGATCATCTGTCAACCGACGTGATCCATCAGGT  >  minE/581808‑581944
                                                                    |                                                                    
tCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCgg                                                                     >  1:2364668/1‑70 (MQ=255)
tCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCgg                                                                     >  1:423341/1‑70 (MQ=255)
        tCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTa                                                            <  1:886140/71‑1 (MQ=255)
        tCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTa                                                            <  1:3329778/71‑1 (MQ=255)
        tCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTa                                                            <  1:1134426/71‑1 (MQ=255)
        tCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTa                                                            <  1:2832762/71‑1 (MQ=255)
          ttATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTAAc                                                          >  1:3784521/1‑71 (MQ=255)
          ttATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTAAc                                                          >  1:2638102/1‑71 (MQ=255)
          ttATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTAAc                                                          >  1:1176213/1‑71 (MQ=255)
           tATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCAGAATTCCGTAAcc                                                         >  1:2523944/1‑71 (MQ=255)
           tATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCAGAATTCCGTAAcc                                                         >  1:2215012/1‑71 (MQ=255)
           tATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCAGAATTCCGTAAcc                                                         >  1:4121645/1‑71 (MQ=255)
           tATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCAGAATTCCGTAAcc                                                         >  1:1476621/1‑71 (MQ=255)
           tATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCAGAATTCCGTAAcc                                                         >  1:841620/1‑71 (MQ=255)
            aTCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTaa                                                           <  1:1323300/68‑1 (MQ=255)
            aTCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTaa                                                           <  1:1597852/68‑1 (MQ=255)
            aTCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTAACCg                                                        >  1:2403995/1‑71 (MQ=255)
            aTCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTAACCg                                                        >  1:788334/1‑71 (MQ=255)
            aTCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTAACCg                                                        >  1:2727127/1‑71 (MQ=255)
            aTCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTAACCg                                                        >  1:271362/1‑71 (MQ=255)
            aTCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTAACCg                                                        >  1:2306650/1‑71 (MQ=255)
                     gATAACAGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTAACCGTCTCGacaa                                               <  1:2775379/71‑1 (MQ=255)
                     gATAACAGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTAACCGTCTCGacaa                                               <  1:3684780/71‑1 (MQ=255)
                     gATAACAGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTAACCGTCTCGacaa                                               <  1:1115385/71‑1 (MQ=255)
                       tAACAGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTAACCGTCTCGacaaca                                             >  1:1372004/1‑71 (MQ=255)
                       tAACAGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTAACCGTCTCGacaaca                                             >  1:2467830/1‑71 (MQ=255)
                           aGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTAACCGTCTCGACAACATTa                                          >  1:870286/1‑70 (MQ=255)
                           aGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCAGAATTCCGTAACCGTCTCGACAACAtt                                           >  1:3039718/1‑69 (MQ=255)
                           aGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCAGAATTCCGTAACCGTCTCGACAACATTa                                          >  1:1283954/1‑70 (MQ=255)
                           aGCACCGATCCGATGGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTAACCGTCTCGACAACATTAt                                         >  1:2509357/1‑71 (MQ=255)
                                  attcgatgGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGttt                                  >  1:3934865/4‑71 (MQ=255)
                                  atgcgatgGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGttt                                  >  1:309094/1‑71 (MQ=255)
                                  atgcgatgGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGtt                                   >  1:1001809/1‑70 (MQ=255)
                                   tgcgatgGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTg                                 <  1:3945602/71‑1 (MQ=255)
                                    gcgatgGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGa                                >  1:2576235/1‑71 (MQ=255)
                                    gcgatgGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGa                                >  1:2612753/1‑71 (MQ=255)
                                    gcgatgGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGa                                >  1:3594868/1‑71 (MQ=255)
                                         ggAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTAAcgg                                                        >  1:1146878/1‑40 (MQ=255)
                                          gAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGatcatc                           <  1:3879333/70‑1 (MQ=255)
                                          gAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGatcatc                           <  1:3884206/70‑1 (MQ=255)
                                          gAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGatcatc                           <  1:3107460/70‑1 (MQ=255)
                                               gaTCAAGAAGATCTTTACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGATCATCTGTCaa                     >  1:3687188/1‑71 (MQ=255)
                                               gaTCAAGAAGATCTTTACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGATCATCTGTCaa                     >  1:3973295/1‑71 (MQ=255)
                                                 tCAAGAAGATCTTTACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGATCATCTGTCAAcc                   <  1:1333867/71‑1 (MQ=255)
                                                 tCAAGAAGATCTTTACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGATCATCTGTCAAcc                   <  1:3265015/71‑1 (MQ=255)
                                                 tCAAGAAGATCTTTACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGATCATCTGTCAAcc                   <  1:1126108/71‑1 (MQ=255)
                                                  cAAGAAGATCTTTACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGATCATCTGTCAAcc                   >  1:2916861/1‑70 (MQ=255)
                                                  cAAGAAGATCTTTACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGATCATCTGTCAAcc                   >  1:2995894/1‑70 (MQ=255)
                                                  cAAGAAGATCTTTACACCAGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGatca                             >  1:1831712/1‑60 (MQ=255)
                                                  cAAGAAGATCTTTACACCAGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGatca                             >  1:3504180/1‑60 (MQ=255)
                                                  cAAGAAGATCTTTACACCAGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGatca                             >  1:2039432/1‑60 (MQ=255)
                                                   aagaagATCTTTTCACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGATCATCTGTCAACCg                  <  1:1322041/70‑1 (MQ=255)
                                                   aagaagATCTTTACACCAGAATTCCGTAACCGTCTCGacaac                                              <  1:2541953/42‑1 (MQ=255)
                                                   aagaagATCTTTACACCAGAATTCCGTAACCGTCTCGacaac                                              <  1:1047795/42‑1 (MQ=255)
                                                   aagaagATCTTTACACCAGAATTCCGTAACCGTCTCGacaac                                              <  1:2116873/42‑1 (MQ=255)
                                                   aagaagATCTTTACACCAGAATTCCGTAACCGTCTCGacaac                                              <  1:2666879/42‑1 (MQ=255)
                                                   aagaagATCTTTACACCAGAATTCCGTAACCGTCTCGacaac                                              <  1:3288200/42‑1 (MQ=255)
                                                   aagaagATCTTTACACCAGAATTCCGTAACCGTCTCGacaac                                              <  1:741881/42‑1 (MQ=255)
                                                   aagaagATCTTTACACCAGAATTCCGTAACCGTCTCGacaac                                              <  1:2952326/42‑1 (MQ=255)
                                                    agaagaTCTTTACACCAGAATTCCGTAACCGTCTCGACAACAt                                            <  1:3647648/43‑1 (MQ=255)
                                                    agaagaTCTTTACACCAGAATTCCGTAACCGTCTCGACAACAt                                            <  1:1707045/43‑1 (MQ=255)
                                                    agaagaTCTTTACACCAGAATTCCGTAACCGTCTCGACAACAt                                            <  1:281539/43‑1 (MQ=255)
                                                     gaagaTCTTTACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGATCATCTGTCAACCGACg               >  1:833878/1‑71 (MQ=255)
                                                     gaagaTCTTTACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGATCATCTGTCAACCGACg               >  1:1014673/1‑71 (MQ=255)
                                                     gaagaTCTTTACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGATCATCTGTCAACCGACg               >  1:1889555/1‑71 (MQ=255)
                                                     gaagaTCTTTACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGATCATCTGTCAACCGACg               >  1:2032407/1‑71 (MQ=255)
                                                     gaagaTCTTTACACCAGAATTCCGTAACCGTCTCGactac                                              <  1:540290/40‑1 (MQ=37)
                                                     gaagaTCTTTACACCAGAATTCCGTAACCGTCTCGACAACAt                                            <  1:3815530/42‑1 (MQ=255)
                                                          tCTTTACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGATCATCTGTCAACCGACGTGAt           <  1:1462229/70‑1 (MQ=255)
                                                          tCTTTACACCGGAATTCCGTAACCGGCTCGACAACATTATCTGGTTTGATCATCTGTCAACCGACGTGAt           <  1:3306131/70‑1 (MQ=255)
                                                           cTTTACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGATCATCTGTCAACCGACGTGATc          >  1:2890856/1‑70 (MQ=255)
                                                           cTTTACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGATCATCTGTCAACCGACGTGATc          >  1:776886/1‑70 (MQ=255)
                                                           cTTTACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGATCATCTCTCAACCGACGTGATc          >  1:1075716/1‑70 (MQ=255)
                                                              tACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGATCATCTGTCAACCGACGTGATc          <  1:2925252/67‑1 (MQ=255)
                                                              tACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGATCATCTGTCAACCGACGTGATc          <  1:841616/67‑1 (MQ=255)
                                                              tACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGATCATCTGTCAACCGACGTGATCCAt       <  1:3586503/70‑1 (MQ=255)
                                                              tACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGATCATCTGTCAACCGACGTGATCCAt       <  1:209646/70‑1 (MQ=255)
                                                              tACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGATCATCTGTCAACCGACGTGATCCAt       >  1:295653/1‑70 (MQ=255)
                                                               acacCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGATCATCTGTCAACCGACGTGATCCATc      <  1:1732472/70‑1 (MQ=255)
                                                               acacCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGATCATCTGTCAACCGACGTGATCCATc      <  1:1083591/70‑1 (MQ=255)
                                                               acacCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGATCATCTGTCAACCGACGTGATCCATc      <  1:887180/70‑1 (MQ=255)
                                                                  ccGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGATCATCTGTCAACCGACGTGATCCATCAggt  <  1:3982788/71‑1 (MQ=255)
                                                                  ccGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGATCATCTGTCAACCGACGTGATCCATCAggt  <  1:3912194/71‑1 (MQ=255)
                                                                    |                                                                    
TCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCAAGAAGATCTTTACACCGGAATTCCGTAACCGTCTCGACAACATTATCTGGTTTGATCATCTGTCAACCGACGTGATCCATCAGGT  >  minE/581808‑581944

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: