Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,494,739 G→A 5.6% P74P (CCG→CCA ibpB → heat shock chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,494,7390GA 5.6% 283.2 / 4.5 89P74P (CCG→CCAibpBheat shock chaperone
Reads supporting (aligned to +/- strand):  ref base G (27/57);  new base A (2/3);  total (29/60)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

GCAGGTTTCCGTCAGGAAGATTTAGAGATTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTT  >  minE/2494671‑2494808
                                                                    |                                                                     
gCAGGTTTCCGTCAGGAAGATTTAGAGATTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCg                                                                       >  1:1915762/1‑69 (MQ=255)
 cAGGTTTCCGTCAGGAAGATTTAGAGATTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGag                                                                    <  1:5170292/71‑1 (MQ=255)
   ggTTTCCGTCAGGAAGATTTAGAGGTTCAACTGGAAGTTACGCGCCTGAGCGTAAAAGGCACGCCGGagca                                                                  <  1:3115060/71‑1 (MQ=255)
   ggTTTCCGTCAGGAAGATTTAGAGATTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGagca                                                                  <  1:2946870/71‑1 (MQ=255)
       tCCGTCAGGAAGATTTAGAGATTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGa                                                                     >  1:1490871/1‑64 (MQ=255)
       tCCGTCAGGAAGATTTAGAGATTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGcc                                                               >  1:221048/1‑70 (MQ=255)
       tCCGTCAGGAAGATTTAGAGATTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGcc                                                               >  1:4657710/1‑70 (MQ=255)
         cGTCAGGAAGATTTAGAGATTCAACTGGAAGGTACGCGTCTGAGCGTAAAAGGCACGCCg                                                                       >  1:2388369/1‑60 (MQ=255)
         cGTCAGGAAGATTTAGAGATTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCg                                                                       >  1:3008449/1‑60 (MQ=255)
         cGTCAGGAAGATTTAGAGATTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCg                                                                       >  1:3443982/1‑60 (MQ=255)
         cGTCAGGAAGATTTAGAGATTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCg                                                                       >  1:356392/1‑60 (MQ=255)
         cGTCAGGAAGATTTAGAGATTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCaa                                                             <  1:3253664/70‑1 (MQ=255)
         cGTCAGGAAGATTTAGAGATTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCaa                                                             <  1:2005413/70‑1 (MQ=255)
             aGGAAGATTTAGAGATTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAaaga                                                         <  1:520547/70‑1 (MQ=255)
              ggAAGATTTAGAGATTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAaagaag                                                       <  1:2371914/71‑1 (MQ=255)
                 aGATTTAGAGATTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAagag                                                     <  1:2523902/70‑1 (MQ=255)
                  gATTTAGAGATTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGaa                                                   <  1:3777467/71‑1 (MQ=255)
                   aTTTAGAGATTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGaaa                                                  <  1:933364/71‑1 (MQ=255)
                    tttAGAGATTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGaaaa                                                 <  1:4155124/71‑1 (MQ=255)
                    tttAGAGATTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGaaa                                                  <  1:341550/70‑1 (MQ=255)
                     ttAGAGATTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGaaaa                                                 <  1:5200199/70‑1 (MQ=255)
                      tAGAGATTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGaaaaaa                                               >  1:1215449/1‑71 (MQ=255)
                      tAGAGATTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGaaaaaa                                               >  1:1486147/1‑71 (MQ=255)
                      tAGAGATTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGaaaaaa                                               >  1:928413/1‑71 (MQ=255)
                      tAGAGATTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGaaaaaa                                               >  1:1719016/1‑71 (MQ=255)
                      tAGAGATTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGaaaaa                                                <  1:3725606/70‑1 (MQ=255)
                       agagATTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGaaaaaa                                               <  1:5185937/70‑1 (MQ=255)
                           aTTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAaag                                                          <  1:3503307/55‑1 (MQ=255)
                           aTTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAaag                                                          <  1:3748126/55‑1 (MQ=255)
                            ttCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGagc                                                                   <  1:2331417/45‑1 (MQ=255)
                            ttCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGagc                                                                   <  1:3223198/45‑1 (MQ=255)
                            ttCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGagc                                                                   <  1:3745951/45‑1 (MQ=255)
                            ttCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGagc                                                                   <  1:3592519/45‑1 (MQ=255)
                            ttCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGagc                                                                   <  1:2018146/45‑1 (MQ=255)
                            ttCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGagc                                                                   <  1:3575128/45‑1 (MQ=255)
                            ttCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGagc                                                                   <  1:62480/45‑1 (MQ=255)
                             tCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCAGAGCAGCCAAAAGAAGAGAAAAAATGGCTg                                         >  1:3238596/1‑70 (MQ=255)
                             tCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCAGAGCAGCCAAAAGAAGAGAAAAAATGGCTg                                         >  1:743853/1‑70 (MQ=255)
                              ccaCTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCaa                                                             <  1:3067773/47‑1 (MQ=255)
                                aCTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATc                                     <  1:3376735/71‑1 (MQ=255)
                                 cTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAaag                                                          <  1:2428532/49‑1 (MQ=255)
                                 cTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAaag                                                          <  1:1306302/49‑1 (MQ=255)
                                 cTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAaag                                                          <  1:3942048/49‑1 (MQ=255)
                                   ggAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAg                                  >  1:3386047/1‑71 (MQ=255)
                                   ggAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAg                                  >  1:1067716/1‑71 (MQ=255)
                                   ggAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAg                                  >  1:4882035/1‑71 (MQ=255)
                                   ggAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAg                                  >  1:780461/1‑71 (MQ=255)
                                   ggAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAg                                  >  1:1284279/1‑71 (MQ=255)
                                         tACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCt                              <  1:249987/69‑1 (MQ=255)
                                         tACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCt                              <  1:243873/69‑1 (MQ=255)
                                           cgcgCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATg                          <  1:3110563/71‑1 (MQ=255)
                                           cgcgCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATg                          <  1:187063/71‑1 (MQ=255)
                                            gcgcCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGa                         <  1:4674491/71‑1 (MQ=255)
                                            gcgcCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGa                         <  1:621070/71‑1 (MQ=255)
                                             cgcCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCaa                                   >  1:280896/1‑60 (MQ=255)
                                                cTGAGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGaaaaa                                                <  1:2570243/44‑1 (MQ=255)
                                                   aGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGaa                                                   >  1:4460866/1‑38 (MQ=255)
                                                   aGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGaa                                                   >  1:279688/1‑38 (MQ=255)
                                                   aGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGaa                                                   >  1:2296435/1‑38 (MQ=255)
                                                   aGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGaa                                                   >  1:1388745/1‑38 (MQ=255)
                                                   aGCGTAAAAGGCACGCCAGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAgt                   <  1:941936/70‑2 (MQ=255)
                                                   aGCGTAAAAGGCACGCCAGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAgt                   <  1:3722277/70‑2 (MQ=255)
                                                   aGCGTAAAAGGCACGCCAGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAgt                   <  1:1117710/70‑2 (MQ=255)
                                                    gCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCa                 >  1:4141386/1‑71 (MQ=255)
                                                       tAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCAt                                      <  1:347436/47‑1 (MQ=255)
                                                       tAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCAt                                      <  1:2837461/47‑1 (MQ=255)
                                                       tAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCAt                                      <  1:4580292/47‑1 (MQ=255)
                                                        aaaaGGCACGCCGGAGCAGCCAAAAGAAGAGaaaa                                                 <  1:3210130/35‑1 (MQ=255)
                                                        aaaaGGCACGCCGGAGCAGCCAAAAGAAGAGaaaa                                                 <  1:4105561/35‑1 (MQ=255)
                                                        aaaaGGCACGCCGGAGCAGCCAAAAGAAGAGaaaa                                                 <  1:66572/35‑1 (MQ=255)
                                                        aaaaGGCACGCCGGAGCAGCCAAAAGAAGAGaaaa                                                 <  1:4037941/35‑1 (MQ=255)
                                                        aaaaGGCACGCCGGAGCAGCCAAAAGAAGAGaaaa                                                 <  1:3233645/35‑1 (MQ=255)
                                                        aaaaGGCACGCCGGAGCAGCCAAAAGAAGAGaaaa                                                 <  1:2426877/35‑1 (MQ=255)
                                                          aaGGCACGCCGGAGCAGCCAAAAGAAGAGaaaaaa                                               >  1:1858710/1‑35 (MQ=255)
                                                          aaGGCACGCCGGAGCAGCCAAAAGAAGAGaaaaaa                                               >  1:2262646/1‑35 (MQ=255)
                                                          aaGGCACGCCGGAGCAGCCAAAAGAAGAGaaaaaa                                               >  1:2382710/1‑35 (MQ=255)
                                                          aaGGCACGCCGGAGCAGCCAAAAGAAGAGaaaaaa                                               >  1:5044150/1‑35 (MQ=255)
                                                          aaGGCACGCCGGAGCAGCCAAAAGAAGAGaaaaaa                                               >  1:502611/1‑35 (MQ=255)
                                                                cGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAg            <  1:896132/64‑1 (MQ=255)
                                                                cGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAg            <  1:453041/64‑1 (MQ=255)
                                                                cGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAg            <  1:1560397/64‑1 (MQ=255)
                                                                cGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAg            <  1:1478859/64‑1 (MQ=255)
                                                                cGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAg            <  1:2321581/64‑1 (MQ=255)
                                                                cGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAg            <  1:4168958/64‑1 (MQ=255)
                                                                 gCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAg            <  1:4363975/63‑1 (MQ=255)
                                                                 gCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCt    >  1:3471453/1‑71 (MQ=255)
                                                                 gCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCt    >  1:4435152/1‑71 (MQ=255)
                                                                  ccGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGc                   <  1:1731759/55‑1 (MQ=255)
                                                                  ccGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGc                   <  1:5147152/55‑1 (MQ=255)
                                                                  ccGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCtt   >  1:2402208/1‑71 (MQ=255)
                                                                  ccGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCtt   >  1:2808587/1‑71 (MQ=255)
                                                                   cGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCttt  <  1:2254699/71‑1 (MQ=255)
                                                                   cGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCttt  <  1:1509346/71‑1 (MQ=255)
                                                                   cGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCtt   <  1:3169558/70‑1 (MQ=255)
                                                                    |                                                                     
GCAGGTTTCCGTCAGGAAGATTTAGAGATTCAACTGGAAGGTACGCGCCTGAGCGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTT  >  minE/2494671‑2494808

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: