Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,889,569 A→G 9.1% R220R (CGT→CGC gudX ← predicted glucarate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,889,5690AG 9.1% 214.1 / 10.3 77R220R (CGT→CGCgudXpredicted glucarate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (26/44);  new base G (3/4);  total (29/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

TCATCAAGCAGCCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAGTTTGAAATCTTTAAAGCCGTAGC  >  minE/1889501‑1889639
                                                                    |                                                                      
tcatcaAGCAGCCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCAc                                                                       >  1:1758505/1‑70 (MQ=255)
 catcaAGCAGCCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCAc                                                                       <  1:2290959/69‑1 (MQ=255)
 catcaAGCAGCCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCAc                                                                       <  1:3461961/69‑1 (MQ=255)
 catcaAGCAGCCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACg                                                                      >  1:950165/1‑70 (MQ=255)
  atcaAGCAGCCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGaa                                                                    >  1:1733211/1‑71 (MQ=255)
     aaGCAGCCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACa                                                                  >  1:530792/1‑70 (MQ=255)
     aaGCAGCCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAg                                                                 >  1:1648406/1‑71 (MQ=255)
       gcagcCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAgtg                                                               >  1:3132193/1‑71 (MQ=255)
       gcagcCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAgtg                                                               <  1:3711579/71‑1 (MQ=255)
       gcagcCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAgtg                                                               <  1:3361247/71‑1 (MQ=255)
        cagcCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAgtgt                                                              <  1:2647745/71‑1 (MQ=255)
        cagcCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAgtgt                                                              <  1:2307342/71‑1 (MQ=255)
              tGCACCGTTGGTATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGAttt                                                        <  1:1846975/71‑1 (MQ=255)
              tGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCa                                                                        >  1:1369443/1‑55 (MQ=255)
              tGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGAttt                                                        >  1:3212975/1‑71 (MQ=255)
               gCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCGCGAACAGTGTCGATTTc                                                       >  1:1426352/1‑71 (MQ=255)
                  ccGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCGCGAACAGTGTCGATTTCTTg                                                    <  1:1146290/71‑1 (MQ=255)
                  ccGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCGCGAACAGTGTCGATTTCTTg                                                    <  1:334975/71‑1 (MQ=255)
                   cGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGa                                                           >  1:3863885/1‑63 (MQ=255)
                   cGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGa                                                           >  1:847229/1‑63 (MQ=255)
                   cGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGATTTCTTGc                                                   >  1:157702/1‑71 (MQ=255)
                     ttGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGaa                                                                    <  1:1558663/52‑1 (MQ=255)
                        ggATCAACGGTAATCCGCGCATCCGTGAAGCGTTTCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGcc                                              <  1:658858/71‑1 (MQ=255)
                        ggATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGcc                                              <  1:611757/71‑1 (MQ=255)
                        ggATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGcc                                              <  1:2271177/71‑1 (MQ=255)
                         gATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGc                                               >  1:3170060/1‑69 (MQ=255)
                          aTCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCGCGAACAGTGTCGATTTc                                                       <  1:1780179/60‑1 (MQ=255)
                           tCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACa                                                                  >  1:3093609/1‑48 (MQ=255)
                           tCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGCCAgg                                           >  1:2614518/1‑71 (MQ=255)
                             aaCGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCGCGAACAGTGTCGATTTCTTGCTCGcc                                              >  1:827065/1‑66 (MQ=255)
                             aaCGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGCCAg                                            <  1:1364481/68‑1 (MQ=255)
                               cGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAAc                                                                   <  1:1068836/43‑1 (MQ=255)
                               cGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTc                                                             >  1:1611609/1‑49 (MQ=255)
                               cGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTaa                                        >  1:1425985/1‑70 (MQ=255)
                               cGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTa                                         >  1:2273621/1‑69 (MQ=255)
                                  tAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACAcg                                    >  1:2505629/1‑71 (MQ=255)
                                     tCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGAtt                                                         <  1:3041518/47‑1 (MQ=255)
                                     tCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGAtt                                                         <  1:2931885/47‑1 (MQ=255)
                                     tCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACAcgc                                   <  1:1262324/69‑1 (MQ=255)
                                     tCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACAcgc                                   <  1:2661255/69‑1 (MQ=255)
                                      ccGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACAcgccg                                 <  1:1860804/70‑1 (MQ=255)
                                        gcgcATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGATTTCTTGCTc                                                 >  1:711628/1‑52 (MQ=255)
                                             tCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGAt                                                          <  1:2085924/38‑1 (MQ=255)
                                             tCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGAt                                                          <  1:2600078/38‑1 (MQ=255)
                                             tCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTa                                         <  1:607122/55‑1 (MQ=255)
                                              ccGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGc                                               <  1:3403939/48‑1 (MQ=255)
                                                  gAAGCGTTTCTTCAATGCACGAACAGTGTCGATTTCtt                                                     >  1:597221/1‑38 (MQ=255)
                                                   aaGCGTTTCTTCAATGCACGAACAGTGTCGATTTCt                                                      <  1:2878002/36‑1 (MQ=255)
                                                   aaGCGTTTCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGt                                          <  1:1379573/48‑1 (MQ=255)
                                                     gcttttCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAGTTTTAAAt                 >  1:1755336/4‑71 (MQ=255)
                                                       gTTTCTTCAATGCGCGAACAGTTTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAGTTTGAAATc                <  1:1821258/70‑1 (MQ=255)
                                                       gTTTCTTCAATGCGCGAACAGAGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAGTTTGAAATCt               >  1:2216941/1‑71 (MQ=255)
                                                       gTTTCTTCAAAGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAGtt                       >  1:3381466/1‑63 (MQ=255)
                                                        tttCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAg                         <  1:2307460/60‑1 (MQ=255)
                                                        tttCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAg                         <  1:2996854/60‑1 (MQ=255)
                                                        tttCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAg                         <  1:520318/60‑1 (MQ=255)
                                                        tttCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAg                         <  1:1511159/60‑1 (MQ=255)
                                                        tttCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAGTTTGAAATCtt              <  1:966073/71‑1 (MQ=255)
                                                        tttCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAGTTTGAAATCtt              <  1:3642418/71‑1 (MQ=255)
                                                         ttcttcAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAGTTTg                     <  1:1641244/63‑1 (MQ=255)
                                                          tcttcAATGCACGAACATTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAGTTTGAAATCtt              >  1:1655052/1‑69 (MQ=255)
                                                          tcttcAATGCACGAACAGTGTCGATTTCTTGCTCGcc                                              >  1:1543770/1‑37 (MQ=255)
                                                          tcttcAATGCACGAACAGTGTCGATTTCTTGCTCGc                                               >  1:1095238/1‑36 (MQ=255)
                                                          tcttcAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAGTTTGAAATCttt             >  1:3467114/1‑70 (MQ=255)
                                                          tcttcAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAGTTTGAAATCttt             <  1:488714/70‑1 (MQ=255)
                                                           cttcAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAGTTTGAAATCTTTa            <  1:184319/70‑1 (MQ=255)
                                                            ttcAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAGTTTGAAATCtt              <  1:791550/67‑1 (MQ=255)
                                                              cAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAg                         <  1:3658869/54‑1 (MQ=255)
                                                              cAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAGTTTGAAATCTTTaa           <  1:1760680/68‑1 (MQ=255)
                                                              cAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAGTTTGAAATCTTTaa           <  1:1536113/68‑1 (MQ=255)
                                                              cAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAGTTTGAAATCTTTAAAg         <  1:3517427/70‑1 (MQ=255)
                                                                aTGCACGAACAGTGTCGATTTCTTGCTCGCCAGGt                                          <  1:2479639/35‑1 (MQ=255)
                                                                aTGCACGAACAGTGTCGATTTCTTGCTCGCCAGGt                                          <  1:20895/35‑1 (MQ=255)
                                                                aTGCACGAACAGTGTCGATTTCTTGCTCGCCAGGt                                          <  1:1653479/35‑1 (MQ=255)
                                                                aTGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAGTTTGAAATCTTTAAAGCCg      <  1:124940/71‑1 (MQ=255)
                                                                 tGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTaa                                        <  1:3685257/36‑1 (MQ=255)
                                                                 tGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTaa                                        <  1:3493471/36‑1 (MQ=255)
                                                                   cACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAGTTTGAAATCTTTAAAGCCGTAg   >  1:3708974/1‑71 (MQ=255)
                                                                    aCGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAGTTTGAAATCTTTAAAGCCGTAGc  <  1:2231568/71‑1 (MQ=255)
                                                                    |                                                                      
TCATCAAGCAGCCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAGTTTGAAATCTTTAAAGCCGTAGC  >  minE/1889501‑1889639

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: