Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 547,835 +G pseudogene (255/261 nt)
pseudogene (255/1259 nt)
ylbE →
ylbE →
ECK0512:JW0507:b4507; hypothetical protein, N‑ter fragment
ECK0512:JW0508+JW0507:b4572; hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb547,8311.G97.2% 124.5 / ‑3.6 36pseudogene (251/261 nt)
pseudogene (251/1259 nt)
ylbE
ylbE
ECK0512:JW0507:b4507; hypothetical protein, N‑ter fragment
ECK0512:JW0508+JW0507:b4572; hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (10/25);  total (10/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

GCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑GGG‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAT  >  W3110S.gb/547771‑547890
                                                             |                                                            
gcgcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑GGGTGAAGTg                                                    <  1:3058/71‑1 (MQ=255)
   cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGG‑GAAGTGAAc                                                 <  1:1388174/71‑1 (MQ=255)
   cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGG‑GAAGTGAAc                                                 <  1:481687/71‑1 (MQ=255)
   cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGG‑GAAGTGAAc                                                 <  1:3449605/71‑1 (MQ=255)
   cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGG‑GAAGTGAAc                                                 <  1:2596054/71‑1 (MQ=255)
   cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGG‑GAAGTGAAc                                                 <  1:674442/71‑1 (MQ=255)
         tGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGG‑GAAGTGAACTTCAtt                                           >  1:502765/1‑71 (MQ=255)
         tGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGG‑GAAGTGAACTTCAtt                                           >  1:19023/1‑71 (MQ=255)
         tGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGG‑GAAGTGAACTTCAtt                                           >  1:2771568/1‑71 (MQ=255)
         tGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGG‑GAAGTGAACTTCAt                                            >  1:1705612/1‑70 (MQ=255)
                 ggTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑g                                                             >  1:2305208/1‑45 (MQ=255)
                   ttGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGG‑GAAGTGAACTTCATTCCTTGTCAcc                                 <  1:3865089/71‑1 (MQ=255)
                   ttGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGG‑GAAGTGAACTTCATTCCTTGTCAcc                                 <  1:713975/71‑1 (MQ=255)
                   ttGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGG‑GAAGTGAACTTCATTCCTTGTCAcc                                 <  1:3471555/71‑1 (MQ=255)
                      ggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGG‑GAAGTGAACTTCATTCCTTGTCACCAt                               <  1:679854/70‑1 (MQ=255)
                      ggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGG‑GAAGTGAACTTCATTCCTTGTCACCAt                               <  1:3713231/70‑1 (MQ=255)
                      ggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGG‑GAAGTGAACTTCATTCCTTGTCACCAt                               <  1:2850138/70‑1 (MQ=255)
                      ggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGG‑GAAGTGAACTTCATTCCTTGTCACCAt                               <  1:2998943/70‑1 (MQ=255)
                      ggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGG‑GAAGTGAACTTCATTCCTTGTCACCAt                               <  1:3155072/70‑1 (MQ=255)
                       gCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGG‑GAAGTGAACTTCATTCCTTGTCACCAt                               <  1:1567788/69‑1 (MQ=255)
                               tGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGG‑GAAGTGAACt                                                >  1:2169057/1‑44 (MQ=38)
                               tGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGG‑GAAGTGAACt                                                >  1:3060881/1‑44 (MQ=38)
                               tGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGG‑GAAGTGAACt                                                >  1:2318940/1‑44 (MQ=38)
                               tGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGG‑GAAGTGAACt                                                >  1:2072930/1‑44 (MQ=38)
                                gAAGCGCAGGCGCTGGCAATACTGGAGCAGGGG‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGt                    >  1:625732/1‑71 (MQ=255)
                                gAAGCGCAGGCGCTGGCAATACTGGAGCAGGGG‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGt                    >  1:250001/1‑71 (MQ=255)
                                    cgcAGGCGCTGGCAATACTGGAGCAGGGG‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCggg                <  1:3650482/71‑1 (MQ=255)
                                     gcAGGCGCTGGCAATACTGGAGCAGGGG‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCggg                <  1:730496/70‑1 (MQ=255)
                                       aGGCTCTGGCAATACTGGAGCAGGGG‑GAAGTGAACTTCATTCCTTGTCAcc                                 <  1:1530093/51‑1 (MQ=255)
                                            cTGGCAATACTGGAGCAGGGG‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc              <  1:565391/65‑1 (MQ=255)
                                            cTGGCAATACTGGAGCAGGGG‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc              <  1:2432292/65‑1 (MQ=255)
                                            cTGGCAATACTGGAGCAGGGG‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc              <  1:1990463/65‑1 (MQ=255)
                                                     cTGGAGCAGGGG‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt  <  1:479268/68‑1 (MQ=255)
                                                     cTGGAGCAGGGG‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt  <  1:415663/68‑1 (MQ=255)
                                                     cTGGAGCAGGGG‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt  <  1:54336/68‑1 (MQ=255)
                                                     cTGGAGCAGGGG‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt  <  1:267293/68‑1 (MQ=255)
                                                     cTGGAGCAGGGG‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt  <  1:2474397/68‑1 (MQ=255)
                                                             |                                                            
GCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑GGG‑GAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAT  >  W3110S.gb/547771‑547890

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: ‑