Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,211,240 C→T 25.4% G38S (GGT→AGT) ‡ cspC ← stress protein, member of the CspA‑family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,211,2400CT25.4% 210.9 / 66.3 114G38S (GGT→AGT) ‡cspCstress protein, member of the CspA‑family
Reads supporting (aligned to +/- strand):  ref base C (37/48);  new base T (8/21);  total (45/69)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 127 total aligned reads displayed.

CCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTT  >  minE/1211173‑1211308
                                                                   |                                                                    
ccTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTg                                                                   <  1:773513/71‑1 (MQ=255)
ccTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTg                                                                   <  1:575558/71‑1 (MQ=255)
    tCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGa                                                                 >  1:626907/1‑69 (MQ=255)
    tCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGa                                                                 <  1:2718993/69‑1 (MQ=255)
    tCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGa                                                                 >  1:2660368/1‑69 (MQ=255)
        gCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCATCCAGAGTTTTGAAGCCATTACCCTGGATAGc                                                             <  1:4264330/69‑1 (MQ=255)
        gCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCg                                                            <  1:1299178/70‑1 (MQ=255)
          cGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGgag                                                         <  1:3148808/71‑1 (MQ=255)
          cGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGgag                                                         <  1:1933932/71‑1 (MQ=255)
           gTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGa                                                                 >  1:63875/1‑62 (MQ=255)
           gTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGgag                                                         >  1:3579978/1‑70 (MQ=255)
           gTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATAACCCTGGa                                                                 >  1:321761/1‑62 (MQ=255)
            tCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTAccc                                                                     <  1:2603352/57‑1 (MQ=255)
            tCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGgag                                                         >  1:854272/1‑69 (MQ=255)
            tCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGga                                                          >  1:1729316/1‑68 (MQ=255)
               tGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAgtg                                                    >  1:57372/1‑71 (MQ=255)
               gtgaTTTCGAACTCTACGCGCTGACCTTCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAgtg                                                    <  1:3133154/68‑1 (MQ=255)
                 aaTTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAgtg                                                    <  1:5304883/69‑1 (MQ=255)
                 aaTTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTa                                                  <  1:1780813/71‑1 (MQ=255)
                 aaTTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTa                                                  >  1:5393203/1‑71 (MQ=255)
                 aaTTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTa                                                  <  1:5534542/71‑1 (MQ=255)
                  aTTTCGAACTCTACGCGCTGACCTTCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTAc                                                 >  1:1514981/1‑71 (MQ=255)
                  aTTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTa                                                  <  1:5208823/70‑1 (MQ=255)
                   tttCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACg                                                <  1:1315390/71‑1 (MQ=255)
                   tttCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACg                                                <  1:394843/71‑1 (MQ=255)
                     tCGAACTCTACGCGCTGACCTTCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGa                                               <  1:2607437/70‑1 (MQ=255)
                     tCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTAcc                                                                      >  1:2481993/1‑47 (MQ=255)
                        aaCTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAgtg                                                    <  1:4910506/62‑1 (MQ=255)
                        aaCTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAaca                                            >  1:202511/1‑70 (MQ=255)
                        aaCTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAaca                                            >  1:2992507/1‑70 (MQ=255)
                         aCTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTAcc                                                                      >  1:5322828/1‑43 (MQ=255)
                         aCTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTAcc                                                                      >  1:3794890/1‑43 (MQ=255)
                             aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTgg                                                                  >  1:532067/1‑43 (MQ=255)
                             aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGa                                                                 >  1:2265028/1‑44 (MQ=255)
                             aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCtt                                      >  1:4738548/1‑71 (MQ=255)
                              aCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGgag                                                         <  1:4306100/51‑1 (MQ=255)
                               ctttCTGACCTTCGTCCAGAGATTTGTAACCATCGTTCTGGATAGCAGAGAAgt                                                     >  1:3540719/3‑54 (MQ=255)
                               cGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTAcc                                                                      >  1:2389265/1‑37 (MQ=255)
                               cGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTAcc                                                                      >  1:148853/1‑37 (MQ=255)
                               cGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTAcc                                                                      >  1:5516915/1‑37 (MQ=255)
                               cGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTAcc                                                                      >  1:1834766/1‑37 (MQ=255)
                               cGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTAcc                                                                      >  1:5339600/1‑37 (MQ=255)
                               cGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtg                                    <  1:1692410/71‑1 (MQ=255)
                                tttCTGACCTTCGTCCAGAGATTTGTAACCATCGTTCTGGATAGCAgaga                                                        <  1:3311892/49‑1 (MQ=255)
                                gTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTAcc                                                                      >  1:2607791/1‑36 (MQ=255)
                                 cgcTGACCTTCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTtgc                                   <  1:5022279/68‑1 (MQ=255)
                                  tCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgctg                                 >  1:4519688/1‑71 (MQ=255)
                                   cTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgctgc                                >  1:1254573/1‑71 (MQ=255)
                                    tGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGAt                                                                <  1:5508843/38‑1 (MQ=255)
                                    tGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgctgc                                >  1:468109/1‑70 (MQ=255)
                                     gACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCa                              >  1:5551196/1‑71 (MQ=255)
                                     gACCTTCAGCAAGAGTTTTAAAACCATTAGTCTGGa                                                                 >  1:2001823/1‑36 (MQ=16)
                                         ttCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTa                                                  <  1:3684160/47‑1 (MQ=255)
                                         ttCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTa                                                  <  1:1258506/47‑1 (MQ=255)
                                         ttCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTa                                                  <  1:138043/47‑1 (MQ=255)
                                         ttCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTa                                                  <  1:4446477/47‑1 (MQ=255)
                                          tCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGc                         >  1:3933927/1‑71 (MQ=255)
                                          tCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGc                         >  1:2276753/1‑71 (MQ=255)
                                          tCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCtt                                      <  1:1788695/58‑1 (MQ=255)
                                          tCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCtt                                      <  1:5328213/58‑1 (MQ=255)
                                          tCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCtt                                      <  1:3821948/58‑1 (MQ=255)
                                           cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGgaga                                                        >  1:1257439/1‑39 (MQ=255)
                                           cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTa                                                  >  1:574708/1‑45 (MQ=255)
                                           cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAg                          <  1:3489174/69‑1 (MQ=255)
                                           cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAg                          <  1:1528557/69‑1 (MQ=255)
                                           cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGc                         <  1:1600824/70‑1 (MQ=255)
                                            aGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCg                       <  1:2632015/71‑1 (MQ=255)
                                              ccAGAGTTTTGAAGCCATTACCCTGGATAGCGgag                                                         <  1:3270136/35‑1 (MQ=255)
                                               aagagTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAg                    >  1:2641007/2‑71 (MQ=255)
                                                agagTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGact                   <  1:2794389/71‑3 (MQ=255)
                                                agagTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAg                    >  1:5312850/1‑70 (MQ=255)
                                                  agTTTTAAATCCATTAGTCTGGATTGCAGAGAAGTGTACGAac                                             <  1:2637403/43‑1 (MQ=16)
                                                   gTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAAt                <  1:5075619/71‑1 (MQ=255)
                                                   gTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAAt                >  1:2661813/1‑71 (MQ=255)
                                                   gTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAAt                >  1:2186337/1‑71 (MQ=255)
                                                   gTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAAt                <  1:2078635/71‑1 (MQ=255)
                                                   gTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAAt                <  1:5394219/71‑1 (MQ=255)
                                                   gTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAAt                <  1:1547175/71‑1 (MQ=16)
                                                     tttAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGccg                              <  1:5545042/55‑2 (MQ=16)
                                                      ttGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGa              <  1:1533115/70‑1 (MQ=255)
                                                      ttGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGa              <  1:2205268/70‑1 (MQ=255)
                                                       tGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAg            >  1:1975074/1‑71 (MQ=255)
                                                        gAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATc                                        >  1:661413/1‑42 (MQ=255)
                                                         aaaccATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAAcc          <  1:663283/68‑1 (MQ=12)
                                                          aGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGcc          >  1:431655/1‑70 (MQ=255)
                                                            ccATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCttt                                     >  1:1058512/1‑41 (MQ=19)
                                                            ccATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCaaa       <  1:2245700/71‑1 (MQ=255)
                                                             cATTGCTCTGGATGGCGGTGAAATGGACGAAAACGTCTTTGCTGCCATCATCAGGa                     >  1:3652710/1‑56 (MQ=2)
                                                             cATTGCTCTGGATGGCGGTGAAATGGACGAAAACGTCTTTGCTGCCATCATCAGGa                     >  1:2481815/1‑56 (MQ=2)
                                                             cATTGCTCTGGATGGCGGTGAAATGGACGAAAACGTCTTTGCTGCCATCATCAGGAGTGATAAAGCCaaa       <  1:951577/70‑1 (MQ=9)
                                                             cATTGCTCTGGATGGCGGTGAAATGGACGAAAACGTCTTTGCTGCCATCATCAGGAGTGATAAAGCCaaa       <  1:5092732/70‑1 (MQ=9)
                                                             cATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGa                     >  1:2581149/1‑56 (MQ=12)
                                                              aTTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAAc     >  1:597120/1‑71 (MQ=255)
                                                              aTTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAAc     >  1:3805245/1‑71 (MQ=255)
                                                               ttACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCaa        <  1:3383326/67‑1 (MQ=255)
                                                                 aCCCTGGATAGCGGAGAAGTGTACGAACACATCtt                                      <  1:1814370/35‑1 (MQ=255)
                                                                 aCCCTGGATAGCGGAGAAGTGTACGAACACATCtt                                      <  1:2290882/35‑1 (MQ=255)
                                                                 aCCCTGGATAGCGGAGAAGTGTACGAACACATCtt                                      <  1:2382160/35‑1 (MQ=255)
                                                                  cccTGGATAGCGGAGAAGTGTACGAACACATCTTTGATGCCATCAGCCGGAGTAATGAAGCCAAAACCtt  >  1:4280542/1‑70 (MQ=255)
                                                                  cccTGGATAGCGGAGAAGTGTAAGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCtt  <  1:4050024/70‑1 (MQ=255)
                                                                   |                                                                    
CCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTT  >  minE/1211173‑1211308

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