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

CACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAA  >  W3110S.gb/547642‑547756
                                                    |                                                              
cacaCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTc                                                >  1:2598876/1‑69 (MQ=255)
    ctggctgtATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGc                                            <  1:3420498/69‑1 (MQ=255)
    ctggctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCtt                                                 <  1:1994752/64‑1 (MQ=255)
    ctggctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCtt                                                 <  1:2026799/64‑1 (MQ=255)
    ctggctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGc                                            <  1:186785/69‑1 (MQ=255)
     tggctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                                         >  1:3579086/1‑55 (MQ=255)
     tggctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                                         >  1:3349067/1‑55 (MQ=255)
     tggctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                                         >  1:217410/1‑55 (MQ=255)
     tggctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCa                                                            <  1:2825882/52‑1 (MQ=255)
        ctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGc                                            <  1:422640/65‑1 (MQ=255)
          ggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCgg                                         >  1:134203/1‑66 (MQ=255)
                       gCTTCTTCACTTATCAGCGAACTAAACGAGCGCAAAACACTGCTTCACGCCGGgccgcc                                   <  1:3141248/59‑1 (MQ=39)
                             tcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGga                    <  1:584183/68‑1 (MQ=255)
                                 ttATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGa               <  1:193057/69‑1 (MQ=255)
                                      aGCCAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCgg           <  1:657546/68‑1 (MQ=255)
                                             tAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGa   <  1:2347694/69‑1 (MQ=255)
                                             tAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGa   <  1:298126/69‑1 (MQ=255)
                                             tAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGa   <  1:1765546/69‑1 (MQ=255)
                                               aaCGAGGGCAAAACACTGCTTCACGCCGGgccgcc                                   <  1:1714016/35‑1 (MQ=255)
                                               aaCGAGGGCAAAACACTGCTTCACGCCGGgccgcc                                   <  1:1044249/35‑1 (MQ=255)
                                               aaCGAGGGCAAAACACTGCTTCACGACGGGCCGCCaa                                 >  1:3477393/1‑37 (MQ=37)
                                                 cGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGaa  <  1:2029447/66‑1 (MQ=255)
                                                    |                                                              
CACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAA  >  W3110S.gb/547642‑547756

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: