Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,826,741 G→A 7.7% intergenic (‑418/‑17) hydN ← / → mutS formate dehydrogenase‑H, [4Fe‑4S] ferredoxin subunit/methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,826,7410GA 7.7% 210.1 / 6.5 65intergenic (‑418/‑17)hydN/mutSformate dehydrogenase‑H, [4Fe‑4S] ferredoxin subunit/methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (28/32);  new base A (3/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

GTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTATCT  >  minE/1826676‑1826810
                                                                 |                                                                     
gTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCg                                                                  <  1:2078746/71‑1 (MQ=255)
  aTCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGa                                                                >  1:1737176/1‑71 (MQ=255)
     aTGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAAc                                                                    <  1:885333/64‑1 (MQ=255)
        tgCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAAc                                                          >  1:1231343/1‑71 (MQ=255)
          cgcCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCC                                                        <  1:1067002‑M1/71‑1 (MQ=255)
          cgcCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCC                                                        <  1:541502‑M1/71‑1 (MQ=255)
               tATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCAt                                                     >  1:2942566/1‑69 (MQ=255)
                aTGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACa                                                              >  1:2801639/1‑59 (MQ=255)
                aTGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCAt                                                     >  1:3684944/1‑68 (MQ=255)
                aTGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATg                                                    <  1:3374537/69‑1 (MQ=255)
                   tgATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGc                                               <  1:3452205/71‑1 (MQ=255)
                     aTTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGAc                                                               >  1:534606/1‑53 (MQ=255)
                     aTTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACa                                                              >  1:3936130/1‑54 (MQ=255)
                     aTTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAg                                                  >  1:2143039/1‑66 (MQ=255)
                     aTTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCaa                                             <  1:3606552/71‑1 (MQ=255)
                      ttACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGaa                                                                     <  1:443471/46‑1 (MQ=255)
                        aCAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACTGGACATAACCCCATGAGTGCaa                                             <  1:1660161/68‑1 (MQ=255)
                        aCAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACTGGACATAACCCCATGAGTGCaa                                             <  1:1655035/68‑1 (MQ=255)
                        aCAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCaa                                             <  1:3615200/68‑1 (MQ=255)
                        aCAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCaa                                             >  1:2055167/1‑68 (MQ=255)
                        aCAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCaa                                             >  1:2796847/1‑68 (MQ=255)
                         cAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCaa                                             <  1:3197353/67‑1 (MQ=255)
                             gAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGa                                                                >  1:4438604/1‑44 (MQ=255)
                             gAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGa                                                                >  1:492451/1‑44 (MQ=255)
                             gAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCa                                              <  1:1663378/62‑1 (MQ=255)
                             gAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCa                                              <  1:3501629/62‑1 (MQ=255)
                              aaaaTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAAtt                                    <  1:629496/71‑1 (MQ=255)
                                aaTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTc                                  <  1:3491459/71‑1 (MQ=255)
                                aaTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTc                                  <  1:1287793/71‑1 (MQ=255)
                                aaTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTc                                  <  1:1214309/71‑1 (MQ=255)
                                 aTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCg                                 <  1:4114011/71‑1 (MQ=255)
                                  tAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGaa                                        <  1:1567152/63‑1 (MQ=255)
                                  tAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGa                                >  1:1345026/1‑71 (MQ=255)
                                    aaaaCCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacg                              >  1:796259/1‑71 (MQ=255)
                                    aaaaCCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacg                              <  1:2704351/71‑1 (MQ=255)
                                      aaCCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgc                             >  1:1185847/1‑70 (MQ=255)
                                       aCCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgcc                            >  1:416315/1‑70 (MQ=255)
                                         cATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGaa                                        <  1:2990908/56‑1 (MQ=255)
                                         cATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgccca                          <  1:3479406/70‑1 (MQ=255)
                                           tCACACCCCATTTAATATCAGGGAACCGGACATAACCC                                                        <  1:1220617‑M1/38‑1 (MQ=255)
                                           tCACACCCCATTTAATATCAGGGAACCGGACATAACCC                                                        <  1:2395566‑M1/38‑1 (MQ=255)
                                           tCACACCCCATTTAATATCAGGGAACCGGACATAACCC                                                        <  1:249151‑M1/38‑1 (MQ=255)
                                             acacCCCATTTAATATCAGGGAAACGGACATAACCCCATGAGTGCAATAGaa                                        >  1:2693011/1‑52 (MQ=255)
                                                  taattTAATATCAGTAAACCGGACATAACCCatgaagaaaaatcgcgcttttttgaag                             >  1:1442103‑M1/3‑31 (MQ=255)
                                                  taattTAATATCAGTAAACCGGACATAACCCatgaagaaaaatcgcgctttt                                   <  1:2518779‑M1/50‑22 (MQ=255)
                                                   cATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCatga               >  1:588066/1‑71 (MQ=255)
                                                   cATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCatga               >  1:1407137/1‑71 (MQ=255)
                                                    aTTTAATATCAGTAAACCGGACATAACCCatgaagaaaaatcgcgcttttttgaagtgggcagggggcaag              >  1:646539‑M1/1‑29 (MQ=255)
                                                    aTTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCatgat              >  1:3719914/1‑71 (MQ=255)
                                                     tttAATATCAGTAAACCGGACATAACCCatgaagaaaaatcgcgcttttttgaagtggg                         <  1:1083581‑M1/59‑32 (MQ=255)
                                                      ttAATATCAGTAAACCGGACATAACCCatgaagaaaaatcgcgcttttttgaagtgggcagggggcaagta            >  1:968414‑M1/1‑27 (MQ=255)
                                                      ttAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGa                                >  1:675173/1‑51 (MQ=255)
                                                       tAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgccca                          >  1:2866164/1‑56 (MQ=255)
                                                       tAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCatgatg             >  1:534663/1‑69 (MQ=255)
                                                       tAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCatgatg             >  1:2517073/1‑69 (MQ=255)
                                                        aaTATCAGGGAACCGGACATAACCCCATGAGTGCaa                                             >  1:1784008/1‑36 (MQ=255)
                                                         atatCAGGGCACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATgca           <  1:2634792/69‑1 (MQ=255)
                                                            tCAGGGAACCGGCCATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGcagcag       <  1:4164600/70‑1 (MQ=255)
                                                            tCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTc                                  <  1:3558507/43‑1 (MQ=255)
                                                            tCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTc                                  <  1:329972/43‑1 (MQ=255)
                                                            tCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCatga               >  1:2720519/1‑62 (MQ=255)
                                                              aGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTc                                  <  1:763525/41‑1 (MQ=255)
                                                              aGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTc                                  <  1:216454/41‑1 (MQ=255)
                                                              aGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTa     >  1:3780696/1‑70 (MQ=255)
                                                              aGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTAt    >  1:1059520/1‑71 (MQ=255)
                                                               gggAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATgcagc         <  1:2586920/65‑1 (MQ=255)
                                                                ggAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGa                                >  1:1788379/1‑41 (MQ=255)
                                                                ggAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGa                                >  1:506071/1‑41 (MQ=255)
                                                                ggAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTAtct  <  1:2471574/71‑1 (MQ=255)
                                                                 gAACCGGACATAACCCCATGAGTGCAATATAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTAtct  >  1:2905408/1‑70 (MQ=255)
                                                                 gAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTAtc   >  1:627229/1‑69 (MQ=255)
                                                                 |                                                                     
GTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTATCT  >  minE/1826676‑1826810

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: