Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 547,694 A→G 100% pseudogene (114/261 nt)
pseudogene (114/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,6940AG100.0% 166.4 / NA 49pseudogene (114/261 nt)
pseudogene (114/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 A (0/0);  new base G (26/23);  total (26/23)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAG  >  W3110S.gb/547628‑547758
                                                                  |                                                                
tcgtcgCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGc                                                               >  1:3087157/1‑70 (MQ=255)
  gtcgCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaa                                                             <  1:2370823/70‑1 (MQ=255)
  gtcgCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaa                                                             <  1:1718869/70‑1 (MQ=255)
           gTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTg                                                      <  1:2312187/68‑1 (MQ=255)
           gTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTg                                                      <  1:1772940/68‑1 (MQ=255)
              cacaCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCAc                                                >  1:3432093/1‑71 (MQ=255)
                 acTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGcc                                             >  1:214803/1‑71 (MQ=255)
                 acTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGc                                              >  1:1297225/1‑70 (MQ=255)
                  ctggctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                                           <  1:744047/56‑1 (MQ=255)
                  ctggctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCg                                            >  1:1407119/1‑71 (MQ=255)
                  ctggctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCg                                            <  1:2644794/71‑1 (MQ=255)
                   tggctggATGTGCAACCTGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCgg                                           >  1:3562054/1‑71 (MQ=255)
                        ggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGc                                                     <  1:731446/56‑1 (MQ=255)
                          aTGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCa                                    <  1:1871577/71‑1 (MQ=255)
                           tgtgCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCaa                                   <  1:1954306/71‑1 (MQ=255)
                           tgtgCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCaa                                   <  1:247940/71‑1 (MQ=255)
                             tgCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAAt                                  <  1:2491058/70‑1 (MQ=255)
                              gCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATgc                                >  1:3231479/1‑71 (MQ=255)
                              gCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATg                                 >  1:376340/1‑70 (MQ=255)
                                  ccGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTg                            >  1:1089694/1‑71 (MQ=255)
                                  ccGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTg                            >  1:3645693/1‑71 (MQ=255)
                                  ccGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTg                            >  1:2746378/1‑71 (MQ=255)
                                     gCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTgg                           >  1:3099278/1‑69 (MQ=255)
                                         cttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGcc                                             >  1:2197788/1‑47 (MQ=255)
                                         cttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGga                      <  1:675717/70‑1 (MQ=255)
                                         cttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGga                      <  1:2251950/70‑1 (MQ=255)
                                         cttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGga                      <  1:2029532/70‑1 (MQ=255)
                                         cttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGga                      <  1:1948037/70‑1 (MQ=255)
                                          ttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAg                        >  1:1526323/1‑67 (MQ=255)
                                           tcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGgaga                    >  1:3626985/1‑70 (MQ=255)
                                            cACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGAt                   >  1:2609132/1‑70 (MQ=255)
                                             aCTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATg                  <  1:223964/70‑1 (MQ=255)
                                                  tCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGTGCCGCCAATGCGCTGGCAGGAGATGACCGGa            >  1:3562089/1‑71 (MQ=255)
                                                    aGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCgg                                           <  1:193983/38‑1 (MQ=255)
                                                         aCTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATg      >  1:2075312/1‑70 (MQ=255)
                                                          cTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTgg                           <  1:351825/48‑1 (MQ=255)
                                                          cTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGaa    <  1:1537834/71‑1 (MQ=255)
                                                          cTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGa     >  1:2241414/1‑70 (MQ=255)
                                                           tAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGga                      <  1:562751/52‑1 (MQ=255)
                                                            aaaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGc                          <  1:2428318/47‑1 (MQ=255)
                                                            aaaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGga                      >  1:3546244/1‑51 (MQ=255)
                                                            aaaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGga                      >  1:1489974/1‑51 (MQ=255)
                                                            aaaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAg  >  1:2653943/1‑71 (MQ=255)
                                                            aaaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAg  >  1:2421765/1‑71 (MQ=255)
                                                             aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCt                             >  1:3421367/1‑43 (MQ=255)
                                                             aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCa                         >  1:2048158/1‑47 (MQ=255)
                                                             aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATg                  >  1:1589327/1‑54 (MQ=255)
                                                             aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGaaa   <  1:1670317/69‑1 (MQ=255)
                                                             aaCGAGGGCAAAACACTGCTACACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGaaa   <  1:1737052/69‑1 (MQ=255)
                                                                  |                                                                
TCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAG  >  W3110S.gb/547628‑547758

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: