Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,493,839 A→C 66.7% 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,8390AC66.7% 80.8 / 40.5 57intergenic (+221/‑154)gyrB/ibpADNA gyrase, subunit B/heat shock chaperone
Reads supporting (aligned to +/- strand):  ref base A (7/12);  new base C (24/14);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 9.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

TTTTATTTCTGCAATCCGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGATGAG  >  minE/2493777‑2493908
                                                              |                                                                     
ttttATTTCTGCAATCCGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGaaaaaaa                                                                 <  1:5524626/69‑1 (MQ=255)
     tttCTGCAATCCGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAg                                                              <  1:304583/67‑1 (MQ=255)
           cAATCCGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGAt                                                     >  1:4766206/1‑70 (MQ=255)
                 gCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGaaaaaaa                                                                 <  1:92813/52‑1 (MQ=255)
                 gCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGaaaaaaa                                                                 <  1:3406183/52‑1 (MQ=255)
                 gCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGaaaaaaa                                                                 <  1:4902129/52‑1 (MQ=255)
                  cAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCtt                                              >  1:3364734/1‑70 (MQ=255)
                  cAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCtt                                              >  1:746169/1‑70 (MQ=255)
                   aaTTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGaaaaaaa                                                                 <  1:5593102/50‑1 (MQ=255)
                   aaTTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGaaaaaaa                                                                 <  1:5646675/50‑1 (MQ=255)
                     ttATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGaaaaaaaaa                                                               >  1:3743267/1‑50 (MQ=255)
                     ttATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGaaaaaaaaa                                                               >  1:4246022/1‑50 (MQ=255)
                     ttATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGaaaaaaa                                                                 <  1:4793362/48‑1 (MQ=255)
                     ttATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGaaa                                          <  1:5258143/71‑1 (MQ=255)
                       atcatcCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAg                                        >  1:2047740/1‑71 (MQ=255)
                       atcatcCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGATGATAAGGCTTGAAAAg                                        >  1:403524/1‑71 (MQ=255)
                               gCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTc                                <  1:3809274/71‑1 (MQ=255)
                                   aTCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGAt                                                     <  1:4985623/46‑1 (MQ=255)
                                      tactaTGCAATAGACCATAAACTGAAAAAAAAAGTCCGCt                                                        <  1:3912234/40‑1 (MQ=255)
                                       tctaTGCAATAGACCATAAACTGCAAAAAAAAGTCCGCTGAt                                                     <  1:755670/41‑1 (MQ=39)
                                       tctaTGCAATAGACCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTg                                             <  1:4627031/49‑1 (MQ=255)
                                        ctaTGCAATAGACCATAAACTGAAAAAAAAAGTcc                                                           >  1:2019454/1‑35 (MQ=255)
                                        ctaTCCAATAGACCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCAtt                                  <  1:2535585/60‑1 (MQ=255)
                                         taTGCAATAGACCATAAACTGCAAAAAAAAGTCCGCt                                                        <  1:2071256/37‑1 (MQ=255)
                                         taTGCAATAGACCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGaaaa                                         >  1:940286/1‑52 (MQ=255)
                                         taTGCAATAGACCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGaaaa                                         >  1:4678241/1‑52 (MQ=255)
                                         taTGCAATAGACCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGaaaa                                         >  1:1630552/1‑52 (MQ=255)
                                         taTGCAATAGACCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCAtt                      >  1:1796049/1‑71 (MQ=255)
                                           tGCAATAGACCATAAACTGCAAAAAAAAGTCCGCTGATAAGGc                                                >  1:4875959/1‑43 (MQ=255)
                                           tGCAATAGACCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGaaa                                          >  1:2337316/1‑49 (MQ=255)
                                           tGCAATAGACCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCAttt                                 >  1:2201665/1‑58 (MQ=255)
                                           tGCAATAGACCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCAttt                                 >  1:2624954/1‑58 (MQ=255)
                                           tGCAATAGACCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCAtttt                    >  1:2403760/1‑71 (MQ=255)
                                           tGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCAtttt                    <  1:1637716/71‑1 (MQ=255)
                                           tGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCAtttt                    <  1:2475266/71‑1 (MQ=255)
                                            gCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCAtttt                    <  1:123907/70‑1 (MQ=255)
                                                tAGACCATAAACTGCAAAAAAAAGTGCGCTGATAAGGCTTGAAAAGTTCATTTcc                               >  1:4811540/1‑55 (MQ=255)
                                                 aGACCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGaaaa                                         >  1:655484/1‑44 (MQ=255)
                                                 aGACCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAg                                        >  1:3376410/1‑45 (MQ=255)
                                                 aGACCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGt                                       >  1:899895/1‑46 (MQ=255)
                                                  gACCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACAt               >  1:2698320/1‑69 (MQ=255)
                                                  gACCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCg             >  1:4049814/1‑71 (MQ=255)
                                                  gACCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCg             <  1:351272/71‑1 (MQ=255)
                                                  gACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTc                                     <  1:3126821/47‑1 (MQ=255)
                                                   aCCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTc                                     <  1:2879620/46‑1 (MQ=255)
                                                    ccATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTcc                               <  1:4765338/51‑1 (MQ=255)
                                                     cATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAg          >  1:2865444/1‑71 (MQ=255)
                                                      aTAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCAt                                   <  1:5141393/45‑1 (MQ=255)
                                                      aTAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTc                                <  1:2597/48‑1 (MQ=255)
                                                       tAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTc                                <  1:5494310/47‑1 (MQ=255)
                                                       tAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGa                            >  1:215284/1‑51 (MQ=255)
                                                       tAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGa                            >  1:544089/1‑51 (MQ=255)
                                                       tAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGcc        >  1:4767103/1‑71 (MQ=255)
                                                       tAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGcc        >  1:2445544/1‑71 (MQ=255)
                                                       tAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGcc        <  1:610029/71‑1 (MQ=255)
                                                       tAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGcc        <  1:5076355/71‑1 (MQ=255)
                                                        aaaCTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGt            >  1:5436634/1‑66 (MQ=255)
                                                        aaaCTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGt            >  1:4879446/1‑66 (MQ=255)
                                                        aaaCTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGt            >  1:1320467/1‑66 (MQ=255)
                                                        aaaCTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGcc        >  1:3485064/1‑70 (MQ=255)
                                                         aaCTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGa      >  1:2564465/1‑71 (MQ=255)
                                                          aCTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGAt     <  1:3815779/71‑1 (MQ=255)
                                                           cTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTc                                     <  1:1344305/38‑1 (MQ=255)
                                                           cTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTc                                     <  1:86518/38‑1 (MQ=255)
                                                           cTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGATg    <  1:5609190/71‑1 (MQ=255)
                                                           cTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGATg    <  1:2717303/71‑1 (MQ=255)
                                                              aaaaaaaaaGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGATGAg  <  1:5433787/70‑1 (MQ=255)
                                                              aaaaaaaaaGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGATGAg  >  1:3093515/1‑70 (MQ=255)
                                                              |                                                                     
TTTTATTTCTGCAATCCGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGATGAG  >  minE/2493777‑2493908

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