Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,211,331 C→G 12.8% Q7H (CAG→CAC) ‡ 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,3310CG12.8% 246.7 / 23.1 94Q7H (CAG→CAC) ‡cspCstress protein, member of the CspA‑family
Reads supporting (aligned to +/- strand):  ref base C (36/45);  new base G (4/8);  total (41/53)
Fisher's exact test for biased strand distribution p-value = 5.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

CGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGA  >  minE/1211261‑1211399
                                                                      |                                                                    
cGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAAcc                                                                      <  1:4801092/71‑1 (MQ=255)
 gAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCt                                                                     >  1:3011678/1‑71 (MQ=255)
 gAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCt                                                                     >  1:4944556/1‑71 (MQ=255)
   acacATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGa                                                                   <  1:3971008/71‑1 (MQ=255)
   acacATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGa                                                                   <  1:3624785/71‑1 (MQ=255)
   acacATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGa                                                                   >  1:2078914/1‑71 (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)
          tttGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTg                                                                    >  1:5654221/1‑63 (MQ=255)
          tttGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTaa                                                            >  1:4453267/1‑71 (MQ=255)
            tgctgcCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTaa                                                            <  1:3695205/69‑1 (MQ=255)
            tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTa                                                             <  1:2164328/68‑1 (MQ=255)
            tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTa                                                             <  1:3678572/68‑1 (MQ=255)
            tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTa                                                             <  1:3454643/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)
                      aGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCttt                                                       >  1:4610099/1‑64 (MQ=255)
                      aGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt                                                 >  1:3981886/1‑70 (MQ=255)
                        ccGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTAATCTTAGACAtcaa                                               <  1:4366361/70‑4 (MQ=255)
                        ccGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTAATCTTAGACAtcaa                                               <  1:3476175/70‑4 (MQ=255)
                        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)
                          ggAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa                                            <  1:3067404/71‑1 (MQ=255)
                            aGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt                                                 >  1:3768443/1‑64 (MQ=255)
                            aGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTAATCttag                                                      >  1:986/1‑57 (MQ=16)
                            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:172183/1‑49 (MQ=255)
                                aTGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTaa                                                            >  1:4980506/1‑49 (MQ=255)
                                aTGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTcc                                      <  1:3240707/71‑1 (MQ=255)
                                aTGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTcc                                      <  1:2504374/71‑1 (MQ=255)
                                aTGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTcc                                      >  1:3877377/1‑71 (MQ=255)
                                aTGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTcc                                      >  1:3116939/1‑71 (MQ=255)
                                aTGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTcc                                      >  1:3092134/1‑71 (MQ=255)
                                 tGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt                                                  <  1:5568382/58‑1 (MQ=255)
                                 tGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt                                                  <  1:4771719/58‑1 (MQ=255)
                                 tGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCt                                     <  1:541882/71‑1 (MQ=255)
                                   aaGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa                                          >  1:560782/1‑64 (MQ=255)
                                   aaGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa                                          >  1:2473115/1‑64 (MQ=255)
                                   aaGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa                                          >  1:2466670/1‑64 (MQ=255)
                                   aaGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa                                          >  1:1909686/1‑64 (MQ=255)
                                   aaGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa                                          >  1:4413746/1‑64 (MQ=255)
                                     gCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCt                                     >  1:2943908/1‑67 (MQ=255)
                                      ccAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTTcc                                                    <  1:2447139/51‑1 (MQ=255)
                                       cAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGAt                                >  1:2456088/1‑70 (MQ=255)
                                        aaaaCCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATc                                                          <  1:2119062/43‑1 (MQ=255)
                                        aaaaCCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGcc                                                    <  1:3424188/49‑1 (MQ=255)
                                        aaaaCCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt                                                  >  1:2056603/1‑51 (MQ=255)
                                        aaaaCCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTg                              <  1:3300760/71‑1 (MQ=255)
                                         aaaCCTTTAGACTCGTTGAACCACTTAACCTGACCTTTa                                                             <  1:4757809/39‑1 (MQ=255)
                                         aaaCCTTTAGACTCGTTGAACCACTTAACCTGACCTTTa                                                             <  1:56089/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)
                                                       gTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCg                 >  1:5641283/1‑69 (MQ=255)
                                                         tGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAgg             >  1:746258/1‑71 (MQ=255)
                                                          gAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAgg             >  1:726041/1‑70 (MQ=255)
                                                          gAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAgg             <  1:4010527/70‑1 (MQ=255)
                                                          gAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAgg             >  1:3792355/1‑70 (MQ=255)
                                                              cACTTAACCTGACCTTTAATCTTTGCCATTTGaaa                                            <  1:2944527/35‑1 (MQ=255)
                                                              cACTTAACCTGACCTTTAATCTTTGCCATTTGaaa                                            <  1:2895290/35‑1 (MQ=255)
                                                                cTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGtt                            <  1:5642193/49‑1 (MQ=255)
                                                                cTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGtt                            <  1:2023398/49‑1 (MQ=255)
                                                                cTTAACATGACCGTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACa      >  1:6902/1‑71 (MQ=255)
                                                                  tAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCg                 >  1:2381631/1‑58 (MQ=255)
                                                                   aaCCTGACCTTTAATCTTTGCCATTTGAAAAATTcc                                      >  1:4456357/1‑36 (MQ=255)
                                                                   aaCCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGtttt                          >  1:2220450/1‑48 (MQ=255)
                                                                    aCCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAg   <  1:2264444/70‑1 (MQ=255)
                                                                    aCCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGa  <  1:3641857/71‑1 (MQ=255)
                                                                      |                                                                    
CGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGA  >  minE/1211261‑1211399

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: