Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,493,839 A→C 16.4% intergenic (+221/‑154) gyrB → / → ibpA DNA gyrase, subunit B/heat shock chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,493,8390AC16.4% 151.4 / 36.1 82intergenic (+221/‑154)gyrB/ibpADNA gyrase, subunit B/heat shock chaperone
Reads supporting (aligned to +/- strand):  ref base A (41/27);  new base C (6/7);  total (48/34)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCCTTTTTATTTCTGCAATCCGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGATGAG  >  minE/2493771‑2493908
                                                                     |                                                                     
cGCCCTTTTTATTTCTGCAATCCGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGaa                                                                      <  1:1831793/70‑1 (MQ=255)
       tttATTTCTGCAATCCGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGaaaaaaaaa                                                               <  1:3770772/70‑1 (MQ=255)
       tttATTTCTGCAATCCGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAg                                                              >  1:426996/1‑71 (MQ=255)
       tttATTTCTGCAATCCGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAg                                                              >  1:4083432/1‑71 (MQ=255)
                 cAATCCGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCTCTGAt                                                     <  1:2842518/70‑1 (MQ=255)
                 cAATCCGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGAt                                                     >  1:3624888/1‑70 (MQ=255)
                      cGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGaaaaaaaaa                                                               >  1:1404930/1‑55 (MQ=255)
                      cGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGt                                                             >  1:2178522/1‑57 (MQ=255)
                      cGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGt                                                             >  1:3492534/1‑57 (MQ=255)
                      cGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCt                                               >  1:3551498/1‑71 (MQ=255)
                      cGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCt                                               >  1:3801443/1‑71 (MQ=255)
                        cAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGc                                                <  1:2114641/68‑1 (MQ=255)
                         aaTTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGc                                                <  1:2360834/67‑1 (MQ=255)
                          aTTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGaa                                           <  1:1235895/71‑1 (MQ=255)
                           ttATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGaaaaaaa                                                                 <  1:1786134/48‑1 (MQ=255)
                           ttATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCt                                                        <  1:738898/57‑1 (MQ=255)
                            tATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGaa                                           >  1:1267758/1‑69 (MQ=255)
                            tATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGaa                                           >  1:2683639/1‑69 (MQ=255)
                             atcatcCCGCAAATCTACTATGCAATAGACCATAAACTGaaaaaaa                                                                 <  1:1821587/46‑1 (MQ=255)
                             atcatcCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAg                                                              >  1:643664/1‑49 (MQ=255)
                             atcatcCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAg                                                              >  1:2207226/1‑49 (MQ=255)
                              tcatcCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAg                                        <  1:810045/70‑1 (MQ=255)
                              tcatcCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAg                                        <  1:4110706/70‑1 (MQ=255)
                              tcatcCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGt                                       <  1:483489/71‑1 (MQ=255)
                               catcCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATaa                                                   >  1:2696926/1‑58 (MQ=255)
                               catcCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGTCTTGAAAAGt                                       >  1:381384/1‑70 (MQ=255)
                                 tcCCGCAAATCTACTATGCAATAGACCATAAACTGaaaaaaaaa                                                               <  1:3008009/44‑1 (MQ=255)
                                 tcCCGCAAATCTACTATGCAATAGACCATAAACTGaaaaaaaaa                                                               <  1:139883/44‑1 (MQ=255)
                                 tcCCGCAAATCTACTATGCAATAGACCATAAACTGaaaaaaaa                                                                <  1:3666528/43‑1 (MQ=255)
                                 tcCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATaa                                                   >  1:273351/1‑56 (MQ=255)
                                  cccGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCAt                                   >  1:356707/1‑71 (MQ=255)
                                   ccGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCAtt                                  >  1:3212490/1‑71 (MQ=255)
                                   ccGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCAtt                                  >  1:2582823/1‑71 (MQ=255)
                                    cGCAAATCTACTATGCAATAGACCATAAACTGaaaaaaa                                                                 <  1:2015466/39‑1 (MQ=255)
                                      cAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAg                                        <  1:2860538/62‑1 (MQ=255)
                                        aaTCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTc                                >  1:3647727/1‑68 (MQ=255)
                                           agtctaTGCAATAGACCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGAc                           >  1:2027875/4‑70 (MQ=255)
                                            gtctaTGCAATAGACCATAAACTGCAAAAAAAAGTCCGCt                                                        <  1:1253055/38‑1 (MQ=39)
                                            gtctaTGCAATAGACCATAAACTGCAAAAAAAAGTCCGCt                                                        <  1:2709656/38‑1 (MQ=39)
                                            gtctaTGCAATAGACCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGt                                       >  1:454810/3‑57 (MQ=255)
                                             tctaTGCAATAGACCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGaaa                                          <  1:2446124/52‑1 (MQ=255)
                                              ctaTGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCAt                       >  1:796032/1‑71 (MQ=255)
                                              ctaTGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCAt                       <  1:1107748/71‑1 (MQ=255)
                                              ctaTGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTCAAAAGTTCATTTCCAGACCCAt                       >  1:181507/1‑71 (MQ=255)
                                               taTGCAATAGACCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTcc                               >  1:565889/1‑62 (MQ=255)
                                                 tGCAATAGACCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCAtt                      <  1:2063051/69‑1 (MQ=255)
                                                 tGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATaa                                                   >  1:3561493/1‑40 (MQ=255)
                                                 tGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCAttt                     >  1:661346/1‑70 (MQ=255)
                                                 tGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCAtt                      <  1:777433/69‑1 (MQ=255)
                                                 tGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCAtt                      <  1:2236560/69‑1 (MQ=255)
                                                 tGCAATAGACCATAAACTGCAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCAttt                     >  1:287905/1‑69 (MQ=255)
                                                  gCAATAGACCATAAACTGAAAAAAAAAGTCCGCTg                                                       <  1:432638/35‑1 (MQ=255)
                                                   cAATAGACCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTa                  <  1:1721096/71‑1 (MQ=255)
                                                   cAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGaaaa                                         >  1:964335/1‑48 (MQ=255)
                                                   cAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAg                                        <  1:395891/49‑1 (MQ=255)
                                                   cAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCAttttt                   <  1:3392764/70‑1 (MQ=255)
                                                   cAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTa                  >  1:4040170/1‑71 (MQ=255)
                                                   cAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTa                  <  1:109150/71‑1 (MQ=255)
                                                     aTAGACCATAAACTGAAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTAc                 <  1:3446644/71‑1 (MQ=255)
                                                     aTAGACCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGaaaa                                         >  1:2314872/1‑46 (MQ=255)
                                                     aTAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCAt                       <  1:709280/64‑1 (MQ=255)
                                                      tAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACa                >  1:999857/1‑70 (MQ=255)
                                                       aGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGaaaa                                         >  1:3839751/1‑44 (MQ=255)
                                                       aGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGaaaa                                         >  1:1051532/1‑44 (MQ=255)
                                                       aGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATc              >  1:1016831/1‑71 (MQ=255)
                                                       aGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATc              >  1:4013650/1‑71 (MQ=255)
                                                        gACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGaaa                                          >  1:1409729/1‑42 (MQ=255)
                                                        gACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCg             >  1:64342/1‑71 (MQ=255)
                                                        gACCATAAACTGAAAAAAAAAGTACGCTGATAAGGCTTGaaaa                                         >  1:493980/1‑43 (MQ=255)
                                                         aCCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTcc                               <  1:1052421/52‑1 (MQ=255)
                                                          ccATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGaaaa                                         >  1:1796290/1‑41 (MQ=255)
                                                           cATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAg          >  1:589389/1‑71 (MQ=255)
                                                           cATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAg          <  1:3766469/71‑1 (MQ=255)
                                                             tAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGaa                                           <  1:1407232/36‑1 (MQ=255)
                                                             tAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGcc        <  1:1406363/71‑1 (MQ=255)
                                                             tAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTc                                     <  1:4008613/42‑1 (MQ=255)
                                                              aaaCTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCAttt                     >  1:3167680/1‑57 (MQ=255)
                                                              aaaCTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGcc        >  1:3242298/1‑70 (MQ=255)
                                                              aaaCTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGcc        >  1:864856/1‑70 (MQ=255)
                                                              aaaCTGAAAAAAAAAGTCCGCTGATAAGGCTTAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGa      >  1:3399266/1‑71 (MQ=255)
                                                               aaCTTAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCg       <  1:134451/70‑1 (MQ=255)
                                                               aaCTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGc         >  1:3693697/1‑68 (MQ=255)
                                                               aaCTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGa      >  1:1380315/1‑71 (MQ=255)
                                                               aaCTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGa      >  1:915461/1‑71 (MQ=255)
                                                                aCTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGAt     <  1:2007521/71‑1 (MQ=255)
                                                                 cTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGAt     >  1:2271814/1‑70 (MQ=255)
                                                                 cTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGATg    <  1:2812071/71‑1 (MQ=255)
                                                                 cTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGATg    <  1:1403937/71‑1 (MQ=255)
                                                                 cTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGATg    <  1:3126892/71‑1 (MQ=255)
                                                                  tGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGAt     >  1:343005/1‑69 (MQ=255)
                                                                   gAAAAAAAAAGTCCGCTTATAAGGCTTGAAAAGTTCAt                                   <  1:2101784/38‑1 (MQ=255)
                                                                     aaaaaaaaaGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGATGAg  >  1:1210364/1‑70 (MQ=255)
                                                                     aaaaaaaaaGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGATGAg  <  1:3059355/70‑1 (MQ=255)
                                                                     |                                                                     
CGCCCTTTTTATTTCTGCAATCCGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGATGAG  >  minE/2493771‑2493908

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: