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

GCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTC  >  W3110S.gb/547761‑547896
                                                                       |                                                                 
gcgtgcgtGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCa                                                                    <  1:408276/71‑1 (MQ=255)
   tgcgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggg                                                                >  1:1213429/1‑71 (MQ=255)
       tgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgggggaa                                                            <  1:1030441/71‑4 (MQ=255)
            gcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTg                                                          >  1:2293097/1‑69 (MQ=255)
             cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAAc                                                       <  1:3142414/71‑1 (MQ=255)
             cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAAc                                                       >  1:2365447/1‑71 (MQ=255)
                gTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                                                     >  1:1270005/1‑70 (MQ=255)
                 tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                                                     <  1:4302682/69‑1 (MQ=255)
                 tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCa                                                   >  1:396450/1‑71 (MQ=255)
                 tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCa                                                   >  1:1352419/1‑71 (MQ=255)
                  cTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGaa                                                        >  1:4474725/1‑65 (MQ=255)
                   tGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTg                                                          <  1:552817/62‑1 (MQ=255)
                      tCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCAt                                                  <  1:4472036/67‑1 (MQ=255)
                        gAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCAt                                                  <  1:2943895/65‑1 (MQ=255)
                        gAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCAt                                                  <  1:2147691/65‑1 (MQ=255)
                        gAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCAt                                                  <  1:2948947/65‑1 (MQ=255)
                        gAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTc                                                <  1:1485490/67‑1 (MQ=255)
                           ggTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGt                                           <  1:2997040/69‑1 (MQ=255)
                           ggTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCa                                         >  1:1147951/1‑71 (MQ=255)
                                ggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAt                                     <  1:2892967/70‑1 (MQ=255)
                                 gCGAAAGAAGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgt                                   <  1:4519925/71‑1 (MQ=255)
                                  cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                                  <  1:384650/71‑1 (MQ=255)
                                   gAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGa                                 <  1:111145/71‑1 (MQ=255)
                                    aaaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGt                                           >  1:2018500/1‑60 (MQ=255)
                                    aaaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGt                                           >  1:1180143/1‑60 (MQ=255)
                                      aGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGaa                                >  1:4355155/1‑69 (MQ=255)
                                      aGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAAt                               >  1:314366/1‑70 (MQ=255)
                                        aTGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGaa                                                        >  1:3265320/1‑43 (MQ=38)
                                        aTGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGa                                 >  1:1133225/1‑66 (MQ=255)
                                         tGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGa                                                         <  1:2224384/41‑1 (MQ=38)
                                             gcgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                                  <  1:2738589/60‑1 (MQ=255)
                                              cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCggg                      >  1:567036/1‑71 (MQ=255)
                                                cAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc                    <  1:2696959/71‑1 (MQ=255)
                                                 aGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAcc                                       >  1:3303609/1‑51 (MQ=39)
                                                 aGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCa                   <  1:4204846/71‑1 (MQ=255)
                                                  ggCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                                                     <  1:3046775/36‑1 (MQ=37)
                                                   gcgcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAt                                     <  1:4330883/51‑1 (MQ=39)
                                                     gcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgt                                   >  1:2265313/1‑51 (MQ=39)
                                                      cTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTg                                            <  1:4016430/41‑1 (MQ=38)
                                                        ggCAATACTGGAGCAGGGGGAAGTGAACTTCATTcc                                               <  1:245317/36‑1 (MQ=37)
                                                         gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg           >  1:2535659/1‑71 (MQ=255)
                                                           aaTACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                                  <  1:1592553/46‑1 (MQ=255)
                                                            aTACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt        <  1:3241796/71‑1 (MQ=255)
                                                             tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAAt                               >  1:2704623/1‑47 (MQ=255)
                                                             tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg           <  1:3445955/67‑1 (MQ=255)
                                                             tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt        >  1:2929985/1‑70 (MQ=255)
                                                               cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg           <  1:2669929/65‑1 (MQ=255)
                                                                 ggAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACtt   >  1:337874/1‑71 (MQ=255)
                                                                 ggAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACt    >  1:2871332/1‑70 (MQ=255)
                                                                 ggAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACt    >  1:3379172/1‑70 (MQ=255)
                                                                  gAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTc  >  1:2705437/1‑71 (MQ=255)
                                                                       |                                                                 
GCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTC  >  W3110S.gb/547761‑547896

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: