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

TCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAGG  >  W3110S.gb/547631‑547759
                                                               |                                                                 
tcgCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                                            <  1:3789215/71‑1 (MQ=255)
tcgCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                                            <  1:1637935/71‑1 (MQ=255)
   cgcTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAacac                                                         >  1:3168886/1‑71 (MQ=255)
    gcTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACt                                                        >  1:1027075/1‑71 (MQ=255)
        gTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCtt                                                    <  1:872679/71‑1 (MQ=255)
         tCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTc                                                   <  1:4037935/71‑1 (MQ=255)
            acacTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACg                                                <  1:238456/71‑1 (MQ=255)
                      gATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgccgcc                                      <  1:3031432/71‑1 (MQ=255)
                       aTGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCa                                     <  1:3151984/71‑1 (MQ=255)
                       aTGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCa                                     <  1:1968286/71‑1 (MQ=255)
                         gtgCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCt                                                     <  1:1303088/53‑1 (MQ=255)
                            cAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATgc                                 <  1:1516887/70‑1 (MQ=255)
                              aCCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGc                                                      <  1:3677785/47‑1 (MQ=255)
                              aCCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCa                                     <  1:1441066/64‑1 (MQ=255)
                              aCCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCt                              <  1:4339522/71‑1 (MQ=255)
                                cGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCaa                                    >  1:820196/1‑63 (MQ=255)
                                cGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTgg                            <  1:3870283/71‑1 (MQ=255)
                                 ggCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCATc                                                   >  1:3133014/1‑47 (MQ=255)
                                  gCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                                            <  1:1898692/37‑1 (MQ=255)
                                  gCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                                            <  1:2231922/37‑1 (MQ=255)
                                   cttcttCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCa                          >  1:903500/1‑70 (MQ=255)
                                   cttcttCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAg                         <  1:3146467/71‑1 (MQ=255)
                                     tcttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCtt                                                    <  1:2714082/42‑1 (MQ=255)
                                     tcttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAg                         >  1:2527582/1‑69 (MQ=255)
                                       ttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAg                         >  1:4584730/1‑67 (MQ=255)
                                        tcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGAt                    >  1:3602064/1‑71 (MQ=255)
                                         cACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCaa                                    >  1:793617/1‑54 (MQ=255)
                                          aCTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATg                   <  1:1558590/70‑1 (MQ=255)
                                          aCTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATg                   <  1:1885145/70‑1 (MQ=255)
                                          aCTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGa                  <  1:2040912/71‑1 (MQ=255)
                                           cTTATCAGCGAACTAAACGAGGGCAAAACACTGCtt                                                    <  1:4537914/36‑1 (MQ=255)
                                           cTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGAt                    <  1:2612791/68‑1 (MQ=255)
                                           cTTATCAGCGAACTAAACGAGGGCAAAACACAGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATg                   <  1:2108636/69‑1 (MQ=255)
                                            ttATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGga                       >  1:3360970/1‑64 (MQ=255)
                                             tATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCaa                                    <  1:2413468/50‑1 (MQ=255)
                                             tATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGa                  >  1:3077253/1‑68 (MQ=255)
                                              aTCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGgaga                     >  1:4618915/1‑64 (MQ=255)
                                              aTCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCg               >  1:3233467/1‑70 (MQ=255)
                                                    gAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCAt        >  1:2821403/1‑71 (MQ=255)
                                                     aaCTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCaa                                    <  1:1559721/42‑1 (MQ=255)
                                                      aCTAAACGAGGGCAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGa      <  1:451366/70‑1 (MQ=255)
                                                       cTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGga                       <  1:1852331/53‑1 (MQ=255)
                                                         aaaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAg   >  1:651419/1‑71 (MQ=255)
                                                          aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTgg                            <  1:2846097/45‑1 (MQ=255)
                                                          aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGaaa    <  1:1868034/69‑1 (MQ=255)
                                                          aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAgg  <  1:2823824/71‑1 (MQ=255)
                                                               |                                                                 
TCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAGG  >  W3110S.gb/547631‑547759

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: