Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,458,571 A→G 37.8% pseudogene (3384/4605 nt)
pseudogene (3384/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,5710AG37.8% 18.8 / 26.9 45pseudogene (3384/4605 nt)
pseudogene (3384/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 A (5/23);  new base G (5/12);  total (10/35)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-01

TTACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGCCATCTTCACCCCACCAGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACG  >  minE/1458507‑1458639
                                                                |                                                                    
ttACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCAt                                                                  <  1:820295/69‑1 (MQ=255)
ttACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCAt                                                                  <  1:720630/69‑1 (MQ=255)
ttACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCAt                                                                  <  1:1640591/69‑1 (MQ=255)
  aCGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCAtt                                                                 >  1:8509/1‑68 (MQ=255)
  aCGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCAtt                                                                 >  1:2592881/1‑68 (MQ=255)
  aCGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCAtt                                                                 >  1:2432846/1‑68 (MQ=255)
  aCGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCAtt                                                                 >  1:892648/1‑68 (MQ=255)
  aCGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCAtt                                                                 >  1:441563/1‑68 (MQ=255)
     ccGATCGGCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGcc                                                              <  1:267714/68‑1 (MQ=255)
     ccGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGcc                                                              <  1:1832888/68‑1 (MQ=255)
     ccGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGcc                                                              <  1:1874599/68‑1 (MQ=255)
     ccGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGcc                                                              <  1:1898754/68‑1 (MQ=255)
     ccGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGcc                                                              <  1:3370739/68‑1 (MQ=255)
     ccGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGcc                                                              <  1:564189/68‑1 (MQ=255)
     ccGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGcc                                                              <  1:560249/68‑1 (MQ=255)
         tCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGCCATCtt                                                         <  1:3134950/69‑1 (MQ=255)
         tCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGCCATCtt                                                         <  1:2599775/69‑1 (MQ=255)
         tCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGCCATCtt                                                         <  1:2603849/69‑1 (MQ=255)
         tCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGCCATCtt                                                         <  1:2643459/69‑1 (MQ=255)
         tCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGCCATCtt                                                         <  1:897328/69‑1 (MQ=255)
         tCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGCCATCtt                                                         <  1:2191595/69‑1 (MQ=255)
         tCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGCCATCtt                                                         <  1:3255560/69‑1 (MQ=255)
         tCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGCCATCtt                                                         <  1:1978303/69‑1 (MQ=255)
          cGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGCCATCtt                                                         <  1:1586199/68‑1 (MQ=255)
          cGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGCCATCtt                                                         <  1:2093484/68‑1 (MQ=255)
              tggctggcCTGCCAGTCGGCGTACCATGCCTATGCCGTGAGGAAGGCGTCGCCATTGCCATCTTCAccc                                                    <  1:1112293/69‑1 (MQ=255)
              tggctggcCTGCCAGTCGGCGTACCATGCCCATGCCGTGAGGAAGGCGTCGCCATTGCCATCTTCAccc                                                    <  1:2487592/69‑1 (MQ=255)
              tggctggcCTGCCAGTCGGCGTACCATGCCCATGCCGTGAGGAAGGCGTCGCCATTGCCATCTTCAccc                                                    <  1:2352661/69‑1 (MQ=255)
              tggctggcCTGCCAGTCGGCGTACCATGCCCATGCCGTGAGGAAGGCGTCGCCATTGCCATCTTCAccc                                                    <  1:1810568/69‑1 (MQ=255)
              tggctggcCTGCCAGTCGGCGTACCATGCCCATGCCGTGAGGAAGGCGTCGCCATTGCCATCTTCAccc                                                    <  1:2565974/69‑1 (MQ=255)
                                       atgcccatgccTTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTCCGGGcc                           <  1:691474/69‑1 (MQ=255)
                                       atgcccatgccGTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTCCGGGcc                           <  1:1945024/69‑1 (MQ=255)
                                       atgcccatgccGTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTCCGGGcc                           <  1:1637416/69‑1 (MQ=255)
                                       atgcccatgccGTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTCCGGGcc                           <  1:1987051/69‑1 (MQ=255)
                                       atgcccatgccGTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTCCGGGcc                           <  1:2869734/69‑1 (MQ=255)
                                       atgcccatgccGTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTCCGGGcc                           <  1:2704278/69‑1 (MQ=255)
                                       atgcccatgccGTAAGGAAGGCGTCACCATTGCCATCtt                                                         <  1:2796468/39‑1 (MQ=255)
                                       atgcccatgccGTAAGGAAGGCGTCACCATTGCCATCTTCACCCCACCAGGTAAAGCGCGCTCCGGGc                            <  1:687114/68‑1 (MQ=255)
                                       atgcccataccGTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTACGGGcc                           <  1:1308139/69‑1 (MQ=255)
                                        tgcccatgccGTAAGGAAGGCGTCACCATTGCCATCtt                                                         <  1:2391691/38‑1 (MQ=255)
                                               gccGTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTCCgg                              >  1:2196791/1‑58 (MQ=255)
                                               gccGTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTCCGGGc                            >  1:3434234/1‑60 (MQ=255)
                                               gccGTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTCCGGGc                            >  1:522200/1‑60 (MQ=255)
                                               gccGTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTCCGGGc                            >  1:1326567/1‑60 (MQ=255)
                                               gccGTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTCCGGGc                            >  1:1214899/1‑60 (MQ=255)
                                                                aCCATTGCCATCTTCACCCCACCAGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAg     >  1:2226407/1‑66 (MQ=255)
                                                                aCCATTGCCATCTTCACCCCACCAGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACg  >  1:922468/1‑69 (MQ=255)
                                                                aCCATTGCCATCTTCACCCCACCAGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACg  >  1:3085351/1‑69 (MQ=255)
                                                                aCCATTGCCATCTTCACCCCACCAGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACg  >  1:2273580/1‑69 (MQ=255)
                                                                aCCATTGCCATCTTCACCCCACCAGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACg  >  1:1133598/1‑69 (MQ=255)
                                                                |                                                                    
TTACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGCCATCTTCACCCCACCAGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACG  >  minE/1458507‑1458639

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: