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,6940AG98.1% 170.9 / ‑4.2 52pseudogene (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 (1/0);  new base G (19/32);  total (20/32)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-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.

CGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAGGGG  >  W3110S.gb/547629‑547761
                                                                 |                                                                   
cgtcgCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCa                                                                 <  1:2362130/70‑1 (MQ=255)
 gtcgCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaa                                                                <  1:3239086/70‑1 (MQ=255)
   cgcgCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                                              >  1:1362769/1‑70 (MQ=255)
      gcTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAaca                                                            <  1:3736428/69‑1 (MQ=255)
       cTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACt                                                          <  1:2458819/70‑1 (MQ=255)
       cTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTg                                                         <  1:914918/71‑1 (MQ=255)
          gTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaa                                                                <  1:3127839/61‑1 (MQ=255)
          gTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaa                                                                <  1:3907292/61‑1 (MQ=255)
           tCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaa                                                                <  1:1391026/60‑1 (MQ=255)
             cacaCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCAc                                                   >  1:3818353/1‑71 (MQ=255)
               cacTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                                              <  1:981179/58‑1 (MQ=255)
                acTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGcc                                                <  1:2479374/71‑1 (MQ=255)
                 ctggctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCg                                               >  1:1998967/1‑71 (MQ=255)
                   ggctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAaca                                                            >  1:2205551/1‑56 (MQ=255)
                    gctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCgg                                              <  1:3743065/69‑1 (MQ=255)
                         aTGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCt                                                       >  1:2569607/1‑55 (MQ=255)
                           gtgCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACTCTGCTTCACGCCGGGCCGCCAAt                                     >  1:2190677/1‑71 (MQ=255)
                             gCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATg                                    >  1:3853805/1‑70 (MQ=255)
                              cAACCGGCTTCTTAACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATgcg                                  >  1:4148326/1‑71 (MQ=255)
                                  cGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgccgc                                         >  1:313993/1‑60 (MQ=255)
                                  cGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCa                                       >  1:4036482/1‑62 (MQ=255)
                                   ggCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTgg                              >  1:2699560/1‑70 (MQ=255)
                                   ggCTTCTTCACTTATCAGCGAACTAAACGAGGGAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTgg                              >  1:308926/1‑69 (MQ=255)
                                    gCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCa                                                    >  1:3608201/1‑47 (MQ=255)
                                     cttcttGACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCg                                               <  1:1553518/51‑1 (MQ=255)
                                        cttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTgg                              <  1:3917132/65‑1 (MQ=255)
                                        cttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAgg                          <  1:363591/69‑1 (MQ=255)
                                        cttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGga                         <  1:1183579/70‑1 (MQ=255)
                                         ttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCg                                               >  1:970001/1‑47 (MQ=255)
                                          tcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTg                               >  1:1934448/1‑62 (MQ=255)
                                          tcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGAt                      >  1:3340120/1‑71 (MQ=255)
                                           cACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGgag                        <  1:1908306/68‑1 (MQ=255)
                                             cTTATCAGCGAACTAAACGAGGGCAAAACACTGCtt                                                      <  1:320541/36‑1 (MQ=255)
                                             cTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGga                         <  1:2507883/65‑1 (MQ=255)
                                             cTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGAc                   <  1:2822802/71‑1 (MQ=255)
                                              ttATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCt                                <  1:216826/57‑1 (MQ=255)
                                                aTCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCg                 <  1:1158561/70‑1 (MQ=255)
                                                     cGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCa           >  1:2730726/1‑71 (MQ=255)
                                                     cGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCa           >  1:2068920/1‑71 (MQ=255)
                                                       aaCTAAACCAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGa                    <  1:1920895/60‑1 (MQ=255)
                                                        aCTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGa        <  1:2154985/71‑1 (MQ=255)
                                                        aCTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGa        <  1:2203356/71‑1 (MQ=255)
                                                         cTAAACGAGGGCAAAACACTGCTTCACGCCGGgccgcc                                        <  1:1305048/38‑1 (MQ=255)
                                                         cTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGaa       <  1:1322879/71‑1 (MQ=255)
                                                         cTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGaa       <  1:58718/71‑1 (MQ=255)
                                                         cTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGa        >  1:406718/1‑70 (MQ=255)
                                                           aaaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGAcc                  >  1:3381617/1‑58 (MQ=255)
                                                            aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGc                             <  1:456805/46‑1 (MQ=255)
                                                            aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGAccc            <  1:2226244/63‑1 (MQ=255)
                                                            aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAg     <  1:3240741/70‑1 (MQ=255)
                                                             aCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAgg    <  1:4240133/70‑1 (MQ=255)
                                                              cGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAgggg  <  1:1979794/71‑1 (MQ=255)
                                                                 |                                                                   
CGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAGGGG  >  W3110S.gb/547629‑547761

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: