Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,425,818 C→T 16.3% D92N (GAC→AAC) ‡ rffG ← dTDP‑glucose 4,6‑dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,425,8180CT16.3% 200.8 / 26.6 86D92N (GAC→AAC) ‡rffGdTDP‑glucose 4,6‑dehydratase
Reads supporting (aligned to +/- strand):  ref base C (30/42);  new base T (4/10);  total (34/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 = 4.66e-01

AAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTGAATACG  >  minE/2425749‑2425888
                                                                     |                                                                      
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)
            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)
                        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:4566631/57‑2 (MQ=255)
                                               tCAATAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGCAGCCAgg                                      <  1:4306628/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:3736104/70‑1 (MQ=255)
                                                 aaTAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCt                       <  1:4878499/70‑1 (MQ=255)
                                                 aaTAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCt                       <  1:2843927/70‑1 (MQ=255)
                                                 aaTAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCt                       <  1:3665213/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:2181131/54‑3 (MQ=12)
                                                   tAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGCAGCCAggt                                     <  1:4471808/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:2204662/1‑71 (MQ=255)
                                                       cGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGAt                <  1:1752427/71‑1 (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:12693/71‑1 (MQ=255)
                                                               ggCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTg        <  1:3978996/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)
                                                                     cAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGATACAGt                         >  1:4380034/1‑48 (MQ=255)
                                                                     cAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGt                         >  1:2226470/1‑48 (MQ=255)
                                                                     cAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGt                         >  1:1421545/1‑48 (MQ=255)
                                                                     cAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTGAATAcg  >  1:488319/1‑71 (MQ=255)
                                                                     cAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTGAATAcg  >  1:2105939/1‑71 (MQ=255)
                                                                     cAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTGAATAcg  >  1:4225083/1‑71 (MQ=255)
                                                                     |                                                                      
AAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTGAATACG  >  minE/2425749‑2425888

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: