Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,458,271 C→T 37.5% pseudogene (3684/4605 nt)
pseudogene (3684/4104 nt)
yfaS ←
yfaS ←
ECK2220:JW2222+JW2221:b4500; hypothetical protein
ECK2220:JW2222:b2228; hypothetical protein, N‑ter fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,458,2710CT37.5% 67.6 / 68.2 109pseudogene (3684/4605 nt)
pseudogene (3684/4104 nt)
yfaS
yfaS
ECK2220:JW2222+JW2221:b4500; hypothetical protein
ECK2220:JW2222:b2228; hypothetical protein, N‑ter fragment
Reads supporting (aligned to +/- strand):  ref base C (46/22);  new base T (10/31);  total (56/53)
Fisher's exact test for biased strand distribution p-value = 1.29e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

GACTTCAACTGCGCTTTTTTCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGATATCGCGGGTGTCTTCCTCAGAGAA  >  minE/1458203‑1458327
                                                                    |                                                        
gACTTCAACTGCGCTTTTTTCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGAc                                                          >  1:42708/1‑69 (MQ=255)
     cAACTGCGCTTTTTTCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCgg                                                     >  1:922066/1‑69 (MQ=255)
     cAACTGCGCTTTTTTCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCgg                                                     >  1:3084550/1‑69 (MQ=255)
         tGCGCTTTTTTCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGCCTCCGGTGTa                                                 >  1:516606/1‑69 (MQ=255)
         tGCGCTTTTTTCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGgtgt                                                  >  1:702639/1‑68 (MQ=255)
         tGCGCTTTTTTCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTa                                                 >  1:693208/1‑69 (MQ=255)
         tGCGCTTTTTTCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTa                                                 >  1:976793/1‑69 (MQ=255)
         tGCGCTTTTTTCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTa                                                 >  1:3017284/1‑69 (MQ=255)
         tGCGCTTTTTTCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTa                                                 >  1:3041354/1‑69 (MQ=255)
         tGCGCTTTTTTCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTa                                                 >  1:680049/1‑69 (MQ=255)
         tGCGCTTTTTTCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTa                                                 >  1:3251620/1‑69 (MQ=255)
         tGCGCTTTTTTCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTa                                                 >  1:1822171/1‑69 (MQ=255)
         tGCGCTTTTTTCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTa                                                 >  1:1769222/1‑69 (MQ=255)
          gggcTTTTTTCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCC‑‑‑GGAGACTCCGGTGTATCg                                              >  1:3330280/3‑68 (MQ=255)
          gcgcTTTTTTCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTa                                                 >  1:2115267/1‑68 (MQ=255)
          gcgcTTTTTTCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTa                                                 >  1:1008737/1‑68 (MQ=255)
          gcgcTTTTTTCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTa                                                 >  1:895237/1‑68 (MQ=255)
          gcgcTTTTTTCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTa                                                 >  1:2423196/1‑68 (MQ=255)
          gcgcTTTTTTCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTa                                                 >  1:906371/1‑68 (MQ=255)
              ttttttCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGa                                             >  1:1800839/1‑68 (MQ=255)
              ttttttCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAg                                            >  1:1504682/1‑69 (MQ=255)
              ttttttCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAg                                            >  1:660181/1‑69 (MQ=255)
              ttttttCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAg                                            >  1:2950625/1‑69 (MQ=255)
              ttttttCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAg                                            >  1:2129613/1‑69 (MQ=255)
              ttttttCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAg                                            >  1:353137/1‑69 (MQ=255)
              ttttttCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAg                                            >  1:2020860/1‑69 (MQ=255)
              ttttttCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAg                                            >  1:804844/1‑69 (MQ=255)
              ttttttCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAg                                            >  1:2100267/1‑69 (MQ=255)
              atttttCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCACTGGAGATTCAGGTGTATCGAg                                            >  1:2102130/2‑69 (MQ=255)
                  ttCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAgt                                           <  1:768784/66‑2 (MQ=255)
                  ttCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAgg                                           <  1:1609600/66‑1 (MQ=255)
                  ttCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAgg                                           <  1:1201662/66‑1 (MQ=255)
                  ttCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAgg                                           <  1:1290162/66‑1 (MQ=255)
                  ttCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAgg                                           <  1:1369298/66‑1 (MQ=255)
                  ttCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAgg                                           <  1:88082/66‑1 (MQ=255)
                  ttCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAgg                                           <  1:3247437/66‑1 (MQ=255)
                  ttCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAGGa                                          <  1:1272016/67‑1 (MQ=255)
                  ttCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAGGa                                          <  1:2686919/67‑1 (MQ=255)
                                ggTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAGGATCAGTCTATCATTGa                           <  1:3026057/68‑1 (MQ=255)
                                ggTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGat                          <  1:1215919/69‑1 (MQ=255)
                                ggTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGat                          <  1:3337815/69‑1 (MQ=255)
                                ggTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGa                           <  1:3203978/68‑1 (MQ=255)
                                       aCGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGATATcgcg                    >  1:380303/1‑68 (MQ=255)
                                       aCGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGATATc                       >  1:3205901/1‑65 (MQ=255)
                                       aCGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGATATCGCgg                   >  1:989331/1‑69 (MQ=255)
                                       aCGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGATATCGCgg                   >  1:2379384/1‑69 (MQ=255)
                                       aCGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGATATCGCgg                   >  1:1492893/1‑69 (MQ=255)
                                       aCGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGATATCGCgg                   >  1:847138/1‑69 (MQ=255)
                                        cGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAGGATCAGGGTGTCATTg                            <  1:145210/59‑1 (MQ=255)
                                        cGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAGGATCAGGCTGTCATTg                            <  1:2515216/59‑1 (MQ=255)
                                        cGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAGGATCAGGCTGTCATTg                            <  1:670823/59‑1 (MQ=255)
                                        cGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAGGATCAGGCTGTCATTg                            <  1:1465286/59‑1 (MQ=255)
                                        cGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAGGATCAGGCTGTCATTg                            <  1:964588/59‑1 (MQ=255)
                                        cGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAGGATCAGGCTGTCATTg                            <  1:3120443/59‑1 (MQ=255)
                                        cGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAGGATCAGGCTGTCATTg                            <  1:3086875/59‑1 (MQ=255)
                                        cGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAGGATCAGGCTGTCATTg                            <  1:1965475/59‑1 (MQ=255)
                                        cGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAGGATCAGGCTGTCATTg                            <  1:3055445/59‑1 (MQ=255)
                                        cGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAGGATCAGGCTGTCATTg                            <  1:808820/59‑1 (MQ=255)
                                        cGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAGGATCAGGCTGTCATTg                            <  1:2999197/59‑1 (MQ=255)
                                        cGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAGGATCAGGCTGTCATTg                            <  1:2984757/59‑1 (MQ=255)
                                        cGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAGGATCAGGCTGTCATTg                            <  1:2174644/59‑1 (MQ=255)
                                        cGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAGGATCAGGCTGTCATTg                            <  1:2631263/59‑1 (MQ=255)
                                        cGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAGGATCAGGCTGTCATTg                            <  1:2696614/59‑1 (MQ=255)
                                        cGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAGGATCAGGCTGTCATTg                            <  1:724519/59‑1 (MQ=255)
                                        cGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAGGATCAGGCTGTCATTg                            <  1:2488394/59‑1 (MQ=255)
                                        cGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAGGATCAGGCTGTCATTg                            <  1:2866065/59‑1 (MQ=255)
                                        cGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAGGATCAGGCTGTCATTg                            <  1:2974012/59‑1 (MQ=255)
                                        cGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAGGATCAGGCTGTCATTg                            <  1:2923628/59‑1 (MQ=255)
                                        cGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAGGATCAGGCTGTCATTg                            <  1:2946980/59‑1 (MQ=255)
                                        cGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAGGATCAGGCTGTCATTg                            <  1:2953405/59‑1 (MQ=255)
                                        cGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAGGATCAGGCTGTCATTg                            <  1:2209349/59‑1 (MQ=255)
                                        cGTTTGCCACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAGGATCACGCTGTCATTg                            <  1:2048959/59‑1 (MQ=255)
                                         gTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAgg                                           <  1:1406791/43‑1 (MQ=255)
                                         gTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAgg                                           <  1:2150977/43‑1 (MQ=255)
                                          tttGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGATATCGCGGgtg                <  1:560441/69‑1 (MQ=255)
                                          tttGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGATATCGCGGgtg                <  1:211146/69‑1 (MQ=255)
                                               cACCGCATCTGCCAGTGGAGATTCAGGTGTATCGAGGATCAGGCTGTCATTg                            <  1:1443781/52‑1 (MQ=255)
                                                      tCTGCCAGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGATATCGCGGGTGTCTTCCTCaga     >  1:2870418/1‑68 (MQ=255)
                                                      tCTGCCAGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGATATCGCGGGTGTCTTCCTCaga     >  1:1491751/1‑68 (MQ=255)
                                                            aGTGGAGACTCCGGTGTATCGAGGTTCAGGCTATCATTGATATCGCGGGTGTCTTCCTCAGAGaa  <  1:1168176/65‑1 (MQ=255)
                                                            aGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCAttta                           >  1:2286665/1‑38 (MQ=255)
                                                            aGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGat                          >  1:812747/1‑41 (MQ=255)
                                                            aGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGat                          >  1:2803802/1‑41 (MQ=255)
                                                            aGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGat                          >  1:1935126/1‑41 (MQ=255)
                                                            aGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGat                          >  1:1133193/1‑41 (MQ=255)
                                                            aGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGat                          >  1:1473270/1‑41 (MQ=255)
                                                            aGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGat                          >  1:1552927/1‑41 (MQ=255)
                                                            aGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGat                          >  1:408948/1‑41 (MQ=255)
                                                            aGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGat                          >  1:2073604/1‑41 (MQ=255)
                                                            aGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGat                          >  1:2129047/1‑41 (MQ=255)
                                                            aGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGat                          >  1:2312904/1‑41 (MQ=255)
                                                            aGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGat                          >  1:2967146/1‑41 (MQ=255)
                                                            aGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGat                          >  1:2973195/1‑41 (MQ=255)
                                                            aGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGat                          >  1:3062875/1‑41 (MQ=255)
                                                            aGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGa                           >  1:2060392/1‑40 (MQ=255)
                                                            aGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGa                           >  1:1883419/1‑40 (MQ=255)
                                                            aGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGATATCGCGGGTGTCTTCCTCAGAGaa  <  1:1901292/65‑1 (MQ=255)
                                                            aGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGATATCGCGGGTGTCTTCCTCAGAGaa  <  1:1599671/65‑1 (MQ=255)
                                                            aGTCGAGACTCCGGTGTATCGAGGATCAGGCTAACATTGat                          >  1:1219453/1‑41 (MQ=255)
                                                               ggAGATTCAGGTGTATCGAGGATCAGGCTGTCATTg                            <  1:1838782/36‑1 (MQ=25)
                                                                    |                                                        
GACTTCAACTGCGCTTTTTTCAGCAACGTCATGGTTAAGACGTTTGCCACCGCATCTGCCAGTGGAGACTCCGGTGTATCGAGGATCAGGCTATCATTGATATCGCGGGTGTCTTCCTCAGAGAA  >  minE/1458203‑1458327

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: