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

CGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAGG  >  W3110S.gb/547634‑547759
                                                            |                                                                 
cgcTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAacac                                                         >  1:858970/1‑71 (MQ=255)
cgcTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAacac                                                         >  1:5212571/1‑71 (MQ=255)
cgcTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAacac                                                         >  1:5415958/1‑71 (MQ=255)
cgcTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACATATCAGCGAACTAAACGAGGGCAAAacac                                                         >  1:2911936/1‑71 (MQ=255)
   tCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGc                                                      >  1:4008834/1‑71 (MQ=255)
      tCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAacac                                                         >  1:4524925/1‑65 (MQ=255)
      tCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAacac                                                         >  1:4653586/1‑65 (MQ=255)
      tCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAacac                                                         >  1:4906391/1‑65 (MQ=255)
      tCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAacac                                                         >  1:549225/1‑65 (MQ=255)
        cagacTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCa                                                  <  1:5599519/67‑1 (MQ=255)
        cacaCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCa                                                  <  1:910755/70‑1 (MQ=255)
          cacTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaa                                                              <  1:2820580/56‑1 (MQ=255)
           acTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGTTTCACGcc                                              >  1:4323705/1‑71 (MQ=255)
           acTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGcc                                              >  1:2957665/1‑71 (MQ=255)
            ctggctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                                            >  1:237482/1‑56 (MQ=255)
            ctggctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                                            >  1:2361864/1‑56 (MQ=255)
            ctggctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCg                                             >  1:637802/1‑71 (MQ=255)
            ctggctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCg                                             >  1:209598/1‑71 (MQ=255)
            ctggctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCg                                             >  1:2062055/1‑71 (MQ=255)
            ctggctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCg                                             >  1:2018647/1‑71 (MQ=255)
                    aTGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCAc                                                 <  1:2740038/59‑1 (MQ=255)
                    aTGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCAc                                                 <  1:2467736/59‑1 (MQ=255)
                    aTGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCAc                                                 <  1:3060643/59‑1 (MQ=255)
                    aTGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCAc                                                 <  1:3107189/59‑1 (MQ=255)
                    aTGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCAc                                                 <  1:3292478/59‑1 (MQ=255)
                    aTGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCAc                                                 <  1:3441702/59‑1 (MQ=255)
                    aTGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCAc                                                 <  1:464762/59‑1 (MQ=255)
                         cAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATgcg                                >  1:3236869/1‑71 (MQ=255)
                           aCCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAac                                                           <  1:985985/42‑1 (MQ=255)
                           aCCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCt                              >  1:3905904/1‑71 (MQ=255)
                           aCCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCt                              >  1:906345/1‑71 (MQ=255)
                             cGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTgg                            >  1:1323046/1‑71 (MQ=255)
                             cGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTgg                            >  1:3103484/1‑71 (MQ=255)
                                  tcttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTctc                                                 <  1:4017490/45‑3 (MQ=39)
                                   cttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAg                         >  1:1440243/1‑68 (MQ=255)
                                       aCTTATCCGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATg                   <  1:3017536/70‑1 (MQ=255)
                                       aCTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAgg                        >  1:467321/1‑65 (MQ=255)
                                       aCTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATg                   <  1:4668615/70‑1 (MQ=255)
                                       aCTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATg                   <  1:3262917/70‑1 (MQ=255)
                                       aCTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATg                   <  1:5282331/70‑1 (MQ=255)
                                        cTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGAc                 <  1:3192164/71‑1 (MQ=255)
                                           aTCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCg               >  1:4026232/1‑70 (MQ=255)
                                           aTCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCg               >  1:1765291/1‑70 (MQ=255)
                                           aTCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCg               >  1:2835468/1‑70 (MQ=255)
                                           aTCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCg               >  1:5537127/1‑70 (MQ=255)
                                           aTCAGCGAACTAAACGAGGGCAAAACACTGCTTAACGCCGGGCCGCCAATGCGCTGGCAGGAGATGAcc                >  1:657147/1‑69 (MQ=255)
                                                cGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCa         >  1:5461853/1‑71 (MQ=255)
                                                cGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGCAGATGACCGGACCCa         >  1:245424/1‑71 (MQ=255)
                                                 gAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCAt        >  1:4775952/1‑71 (MQ=255)
                                                   aCTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGa      <  1:884909/71‑1 (MQ=255)
                                                       aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAgg  >  1:4488380/1‑71 (MQ=255)
                                                            |                                                                 
CGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAGG  >  W3110S.gb/547634‑547759

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