Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,425,817 T→G 16.5% D92A (GAC→GCC) ‡ rffG ← dTDP‑glucose 4,6‑dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,425,8170TG16.5% 155.9 / 33.7 85D92A (GAC→GCC) ‡rffGdTDP‑glucose 4,6‑dehydratase
Reads supporting (aligned to +/- strand):  ref base T (29/42);  new base G (4/10);  total (33/52)
Fisher's exact test for biased strand distribution p-value = 5.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.59e-01

AAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTGAATA  >  minE/2425749‑2425886
                                                                    |                                                                     
aaGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTc                                                                      <  1:1561731/70‑1 (MQ=255)
aaGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCCTc                                                                      >  1:4848317/1‑70 (MQ=255)
 aGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCaa                                                                    <  1:5107199/71‑1 (MQ=255)
 aGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCa                                                                     <  1:1187261/70‑1 (MQ=255)
    ccGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAg                                                                 <  1:3373636/71‑1 (MQ=255)
     cgcgCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAg                                                                 <  1:1520510/70‑1 (MQ=255)
      gcgcCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGaa                                                               >  1:2461800/1‑71 (MQ=255)
      gcgcCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGaa                                                               >  1:672439/1‑71 (MQ=255)
            cTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGt                                                                       <  1:4145212/57‑1 (MQ=255)
            cTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTc                                                          >  1:2357644/1‑70 (MQ=255)
            cTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTc                                                          >  1:4552253/1‑70 (MQ=255)
               cAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAAc                                                       <  1:4946705/70‑1 (MQ=255)
               cAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAAc                                                       <  1:2728608/70‑1 (MQ=255)
                     aTGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTc                                                                      <  1:4551543/49‑1 (MQ=255)
                      tGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCaa                                                                    <  1:3042882/50‑1 (MQ=255)
                      tGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCttt                                               >  1:2624697/1‑71 (MQ=255)
                        tataAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGt                                                                       >  1:1200551/1‑45 (MQ=255)
                        tataAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCt                                                 >  1:977172/1‑67 (MQ=255)
                        tataAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCg                                             >  1:1051877/1‑71 (MQ=255)
                        tataAGTCCCGACAATGTAGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCt                                                 >  1:2583339/1‑67 (MQ=255)
                         ataAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCgg                                            >  1:2868128/1‑71 (MQ=255)
                          taAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATg                                                    <  1:2085425/62‑1 (MQ=255)
                           aaGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCaa                                                        <  1:2888384/57‑1 (MQ=255)
                           aaGTACCAACAATATTGGTTTCAATAAATGCCGCAGGGCCTGTAATTGAACGGTCAAc                                                       <  1:827688/58‑1 (MQ=255)
                            aGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTg                                         >  1:1189486/1‑71 (MQ=255)
                            aGTACCAACAATATTGGTTTCAATAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGcag                                         >  1:3274020/1‑69 (MQ=255)
                            aGTACCAACAATATTGGTTTCAATAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGca                                          <  1:2768039/70‑2 (MQ=255)
                            aGTACCAACAATATTGGTTTCAATAAATGCCGCACGGCCTGTAATTGAACGGTCAACATGGCTTTCAGca                                          <  1:1468395/70‑2 (MQ=255)
                              tCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGc                                        >  1:2915128/1‑70 (MQ=255)
                              tCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGc                                        >  1:2165144/1‑70 (MQ=255)
                                ccaaCAATATTGGTTTCAATAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGCAGCCAg                                     <  1:4016703/68‑1 (MQ=255)
                                ccGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAg                                     <  1:1554660/71‑1 (MQ=255)
                                  gACAATGTTGGTTTCAATGAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGa                                    <  1:2008773/70‑1 (MQ=255)
                                  gACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGAt                                   >  1:2648124/1‑71 (MQ=255)
                                   aCAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGcc                                       <  1:423368/66‑1 (MQ=255)
                                   aCAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGAt                                   >  1:4480877/1‑70 (MQ=255)
                                    cAATATTGGTTTCAATAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGCAGCCAGGTGc                                 <  1:668410/71‑1 (MQ=255)
                                     aaTGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCttt                                               >  1:4805586/1‑56 (MQ=255)
                                     aaTGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCa                                >  1:2081447/1‑71 (MQ=255)
                                     aaTGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCa                                >  1:3004546/1‑71 (MQ=255)
                                     aaTGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGGCAACATGGCTTTCGGCTGCCAGATGCa                                <  1:4994585/71‑1 (MQ=255)
                                       tGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCAt                               >  1:1766674/1‑70 (MQ=255)
                                       tGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATg                              <  1:1381664/71‑1 (MQ=255)
                                       tGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATg                              <  1:1306547/71‑1 (MQ=255)
                                       tGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATg                              <  1:5304892/71‑1 (MQ=255)
                                       tGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATg                              <  1:1075914/71‑1 (MQ=255)
                                         ttGGTTTCAATAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGCAGCCAGGTGCATCac                            >  1:4258597/1‑71 (MQ=255)
                                         ttGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGc                                           >  1:199525/1‑56 (MQ=255)
                                         ttGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCa                                      <  1:2187802/61‑1 (MQ=255)
                                           ggTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCa                                >  1:5260986/1‑65 (MQ=255)
                                            gTTTCAATAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGcag                                         >  1:834593/1‑53 (MQ=255)
                                             tttCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCt                                                 <  1:2382549/46‑1 (MQ=255)
                                             tttCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCt                                                 <  1:4754659/46‑1 (MQ=255)
                                             tttCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAg                        <  1:846269/71‑1 (MQ=255)
                                             tttCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAg                        <  1:4314033/71‑1 (MQ=255)
                                             tttCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAg                        <  1:3643897/71‑1 (MQ=255)
                                               tCAATAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGCAGCCAgg                                    <  1:4306628/57‑2 (MQ=255)
                                               tCAATAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGCAGCCAgg                                    <  1:4566631/57‑2 (MQ=255)
                                               tCAATAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGCAGCCAg                                     >  1:5231748/1‑56 (MQ=255)
                                               tCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACa                                                      >  1:3626959/1‑39 (MQ=255)
                                                cAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCt                     >  1:4457491/1‑71 (MQ=255)
                                                 aaTAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCt                     <  1:3665213/70‑1 (MQ=255)
                                                 aaTAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCt                     <  1:3736104/70‑1 (MQ=255)
                                                 aaTAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCt                     <  1:4878499/70‑1 (MQ=255)
                                                 aaTAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCt                     <  1:2843927/70‑1 (MQ=255)
                                                 aaTAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTg                    >  1:4198487/1‑71 (MQ=255)
                                                  aTAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGcag                                         <  1:3045889/49‑3 (MQ=255)
                                                  aTAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTgg                   <  1:630480/71‑1 (MQ=255)
                                                   tAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGCAGCCAggt                                   <  1:4471808/54‑3 (MQ=12)
                                                   tAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGCAGCCAggt                                   <  1:2181131/54‑3 (MQ=12)
                                                    aaaCGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCa                                      >  1:1235619/1‑50 (MQ=255)
                                                    aaaCGCTGCCGGGCCGTCAATAGAACGGTAAACATGGCTTTCGGCTGCCa                                      >  1:5329179/1‑50 (MQ=255)
                                                     aaCGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCt                 <  1:1390005/70‑1 (MQ=255)
                                                      aCGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGa               <  1:2618622/71‑1 (MQ=255)
                                                       cGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGAt              <  1:1752427/71‑1 (MQ=255)
                                                       cGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGAt              >  1:2204662/1‑71 (MQ=255)
                                                          tGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCAt                               <  1:2917800/51‑1 (MQ=255)
                                                          tGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCAt                               <  1:2851498/51‑1 (MQ=255)
                                                              gggCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGt       <  1:4754421/71‑1 (MQ=255)
                                                              gggCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGGCTGGCTGAt              <  1:1568776/64‑1 (MQ=255)
                                                               ggCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCt                     <  1:2145930/56‑1 (MQ=255)
                                                               ggCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTg      <  1:3978996/71‑1 (MQ=255)
                                                               ggCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTg      <  1:12693/71‑1 (MQ=255)
                                                                gCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAg        >  1:2883026/1‑68 (MQ=255)
                                                                   gTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTc          <  1:2429509/63‑1 (MQ=255)
                                                                   gTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTc          <  1:1129379/63‑1 (MQ=255)
                                                                   gTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTGAATa  >  1:3336843/1‑71 (MQ=255)
                                                                    tCAATAGTACGGTCAACATGGCTTTCGGCTGCCAGATGCATg                              <  1:2376346/42‑1 (MQ=255)
                                                                    tCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGt       <  1:2928154/65‑1 (MQ=255)
                                                                    tCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTg      >  1:2307300/1‑66 (MQ=255)
                                                                    tCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTGa     >  1:1388735/1‑67 (MQ=255)
                                                                    |                                                                     
AAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTGAATA  >  minE/2425749‑2425886

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: