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.G100.0% 159.7 / NA 48pseudogene (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);  new base G (29/19);  total (29/19)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTT  >  W3110S.gb/547771‑547895
                                                             |                                                                
gcgcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTg                                                         <  1:1889767/71‑1 (MQ=255)
    aTGTCTGTTCGAATGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggggga                                                            <  1:1026824/63‑3 (MQ=255)
     tGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                                                    <  1:2916959/71‑1 (MQ=255)
     tGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                                                    >  1:4050690/1‑71 (MQ=255)
       tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTc                                                   <  1:2681623/70‑1 (MQ=255)
       tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTc                                                   <  1:362591/70‑1 (MQ=255)
       tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCa                                                  >  1:1644524/1‑71 (MQ=255)
       tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAg                                                                 >  1:496654/1‑55 (MQ=255)
         tGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCAt                                                 >  1:1530588/1‑70 (MQ=255)
           ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCa                                                  >  1:2276753/1‑67 (MQ=255)
             cGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCt                                             >  1:2146489/1‑70 (MQ=255)
                  gTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAc                                       <  1:3706382/71‑1 (MQ=255)
                   ttGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAcc                                      <  1:3365364/71‑1 (MQ=255)
                   ttGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAc                                       >  1:1345045/1‑70 (MQ=255)
                    tGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCa                                     >  1:3626417/1‑71 (MQ=255)
                    tGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCa                                     >  1:3056327/1‑71 (MQ=255)
                     gggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCa                                     <  1:1052265/70‑1 (MQ=255)
                        cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                                 <  1:3835874/71‑1 (MQ=255)
                           aaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAAt                              >  1:2433357/1‑71 (MQ=255)
                           aaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAAt                              >  1:928564/1‑71 (MQ=255)
                            aGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATg                             >  1:2579763/1‑71 (MQ=255)
                             gATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGc                            <  1:2719072/71‑1 (MQ=255)
                              aTGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc                           <  1:160240/71‑1 (MQ=255)
                              aTGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc                           <  1:3556361/71‑1 (MQ=255)
                               tGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGc                            >  1:1020238/1‑69 (MQ=255)
                               tGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCg                          >  1:115634/1‑71 (MQ=255)
                               tGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCg                          >  1:3066216/1‑71 (MQ=255)
                                       aGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                                 <  1:952368/56‑1 (MQ=255)
                                       aGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                                 <  1:1350298/56‑1 (MQ=255)
                                        tgcgcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCaa                 <  1:2744679/70‑1 (MQ=255)
                                        ggCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCaa                 <  1:1916613/71‑1 (MQ=255)
                                               gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg          <  1:234740/71‑1 (MQ=255)
                                               gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTAGTCACCATGTGAATGCCGTCGGGcc                   <  1:154237/62‑1 (MQ=255)
                                                cAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg          >  1:2791582/1‑70 (MQ=255)
                                                 aaTACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAt                                    >  1:718151/1‑43 (MQ=38)
                                                  aTACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTCAATGCCGTCGGGcc                   <  1:3292033/59‑1 (MQ=255)
                                                    aCTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGt                         >  1:2811722/1‑51 (MQ=255)
                                                    aCTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGt                         >  1:1935752/1‑51 (MQ=255)
                                                     cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGaa                               >  1:2137508/1‑44 (MQ=38)
                                                     cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGaa                               >  1:322538/1‑44 (MQ=38)
                                                     cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGaa                               >  1:705431/1‑44 (MQ=38)
                                                     cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAc    >  1:813744/1‑71 (MQ=255)
                                                     cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAc    >  1:2836595/1‑71 (MQ=255)
                                                     cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAc    <  1:3472854/71‑1 (MQ=255)
                                                     cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAc    >  1:2942133/1‑71 (MQ=255)
                                                     cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGACACCATGTGAAt                              >  1:574122/1‑45 (MQ=255)
                                                       ggAGCAGGGGGAAGTGAACTTCATTCGTTGTCACCATGTGAATGCCGTCGGGCCAAt                >  1:228001/1‑57 (MQ=255)
                                                       ggAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAAt                >  1:2741878/1‑57 (MQ=255)
                                                       ggAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACtt  >  1:2618544/1‑71 (MQ=255)
                                                       ggAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCAGGCCAATGGGCGGTATTACtt  >  1:1063894/1‑71 (MQ=255)
                                                             |                                                                
GCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTT  >  W3110S.gb/547771‑547895

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: