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.

CGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAT  >  W3110S.gb/547766‑547890
                                                                  |                                                            
cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggggg                                                        <  1:341211/71‑2 (MQ=255)
       gcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGaa                                                  >  1:2724306/1‑71 (MQ=255)
        cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAAc                                                 <  1:2379408/71‑1 (MQ=255)
        cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAAc                                                 <  1:3532245/71‑1 (MQ=255)
          tGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                                               <  1:2399996/71‑1 (MQ=255)
            tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGTGGAAGTGAACTTCa                                             <  1:2074623/71‑1 (MQ=255)
              tGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTc                                              <  1:1789727/68‑1 (MQ=255)
                  cGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCa                                             >  1:173823/1‑65 (MQ=255)
                  cGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCa                                             >  1:3277348/1‑65 (MQ=255)
                  cGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCtt                                       >  1:3110318/1‑71 (MQ=255)
                  cGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCt                                        >  1:3334044/1‑70 (MQ=255)
                   gAAGGTTGGACGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCtt                                       <  1:1345562/70‑1 (MQ=255)
                    aaGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTg                                      <  1:1617309/70‑1 (MQ=255)
                      ggTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCa                                   >  1:2234629/1‑71 (MQ=255)
                            gCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgt                             <  1:1718321/71‑1 (MQ=255)
                              gAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGa                           <  1:1679413/71‑1 (MQ=255)
                               aaaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgg                                                            >  1:429073/1‑37 (MQ=255)
                                aaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAAt                         >  1:1735756/1‑71 (MQ=255)
                                 aGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATg                        >  1:2885573/1‑71 (MQ=255)
                                 aGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATg                        >  1:1175608/1‑71 (MQ=255)
                                   aTGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc                      >  1:2669753/1‑71 (MQ=255)
                                     gAAGCGCAGGCGGTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCg                     >  1:3324299/1‑70 (MQ=255)
                                      aaGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTc                   >  1:1233226/1‑71 (MQ=255)
                                      aaGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTc                   >  1:1256434/1‑71 (MQ=255)
                                        gcgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCa                                             >  1:226216/1‑43 (MQ=38)
                                         cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCggg                >  1:2695646/1‑71 (MQ=255)
                                            aGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCa             <  1:2125620/71‑1 (MQ=255)
                                            aGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCa             <  1:587829/71‑1 (MQ=255)
                                              gcgcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAAt           <  1:219691/71‑1 (MQ=255)
                                               cgcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATg          <  1:2926476/71‑1 (MQ=255)
                                                 cTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgt                             <  1:2962168/50‑1 (MQ=255)
                                                    gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCg      <  1:3175947/70‑1 (MQ=255)
                                                    gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCg      <  1:3526457/70‑1 (MQ=255)
                                                     cAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGt    <  1:234060/71‑1 (MQ=255)
                                                     cAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGt    <  1:1240173/71‑1 (MQ=255)
                                                      aaTACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATgg         <  1:2259867/65‑1 (MQ=255)
                                                       aTACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt  >  1:575051/1‑71 (MQ=255)
                                                       aTACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt  <  1:2848258/71‑1 (MQ=255)
                                                          cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCa                                >  1:1953765/1‑38 (MQ=38)
                                                          cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATg                        <  1:1418880/46‑1 (MQ=255)
                                                          cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATg                        <  1:3587424/46‑1 (MQ=255)
                                                           tGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCggg                >  1:2910708/1‑53 (MQ=255)
                                                           tGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCggg                >  1:829417/1‑53 (MQ=255)
                                                             gAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGc       >  1:1265204/1‑60 (MQ=255)
                                                                  |                                                            
CGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAT  >  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: