Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,458,559 A→G 27.8% pseudogene (3396/4605 nt)
pseudogene (3396/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,5590AG27.8% 62.4 / 23.9 61pseudogene (3396/4605 nt)
pseudogene (3396/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/39);  new base G (5/12);  total (10/51)
Fisher's exact test for biased strand distribution p-value = 1.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTCCGGTTGTTGCGTTACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGCCATCTTCACCCCACCAGGTAAAGCGCGCTCCGGGCC  >  minE/1458492‑1458614
                                                                   |                                                       
aTTCCGGTTGTTGCGTTACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTa                                                         >  1:2641516/1‑68 (MQ=255)
   ccGGTTGTTGCGTTACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTaagga                                                     >  1:2008942/1‑69 (MQ=255)
   ccGGTTGTTGCGTTACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTaagga                                                     >  1:1563433/1‑69 (MQ=255)
   ccGGTTGTTGCGTTACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTaagga                                                     >  1:3136022/1‑69 (MQ=255)
       ttgttgCGTTACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTaaggaagac                                                 <  1:1906737/69‑3 (MQ=255)
       ttgttgCGTTACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGc                                                 <  1:1136176/69‑1 (MQ=255)
       ttgttgCGTTACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGc                                                 <  1:2328678/69‑1 (MQ=255)
       ttgttgCGTTACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGc                                                 <  1:255347/69‑1 (MQ=255)
       ttgttgCGTTACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGc                                                 <  1:2014820/69‑1 (MQ=255)
       ttgttgCGTTACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGc                                                 <  1:2945291/69‑1 (MQ=255)
       ttgttgCGTTACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGc                                                 <  1:1940278/69‑1 (MQ=255)
       ttgttgCGTTACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGc                                                 <  1:1882943/69‑1 (MQ=255)
       ttgttgCGTTACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGc                                                 <  1:51350/69‑1 (MQ=255)
       ttgttgCGTTACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGc                                                 <  1:1392169/69‑1 (MQ=255)
       ttgttgCGTTACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGc                                                 <  1:84181/69‑1 (MQ=255)
       ttgttgCGTTACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGc                                                 <  1:1240638/69‑1 (MQ=255)
       ttgttgCGTTACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGc                                                 <  1:1196825/69‑1 (MQ=255)
       ttgttgCGTTACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAACGc                                                 <  1:1465559/69‑1 (MQ=255)
       ttgttgCGTTACGCCGATCGCCTGGCTGGCCTGCCAGGCGGCGTACCATGCCCATGCCGTAAGGAAGGc                                                 <  1:1681379/69‑1 (MQ=255)
        tgttgCGTTACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGc                                                 <  1:416442/68‑1 (MQ=255)
              gTTACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGc                                                 <  1:1681382/62‑1 (MQ=255)
               ttACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCAt                                         <  1:1640591/69‑1 (MQ=255)
               ttACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCAt                                         <  1:720630/69‑1 (MQ=255)
               ttACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCAt                                         <  1:820295/69‑1 (MQ=255)
                 aCGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCAtt                                        >  1:441563/1‑68 (MQ=255)
                 aCGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCAtt                                        >  1:8509/1‑68 (MQ=255)
                 aCGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCAtt                                        >  1:2592881/1‑68 (MQ=255)
                 aCGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCAtt                                        >  1:892648/1‑68 (MQ=255)
                 aCGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCAtt                                        >  1:2432846/1‑68 (MQ=255)
                    ccGATCGGCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGcc                                     <  1:267714/68‑1 (MQ=255)
                    ccGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGcc                                     <  1:564189/68‑1 (MQ=255)
                    ccGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGcc                                     <  1:560249/68‑1 (MQ=255)
                    ccGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGcc                                     <  1:1832888/68‑1 (MQ=255)
                    ccGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGcc                                     <  1:1874599/68‑1 (MQ=255)
                    ccGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGcc                                     <  1:3370739/68‑1 (MQ=255)
                    ccGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGcc                                     <  1:1898754/68‑1 (MQ=255)
                        tCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGCCATCtt                                <  1:3134950/69‑1 (MQ=255)
                        tCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGCCATCtt                                <  1:897328/69‑1 (MQ=255)
                        tCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGCCATCtt                                <  1:3255560/69‑1 (MQ=255)
                        tCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGCCATCtt                                <  1:2643459/69‑1 (MQ=255)
                        tCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGCCATCtt                                <  1:2603849/69‑1 (MQ=255)
                        tCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGCCATCtt                                <  1:2191595/69‑1 (MQ=255)
                        tCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGCCATCtt                                <  1:1978303/69‑1 (MQ=255)
                        tCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGCCATCtt                                <  1:2599775/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:1810568/69‑1 (MQ=255)
                             tggctggcCTGCCAGTCGGCGTACCATGCCCATGCCGTGAGGAAGGCGTCGCCATTGCCATCTTCAccc                           <  1:2352661/69‑1 (MQ=255)
                             tggctggcCTGCCAGTCGGCGTACCATGCCCATGCCGTGAGGAAGGCGTCGCCATTGCCATCTTCAccc                           <  1:2487592/69‑1 (MQ=255)
                             tggctggcCTGCCAGTCGGCGTACCATGCCCATGCCGTGAGGAAGGCGTCGCCATTGCCATCTTCAccc                           <  1:2565974/69‑1 (MQ=255)
                                                      atgcccatgccTTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTCCGGGcc  <  1:691474/69‑1 (MQ=255)
                                                      atgcccatgccGTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTCCGGGcc  <  1:1987051/69‑1 (MQ=255)
                                                      atgcccatgccGTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTCCGGGcc  <  1:1945024/69‑1 (MQ=255)
                                                      atgcccatgccGTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTCCGGGcc  <  1:2869734/69‑1 (MQ=255)
                                                      atgcccatgccGTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTCCGGGcc  <  1:1637416/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)
                                                                   |                                                       
ATTCCGGTTGTTGCGTTACGCCGATCGCCTGGCTGGCCTGCCAGTCGGCGTACCATGCCCATGCCGTAAGGAAGGCGTCACCATTGCCATCTTCACCCCACCAGGTAAAGCGCGCTCCGGGCC  >  minE/1458492‑1458614

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: