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

CGCCCTTTTTATTTCTGCAATCCGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGATGAG  >  minE/2493771‑2493908
                                                                     |                                                                     
cGCCCTTTTTATTTCTGCAATCCGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTG‑aa                                                                      <  1:1831793/70‑1 (MQ=255)
       tttATTTCTGCAATCCGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTG‑aaaaaaaaa                                                               <  1:3770772/70‑1 (MQ=255)
       tttATTTCTGCAATCCGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAg                                                              >  1:426996/1‑71 (MQ=255)
       tttATTTCTGCAATCCGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAg                                                              >  1:4083432/1‑71 (MQ=255)
                 cAATCCGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAGTCCTCTGAt                                                     <  1:2842518/70‑1 (MQ=255)
                 cAATCCGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGAt                                                     >  1:3624888/1‑70 (MQ=255)
                      cGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTG‑aaaaaaaaa                                                               >  1:1404930/1‑55 (MQ=255)
                      cGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAGt                                                             >  1:2178522/1‑57 (MQ=255)
                      cGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAGt                                                             >  1:3492534/1‑57 (MQ=255)
                      cGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCt                                               >  1:3551498/1‑71 (MQ=255)
                      cGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCt                                               >  1:3801443/1‑71 (MQ=255)
                        cAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGc                                                <  1:2114641/68‑1 (MQ=255)
                         aaTTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGc                                                <  1:2360834/67‑1 (MQ=255)
                          aTTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGaa                                           <  1:1235895/71‑1 (MQ=255)
                           ttATCATCCCGCAAATCTACTATGCAATAGACCATAAACTG‑aaaaaaa                                                                 <  1:1786134/48‑1 (MQ=255)
                           ttATCATCCCGCAAATCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCt                                                        <  1:738898/57‑1 (MQ=255)
                            tATCATCCCGCAAATCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGaa                                           >  1:1267758/1‑69 (MQ=255)
                            tATCATCCCGCAAATCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGaa                                           >  1:2683639/1‑69 (MQ=255)
                             atcatcCCGCAAATCTACTATGCAATAGACCATAAACTG‑aaaaaaa                                                                 <  1:1821587/46‑1 (MQ=255)
                             atcatcCCGCAAATCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAg                                                              >  1:643664/1‑49 (MQ=255)
                             atcatcCCGCAAATCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAg                                                              >  1:2207226/1‑49 (MQ=255)
                              tcatcCCGCAAATCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAg                                        <  1:810045/70‑1 (MQ=255)
                              tcatcCCGCAAATCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAg                                        <  1:4110706/70‑1 (MQ=255)
                              tcatcCCGCAAATCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGt                                       <  1:483489/71‑1 (MQ=255)
                               catcCCGCAAATCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATaa                                                   >  1:2696926/1‑58 (MQ=255)
                               catcCCGCAAATCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGTCTTGAAAAGt                                       >  1:381384/1‑70 (MQ=255)
                                 tcCCGCAAATCTACTATGCAATAGACCATAAACTG‑aaaaaaaaa                                                               <  1:3008009/44‑1 (MQ=255)
                                 tcCCGCAAATCTACTATGCAATAGACCATAAACTG‑aaaaaaaaa                                                               <  1:139883/44‑1 (MQ=255)
                                 tcCCGCAAATCTACTATGCAATAGACCATAAACTG‑aaaaaaaa                                                                <  1:3666528/43‑1 (MQ=255)
                                 tcCCGCAAATCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATaa                                                   >  1:273351/1‑56 (MQ=255)
                                  cccGCAAATCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCAt                                   >  1:356707/1‑71 (MQ=255)
                                   ccGCAAATCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCAtt                                  >  1:3212490/1‑71 (MQ=255)
                                   ccGCAAATCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCAtt                                  >  1:2582823/1‑71 (MQ=255)
                                    cGCAAATCTACTATGCAATAGACCATAAACTG‑aaaaaaa                                                                 <  1:2015466/39‑1 (MQ=255)
                                      cAAATCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAg                                        <  1:2860538/62‑1 (MQ=255)
                                        aaTCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTc                                >  1:3647727/1‑68 (MQ=255)
                                           agtctaTGCAATAGACCATAAACTG‑CAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGAc                           >  1:2027875/4‑70 (MQ=255)
                                            gtctaTGCAATAGACCATAAACTG‑CAAAAAAAAGTCCGCt                                                        <  1:1253055/38‑1 (MQ=39)
                                            gtctaTGCAATAGACCATAAACTG‑CAAAAAAAAGTCCGCt                                                        <  1:2709656/38‑1 (MQ=39)
                                            gtctaTGCAATAGACCATAAACTG‑CAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGt                                       >  1:454810/3‑57 (MQ=255)
                                             tctaTGCAATAGACCATAAACTG‑CAAAAAAAAGTCCGCTGATAAGGCTTGaaa                                          <  1:2446124/52‑1 (MQ=255)
                                              ctaTGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCAt                       >  1:796032/1‑71 (MQ=255)
                                              ctaTGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCAt                       <  1:1107748/71‑1 (MQ=255)
                                              ctaTGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTCAAAAGTTCATTTCCAGACCCAt                       >  1:181507/1‑71 (MQ=255)
                                               taTGCAATAGACCATAAACTG‑CAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTcc                               >  1:565889/1‑62 (MQ=255)
                                                 tGCAATAGACCATAAACTG‑CAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCAtt                      <  1:2063051/69‑1 (MQ=255)
                                                 tGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATaa                                                   >  1:3561493/1‑40 (MQ=255)
                                                 tGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCAttt                     >  1:661346/1‑70 (MQ=255)
                                                 tGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCAtt                      <  1:777433/69‑1 (MQ=255)
                                                 tGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCAtt                      <  1:2236560/69‑1 (MQ=255)
                                                 tGCAATAGACCATAAACTG‑‑CAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCAttt                     >  1:287905/1‑69 (MQ=255)
                                                  gCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTg                                                       <  1:432638/35‑1 (MQ=255)
                                                   cAATAGACCATAAACTG‑CAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTa                  <  1:1721096/71‑1 (MQ=255)
                                                   cAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGaaaa                                         >  1:964335/1‑48 (MQ=255)
                                                   cAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAg                                        <  1:395891/49‑1 (MQ=255)
                                                   cAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCAttttt                   <  1:3392764/70‑1 (MQ=255)
                                                   cAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTa                  >  1:4040170/1‑71 (MQ=255)
                                                   cAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTa                  <  1:109150/71‑1 (MQ=255)
                                                     aTAGACCATAAACTGAAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTAc                 <  1:3446644/71‑1 (MQ=255)
                                                     aTAGACCATAAACTG‑CAAAAAAAAGTCCGCTGATAAGGCTTGaaaa                                         >  1:2314872/1‑46 (MQ=255)
                                                     aTAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCAt                       <  1:709280/64‑1 (MQ=255)
                                                      tAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACa                >  1:999857/1‑70 (MQ=255)
                                                       aGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGaaaa                                         >  1:3839751/1‑44 (MQ=255)
                                                       aGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGaaaa                                         >  1:1051532/1‑44 (MQ=255)
                                                       aGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATc              >  1:1016831/1‑71 (MQ=255)
                                                       aGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATc              >  1:4013650/1‑71 (MQ=255)
                                                        gACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGaaa                                          >  1:1409729/1‑42 (MQ=255)
                                                        gACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCg             >  1:64342/1‑71 (MQ=255)
                                                        gACCATAAACTG‑AAAAAAAAAGTACGCTGATAAGGCTTGaaaa                                         >  1:493980/1‑43 (MQ=255)
                                                         aCCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTcc                               <  1:1052421/52‑1 (MQ=255)
                                                          ccATAAACTG‑CAAAAAAAAGTCCGCTGATAAGGCTTGaaaa                                         >  1:1796290/1‑41 (MQ=255)
                                                           cATAAACTG‑CAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAg          >  1:589389/1‑71 (MQ=255)
                                                           cATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAg          <  1:3766469/71‑1 (MQ=255)
                                                             tAAACTG‑CAAAAAAAAGTCCGCTGATAAGGCTTGaa                                           <  1:1407232/36‑1 (MQ=255)
                                                             tAAACTG‑CAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGcc        <  1:1406363/71‑1 (MQ=255)
                                                             tAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTc                                     <  1:4008613/42‑1 (MQ=255)
                                                              aaaCTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCAttt                     >  1:3167680/1‑57 (MQ=255)
                                                              aaaCTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGcc        >  1:3242298/1‑70 (MQ=255)
                                                              aaaCTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGcc        >  1:864856/1‑70 (MQ=255)
                                                              aaaCTG‑AAAAAAAAAGTCCGCTGATAAGGCTT‑AAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGa      >  1:3399266/1‑71 (MQ=255)
                                                               aaCTT‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCg       <  1:134451/70‑1 (MQ=255)
                                                               aaCTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGc         >  1:3693697/1‑68 (MQ=255)
                                                               aaCTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGa      >  1:1380315/1‑71 (MQ=255)
                                                               aaCTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGa      >  1:915461/1‑71 (MQ=255)
                                                                aCTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGAt     <  1:2007521/71‑1 (MQ=255)
                                                                 cTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGAt     >  1:2271814/1‑70 (MQ=255)
                                                                 cTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGATg    <  1:2812071/71‑1 (MQ=255)
                                                                 cTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGATg    <  1:1403937/71‑1 (MQ=255)
                                                                 cTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGATg    <  1:3126892/71‑1 (MQ=255)
                                                                  tG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGAt     >  1:343005/1‑69 (MQ=255)
                                                                   g‑AAAAAAAAAGTCCGCTTATAAGGCTTGAAAAGTTCAt                                   <  1:2101784/38‑1 (MQ=255)
                                                                     aaaaaaaaaGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGATGAg  >  1:1210364/1‑70 (MQ=255)
                                                                     aaaaaaaaaGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGATGAg  <  1:3059355/70‑1 (MQ=255)
                                                                     |                                                                     
CGCCCTTTTTATTTCTGCAATCCGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTG‑AAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGATGAG  >  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: ‑