Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 547,831:1 +G 100% pseudogene (251/261 nt)
pseudogene (251/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% 173.8 / ‑3.8 51pseudogene (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 (34/16);  total (34/17)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

CGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTT  >  W3110S.gb/547766‑547895
                                                                  |                                                                 
cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggggg                                                              <  1:364256/71‑2 (MQ=255)
cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggggg                                                              <  1:785637/71‑2 (MQ=255)
cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgggg                                                               >  1:1115433/1‑70 (MQ=255)
cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgggg                                                               >  1:4877629/1‑70 (MQ=255)
cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgggg                                                               >  1:1953576/1‑70 (MQ=255)
  tgGGCGCATGTCTGTTCGAAGTTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggggg                                                              <  1:4713806/69‑2 (MQ=255)
    ggCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCa                                                                    <  1:4576640/62‑1 (MQ=255)
    ggCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCa                                                                    <  1:2630711/62‑1 (MQ=255)
    ggCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCa                                                                    <  1:3225807/62‑1 (MQ=255)
    ggCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCa                                                                    <  1:3264972/62‑1 (MQ=255)
    ggCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCa                                                                    <  1:3335550/62‑1 (MQ=255)
    ggCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCa                                                                    <  1:2312093/62‑1 (MQ=255)
    ggCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCa                                                                    <  1:2304015/62‑1 (MQ=255)
    ggCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCa                                                                    <  1:1983754/62‑1 (MQ=255)
     gcgcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCa                                                                    <  1:489940/61‑1 (MQ=255)
     gcgcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTg                                                         <  1:1340552/71‑1 (MQ=255)
       gcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGaa                                                       >  1:1878517/1‑71 (MQ=255)
       gcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGaa                                                       >  1:740414/1‑71 (MQ=255)
           gTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                                                    >  1:1121201/1‑70 (MQ=255)
           gTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                                                    >  1:4491733/1‑70 (MQ=255)
                          gggCGAAATATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGTAAGTGAACTTCATTCCTTGTCACCa                                     <  1:5148793/70‑1 (MQ=255)
                          gggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCa                                     <  1:2358198/70‑1 (MQ=255)
                             cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                                 <  1:823455/71‑1 (MQ=255)
                             cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                                 <  1:923411/71‑1 (MQ=255)
                               aaaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGaa                               >  1:4781397/1‑71 (MQ=255)
                               aaaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGaa                               >  1:2303645/1‑71 (MQ=255)
                               aaaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGaa                               >  1:1535404/1‑71 (MQ=255)
                                   aTGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                                 <  1:2640893/65‑1 (MQ=255)
                                    tGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                                 <  1:2882416/64‑1 (MQ=255)
                                     gAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGGGAACTTCATTCCTTg                                           <  1:2815633/53‑1 (MQ=38)
                                       aGCGCCGGCGCTGGCAATACTGGAGAAGGGGGAAGTGAACTTCATTCCTTGTc                                         >  1:1713944/1‑53 (MQ=37)
                                       aGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTc                                         >  1:2054849/1‑53 (MQ=39)
                                       aGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTc                                         >  1:2466417/1‑53 (MQ=39)
                                       aGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTc                                         >  1:3852680/1‑53 (MQ=39)
                                       aGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTc                                         >  1:4086437/1‑53 (MQ=39)
                                           cAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc                   <  1:3849696/71‑1 (MQ=255)
                                           cAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc                   <  1:4348328/71‑1 (MQ=255)
                                                gcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATgg              <  1:2939995/71‑1 (MQ=255)
                                                gcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATgg              >  1:5121628/1‑71 (MQ=255)
                                                gcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATgg              <  1:856327/71‑1 (MQ=255)
                                                 cTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATgg              <  1:5187924/70‑1 (MQ=255)
                                                    gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg          >  1:3083145/1‑71 (MQ=255)
                                                    gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg          <  1:3259340/71‑1 (MQ=255)
                                                    gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg          <  1:5116013/71‑1 (MQ=255)
                                                    gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg          <  1:3303026/71‑1 (MQ=255)
                                                     cAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg          >  1:602836/1‑70 (MQ=255)
                                                       aTACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt       <  1:1062820/71‑1 (MQ=255)
                                                        tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt       >  1:599083/1‑70 (MQ=255)
                                                        tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt       >  1:4238031/1‑70 (MQ=255)
                                                        tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt       >  1:1676236/1‑70 (MQ=255)
                                                        tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt       >  1:1820202/1‑70 (MQ=255)
                                                        tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt       >  1:1442373/1‑70 (MQ=255)
                                                        tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt       >  1:2093612/1‑70 (MQ=255)
                                                        tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt       >  1:4477088/1‑70 (MQ=255)
                                                        tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt       >  1:2379685/1‑70 (MQ=255)
                                                        tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt       >  1:242241/1‑70 (MQ=255)
                                                        tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt       >  1:2660797/1‑70 (MQ=255)
                                                         aCTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGc            >  1:3183817/1‑64 (MQ=255)
                                                         aCTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTa     >  1:2207107/1‑71 (MQ=255)
                                                         aCTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTa     >  1:3854028/1‑71 (MQ=255)
                                                         aCTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTa     >  1:25352/1‑71 (MQ=255)
                                                         aCTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATCCCGt                         >  1:5001129/1‑51 (MQ=255)
                                                            ggAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACtt  >  1:2265482/1‑71 (MQ=255)
                                                            ggAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACtt  >  1:2557672/1‑71 (MQ=255)
                                                            ggAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACtt  >  1:399415/1‑71 (MQ=255)
                                                            ggAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACtt  >  1:3201803/1‑71 (MQ=255)
                                                                  |                                                                 
CGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTT  >  W3110S.gb/547766‑547895

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