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.G98.0% 176.5 / ‑4.3 49pseudogene (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/0);  major base G (42/6);  minor base T (1/0);  total (43/6)
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.92e-01

CATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCG  >  W3110S.gb/547774‑547874
                                                          |                                           
cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgggggaa                                      >  1:788483/1‑62 (MQ=255)
cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgggggaa                                      >  1:1121833/1‑62 (MQ=255)
cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgggggaa                                      >  1:1152056/1‑62 (MQ=255)
cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgggggaa                                      >  1:2769414/1‑62 (MQ=255)
    tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgggggt                                       <  1:2452340/60‑3 (MQ=255)
    tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggggga                                       <  1:85868/60‑3 (MQ=255)
    tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggggga                                       <  1:2463480/60‑3 (MQ=255)
      tGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggggga                                       <  1:697781/58‑3 (MQ=255)
        ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                               >  1:2398875/1‑65 (MQ=255)
        ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                               >  1:2437900/1‑65 (MQ=255)
        ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                               >  1:2275332/1‑65 (MQ=255)
        ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                               >  1:2517181/1‑65 (MQ=255)
        ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                               >  1:2616803/1‑65 (MQ=255)
        ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                               >  1:2651104/1‑65 (MQ=255)
        ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                               >  1:1777332/1‑65 (MQ=255)
        ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                               >  1:167018/1‑65 (MQ=255)
        ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                               >  1:1514984/1‑65 (MQ=255)
        ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                               >  1:625768/1‑65 (MQ=255)
        ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                               >  1:1287331/1‑65 (MQ=255)
        ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                               >  1:1146400/1‑65 (MQ=255)
        ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                               >  1:1114994/1‑65 (MQ=255)
        ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                               >  1:798474/1‑65 (MQ=255)
                     cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAc                  >  1:1599088/1‑65 (MQ=255)
                     cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAc                  >  1:3012773/1‑65 (MQ=255)
                     cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAc                  >  1:1738746/1‑65 (MQ=255)
                     cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAc                  >  1:1406427/1‑65 (MQ=255)
                     cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAc                  >  1:2563755/1‑65 (MQ=255)
                     cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAc                  >  1:1260092/1‑65 (MQ=255)
                       aaaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCa                <  1:473116/65‑1 (MQ=255)
                       aaaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCa                <  1:1030965/65‑1 (MQ=255)
                       aaaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCa                <  1:2146212/65‑1 (MQ=255)
                       aaaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCa                <  1:1683153/65‑1 (MQ=255)
                                 cgcAGGCGCTGGCAATACTGGAGCATGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc      >  1:2607953/1‑65 (MQ=255)
                                 cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTc                              >  1:2851079/1‑41 (MQ=38)
                                 cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc      >  1:3054170/1‑65 (MQ=255)
                                 cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc      >  1:1024705/1‑65 (MQ=255)
                                 cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc      >  1:842609/1‑65 (MQ=255)
                                 cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc      >  1:107138/1‑65 (MQ=255)
                                 cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc      >  1:1075730/1‑65 (MQ=255)
                                 cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc      >  1:1221905/1‑65 (MQ=255)
                                 cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc      >  1:1402212/1‑65 (MQ=255)
                                 cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc      >  1:588193/1‑65 (MQ=255)
                                 cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc      >  1:558036/1‑65 (MQ=255)
                                 cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc      >  1:1471978/1‑65 (MQ=255)
                                 cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc      >  1:43904/1‑65 (MQ=255)
                                 cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc      >  1:3087470/1‑65 (MQ=255)
                                 cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc      >  1:2289757/1‑65 (MQ=255)
                                 cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc      >  1:1580687/1‑65 (MQ=255)
                                 cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc      >  1:3011309/1‑65 (MQ=255)
                                 cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc      >  1:2729823/1‑65 (MQ=255)
                                 cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc      >  1:1828595/1‑65 (MQ=255)
                                 cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc      >  1:1908030/1‑65 (MQ=255)
                                 cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc      >  1:2014567/1‑65 (MQ=255)
                                 cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc      >  1:202613/1‑65 (MQ=255)
                                 cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc      >  1:2201630/1‑65 (MQ=255)
                                     ggCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCg  <  1:3084730/65‑1 (MQ=255)
                                     ggCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCg  <  1:497531/65‑1 (MQ=255)
                                                          |                                           
CATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCG  >  W3110S.gb/547774‑547874

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