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.G94.5% 111.8 / ‑0.7 36pseudogene (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 . (1/1);  new base G (17/17);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
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.

TGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTC  >  W3110S.gb/547764‑547896
                                                                    |                                                                   
tgcgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggg                                                                 <  1:1939924/71‑1 (MQ=255)
  cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgggg                                                               >  1:3028596/1‑70 (MQ=255)
        cgcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGa                                                         <  1:440141/71‑1 (MQ=255)
        cgcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGa                                                         <  1:1283608/71‑1 (MQ=255)
              tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCa                                                   >  1:2948435/1‑71 (MQ=255)
                  ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGCGAAGTGAACTTCATTc                                                <  1:1625386/70‑1 (MQ=255)
                    cGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCtt                                             >  1:2511483/1‑71 (MQ=255)
                      aaGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGt                                           <  1:918651/71‑1 (MQ=255)
                            gggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAt                                     >  1:348079/1‑71 (MQ=255)
                            gggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAt                                     >  1:2513037/1‑71 (MQ=255)
                             ggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAc                                        >  1:3579193/1‑67 (MQ=255)
                               cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                                  <  1:1247154/71‑1 (MQ=255)
                               cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                                  <  1:2541513/71‑1 (MQ=255)
                                gAAAGATGAAGCGCATGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgt                                   <  1:2842004/69‑1 (MQ=255)
                                        aaGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGt                                           >  1:292177/1‑53 (MQ=39)
                                        aaGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTc                         <  1:3007304/71‑1 (MQ=255)
                                                gcgcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                                  <  1:11380/54‑1 (MQ=39)
                                                gcgcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                                  <  1:2548129/54‑1 (MQ=39)
                                                 cgcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                                  <  1:420687/53‑1 (MQ=39)
                                                    tGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                                  <  1:1209504/50‑1 (MQ=255)
                                                    tGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                                  <  1:3012554/50‑1 (MQ=255)
                                                    tGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc                            <  1:2171679/56‑1 (MQ=255)
                                                     ggCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtg                                    >  1:1116702/1‑47 (MQ=255)
                                                     ggCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtg                                    >  1:2425221/1‑47 (MQ=255)
                                                      gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGaa                                >  1:1882796/1‑50 (MQ=255)
                                                      gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg           >  1:2191164/1‑71 (MQ=255)
                                                         aTACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCaa                  <  1:49102/61‑1 (MQ=255)
                                                          tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAc                                        <  1:51615/38‑1 (MQ=38)
                                                          tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                                  <  1:449137/44‑1 (MQ=38)
                                                          tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                                  <  1:814000/44‑1 (MQ=38)
                                                          tACTGGAGCAGGTGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt        >  1:2876487/1‑70 (MQ=255)
                                                            cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCa                                         >  1:3195245/1‑35 (MQ=37)
                                                            cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAc     >  1:1782034/1‑71 (MQ=255)
                                                            cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAc     >  1:3682279/1‑71 (MQ=255)
                                                            cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAc     >  1:2490106/1‑71 (MQ=255)
                                                              ggAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACt    >  1:459259/1‑70 (MQ=255)
                                                               gAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTc  >  1:2946143/1‑71 (MQ=255)
                                                               gAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTc  >  1:2548185/1‑71 (MQ=255)
                                                                    |                                                                   
TGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTC  >  W3110S.gb/547764‑547896

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: