Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,211,310 A→G 29.0% S14S (TCT→TCC 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,3100AG29.0% 199.7 / 89.8 124S14S (TCT→TCCcspCstress protein, member of the CspA‑family
Reads supporting (aligned to +/- strand):  ref base A (37/51);  new base G (22/14);  total (59/65)
Fisher's exact test for biased strand distribution p-value = 7.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 128 total aligned reads displayed.

CTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCT  >  minE/1211241‑1211377
                                                                     |                                                                   
cTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAg                                                                    <  1:3836516/71‑1 (MQ=255)
 tGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGa                                                                   >  1:290164/1‑71 (MQ=255)
  ggATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGAc                                                                  >  1:1426323/1‑71 (MQ=255)
   gATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGa                                                                   >  1:4172713/1‑69 (MQ=255)
       gCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTa                                                                     <  1:3655965/63‑1 (MQ=255)
        cagagaAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCAtta                                                            <  1:1167976/69‑2 (MQ=255)
        cagagaAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCAtt                                                             >  1:4202839/3‑70 (MQ=12)
        cagagaAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCAtt                                                             <  1:3159219/68‑1 (MQ=12)
         ggAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGa                                                           <  1:1613290/71‑1 (MQ=255)
         ggAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGa                                                           <  1:303111/71‑1 (MQ=255)
          gagaAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGAc                                                                  >  1:3861295/1‑63 (MQ=12)
            gaAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAAc                                                         >  1:5097356/1‑70 (MQ=255)
                  tACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTaa                                                  >  1:702603/1‑71 (MQ=255)
                  tACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTa                                                                     >  1:4843779/1‑52 (MQ=255)
                  tACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGt                                                              >  1:2812273/1‑59 (MQ=255)
                  tACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTa                                                   <  1:2181620/70‑1 (MQ=255)
                  tACGAACACATCTTTAGAGCCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCAtttta                                                  >  1:4551257/1‑69 (MQ=255)
                    cGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTaaa                                                          >  1:3629146/1‑61 (MQ=9)
                    cGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAAc                                                         >  1:2622983/1‑62 (MQ=9)
                    cGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAAc                                                         >  1:2731558/1‑62 (MQ=9)
                    cGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAAcg                                                >  1:192844/1‑70 (MQ=12)
                    cGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAAcg                                                >  1:4127071/1‑70 (MQ=12)
                    cGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAAc                                                 <  1:3081630/70‑1 (MQ=12)
                     gAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCt                                               >  1:3011678/1‑71 (MQ=255)
                       acacATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGa                                             <  1:3971008/71‑1 (MQ=255)
                       acacATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGa                                             >  1:2078914/1‑71 (MQ=255)
                        cacaTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCAc                                                      <  1:3855277/61‑1 (MQ=255)
                        cacaTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCAc                                                      <  1:476376/61‑1 (MQ=255)
                          cgtCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCt                                          >  1:3944241/3‑71 (MQ=12)
                          caTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCt                                          <  1:1152790/71‑1 (MQ=255)
                          caTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCt                                          <  1:4781149/71‑1 (MQ=255)
                           aTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCtt                                         <  1:4970200/71‑1 (MQ=255)
                            tCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAg                                                                    >  1:4588488/1‑43 (MQ=255)
                            tCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAg                                                                    >  1:1835602/1‑43 (MQ=255)
                            tCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTg                                                            <  1:5642545/51‑1 (MQ=255)
                            tCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTa                                                   <  1:4881688/60‑1 (MQ=255)
                            tCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTa                                                   <  1:5158342/60‑1 (MQ=255)
                             cTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTaa                                                  >  1:5689306/1‑60 (MQ=255)
                              tttGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTaa                                      >  1:4453267/1‑71 (MQ=255)
                                tgctgcCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTaa                                      <  1:3695205/69‑1 (MQ=255)
                                tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTa                                       <  1:3454643/68‑1 (MQ=255)
                                tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTa                                       <  1:2164328/68‑1 (MQ=255)
                                 gctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCt                                   >  1:3595011/1‑71 (MQ=255)
                                 gctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCt                                   >  1:5137044/1‑71 (MQ=255)
                                  ctgcCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTAATCtt                                  <  1:5041135/71‑1 (MQ=255)
                                  ctgcCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTAATCt                                   <  1:3743808/70‑1 (MQ=255)
                                   tgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTa                                       <  1:402581/65‑1 (MQ=255)
                                   tgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTa                                       <  1:1711364/65‑1 (MQ=255)
                                     ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGc                               >  1:1797797/1‑71 (MQ=255)
                                       aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa                             <  1:5384946/71‑1 (MQ=255)
                                        tCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCt                                   <  1:2535722/64‑1 (MQ=255)
                                        tCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCt                                   <  1:3575/64‑1 (MQ=255)
                                        tCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt                            <  1:1346808/71‑1 (MQ=255)
                                        tCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt                            >  1:5457350/1‑71 (MQ=255)
                                          ttccGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCa                                                       >  1:3968665/3‑42 (MQ=255)
                                          ttccGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCa                                                       >  1:4018360/3‑42 (MQ=255)
                                          aGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTaa                                                  >  1:5227390/1‑47 (MQ=255)
                                            ccGTAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAAc                                                 >  1:3683862/1‑46 (MQ=12)
                                            ccGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTa                                                   <  1:3480573/44‑1 (MQ=255)
                                            ccGGAGTAATGAAACCGAATCCTTTGGACTCATTaaa                                                          >  1:3854253/1‑37 (MQ=12)
                                            ccGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAAc                                                 >  1:2053871/1‑46 (MQ=16)
                                            ccGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTAATCTTAGACAtcaa                         <  1:3496766/70‑4 (MQ=255)
                                            ccGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTAATCTTAGACAtca                          <  1:3627065/69‑3 (MQ=255)
                                            ccGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTAATCTTAGACAtca                          <  1:3136418/69‑3 (MQ=255)
                                              ggAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa                             <  1:2943115/64‑1 (MQ=255)
                                              ggAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa                      <  1:3113148/71‑1 (MQ=255)
                                                aGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt                           >  1:3768443/1‑64 (MQ=255)
                                                aGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTAATCtta                                 >  1:810974/1‑57 (MQ=16)
                                                 gTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTAcc                                           >  1:3366984/1‑47 (MQ=9)
                                                  tAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCtt                                         <  1:1822953/48‑1 (MQ=255)
                                                  tAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCtt                                         <  1:4339395/48‑1 (MQ=255)
                                                  tAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa                      <  1:4087591/67‑1 (MQ=255)
                                                   aaTGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCttt                                        <  1:3264398/48‑1 (MQ=255)
                                                    aTGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTaa                                      >  1:4980506/1‑49 (MQ=255)
                                                    aTGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTcc                >  1:3092134/1‑71 (MQ=255)
                                                    aTGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTcc                <  1:3240707/71‑1 (MQ=255)
                                                    aTGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTcc                <  1:2504374/71‑1 (MQ=255)
                                                     tGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt                            <  1:4771719/58‑1 (MQ=255)
                                                     tGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt                            <  1:5568382/58‑1 (MQ=255)
                                                     tGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCt               <  1:541882/71‑1 (MQ=255)
                                                       aaGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa                    >  1:2473115/1‑64 (MQ=255)
                                                       aaGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa                    >  1:4413746/1‑64 (MQ=255)
                                                       aaGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa                    >  1:560782/1‑64 (MQ=255)
                                                       aaGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa                    >  1:2466670/1‑64 (MQ=255)
                                                         gCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCt               >  1:2943908/1‑67 (MQ=255)
                                                          ccAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTTcc                              <  1:2447139/51‑1 (MQ=255)
                                                            aaaaCCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGcc                              <  1:3424188/49‑1 (MQ=255)
                                                            aaaaCCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt                            >  1:2056603/1‑51 (MQ=255)
                                                             aaaCCTTTAGACTCGTTGAACCACTTAACCTGACCTTTa                                       <  1:56089/39‑1 (MQ=255)
                                                             aaaCCTTTAGACTCGTTGAACCACTTAACCTGACCTTTa                                       <  1:4757809/39‑1 (MQ=255)
                                                              aaCCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGtt      >  1:292420/1‑71 (MQ=255)
                                                              aaCCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGtt      <  1:4897250/71‑1 (MQ=255)
                                                              aaCCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGtt      <  1:2991673/71‑1 (MQ=255)
                                                               aCCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAAt                                     >  1:3057925/1‑39 (MQ=255)
                                                                 cTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTa             <  1:4466759/61‑1 (MQ=255)
                                                                 cTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTa             <  1:455692/61‑1 (MQ=255)
                                                                 cTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTttc   <  1:1600336/71‑1 (MQ=255)
                                                                 cTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTttc   <  1:1593635/71‑1 (MQ=255)
                                                                  tttAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTttct  >  1:492893/1‑71 (MQ=255)
                                                                   ttAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTttc   >  1:4366700/1‑69 (MQ=255)
                                                                     |                                                                   
CTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCT  >  minE/1211241‑1211377

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: