Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,332,248 G→C 32.2% V37V (GTC→GTG fdhD ← formate dehydrogenase formation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,332,2480GC32.2% 122.6 / 82.8 112V37V (GTC→GTGfdhDformate dehydrogenase formation protein
Reads supporting (aligned to +/- strand):  ref base G (39/36);  new base C (15/21);  total (54/58)
Fisher's exact test for biased strand distribution p-value = 3.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 120 total aligned reads displayed.

GGCGACGCCATCATCACCACATGCGAAATGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCATAACTCAATTTGACG  >  minE/2332183‑2332304
                                                                 |                                                        
ggCGACGCCATCATCACCACATGCGAAATGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCcgcg                                                          <  1:1055802/66‑1 (MQ=255)
ggCGACGCCATCATCACCACATGCGAAATGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCcgcg                                                          <  1:3203175/66‑1 (MQ=255)
ggCGACGCCATCATCACCACATGCGAAATGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCcgcg                                                          <  1:3175849/66‑1 (MQ=255)
ggCGACGCCATCATCACCACATGCGAAATGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCcgcg                                                          <  1:3104880/66‑1 (MQ=255)
ggCGACGCCATCATCACCACATGCGAAATGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCcgcg                                                          <  1:1327199/66‑1 (MQ=255)
ggCGACGCCATCATCACCACATGCGAAATGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCcgcg                                                          <  1:2881504/66‑1 (MQ=255)
ggCGACGCCATCATCACCACATGCGAAATGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCcgcg                                                          <  1:17821/66‑1 (MQ=255)
ggCGACGCCATCATCACCACATGCGAAATGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCcgcg                                                          <  1:1725919/66‑1 (MQ=255)
      gCCATCATCACCACATGCGAAATGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCa                                                    <  1:973999/66‑1 (MQ=255)
      gCCATCATCACCACATGCGAAATGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCa                                                    <  1:949410/66‑1 (MQ=255)
      gCCATCATCACCACATGCGAAATGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCa                                                    <  1:2615626/66‑1 (MQ=255)
      gCCATCATCACCACATGCGAAATGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCa                                                    <  1:216591/66‑1 (MQ=255)
      gCCATCATCACCACATGCGAAATGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCa                                                    <  1:1847524/66‑1 (MQ=255)
      gCCATCATCACCACATGCGAAATGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCa                                                    <  1:1734413/66‑1 (MQ=255)
      gCCATCATCACCACATGCGAAATGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCATCt                                                 >  1:977555/1‑69 (MQ=255)
      gCCATCATCACCACATGCGAAATGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCATCt                                                 >  1:777597/1‑69 (MQ=255)
      gCCATCATCACCACATGCGAAATGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCATCt                                                 >  1:1078747/1‑69 (MQ=255)
       ccATCATCACCACATGCGAAATGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCATCTa                                                <  1:3035243/69‑1 (MQ=255)
       ccATCATCACCACATGCGAAATGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCATCTa                                                <  1:391185/69‑1 (MQ=255)
       ccATCATCACCACATGCGAAATGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCATCTa                                                <  1:478871/69‑1 (MQ=255)
           catcaCCACATGCGAAATGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCg                                            >  1:710252/1‑69 (MQ=255)
            atcaCCACATGCGAAATGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGt                                           >  1:3022054/1‑69 (MQ=255)
            atcaCCACATGCGAAATGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGt                                           >  1:610417/1‑69 (MQ=255)
                          aaTGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAatc                             <  1:974873/69‑1 (MQ=255)
                          aaTGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAatc                             <  1:2850164/69‑1 (MQ=255)
                          aaTGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAatc                             <  1:385291/69‑1 (MQ=255)
                            tGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAatcat                           <  1:3287809/69‑1 (MQ=255)
                            tGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAatcat                           <  1:808653/69‑1 (MQ=255)
                            tGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAatcat                           <  1:680257/69‑1 (MQ=255)
                                 ttGTAGACCAGCGCAACGGGAACTTCTTCGGCCACCTCATCTAACCGTGGGTGTTGTaa                                <  1:539990/59‑1 (MQ=255)
                                 ttGTAGACCAGCGCAACGGGAACTTCTTCCGCCACCTCATCTAACCGTGGGTGTTGTaa                                <  1:3134569/59‑1 (MQ=255)
                                 ttGTAGACCAGCGCAACGGGAACTTCTTCCGCCACCTCATCTAACCGTGGGTGTTGTaa                                <  1:3191305/59‑1 (MQ=255)
                                 ttGTAGACCAGCGCAACGGGAACTTCTTCCGCCACCTCATCTAACCGTGGGTGTTGTaa                                <  1:2268757/59‑1 (MQ=255)
                                 ttGTAGACCAGCGCAACGGGAACTTCTTCCGCCACCTCATCTAACCGTGGGTGTTGTaa                                <  1:1006601/59‑1 (MQ=255)
                                 ttGTAGACCAGCGCAACGGGAACTTCTTCCGCCACCTCATCTAACCGTGGGTGTTGTaa                                <  1:3098506/59‑1 (MQ=255)
                                 ttGTAGACCAGCGCAACGGGAACTTCTTCCGCCACCTCATCTAACCGTGGGTGTTGTaa                                <  1:2967941/59‑1 (MQ=255)
                                 ttGTAGACCAGCGCAACGGGAACTTCTTCCGCCACCTCATCTAACCGTGGGTGTTGTaa                                <  1:2658574/59‑1 (MQ=255)
                                 ttGTAGACCAGCGCAACGGGAACTTCTTCCGCCACCTCATCTAACCGTGGGTGTTGTaa                                <  1:2632941/59‑1 (MQ=255)
                                 ttGTAGACCAGCGCAACGGGAACTTCTTCCGCCACCTCATCTAACCGTGGGTGTTGTaa                                <  1:225588/59‑1 (MQ=255)
                                 ttGTAGACCAGCGCAACGGGAACTTCTTCCGCCACCTCATCTAACCGTGGGTGTTGTaa                                <  1:705871/59‑1 (MQ=255)
                                 ttGTAGACCAGCGCAACGGGAACTTCTTCCGCCACCTCATCTAACCGTGGGTGTTGTaa                                <  1:771072/59‑1 (MQ=255)
                                 ttGTAGACCAGCGCAACGGGAACTTCTTCCGCCACCTCATCTAACCGTGGGTGTTGTaa                                <  1:1983173/59‑1 (MQ=255)
                                 ttGTAGACCAGCGCAACGGGAACTTCTTCCGCCACCTCATCTAACCGTGGGTGTTGTaa                                <  1:811280/59‑1 (MQ=255)
                                 ttGTAGACCAGCGCAACGGGAACTTCTTCCGCCACCTCATCTAACCGTGGGTGTTGTaa                                <  1:1888332/59‑1 (MQ=255)
                                 ttGTAGACCAGCGCAACGGGAACTTCTTCCGCCACCTCATCTAACCGTGGGTGTTGTaa                                <  1:1884677/59‑1 (MQ=255)
                                 ttGTAGACCAGCGCAACGGGAACTTCTTCCGCCACCTCATCTAACCGTGGGTGTTGTaa                                <  1:1870171/59‑1 (MQ=255)
                                 ttGTAGACCAGCGCAACGGGAACTTCTTCCGCCACCTCATCTAACCGTGGGTGTTGTaa                                <  1:1560161/59‑1 (MQ=255)
                                   gTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAatc                             <  1:1316242/60‑1 (MQ=255)
                                   gTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAatc                             <  1:1696003/60‑1 (MQ=255)
                                   gTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAatc                             <  1:136695/60‑1 (MQ=255)
                                   gTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAatc                             <  1:1387385/60‑1 (MQ=255)
                                   gTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAatc                             <  1:2559635/60‑1 (MQ=255)
                                   gTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAatc                             <  1:2255143/60‑1 (MQ=255)
                                   gTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAatc                             <  1:1985545/60‑1 (MQ=255)
                                   gTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAATCATCAcgcc                     >  1:2517550/1‑68 (MQ=255)
                                   gTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAATCATCAcgcc                     >  1:3017409/1‑68 (MQ=255)
                                   gTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAATCATCAcgcc                     >  1:812522/1‑68 (MQ=255)
                                   gTAGACCAGCGCAACGGGAACTTCTTCCGCAACCTCATCTAACCGTGGGTGTTGTAAatc                             <  1:2008070/60‑1 (MQ=255)
                                    tAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAatc                             <  1:2613881/59‑1 (MQ=255)
                                       aCCAGCGCAACGGGAACTTCTTCCGCCACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCa                 >  1:1949827/1‑68 (MQ=255)
                                       aCCAGCGCAACGGGAACTTCTTCCGCCACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCAt                >  1:464693/1‑69 (MQ=255)
                                       aCCAGCGCAACGGGAACTTCTTCCGCCACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCAt                >  1:512428/1‑69 (MQ=255)
                                       aCCAGCGCAACGGGAACTTCTTCCGCCACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCAt                >  1:3428126/1‑69 (MQ=255)
                                       aCCAGCGCAACGGGAACTTCTTCCGCCACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCAt                >  1:2265392/1‑69 (MQ=255)
                                       aCCAGCGCAACGGGAACTTCTTCCGCCACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCAt                >  1:721774/1‑69 (MQ=255)
                                       aCCAGCGCAACGGGAACTTCTTCCGCCACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCAt                >  1:3169403/1‑69 (MQ=255)
                                       aCCAGCGCAACGGGAACTTCTTCCGCCACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCAt                >  1:2372876/1‑69 (MQ=255)
                                       aCCAGCGCAACGGGAACTTCTTCCGCCACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCAt                >  1:922837/1‑69 (MQ=255)
                                       aCCAGCGCAACGGGAACTTCTTCCGCCACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCAt                >  1:1077253/1‑69 (MQ=255)
                                       aCCAGCGCAACGGGAACTTCTTCCGCCACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCAt                >  1:2664817/1‑69 (MQ=255)
                                        ccAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAatcat                           <  1:2429751/57‑1 (MQ=255)
                                         cAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTg                                          >  1:2114564/1‑41 (MQ=255)
                                         cAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAatcatc                          >  1:2508973/1‑57 (MQ=255)
                                         cAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCATaa              >  1:1099102/1‑69 (MQ=255)
                                         cAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCATaa              >  1:1531231/1‑69 (MQ=255)
                                         cAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCATaa              >  1:689395/1‑69 (MQ=255)
                                         cAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCATaa              >  1:155075/1‑69 (MQ=255)
                                         cAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCATaa              >  1:2058070/1‑69 (MQ=255)
                                         cAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCATaa              >  1:1655956/1‑69 (MQ=255)
                                         cAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCATaa              >  1:2561000/1‑69 (MQ=255)
                                         cAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCATaa              >  1:2430316/1‑69 (MQ=255)
                                         cAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCATaa              >  1:362271/1‑69 (MQ=255)
                                         cAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCATaa              >  1:2334305/1‑69 (MQ=255)
                                         cAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCATaa              >  1:3286630/1‑69 (MQ=255)
                                         cAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCATaa              >  1:3180180/1‑69 (MQ=255)
                                         cAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCATaa              >  1:2932071/1‑69 (MQ=255)
                                                aCGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCATAACTCAAt        >  1:2478045/1‑68 (MQ=255)
                                                aCGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCATAACTCAAt        >  1:3025025/1‑68 (MQ=255)
                                                aCGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCATAACTCAAt        >  1:2360431/1‑68 (MQ=255)
                                                aCGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCATAACTCAAt        >  1:3420725/1‑68 (MQ=255)
                                                aCGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCATAACTCAAt        >  1:2865729/1‑68 (MQ=255)
                                                aCGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCATAACTCAAt        >  1:2117808/1‑68 (MQ=255)
                                                     aaCTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCATAACTCAATTTGACg  <  1:1170112/69‑1 (MQ=255)
                                                     aaCTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCATAACTCAATTTGACg  <  1:1311749/69‑1 (MQ=255)
                                                     aaCTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCATAACTCAATTTGACg  <  1:1976033/69‑1 (MQ=255)
                                                     aaCTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCATAACTCAATTTGACg  <  1:2098863/69‑1 (MQ=255)
                                                       cttcttCCGCGACCTCATCTAACCGTGGGTGTTGTAAatcat                           >  1:279852/1‑42 (MQ=255)
                                                       cttcttCCGCGACCTCATCTAACCGTGGGTGTTGTAAatcat                           >  1:2004453/1‑42 (MQ=255)
                                                       cttcttCCGCGACCTCATCTAACCGTGGGTGTTGTAAatcat                           >  1:1413409/1‑42 (MQ=255)
                                                       cttcttCCGCGACCTCATCTAACCGTGGGTGTTGTAAatcat                           >  1:2114733/1‑42 (MQ=255)
                                                                 |                                                        
GGCGACGCCATCATCACCACATGCGAAATGCCGTTGTAGACCAGCGCAACGGGAACTTCTTCCGCGACCTCATCTAACCGTGGGTGTTGTAAATCATCACGCCGCCATAACTCAATTTGACG  >  minE/2332183‑2332304

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: