Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,493,676 T→C 5.6% intergenic (+58/‑317) 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,6760TC 5.6% 205.7 / 3.5 71intergenic (+58/‑317)gyrB/ibpADNA gyrase, subunit B/heat shock chaperone
Reads supporting (aligned to +/- strand):  ref base T (27/40);  new base C (2/2);  total (29/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

GATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAGC  >  minE/2493610‑2493744
                                                                  |                                                                    
gATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGc                                                                   <  1:541102/70‑1 (MQ=255)
  tatTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCt                                                                      >  1:824598/1‑65 (MQ=255)
   atTTAATGGCGCTAACCATGCGAGCATGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATg                                                               >  1:2650945/1‑71 (MQ=255)
    tttAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTg                                                                    <  1:1124758/65‑1 (MQ=255)
    tttAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTg                                                                    <  1:2866108/65‑1 (MQ=255)
    tttAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATg                                                               <  1:2886743/70‑1 (MQ=255)
      tAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTa                                                             <  1:708616/70‑1 (MQ=255)
      tAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTAt                                                            >  1:3550996/1‑71 (MQ=255)
        aTGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCt                                                                      <  1:2867088/59‑1 (MQ=255)
              cTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAAtttt                                                    <  1:1960570/71‑1 (MQ=255)
                 aCCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTc                                                   <  1:1710278/69‑1 (MQ=255)
                   cATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAAt                                                                <  1:1362020/54‑1 (MQ=255)
                    aTGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGAt                                              <  1:988532/71‑1 (MQ=255)
                     tGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTg                                                          <  1:3020651/58‑1 (MQ=255)
                        gAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTg                                          <  1:1095152/71‑1 (MQ=255)
                        gAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTg                                          <  1:3783501/71‑1 (MQ=255)
                          gCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTg                                          <  1:72727/69‑1 (MQ=255)
                           aatGCCTGATGCGCTCCGCTTATCAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTAg                                       >  1:3604071/3‑71 (MQ=14)
                              gCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGcc                                    >  1:487445/1‑71 (MQ=255)
                               ccTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCg                                   <  1:254587/71‑1 (MQ=255)
                                 tGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGa                                 >  1:2416568/1‑71 (MQ=255)
                                  gATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCaa                                                                 <  1:3937020/38‑1 (MQ=255)
                                   aTGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGAt                                <  1:2665850/70‑1 (MQ=255)
                                   aTGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGAt                                <  1:2749763/70‑1 (MQ=255)
                                   aTGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGAt                                <  1:727083/70‑1 (MQ=255)
                                   aTGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATa                               <  1:880146/71‑1 (MQ=255)
                                   aTGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGCAGGCCGGAt                                <  1:3768838/70‑1 (MQ=255)
                                    tGCGCTACGCTTATCAGGCCTACGAAAATCTTGCaa                                                                 >  1:2286396/1‑36 (MQ=255)
                                    tGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGAt                                >  1:3493317/1‑69 (MQ=255)
                                     gcgcTCCGCTTATCAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGAtt                                             >  1:2329515/1‑55 (MQ=14)
                                      cgcTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAg                             <  1:3378400/70‑1 (MQ=255)
                                      cgcTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAg                             <  1:398353/70‑1 (MQ=255)
                                      cgcTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAg                             <  1:2201592/70‑1 (MQ=255)
                                        cTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGa                                               <  1:1335706/50‑1 (MQ=255)
                                        cTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTg                                          >  1:3079603/1‑55 (MQ=255)
                                        cTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTg                                          >  1:3047426/1‑55 (MQ=255)
                                          aCGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGAttt                                            >  1:1628117/1‑51 (MQ=255)
                                          aCGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGAt                                              <  1:3073025/49‑1 (MQ=255)
                                           cGCTTATCAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTAGGCCGGATAAGGCGtt                        <  1:667801/70‑1 (MQ=14)
                                              ttATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGc                                     <  1:3003743/54‑1 (MQ=255)
                                                aTCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTc                                                   <  1:400982/38‑1 (MQ=255)
                                                 tCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCa                      >  1:3205847/1‑66 (MQ=255)
                                                 tCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCa                      >  1:1305930/1‑66 (MQ=255)
                                                  cAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTAGGCCGGATAAGGCGtt                        <  1:1900159/63‑1 (MQ=14)
                                                   aGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGAtt                                             >  1:1809563/1‑41 (MQ=255)
                                                   aGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTg                                          <  1:4085840/44‑1 (MQ=255)
                                                     gCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTa                                        >  1:3320923/1‑44 (MQ=255)
                                                       cTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCg           >  1:3170884/1‑71 (MQ=255)
                                                        tACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTg                                          <  1:2933856/39‑1 (MQ=255)
                                                        tACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCg                                   <  1:64223/46‑1 (MQ=255)
                                                        tACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGAt                                <  1:1630983/49‑1 (MQ=255)
                                                        tACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATaa                              <  1:927463/51‑1 (MQ=255)
                                                        tACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGt                         <  1:2464765/56‑1 (MQ=255)
                                                         aCGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCg                                   <  1:2654890/45‑1 (MQ=255)
                                                         aCGAAAATCTTGAAATGTATTGAATTTTCATGAttt                                            >  1:3405042/1‑36 (MQ=255)
                                                          cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCg                                   <  1:1911285/44‑1 (MQ=255)
                                                          cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc         >  1:2445799/1‑70 (MQ=255)
                                                          cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc         >  1:2406494/1‑70 (MQ=255)
                                                          cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc         >  1:3665950/1‑70 (MQ=255)
                                                          cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc         <  1:1365970/70‑1 (MQ=255)
                                                          cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc         <  1:3234391/70‑1 (MQ=255)
                                                           gAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc         >  1:2276484/1‑69 (MQ=255)
                                                           gAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc         >  1:2897441/1‑69 (MQ=255)
                                                           gAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc         >  1:3433086/1‑69 (MQ=255)
                                                            aaaaTCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATaa                              >  1:4049445/1‑47 (MQ=255)
                                                            aaaaTCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc         >  1:1692379/1‑68 (MQ=255)
                                                             aaaTCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaaga     <  1:965684/71‑1 (MQ=255)
                                                             aaaTCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaaga     <  1:3028915/71‑1 (MQ=255)
                                                             aaaTCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaaga     <  1:3023940/71‑1 (MQ=255)
                                                             aaaTCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaag      >  1:1688900/1‑70 (MQ=255)
                                                             aaaTCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaag      >  1:3892193/1‑70 (MQ=255)
                                                                tCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAGc  >  1:1268453/1‑71 (MQ=255)
                                                                tCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAGc  >  1:2905031/1‑71 (MQ=255)
                                                                  ttGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGa                                 >  1:2690652/1‑38 (MQ=255)
                                                                  ttGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCAcgcc                  >  1:3996887/1‑53 (MQ=255)
                                                                  ttGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATc             >  1:1854759/1‑58 (MQ=255)
                                                                  |                                                                    
GATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAGC  >  minE/2493610‑2493744

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: