Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,285,590 T→A 33.0% pseudogene (324/327 nt)
pseudogene (1050/1053 nt)
yeeL ←
yeeL ←
ECK1975:JW5325:b1979; hypothetical protein, C‑ter fragment
ECK1975:JW1961+JW5325:b4497; hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,285,5900TA33.0% 125.7 / 92.3 121pseudogene (324/327 nt)
pseudogene (1050/1053 nt)
yeeL
yeeL
ECK1975:JW5325:b1979; hypothetical protein, C‑ter fragment
ECK1975:JW1961+JW5325:b4497; hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (61/20);  new base A (16/24);  total (77/44)
Fisher's exact test for biased strand distribution p-value = 2.49e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

Only 100 of 123 total aligned reads displayed.

TTGATGACACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGC  >  minE/1285534‑1285656
                                                        |                                                                  
ttGATGACACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCaa                                                         >  1:3194387/1‑68 (MQ=255)
     gaCACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAacac                                                    <  1:1005506/68‑1 (MQ=255)
     gaCACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAacac                                                    <  1:2375522/68‑1 (MQ=255)
     gaCACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAacac                                                    >  1:2430781/1‑68 (MQ=255)
     gaCACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAacac                                                    <  1:2706173/68‑1 (MQ=255)
     gaCACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAacac                                                    <  1:1981090/68‑1 (MQ=255)
     gaCACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAacac                                                    >  1:1902164/1‑68 (MQ=255)
     gaCACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAacac                                                    <  1:1504442/68‑1 (MQ=255)
     gaCACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAacac                                                    <  1:1499557/68‑1 (MQ=255)
     gaCACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACt                                                   >  1:3424215/1‑69 (MQ=255)
     gaCACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACt                                                   >  1:2052393/1‑69 (MQ=255)
     gaCACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACt                                                   >  1:2044575/1‑69 (MQ=255)
     gaCACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACt                                                   >  1:3013479/1‑69 (MQ=255)
     gaCACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACt                                                   >  1:1865610/1‑69 (MQ=255)
     gaCACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACt                                                   >  1:3177385/1‑69 (MQ=255)
     gaCACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACt                                                   >  1:175306/1‑69 (MQ=255)
     gaCACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACt                                                   >  1:1751177/1‑69 (MQ=255)
     gaCACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACt                                                   >  1:32118/1‑69 (MQ=255)
     gaCACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACt                                                   >  1:3286599/1‑69 (MQ=255)
     gaCACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACt                                                   >  1:3301558/1‑69 (MQ=255)
     gaCACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACt                                                   >  1:434823/1‑69 (MQ=255)
     gaCACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACt                                                   >  1:1111805/1‑69 (MQ=255)
     gaCACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACt                                                   >  1:3311707/1‑69 (MQ=255)
     gaCACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACt                                                   >  1:1353961/1‑69 (MQ=255)
     gaCACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAAAACt                                                   >  1:1446229/1‑69 (MQ=255)
           accGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTaa                                             >  1:3347332/1‑69 (MQ=255)
           accGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTaa                                             >  1:2006617/1‑69 (MQ=255)
           accGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTaa                                             >  1:1258344/1‑69 (MQ=255)
           accGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTaa                                             >  1:1642281/1‑69 (MQ=255)
                     aCAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCg                                   <  1:1436269/69‑1 (MQ=255)
                     aCAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCg                                   <  1:2630388/69‑1 (MQ=255)
                     aCAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCg                                   <  1:857851/69‑1 (MQ=255)
                     aCAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCg                                   <  1:3129773/69‑1 (MQ=255)
                     aCAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCg                                   <  1:2146553/69‑1 (MQ=255)
                     aCAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCg                                   <  1:588502/69‑1 (MQ=255)
                     aCAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCg                                   <  1:1952227/69‑1 (MQ=255)
                     aCAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCg                                   <  1:434784/69‑1 (MQ=255)
                     aCAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCg                                   <  1:1725135/69‑1 (MQ=255)
                     aCAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCg                                   <  1:79537/69‑1 (MQ=255)
                            ggCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCa                             >  1:1499118/1‑68 (MQ=255)
                            ggCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCa                             >  1:2777931/1‑68 (MQ=255)
                            ggCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCa                             >  1:1264106/1‑68 (MQ=255)
                            ggCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCa                             >  1:1416886/1‑68 (MQ=255)
                                  tcTATGCTTAAGTACTGCCTCAAATTCTCGGCAAAACACttt                                                 <  1:553176/41‑1 (MQ=255)
                                  aCTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGa                       >  1:2938223/1‑68 (MQ=255)
                                  aCTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGa                       >  1:709047/1‑68 (MQ=255)
                                  aCTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGa                       >  1:2690665/1‑68 (MQ=255)
                                  aCTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGa                       >  1:2584175/1‑68 (MQ=255)
                                  aCTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGa                       >  1:1006568/1‑68 (MQ=255)
                                  aCTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGa                       >  1:1068544/1‑68 (MQ=255)
                                  aCTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGa                       >  1:168573/1‑68 (MQ=255)
                                  aCTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGa                       >  1:1817218/1‑68 (MQ=255)
                                  aCTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCANCCAg                        >  1:964914/1‑67 (MQ=255)
                                  aCTATGCATAAGTACTGCCTCAAATTCTCGCCAAAACACttt                                                 <  1:1636252/42‑1 (MQ=255)
                                  aCTATGCATAAGTACTGCCTCAAATTCTCGCCAAAACACttt                                                 <  1:232633/42‑1 (MQ=255)
                                  aCTATGCATAAGTACTGCCTCAAATTCTCGCCAAAACACttt                                                 <  1:941350/42‑1 (MQ=255)
                                  aCTATGCATAAGTACTGCCTCAAATTCTCGCCAAAACACttt                                                 <  1:939344/42‑1 (MQ=255)
                                  aCTATGCATAAGTACTGCCTCAAATTCTCGCCAAAACACttt                                                 <  1:866138/42‑1 (MQ=255)
                                  aCTATGCATAAGTACTGCCTCAAATTCTCGCCAAAACACttt                                                 <  1:1031957/42‑1 (MQ=255)
                                  aCTATGCATAAGTACTGCCTCAAATTCTCGCCAAAACACttt                                                 <  1:1035840/42‑1 (MQ=255)
                                  aCTATGCATAAGTACTGCCTCAAATTCTCGCCAAAACACttt                                                 <  1:571946/42‑1 (MQ=255)
                                  aCTATGCATAAGTACTGCCTCAAATTCTCGCCAAAACACttt                                                 <  1:446997/42‑1 (MQ=255)
                                  aCTATGCATAAGTACTGCCTCAAATTCTCGCCAAAACACttt                                                 <  1:1633257/42‑1 (MQ=255)
                                  aCTATGCATAAGTACTGCCTCAAATTCTCGCCAAAACACttt                                                 <  1:2429514/42‑1 (MQ=255)
                                  aCTATGCATAAGTACTGCCTCAAATTCTCGCCAAAACACttt                                                 <  1:319498/42‑1 (MQ=255)
                                  aCTATGCATAAGTACTGCCTCAAATTCTCGCCAAAACACttt                                                 <  1:1868329/42‑1 (MQ=255)
                                  aCTATGCATAAGTACTGCCTCAAATTCTCGCCAAAACACttt                                                 <  1:1991665/42‑1 (MQ=255)
                                  aCTATGCATAAGTACTGCCTCAAATTCTCGCCAAAACACttt                                                 <  1:2638957/42‑1 (MQ=255)
                                  aCTATGCATAAGTACTGCCTCAAATTCTCGCCAAAACACttt                                                 <  1:3004902/42‑1 (MQ=255)
                                  aCTATGCATAAGTACTGCCTCAAATTCTCGCCAAAACACttt                                                 <  1:1884489/42‑1 (MQ=255)
                                   cTATGCATAAGTACTGCCTCAAATTCTCGCCAAAACACttt                                                 <  1:3172507/41‑1 (MQ=255)
                                   cTATGCATAAGTACTGCCTCAAATTCTCGCCAAAACACttt                                                 <  1:1484284/41‑1 (MQ=255)
                                     aTGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTta                   <  1:2997255/69‑1 (MQ=255)
                                       gCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTtata                 >  1:319551/1‑69 (MQ=255)
                                       gCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTtata                 >  1:951312/1‑69 (MQ=255)
                                       gCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTtata                 >  1:2221064/1‑69 (MQ=255)
                                       gCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTtata                 >  1:3187560/1‑69 (MQ=255)
                                       gCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTtata                 >  1:3042665/1‑69 (MQ=255)
                                       gCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTtat                  >  1:530838/1‑68 (MQ=255)
                                            aGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTc            >  1:2259341/1‑69 (MQ=255)
                                            aGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTc            >  1:3377717/1‑69 (MQ=255)
                                            aGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTc            >  1:3381242/1‑69 (MQ=255)
                                            aGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTc            >  1:2157412/1‑69 (MQ=255)
                                            aGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTc            >  1:2031532/1‑69 (MQ=255)
                                                  gcctcAAATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAAtttt      >  1:2160186/3‑69 (MQ=255)
                                                  gcctcAAATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAAtttt      >  1:2629308/3‑69 (MQ=255)
                                                  gcctcAAATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAAtttt      >  1:1021365/3‑69 (MQ=255)
                                                  gcctcAAATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAAtttt      >  1:658980/3‑69 (MQ=255)
                                                  gcctcAAATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAAtttt      >  1:2011214/3‑69 (MQ=255)
                                                  gcctcAAATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAAtttt      >  1:106879/3‑69 (MQ=255)
                                                  gcctcAAATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAAtttt      >  1:449717/3‑69 (MQ=255)
                                                  gcctcAAATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAAtttt      >  1:1093493/3‑69 (MQ=255)
                                                  gcctcAAATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAAtttt      >  1:117249/3‑69 (MQ=255)
                                                  gcctcAAATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAAtttt      >  1:1386110/3‑69 (MQ=255)
                                                  gcctcAAATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAAtttt      >  1:3224519/3‑69 (MQ=255)
                                                  gcctcAAATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAAtttt      >  1:1741052/3‑69 (MQ=255)
                                                  gcctcAAATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCAAGACTTATATTGTCAAtttt      >  1:1228187/3‑69 (MQ=255)
                                                  gcctcAAATTCTCGCCAAAACACTTTTTAATGCCTTTCCGAATTCATCCAGACTTATATTGTCAAtttt      >  1:2299907/3‑69 (MQ=255)
                                                       atatTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGc  >  1:2769845/1‑68 (MQ=255)
                                                       atatTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGc  >  1:3172645/1‑68 (MQ=255)
                                                        |                                                                  
TTGATGACACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGC  >  minE/1285534‑1285656

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: