Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,910,924 T→C 42.5% P97P (CCA→CCG ppdA ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,910,9240TC42.5% 26.2 / 72.3 94P97P (CCA→CCGppdAconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (21/33);  new base C (16/24);  total (37/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.38e-01

CGGTATTGCGCAGGCCAAAGAAAGCAAGCGAAGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTGGCCCCAGCAGCGGAACTCACAAGGCAC  >  minE/1910857‑1910990
                                                                   |                                                                  
cGGTATTGCGCAGGCCAAAGAAAGCAAGCGAAGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGt                                                                    >  1:2936854/1‑68 (MQ=255)
cGGTATTGCGCAGGCCAAAGAAAGCAAGCGAAGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGt                                                                    >  1:2550867/1‑68 (MQ=255)
cGGTATTGCGCAGGCCAAAGAAAGCAAGCGAAGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGt                                                                    >  1:335567/1‑68 (MQ=255)
       gcgcAGGCCAAAGAAAGCAAGCGAAGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACaaa                                                            >  1:2820061/1‑69 (MQ=255)
       gcgcAGGCCAAAGAAAGCAAGCGAAGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACaaa                                                            >  1:2735427/1‑69 (MQ=255)
       gcgcAGGCCAAAGAAAGCAAGCGAAGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACaaa                                                            >  1:2551054/1‑69 (MQ=255)
       gcgcAGGCCAAAGAAAGCAAGCGAAGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACaaa                                                            >  1:1705980/1‑69 (MQ=255)
           aGGCCAAAGAAAGCAAGCGAAGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGAcc                                                        <  1:2021236/69‑1 (MQ=255)
             gCCAAAGAAAGCAAGCGAAGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCaa                                                      >  1:112782/1‑69 (MQ=255)
             gCCAAAGAAAGCAAGCGAAGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCaa                                                      >  1:3055370/1‑69 (MQ=255)
                aaagaaagCAAGCGAAGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATg                                                    >  1:973064/1‑68 (MQ=255)
                aaagaaagCAAGCGAAGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATg                                                    >  1:2303311/1‑68 (MQ=255)
                   gaaagCAAGCGAAGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGa                                                <  1:2764903/69‑1 (MQ=255)
                   gaaagCAAGCGAAGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGa                                                <  1:2692431/69‑1 (MQ=255)
                   gaaagCAAGCGAAGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGa                                                <  1:3371120/69‑1 (MQ=255)
                   gaaagCAAGCGAAGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGa                                                <  1:465000/69‑1 (MQ=255)
                   gaaagCAAGCGAAGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGa                                                <  1:1625426/69‑1 (MQ=255)
                   gaaagCAAGCGAAGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGa                                                <  1:1516248/69‑1 (MQ=255)
                           gcGAAGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGCCAt                                        <  1:1808215/69‑1 (MQ=255)
                           gcGAAGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGCCAt                                        <  1:934098/69‑1 (MQ=255)
                           gcGAAGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGCCAt                                        <  1:1371946/69‑1 (MQ=255)
                           gcGAAGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGCCAt                                        <  1:87363/69‑1 (MQ=255)
                           gcGAAGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGCCAt                                        <  1:2349939/69‑1 (MQ=255)
                           gcGAAGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGCCAt                                        <  1:3204884/69‑1 (MQ=255)
                           gcGAAGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGCCAt                                        <  1:806485/69‑1 (MQ=255)
                           gcGAAGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGCCAt                                        <  1:3415834/69‑1 (MQ=255)
                               aGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGCGTCACAAAGACCAATGGTGAACTg                                            <  1:916835/61‑1 (MQ=255)
                               aGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTg                                            <  1:3131297/61‑1 (MQ=255)
                               aGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTg                                            <  1:1402965/61‑1 (MQ=255)
                               aGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTg                                            <  1:2695357/61‑1 (MQ=255)
                               aGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTg                                            <  1:1413104/61‑1 (MQ=255)
                               aGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTg                                            <  1:2670677/61‑1 (MQ=255)
                               aGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTg                                            <  1:2655279/61‑1 (MQ=255)
                               aGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTg                                            <  1:2962210/61‑1 (MQ=255)
                               aGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTg                                            <  1:3060780/61‑1 (MQ=255)
                               aGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTg                                            <  1:1487689/61‑1 (MQ=255)
                               aGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTg                                            <  1:192793/61‑1 (MQ=255)
                               aGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTg                                            <  1:800105/61‑1 (MQ=255)
                               aGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTg                                            <  1:3176602/61‑1 (MQ=255)
                               aGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTg                                            <  1:2431483/61‑1 (MQ=255)
                               aGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTg                                            <  1:1935261/61‑1 (MQ=255)
                               aGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTg                                            <  1:639220/61‑1 (MQ=255)
                               aGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTg                                            <  1:3199897/61‑1 (MQ=255)
                               aGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTg                                            <  1:733690/61‑1 (MQ=255)
                               aGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTg                                            <  1:3319926/61‑1 (MQ=255)
                               aGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTg                                            <  1:3323984/61‑1 (MQ=255)
                               aGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTg                                            <  1:3324080/61‑1 (MQ=255)
                               aGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTg                                            <  1:204776/61‑1 (MQ=255)
                               aGGTGGCAGGTCGCTCATTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTg                                            <  1:634980/61‑1 (MQ=255)
                                  tgtCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGCCATGAGAGGt                                 <  1:1828081/69‑1 (MQ=255)
                                  tgtCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGCCATGACAGGt                                 <  1:832248/69‑1 (MQ=255)
                                  tgtCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGCCATGACAGGt                                 <  1:12538/69‑1 (MQ=255)
                                  tgtCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGCCATGACAGGt                                 <  1:1992735/69‑1 (MQ=255)
                                  tgtCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGCCATGACAGGt                                 <  1:210023/69‑1 (MQ=255)
                                  tgtCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGCCATGACAGGt                                 <  1:2554889/69‑1 (MQ=255)
                                  tgtCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGCCATGACAGGt                                 <  1:2557304/69‑1 (MQ=255)
                                   gtCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGCCATGACAGGt                                 <  1:2295838/68‑1 (MQ=255)
                                      aGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGa                                                <  1:1351542/50‑1 (MQ=255)
                                          cGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGc                                           >  1:2637591/1‑51 (MQ=255)
                                          cGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGc                                           <  1:2212556/51‑1 (MQ=255)
                                          cGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGc                                           >  1:69525/1‑51 (MQ=255)
                                          cGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGc                                           >  1:784382/1‑51 (MQ=255)
                                          cGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGc                                           >  1:1604195/1‑51 (MQ=255)
                                          cGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGc                                           >  1:2498016/1‑51 (MQ=255)
                                          cGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGc                                           >  1:1540468/1‑51 (MQ=255)
                                          cGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGc                                           >  1:1482033/1‑51 (MQ=255)
                                          cGCTCATTTCGACTTCGGGCCAGAGTGGCACAAAGACCAATGGTGAACTg                                            >  1:890950/1‑50 (MQ=255)
                                           gCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGCCCAATGGTGAACTGCCATGACAGGt                                 <  1:2811041/60‑1 (MQ=255)
                                              cATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTGGCCCcagc                     <  1:2475507/69‑1 (MQ=255)
                                              cATTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTg                                            <  1:2232081/46‑1 (MQ=255)
                                               aattCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTCGcc                          >  1:1492638/3‑63 (MQ=255)
                                               aTTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTGGCCCcagc                     <  1:331873/68‑1 (MQ=255)
                                               aTTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTCGccc                         >  1:2739452/1‑64 (MQ=255)
                                               aTTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTCGccc                         >  1:725684/1‑64 (MQ=255)
                                               aTTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTCGccc                         >  1:1685608/1‑64 (MQ=255)
                                               aTTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTCGccc                         >  1:1791411/1‑64 (MQ=255)
                                               aTTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTCGccc                         >  1:1348940/1‑64 (MQ=255)
                                               aTTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTCGccc                         >  1:1951906/1‑64 (MQ=255)
                                               aTTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTCGccc                         >  1:2851086/1‑64 (MQ=255)
                                               aTTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTCGccc                         >  1:29239/1‑64 (MQ=255)
                                               aTTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTCGccc                         >  1:2731172/1‑64 (MQ=255)
                                               aTTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTCGccc                         >  1:3128661/1‑64 (MQ=255)
                                               aTTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTCGccc                         >  1:3242617/1‑64 (MQ=255)
                                               aTTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTCGccc                         >  1:3259230/1‑64 (MQ=255)
                                               aTTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTCGcc                          >  1:1616782/1‑63 (MQ=255)
                                               aTTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTCGcc                          >  1:2025100/1‑63 (MQ=255)
                                               aTTTCGACTTCGGGCCAGCGCGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTCGcc                          >  1:328058/1‑63 (MQ=255)
                                                  tCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTGGCCCcagcag                   >  1:2916947/1‑67 (MQ=255)
                                                  tCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTGGCCCCagcagc                  >  1:2036855/1‑68 (MQ=255)
                                                       ttCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTGGCCCCAGCAGCGGAACt            >  1:724978/1‑69 (MQ=255)
                                                       ttCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTGGCCCCAGCAGCGGAACt            >  1:3270773/1‑69 (MQ=255)
                                                       ttCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTGGCCCCAGCAGCGGAACt            >  1:154805/1‑69 (MQ=255)
                                                         cGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGCCATGACAGGt                                 <  1:550278/46‑1 (MQ=255)
                                                           ggCCAGCGTGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTGGCCCCAGCAGCGGAACTcac         <  1:3002598/68‑1 (MQ=255)
                                                                 cGTGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTGGCCCCAGCAGCGGAACTCACAAGGCAc  <  1:3218381/69‑1 (MQ=255)
                                                                 cGTGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTGGCCCCAGCAGCGGAACTCACAAGGCAc  <  1:2119394/69‑1 (MQ=255)
                                                                 cGTGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTGGCCCCAGCAGCGGAACTCACAAGGCAc  <  1:3357710/69‑1 (MQ=255)
                                                                   |                                                                  
CGGTATTGCGCAGGCCAAAGAAAGCAAGCGAAGGTGTCAGGTCGCTCATTTCGACTTCGGGCCAGCGTGGCACAAAGACCAATGGTGAACTGCCATGACAGGTATTGGCCCCAGCAGCGGAACTCACAAGGCAC  >  minE/1910857‑1910990

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: