Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,211,250 G→A 9.4% S34S (TCC→TCT 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,2500GA 9.4% 248.0 / 13.0 85S34S (TCC→TCTcspCstress protein, member of the CspA‑family
Reads supporting (aligned to +/- strand):  ref base G (38/39);  new base A (2/6);  total (40/45)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01

CCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTG  >  minE/1211182‑1211319
                                                                    |                                                                     
ccGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGga                                                                     <  1:3319079/71‑1 (MQ=255)
ccGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGga                                                                     <  1:408955/71‑1 (MQ=255)
 cGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGgag                                                                    <  1:832819/71‑1 (MQ=255)
   tCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGgag                                                                    >  1:3493370/1‑69 (MQ=255)
     cTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAgt                                                                <  1:3163020/71‑1 (MQ=255)
      gtgaTTTCGAACTCTACGCGCTGACCTTCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAgtg                                                               >  1:591151/4‑71 (MQ=255)
        aaTTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAgtg                                                               <  1:2205189/69‑1 (MQ=255)
         aTTTCGAACTCTACGCGCTGACCTTCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCagag                                                                    <  1:3136989/63‑1 (MQ=255)
          tttCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACg                                                           <  1:2421040/71‑1 (MQ=255)
            tCGAACTCTACGCGCTGACCTTCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGa                                                          <  1:803385/70‑1 (MQ=255)
            tCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCATGGATAGCGGAGAAGTGTACGaa                                                         >  1:629314/1‑71 (MQ=255)
               aaCTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAacac                                                      <  1:3285844/71‑1 (MQ=255)
               aaCTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAaca                                                       >  1:1987656/1‑70 (MQ=255)
                  tCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAacac                                                      <  1:2279969/68‑1 (MQ=255)
                   cAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCt                                                  >  1:149205/1‑71 (MQ=255)
                    aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGaa                                                                  >  1:1486940/1‑54 (MQ=255)
                    aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAacac                                                      >  1:2877128/1‑66 (MQ=255)
                    aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAaca                                                       <  1:2092200/65‑1 (MQ=255)
                    aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCtt                                                 >  1:2726902/1‑71 (MQ=255)
                    aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCtt                                                 >  1:3693883/1‑71 (MQ=255)
                    aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCtt                                                 <  1:1291169/71‑1 (MQ=255)
                     aCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCttt                                                <  1:3964724/71‑1 (MQ=255)
                     aCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCttt                                                <  1:818930/71‑1 (MQ=255)
                       gTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTa                                                             <  1:3536567/56‑1 (MQ=255)
                       gTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgc                                              >  1:3809707/1‑71 (MQ=255)
                        ttCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgct                                             >  1:2979536/1‑71 (MQ=255)
                         tCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAacac                                                      >  1:770969/1‑61 (MQ=255)
                          cTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgctgc                                           >  1:941143/1‑71 (MQ=255)
                          cTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgctg                                            >  1:3447212/1‑70 (MQ=255)
                          cTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgctg                                            <  1:3052052/70‑1 (MQ=255)
                           tGACCTTCGTCCAGAGATTTGTAACCATCGTTCTGGATAGCAGAGAAGTGTACGaa                                                         <  1:1301994/56‑1 (MQ=255)
                           tGACCTTCGTCCAGAGATTTGTAACCATCGTTCTGGATAGCAGAGAAGTGTACGaa                                                         <  1:1278779/56‑1 (MQ=255)
                             aCCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCa                                         <  1:3725181/70‑1 (MQ=255)
                             aCCTTCAGCCAGAGTTTTAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCa                                         <  1:2341625/69‑1 (MQ=255)
                                 tCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGc                                    >  1:3223482/1‑71 (MQ=255)
                                 tCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGc                                    >  1:81415/1‑71 (MQ=255)
                                  cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATc                                                   >  1:1257206/1‑55 (MQ=255)
                                  cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgc                                              >  1:191753/1‑60 (MQ=255)
                                  cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGcc                                   >  1:2992955/1‑71 (MQ=255)
                                  cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGcc                                   >  1:3032835/1‑71 (MQ=255)
                                  cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGcc                                   >  1:3071137/1‑71 (MQ=255)
                                  cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGcc                                   >  1:2272288/1‑71 (MQ=255)
                                  cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGcc                                   >  1:2591574/1‑71 (MQ=255)
                                  cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGc                                    <  1:1122496/70‑1 (MQ=255)
                                  cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGc                                    <  1:864049/70‑1 (MQ=255)
                                    gCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAgtg                                                               <  1:873024/41‑1 (MQ=255)
                                      cAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCt                                                  >  1:1740751/1‑52 (MQ=255)
                                       agagTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGt                              <  1:2435086/71‑1 (MQ=255)
                                       agagTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGt                              <  1:436698/71‑1 (MQ=255)
                                        gagTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGt                              >  1:3538257/1‑70 (MQ=255)
                                           ttttGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATc                                                   <  1:1960922/46‑1 (MQ=255)
                                           ttttAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAacacgt                                                    <  1:477609/45‑3 (MQ=255)
                                           ttttAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCtt                                                 <  1:2600963/48‑1 (MQ=255)
                                             ttGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgct                                             >  1:2327533/1‑50 (MQ=255)
                                             ttGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgc                                              >  1:2066364/1‑49 (MQ=255)
                                             ttGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgc                                              >  1:1325835/1‑49 (MQ=255)
                                              tGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGaa                        <  1:4144428/70‑1 (MQ=255)
                                              tGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAg                       >  1:1055842/1‑71 (MQ=255)
                                               gAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCa                                      <  1:294587/55‑1 (MQ=255)
                                                      ttACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCAt                                        <  1:4052745/46‑1 (MQ=255)
                                                      ttACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGa                         <  1:3547024/61‑1 (MQ=255)
                                                       tACCCTGGATAGCGGAGAAGTGTACGAACACATCtt                                                 <  1:3053343/36‑1 (MQ=255)
                                                        aCCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgc                                              >  1:2597496/1‑38 (MQ=255)
                                                        aCCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAg                                     >  1:2499169/1‑47 (MQ=255)
                                                         gtcTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCtt             <  1:176250/68‑1 (MQ=255)
                                                         ctcTGGATGGCGGTGAAATGGACGAAAACGTCTTTGCTGCCATCATCAGGAGTGATAAAGCCAAAACCttt            <  1:998392/69‑1 (MQ=9)
                                                         ctcTGGATGGCGGTGAAATGGACGAAAACGTCTTTGCTGCCATCATCAGGAGTGATAAAGCCAAAACCttt            <  1:3928632/69‑1 (MQ=9)
                                                         ctcTGGATGGCGGTGAAATGGACGAAAACGTCTTTGCTGCCATCATCAGGAGTGATAAAGCCAAAACCttt            <  1:3281052/69‑1 (MQ=9)
                                                         cccTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCtt             <  1:519036/70‑1 (MQ=255)
                                                          ccTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCTTCAGCCGGAGTAATGAAGCCAAAACCTTTa           <  1:1024990/71‑1 (MQ=255)
                                                          ccTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTa           <  1:174806/71‑1 (MQ=255)
                                                           cTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATc                                       >  1:2265560/1‑42 (MQ=255)
                                                            tGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGt                              <  1:2012236/50‑1 (MQ=255)
                                                            tGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGa         >  1:1306103/1‑71 (MQ=255)
                                                            tGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGACTAATGAAg                       >  1:2273317/1‑57 (MQ=255)
                                                              gATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTaa                            >  1:3880776/1‑50 (MQ=255)
                                                              gATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTaa                            >  1:1361331/1‑50 (MQ=255)
                                                              gATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGa                         >  1:2403135/1‑53 (MQ=255)
                                                               aTAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTc      <  1:987740/71‑1 (MQ=255)
                                                               aTAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTc      <  1:1077156/71‑1 (MQ=255)
                                                                tAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCg                                  >  1:2372705/1‑42 (MQ=255)
                                                                tAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCg                                  >  1:3045156/1‑42 (MQ=255)
                                                                tAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCg                                  >  1:1498796/1‑42 (MQ=255)
                                                                  gCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCtt             <  1:1532013/61‑1 (MQ=255)
                                                                   cGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTg  <  1:2085230/71‑1 (MQ=255)
                                                                   cGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTg  <  1:3129450/71‑1 (MQ=255)
                                                                    ggAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCaa                   <  1:920578/53‑1 (MQ=255)
                                                                    ggAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCaa                   <  1:1566687/53‑1 (MQ=255)
                                                                    ggAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTg  <  1:1812119/70‑1 (MQ=255)
                                                                    |                                                                     
CCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTG  >  minE/1211182‑1211319

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: