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.G97.6% 137.6 / ‑3.8 42pseudogene (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/1);  new base G (20/21);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-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.

CGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑GG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAT  >  W3110S.gb/547766‑547890
                                                                  |                                                            
cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑gg‑ggg                                                        <  1:341211/71‑2 (MQ=255)
       gcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGaa                                                  >  1:2724306/1‑71 (MQ=255)
        cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAAc                                                 <  1:2379408/71‑1 (MQ=255)
        cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAAc                                                 <  1:3532245/71‑1 (MQ=255)
          tGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACtt                                               <  1:2399996/71‑1 (MQ=255)
            tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑GGTGGAAGTGAACTTCa                                             <  1:2074623/71‑1 (MQ=255)
              tGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTc                                              <  1:1789727/68‑1 (MQ=255)
                  cGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCa                                             >  1:173823/1‑65 (MQ=255)
                  cGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCa                                             >  1:3277348/1‑65 (MQ=255)
                  cGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCtt                                       >  1:3110318/1‑71 (MQ=255)
                  cGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCt                                        >  1:3334044/1‑70 (MQ=255)
                   gAAGGTTGGACGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCtt                                       <  1:1345562/70‑1 (MQ=255)
                    aaGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTg                                      <  1:1617309/70‑1 (MQ=255)
                      ggTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCa                                   >  1:2234629/1‑71 (MQ=255)
                            gCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCAtgt                             <  1:1718321/71‑1 (MQ=255)
                              gAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGa                           <  1:1679413/71‑1 (MQ=255)
                               aaaGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑gg                                                            >  1:429073/1‑37 (MQ=255)
                                aaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAAt                         >  1:1735756/1‑71 (MQ=255)
                                 aGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATg                        >  1:2885573/1‑71 (MQ=255)
                                 aGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATg                        >  1:1175608/1‑71 (MQ=255)
                                   aTGAAGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc                      >  1:2669753/1‑71 (MQ=255)
                                     gAAGCGCAGGCGGTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCg                     >  1:3324299/1‑70 (MQ=255)
                                      aaGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTc                   >  1:1233226/1‑71 (MQ=255)
                                      aaGCGCAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTc                   >  1:1256434/1‑71 (MQ=255)
                                        gcgcAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCa                                             >  1:226216/1‑43 (MQ=38)
                                         cgcAGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCggg                >  1:2695646/1‑71 (MQ=255)
                                            aGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCa             <  1:2125620/71‑1 (MQ=255)
                                            aGGCGCTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCa             <  1:587829/71‑1 (MQ=255)
                                              gcgcTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAAt           <  1:219691/71‑1 (MQ=255)
                                               cgcTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATg          <  1:2926476/71‑1 (MQ=255)
                                                 cTGGCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCAtgt                             <  1:2962168/50‑1 (MQ=255)
                                                    gCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCg      <  1:3175947/70‑1 (MQ=255)
                                                    gCAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCg      <  1:3526457/70‑1 (MQ=255)
                                                     cAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGt    <  1:234060/71‑1 (MQ=255)
                                                     cAATACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGt    <  1:1240173/71‑1 (MQ=255)
                                                      aaTACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATgg         <  1:2259867/65‑1 (MQ=255)
                                                       aTACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt  >  1:575051/1‑71 (MQ=255)
                                                       aTACTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt  <  1:2848258/71‑1 (MQ=255)
                                                          cTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCa                                >  1:1953765/1‑38 (MQ=38)
                                                          cTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATg                        <  1:1418880/46‑1 (MQ=255)
                                                          cTGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATg                        <  1:3587424/46‑1 (MQ=255)
                                                           tGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCggg                >  1:2910708/1‑53 (MQ=255)
                                                           tGGAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCggg                >  1:829417/1‑53 (MQ=255)
                                                             gAGCAGGG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGc       >  1:1265204/1‑60 (MQ=255)
                                                                  |                                                            
CGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑GG‑GGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAT  >  W3110S.gb/547766‑547890

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: ‑