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% 129.9 / NA 40pseudogene (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 (23/17);  total (23/17)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAC  >  W3110S.gb/547764‑547893
                                                                    |                                                              
tgcgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggg                                                             <  1:2555290/71‑1 (MQ=255)
  cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggggg                                                           <  1:4247034/71‑2 (MQ=255)
  cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggggg                                                           <  1:1585715/71‑2 (MQ=255)
  cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggggg                                                           <  1:333480/71‑2 (MQ=255)
  cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggggg                                                           <  1:3321898/71‑2 (MQ=255)
   gtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggggga                                                          <  1:2463724/71‑3 (MQ=255)
    tgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCTCTGGCAATACTGGAGCAggggga                                                          <  1:2879795/70‑3 (MQ=255)
           aTGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACt                                                   >  1:1791195/1‑71 (MQ=255)
              tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCa                                                >  1:4241554/1‑71 (MQ=255)
              tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCa                                                >  1:3193156/1‑71 (MQ=255)
                  ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTcc                                            >  1:3756451/1‑71 (MQ=255)
                  ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggggg                                                           <  1:1400401/55‑2 (MQ=255)
                  ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggggg                                                           <  1:2914977/55‑2 (MQ=255)
                    cGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGa                                                      >  1:1811467/1‑59 (MQ=39)
                    cGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCa                                                >  1:2141372/1‑65 (MQ=255)
                     gAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGAAGGGGGAAGTGAACTTCATTCCtt                                          <  1:2719901/70‑1 (MQ=255)
                       aGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTc                                       <  1:4366979/71‑1 (MQ=255)
                          ttGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAcc                                    <  1:4134574/71‑1 (MQ=255)
                           tGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTcc                                            <  1:249441/62‑1 (MQ=255)
                              gCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACAAtgt                                <  1:2807485/71‑1 (MQ=255)
                               cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgt                                <  1:4259574/70‑1 (MQ=255)
                                   aGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATg                           >  1:933069/1‑71 (MQ=255)
                                      tGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCg                        <  1:3530703/71‑1 (MQ=255)
                                       gAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCtt                                          >  1:2066933/1‑52 (MQ=39)
                                       gAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTAATTCCtt                                          >  1:485202/1‑52 (MQ=38)
                                          gcgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCgg                    <  1:1843260/71‑1 (MQ=255)
                                           cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCggg                   <  1:2549932/71‑1 (MQ=255)
                                           cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCggg                   >  1:1942540/1‑71 (MQ=255)
                                             cAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc                 <  1:2265912/71‑1 (MQ=255)
                                             cAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc                 <  1:2988015/71‑1 (MQ=255)
                                               ggCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCaa               <  1:2599419/71‑1 (MQ=255)
                                                gcgcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                               <  1:3718132/54‑1 (MQ=39)
                                                  gcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATgg            >  1:3763726/1‑71 (MQ=255)
                                                   cTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGTGCCAATggg           >  1:973556/1‑71 (MQ=255)
                                                   cTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAAt              >  1:2360195/1‑68 (MQ=255)
                                                   cTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAAt              >  1:2207148/1‑68 (MQ=255)
                                                    tGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGc          >  1:554394/1‑71 (MQ=255)
                                                     ggCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGaa                             >  1:2986957/1‑51 (MQ=255)
                                                      gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGc          <  1:4065579/69‑1 (MQ=255)
                                                      gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg        >  1:2037464/1‑71 (MQ=255)
                                                       cAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAc                                     <  1:2264033/41‑1 (MQ=38)
                                                        aaTACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCa                >  1:4245002/1‑61 (MQ=255)
                                                        aaTACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTa      >  1:1957793/1‑71 (MQ=255)
                                                          tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAAt                            <  1:146583/47‑1 (MQ=255)
                                                           aCTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCa                                   >  1:1482318/1‑39 (MQ=38)
                                                           aCTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAtt    >  1:79956/1‑70 (MQ=255)
                                                            cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTa      >  1:611673/1‑67 (MQ=255)
                                                            cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAc  <  1:4225481/71‑1 (MQ=255)
                                                               gAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCa                >  1:139296/1‑54 (MQ=255)
                                                                    |                                                              
TGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAC  >  W3110S.gb/547764‑547893

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