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.2% 169.3 / ‑4.5 54pseudogene (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/1);  new base G (30/23);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-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.

CGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAG  >  W3110S.gb/547629‑547758
                                                                 |                                                                
cgtcgCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCa                                                              >  1:1066540/1‑70 (MQ=255)
     cgcTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                                           <  1:1425343/68‑1 (MQ=255)
          gTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCtt                                                   <  1:694392/71‑1 (MQ=255)
          gTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCtt                                                   <  1:4332163/71‑1 (MQ=255)
          gTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCtt                                                   <  1:3070827/71‑1 (MQ=255)
          gTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCt                                                    >  1:1020480/1‑70 (MQ=255)
           tCCACACTGGCTGGATGTGAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCa                                                              <  1:2504466/58‑1 (MQ=255)
              acacTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaa                                                             <  1:2511527/57‑1 (MQ=255)
              acacTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAacac                                                        >  1:2077924/1‑62 (MQ=255)
               cacTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                                           >  1:1772818/1‑58 (MQ=255)
                    gctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgc                                         <  1:4403188/71‑1 (MQ=255)
                         aTGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCa                                                 >  1:932031/1‑58 (MQ=255)
                         aTGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCa                                                 >  1:146131/1‑58 (MQ=255)
                               aaCCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATgcgc                              <  1:3685361/71‑1 (MQ=255)
                                      ttcttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAg                        >  1:1437526/1‑70 (MQ=255)
                                       tcttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAg                        >  1:527570/1‑69 (MQ=255)
                                        cttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGga                      <  1:2011279/70‑1 (MQ=255)
                                        cttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGga                      <  1:2404828/70‑1 (MQ=255)
                                         ttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAg                        >  1:2738101/1‑67 (MQ=255)
                                         ttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAg                        >  1:3202323/1‑67 (MQ=255)
                                          tcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTc                                                  >  1:2821063/1‑40 (MQ=255)
                                           cACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATg                                 >  1:1996678/1‑56 (MQ=255)
                                           cACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATg                                 >  1:3508330/1‑56 (MQ=255)
                                           cACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGgaga                    >  1:2869253/1‑69 (MQ=255)
                                           cACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGgag                     >  1:3816360/1‑68 (MQ=255)
                                            aCTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGGGATGa                 <  1:3511806/71‑1 (MQ=255)
                                             cTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGAc                >  1:1432348/1‑71 (MQ=255)
                                              ttATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGgag                     >  1:3105918/1‑65 (MQ=255)
                                              ttATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGgag                     >  1:1739274/1‑65 (MQ=255)
                                              ttATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGga                      <  1:4255018/64‑1 (MQ=255)
                                              ttATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGga                      <  1:2956590/64‑1 (MQ=255)
                                                aTCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCg              >  1:1258035/1‑70 (MQ=255)
                                                aTCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCg              >  1:116137/1‑70 (MQ=255)
                                                aTCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATCACCg              >  1:3962854/1‑70 (MQ=255)
                                                 tCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCgg             <  1:165868/70‑1 (MQ=255)
                                                   aGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgccgc                                      <  1:2658802/43‑1 (MQ=255)
                                                   aGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgccgc                                      <  1:3818894/43‑1 (MQ=255)
                                                   aGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGAcc          <  1:1758310/71‑1 (MQ=255)
                                                   aGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGAcc          >  1:2129417/1‑71 (MQ=255)
                                                      gAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAAt                                  >  1:2369819/1‑44 (MQ=255)
                                                      gAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAAt                                  >  1:2364449/1‑44 (MQ=255)
                                                      gAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAAt                                  >  1:3364247/1‑44 (MQ=255)
                                                      gAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGGGCTGGCAGGAGATGACCGGACCCAt       >  1:2993446/1‑71 (MQ=255)
                                                       aaCTAAACGCGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCAt       <  1:2524165/70‑1 (MQ=255)
                                                       aaCTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGAt                   <  1:2272784/58‑1 (MQ=255)
                                                       aaCTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGAt                   <  1:1249499/58‑1 (MQ=255)
                                                       aaCTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCAt       <  1:2529383/70‑1 (MQ=255)
                                                         cTAAACGAGGGCAAAACACTGCTTCACGCCGGgccgcc                                     <  1:499434/38‑1 (MQ=255)
                                                         cTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGa     >  1:2508444/1‑70 (MQ=255)
                                                          tAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGga                      <  1:1550048/52‑1 (MQ=255)
                                                          tAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGga                      <  1:2621064/52‑1 (MQ=255)
                                                           aaaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAg  >  1:3914298/1‑71 (MQ=255)
                                                            aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCaa                                   >  1:1303054/1‑37 (MQ=255)
                                                            aaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCaa                                   >  1:661016/1‑37 (MQ=255)
                                                                 |                                                                
CGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAG  >  W3110S.gb/547629‑547758

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