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.G95.0% 126.6 / ‑1.0 40pseudogene (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/2);  new base G (17/21);  total (17/23)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

GCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAC  >  W3110S.gb/547761‑547893
                                                                       |                                                                
gcgtgcgtGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCa                                                                   <  1:2831291/71‑1 (MQ=255)
  gtgcgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgg                                                                >  1:3566459/1‑71 (MQ=255)
     cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCTGGCGCTGGCAATACTGGAGCAgggg                                                            >  1:859205/1‑70 (MQ=255)
     cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgggg                                                            >  1:3875944/1‑70 (MQ=255)
     cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgggg                                                            >  1:2274222/1‑70 (MQ=255)
      gtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggggg                                                           >  1:2454616/1‑69 (MQ=255)
      gtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggggg                                                           >  1:2449164/1‑69 (MQ=255)
               tGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTg                                                       <  1:876204/66‑1 (MQ=255)
               tGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                                                  <  1:3281078/71‑1 (MQ=255)
               tGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                                                  >  1:2928088/1‑71 (MQ=255)
                 tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCa                                                <  1:3660080/71‑1 (MQ=255)
                 tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCa                                                >  1:1316543/1‑71 (MQ=255)
                   tGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCAt                                               >  1:465592/1‑70 (MQ=255)
                     ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTc                                             <  1:611123/70‑1 (MQ=255)
                       cGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCtt                                          >  1:3080263/1‑71 (MQ=255)
                            gTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAc                                     <  1:1042782/71‑1 (MQ=255)
                              tGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAg                                                                 >  1:3035499/1‑42 (MQ=255)
                                ggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGTGCAggggg                                                           <  1:237394/44‑2 (MQ=255)
                                  cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGTGGAAGTGAACTTCATTCCTTGTCACCAtgtg                               <  1:1185101/71‑1 (MQ=255)
                                   gAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                               <  1:4322588/70‑1 (MQ=255)
                                     aaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAAt                            <  1:2106952/71‑1 (MQ=255)
                                        aTGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc                         <  1:2057386/71‑1 (MQ=255)
                                        aTGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc                         >  1:3012564/1‑71 (MQ=255)
                                          gAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGt                       >  1:1930162/1‑71 (MQ=255)
                                          gAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGt                       >  1:2899962/1‑71 (MQ=255)
                                           aaGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTc                      >  1:2589897/1‑71 (MQ=255)
                                           aaGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAAATTCATTCCTTGt                                        >  1:442320/1‑53 (MQ=38)
                                            aGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCg                     <  1:1205771/71‑1 (MQ=255)
                                                cAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc                 <  1:2820154/71‑1 (MQ=255)
                                                  ggCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCaa               <  1:1547981/71‑1 (MQ=255)
                                                    cgcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATg             <  1:1955278/71‑1 (MQ=255)
                                                    cgcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATg             <  1:3352112/71‑1 (MQ=255)
                                                     gcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtg                                 >  1:3189984/1‑50 (MQ=39)
                                                      cTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAcc                                    <  1:2517227/46‑1 (MQ=39)
                                                        ggCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCg         >  1:4204479/1‑71 (MQ=255)
                                                         gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg        >  1:3402093/1‑71 (MQ=255)
                                                         gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg        >  1:2875596/1‑71 (MQ=255)
                                                            aTACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt     <  1:4287844/71‑1 (MQ=255)
                                                             tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                               <  1:1209432/44‑1 (MQ=38)
                                                             tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                               <  1:4419894/44‑1 (MQ=38)
                                                             tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTc                      <  1:700613/53‑1 (MQ=255)
                                                             tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTc                      <  1:4473427/53‑1 (MQ=255)
                                                              aCTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGc          >  1:3917348/1‑64 (MQ=255)
                                                              aCTGGAGCAGGGTGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTc                      <  1:2974683/52‑1 (MQ=255)
                                                               cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGaa                             >  1:1136175/1‑44 (MQ=38)
                                                               cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGt                       >  1:2929024/1‑50 (MQ=255)
                                                               cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAc  <  1:3096003/71‑1 (MQ=255)
                                                               cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAc  <  1:2936806/71‑1 (MQ=255)
                                                               cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAc  >  1:1198626/1‑71 (MQ=255)
                                                                       |                                                                
GCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAC  >  W3110S.gb/547761‑547893

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: