Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,211,319 G→A 23.1% F11F (TTC→TTT 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,3190GA23.1% 217.2 / 56.1 108F11F (TTC→TTTcspCstress protein, member of the CspA‑family
Reads supporting (aligned to +/- strand):  ref base G (35/48);  new base A (14/11);  total (49/59)
Fisher's exact test for biased strand distribution p-value = 2.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 114 total aligned reads displayed.

GGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGG  >  minE/1211250‑1211388
                                                                     |                                                                     
ggAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGa                                                                      <  1:1613290/71‑1 (MQ=255)
ggAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGa                                                                      <  1:303111/71‑1 (MQ=255)
ggAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGa                                                                      <  1:5539536/71‑1 (MQ=255)
   gaAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAAc                                                                    >  1:5097356/1‑70 (MQ=255)
     aGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCAc                                                                 <  1:4487975/71‑1 (MQ=255)
         tACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTaa                                                             >  1:2964801/1‑71 (MQ=255)
         tACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTaa                                                             >  1:702603/1‑71 (MQ=255)
         tACGAACACATCTTTAGAGCCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCAtttta                                                             >  1:4551257/1‑69 (MQ=255)
           cGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTaaa                                                                     >  1:3629146/1‑61 (MQ=9)
           cGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAAc                                                                    >  1:2622983/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)
           cGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAAcc                                                           <  1:4801092/71‑1 (MQ=255)
            gAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCt                                                          >  1:3011678/1‑71 (MQ=255)
            gAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCt                                                          >  1:4944556/1‑71 (MQ=255)
             aaCACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAAcgtt                                                         <  1:4154567/71‑4 (MQ=12)
              acacATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGa                                                        <  1:3624785/71‑1 (MQ=255)
              acacATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGa                                                        >  1:2078914/1‑71 (MQ=255)
              acacATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGa                                                        <  1:3971008/71‑1 (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:4781149/71‑1 (MQ=255)
                 caTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCt                                                     <  1:1152790/71‑1 (MQ=255)
                  aTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCtt                                                    <  1:4970200/71‑1 (MQ=255)
                   tCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTg                                                                       <  1:5642545/51‑1 (MQ=255)
                   tCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTa                                                              <  1:5158342/60‑1 (MQ=255)
                   tCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTa                                                              <  1:4881688/60‑1 (MQ=255)
                    cTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTaa                                                             >  1:5689306/1‑60 (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:3454643/68‑1 (MQ=255)
                       tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTa                                                  <  1:3678572/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)
                          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)
                               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)
                                   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)
                                   ccGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTAATCTTAGACAtcaa                                    <  1:3476175/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)
                                       aGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTAATCtta                                            >  1:810974/1‑57 (MQ=16)
                                        gTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTAcc                                                      >  1:3366984/1‑47 (MQ=9)
                                         tAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCtt                                                    <  1:4339395/48‑1 (MQ=255)
                                         tAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa                                 <  1:4087591/67‑1 (MQ=255)
                                           aTGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTaa                                                 >  1:172183/1‑49 (MQ=255)
                                           aTGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTaa                                                 >  1:4980506/1‑49 (MQ=255)
                                           aTGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTcc                           >  1:3877377/1‑71 (MQ=255)
                                           aTGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTcc                           <  1:2504374/71‑1 (MQ=255)
                                           aTGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTcc                           <  1:3240707/71‑1 (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:1909686/1‑64 (MQ=255)
                                              aaGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa                               >  1:560782/1‑64 (MQ=255)
                                              aaGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa                               >  1:2466670/1‑64 (MQ=255)
                                              aaGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa                               >  1:2473115/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:4897250/71‑1 (MQ=255)
                                                     aaCCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGtt                 >  1:292420/1‑71 (MQ=255)
                                                     aaCCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGtt                 <  1:2991673/71‑1 (MQ=255)
                                                        cTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTa                        <  1:455692/61‑1 (MQ=255)
                                                        cTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTa                        <  1:4466759/61‑1 (MQ=255)
                                                        cTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTttc              <  1:1600336/71‑1 (MQ=255)
                                                        cTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTttc              <  1:1593635/71‑1 (MQ=255)
                                                          ttAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTttc              >  1:4366700/1‑69 (MQ=255)
                                                                    tGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAgg  >  1:746258/1‑71 (MQ=255)
                                                                     gAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAgg  >  1:3792355/1‑70 (MQ=255)
                                                                     gAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAgg  >  1:726041/1‑70 (MQ=255)
                                                                     gAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAgg  <  1:4010527/70‑1 (MQ=255)
                                                                     |                                                                     
GGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGG  >  minE/1211250‑1211388

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: