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.G92.5% 128.5 / 2.7 40pseudogene (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/2);  new base G (18/19);  total (19/21)
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: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

GCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACT  >  W3110S.gb/547761‑547894
                                                                       |                                                                  
gcgtgcgtGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCa                                                                     <  1:2513400/71‑1 (MQ=255)
gcgtgcgtGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGAAATACTGGAGCa                                                                     <  1:2590288/71‑1 (MQ=255)
  gtgcgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgg                                                                 >  1:4569210/1‑71 (MQ=255)
          gcgcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTg                                                        <  1:4363194/71‑1 (MQ=255)
           cgcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGa                                                       <  1:4290344/71‑1 (MQ=255)
           cgcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGTAAGTGa                                                       <  1:4597/71‑1 (MQ=255)
                 tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCa                                                 >  1:2265851/1‑71 (MQ=255)
                    gTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAg                                                                   <  1:4251300/52‑1 (MQ=255)
                     ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTc                                              <  1:1776418/70‑1 (MQ=255)
                     ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgg                                                                 <  1:1685223/52‑1 (MQ=255)
                        gAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTg                                          >  1:2014374/1‑71 (MQ=255)
                         aaGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGt                                         >  1:2027659/1‑71 (MQ=255)
                          aGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTc                                        <  1:238300/71‑1 (MQ=255)
                            gTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCa                                                                     <  1:278030/43‑1 (MQ=255)
                            gTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCa                                                                     <  1:599407/43‑1 (MQ=255)
                            gTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAc                                      <  1:1734562/71‑1 (MQ=255)
                            gTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAc                                      <  1:744387/71‑1 (MQ=255)
                               gggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCa                                    <  1:69128/70‑1 (MQ=255)
                               gggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCa                                    <  1:2825504/70‑1 (MQ=255)
                               gggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAt                                   >  1:602803/1‑71 (MQ=255)
                                  cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                                <  1:4069224/71‑1 (MQ=255)
                                  cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                                <  1:3441416/71‑1 (MQ=255)
                                    aaaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCtt                                           <  1:4622680/58‑1 (MQ=255)
                                    aaaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGCGGAAGTGAACTTCATTCCTTGTCACCATGTGaa                              >  1:3754836/1‑71 (MQ=255)
                                      aGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAAt                             >  1:348707/1‑70 (MQ=255)
                                      aGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTCAATg                            >  1:3109136/1‑71 (MQ=255)
                                       gATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTc                                              >  1:2541057/1‑52 (MQ=39)
                                        aTGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTc                                        <  1:2396940/57‑1 (MQ=255)
                                        aTGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc                          <  1:433775/71‑1 (MQ=255)
                                         tGAAGCGCAGGCGCTGGCTATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc                          >  1:2707362/1‑70 (MQ=255)
                                          gAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGaa                              >  1:4870622/1‑65 (MQ=255)
                                           aaGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTc                       >  1:4216309/1‑71 (MQ=255)
                                            aGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCg                      <  1:2056170/71‑1 (MQ=255)
                                                     gcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCa                                       >  1:5094776/1‑45 (MQ=38)
                                                         gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg         >  1:2640908/1‑71 (MQ=255)
                                                            aTACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt      <  1:2753406/71‑1 (MQ=255)
                                                             tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCa                 <  1:294237/59‑1 (MQ=255)
                                                             tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg         <  1:5140333/67‑1 (MQ=255)
                                                             tACTGGAGCAGTGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg         <  1:4294683/67‑1 (MQ=255)
                                                               cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc                          >  1:2810891/1‑48 (MQ=255)
                                                               cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGcc                          >  1:175565/1‑48 (MQ=255)
                                                               cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAc   <  1:4678957/71‑1 (MQ=255)
                                                                tGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg         <  1:3905491/64‑1 (MQ=255)
                                                                tGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACt  >  1:963000/1‑71 (MQ=255)
                                                                 ggAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCaa                >  1:3886961/1‑56 (MQ=255)
                                                                  gAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCaa                >  1:5106465/1‑55 (MQ=255)
                                                                  gAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCaa                >  1:1763327/1‑55 (MQ=255)
                                                                       |                                                                  
GCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACT  >  W3110S.gb/547761‑547894

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: