Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,211,247 A→T 8.3% A35A (GCT→GCA 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,2470AT 8.3% 230.4 / 11.9 85A35A (GCT→GCAcspCstress protein, member of the CspA‑family
Reads supporting (aligned to +/- strand):  ref base A (39/36);  new base T (2/5);  total (41/44)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

GCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTC  >  minE/1211181‑1211315
                                                                  |                                                                    
gCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAg                                                                     <  1:2579163/68‑1 (MQ=255)
 ccGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGga                                                                 <  1:408955/71‑1 (MQ=255)
 ccGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGga                                                                 <  1:3319079/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:3693883/1‑71 (MQ=255)
                     aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCtt                                             >  1:2726902/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:3052052/70‑1 (MQ=255)
                           cTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgctg                                        >  1:3447212/1‑70 (MQ=255)
                            tGACCTTCGTCCAGAGATTTGTAACCATCGTTCTGGATAGCAGAGAAGTGTACGaa                                                     <  1:1278779/56‑1 (MQ=255)
                            tGACCTTCGTCCAGAGATTTGTAACCATCGTTCTGGATAGCAGAGAAGTGTACGaa                                                     <  1:1301994/56‑1 (MQ=255)
                             gACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAg                                                                     >  1:1740228/1‑39 (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:3071137/1‑71 (MQ=255)
                                   cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGcc                               >  1:2272288/1‑71 (MQ=255)
                                   cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGcc                               >  1:3032835/1‑71 (MQ=255)
                                   cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGcc                               >  1:2591574/1‑71 (MQ=255)
                                   cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGc                                <  1:864049/70‑1 (MQ=255)
                                   cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGc                                <  1:1122496/70‑1 (MQ=255)
                                     gCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAgtg                                                           <  1:873024/41‑1 (MQ=255)
                                       cAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCt                                              >  1:1740751/1‑52 (MQ=255)
                                        agagTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGt                          <  1:436698/71‑1 (MQ=255)
                                        agagTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGt                          <  1:2435086/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:3281052/69‑1 (MQ=9)
                                                          ctcTGGATGGCGGTGAAATGGACGAAAACGTCTTTGCTGCCATCATCAGGAGTGATAAAGCCAAAACCttt        <  1:998392/69‑1 (MQ=9)
                                                          ctcTGGATGGCGGTGAAATGGACGAAAACGTCTTTGCTGCCATCATCAGGAGTGATAAAGCCAAAACCttt        <  1:3928632/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:1361331/1‑50 (MQ=255)
                                                               gATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTaa                        >  1:3880776/1‑50 (MQ=255)
                                                               gATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGa                     >  1:2403135/1‑53 (MQ=255)
                                                                aTAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTc  <  1:1077156/71‑1 (MQ=255)
                                                                aTAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTc  <  1:987740/71‑1 (MQ=255)
                                                                 tAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCg                              >  1:3045156/1‑42 (MQ=255)
                                                                 tAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCg                              >  1:2372705/1‑42 (MQ=255)
                                                                 tAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCg                              >  1:1498796/1‑42 (MQ=255)
                                                                  |                                                                    
GCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTC  >  minE/1211181‑1211315

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: