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

GCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAT  >  W3110S.gb/547771‑547890
                                                             |                                                            
gcgcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGTGAAGTg                                                    <  1:3058/71‑1 (MQ=255)
   cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAAc                                                 <  1:1388174/71‑1 (MQ=255)
   cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAAc                                                 <  1:481687/71‑1 (MQ=255)
   cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAAc                                                 <  1:3449605/71‑1 (MQ=255)
   cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAAc                                                 <  1:2596054/71‑1 (MQ=255)
   cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAAc                                                 <  1:674442/71‑1 (MQ=255)
         tGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCAtt                                           >  1:502765/1‑71 (MQ=255)
         tGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCAtt                                           >  1:19023/1‑71 (MQ=255)
         tGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCAtt                                           >  1:2771568/1‑71 (MQ=255)
         tGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCAt                                            >  1:1705612/1‑70 (MQ=255)
                 ggTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAg                                                             >  1:2305208/1‑45 (MQ=255)
                   ttGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAcc                                 <  1:3865089/71‑1 (MQ=255)
                   ttGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAcc                                 <  1:713975/71‑1 (MQ=255)
                   ttGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAcc                                 <  1:3471555/71‑1 (MQ=255)
                      ggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAt                               <  1:679854/70‑1 (MQ=255)
                      ggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAt                               <  1:3713231/70‑1 (MQ=255)
                      ggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAt                               <  1:2850138/70‑1 (MQ=255)
                      ggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAt                               <  1:2998943/70‑1 (MQ=255)
                      ggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAt                               <  1:3155072/70‑1 (MQ=255)
                       gCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAt                               <  1:1567788/69‑1 (MQ=255)
                               tGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACt                                                >  1:2169057/1‑44 (MQ=38)
                               tGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACt                                                >  1:3060881/1‑44 (MQ=38)
                               tGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACt                                                >  1:2318940/1‑44 (MQ=38)
                               tGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACt                                                >  1:2072930/1‑44 (MQ=38)
                                gAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGt                    >  1:625732/1‑71 (MQ=255)
                                gAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGt                    >  1:250001/1‑71 (MQ=255)
                                    cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCggg                <  1:3650482/71‑1 (MQ=255)
                                     gcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCggg                <  1:730496/70‑1 (MQ=255)
                                       aGGCTCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAcc                                 <  1:1530093/51‑1 (MQ=255)
                                            cTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc              <  1:565391/65‑1 (MQ=255)
                                            cTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc              <  1:2432292/65‑1 (MQ=255)
                                            cTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc              <  1:1990463/65‑1 (MQ=255)
                                                     cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt  <  1:479268/68‑1 (MQ=255)
                                                     cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt  <  1:415663/68‑1 (MQ=255)
                                                     cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt  <  1:54336/68‑1 (MQ=255)
                                                     cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt  <  1:267293/68‑1 (MQ=255)
                                                     cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt  <  1:2474397/68‑1 (MQ=255)
                                                             |                                                            
GCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAT  >  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: