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

GCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATT  >  W3110S.gb/547771‑547891
                                                             |                                                            
gcgcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTg                                                     <  1:1655641/71‑1 (MQ=255)
   cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAAc                                                  <  1:255259/71‑1 (MQ=255)
    aTGTCTGTTTGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACt                                                 >  1:5175174/1‑71 (MQ=255)
    aTGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACt                                                 <  1:57073/71‑1 (MQ=255)
    aTGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACt                                                 <  1:4500985/71‑1 (MQ=255)
       tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTAGAGCAGGGGGAAGTGAACTTCa                                              >  1:4192738/1‑71 (MQ=255)
         tGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTg                                                     <  1:618204/62‑1 (MQ=255)
           ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCAt                                             <  1:3154971/68‑1 (MQ=255)
           ttCGAAGGTTGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTcc                                          <  1:708090/70‑1 (MQ=255)
              gAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCtt                                        <  1:4430239/70‑1 (MQ=255)
                aGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGt                                      >  1:1548497/1‑70 (MQ=255)
                    tGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCa                                 >  1:871542/1‑71 (MQ=255)
                       gCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgt                              >  1:5110920/1‑71 (MQ=255)
                               tGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCg                      <  1:4646260/71‑1 (MQ=255)
                                gAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGt                     >  1:4539634/1‑71 (MQ=255)
                                  aGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgt                              >  1:891942/1‑60 (MQ=255)
                                        ggCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCaa             <  1:4380163/71‑1 (MQ=255)
                                        ggCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCaa             <  1:4271594/71‑1 (MQ=255)
                                        ggCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCaa             <  1:1116333/71‑1 (MQ=255)
                                         gcgcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgt                              <  1:1200854/53‑1 (MQ=39)
                                          cgcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATg           >  1:3262809/1‑71 (MQ=255)
                                           gcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATgg          >  1:4236602/1‑71 (MQ=255)
                                             tGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCt                                         >  1:3552260/1‑38 (MQ=37)
                                              ggCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCa                                 >  1:3606510/1‑45 (MQ=255)
                                              ggCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCa                                 >  1:2829437/1‑45 (MQ=255)
                                              ggCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCa                                 >  1:2788080/1‑45 (MQ=255)
                                               gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCa              <  1:1881434/63‑1 (MQ=255)
                                               gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg      >  1:1686483/1‑71 (MQ=255)
                                               gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg      <  1:2250230/71‑1 (MQ=255)
                                                cAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGt     <  1:147970/71‑1 (MQ=255)
                                                 aaTACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAAt            >  1:2110589/1‑63 (MQ=255)
                                                 aaTACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAAt            >  1:3067996/1‑63 (MQ=255)
                                                   tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                             <  1:5090352/44‑1 (MQ=38)
                                                   tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc               <  1:2446440/58‑1 (MQ=255)
                                                   tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCaa             >  1:1987297/1‑60 (MQ=255)
                                                   tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAtt  <  1:4445467/71‑1 (MQ=255)
                                                   tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAtt  <  1:2154598/71‑1 (MQ=255)
                                                   tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt   >  1:4548012/1‑70 (MQ=255)
                                                    aCTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAtt  >  1:3851468/1‑70 (MQ=255)
                                                     cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                             >  1:472513/1‑42 (MQ=38)
                                                     cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCg                      <  1:234537/49‑1 (MQ=255)
                                                     cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg      <  1:336141/65‑1 (MQ=255)
                                                             |                                                            
GCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATT  >  W3110S.gb/547771‑547891

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: