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.

GCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑G‑G‑GG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACT  >  W3110S.gb/547761‑547894
                                                                       |                                                                  
gcgtgcgtGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCa                                                                     <  1:2513400/71‑1 (MQ=255)
gcgtgcgtGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGAAATACTGGAGCa                                                                     <  1:2590288/71‑1 (MQ=255)
  gtgcgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑g‑g                                                                 >  1:4569210/1‑71 (MQ=255)
          gcgcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGG‑G‑GG‑AAGTg                                                        <  1:4363194/71‑1 (MQ=255)
           cgcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGG‑G‑GG‑AAGTGa                                                       <  1:4290344/71‑1 (MQ=255)
           cgcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑G‑G‑GGTAAGTGa                                                       <  1:4597/71‑1 (MQ=255)
                 tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGG‑G‑GG‑AAGTGAACTTCa                                                 >  1:2265851/1‑71 (MQ=255)
                    gTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑g                                                                   <  1:4251300/52‑1 (MQ=255)
                     ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTc                                              <  1:1776418/70‑1 (MQ=255)
                     ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑g‑g                                                                 <  1:1685223/52‑1 (MQ=255)
                        gAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTg                                          >  1:2014374/1‑71 (MQ=255)
                         aaGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGt                                         >  1:2027659/1‑71 (MQ=255)
                          aGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTc                                        <  1:238300/71‑1 (MQ=255)
                            gTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCa                                                                     <  1:278030/43‑1 (MQ=255)
                            gTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCa                                                                     <  1:599407/43‑1 (MQ=255)
                            gTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTCAc                                      <  1:1734562/71‑1 (MQ=255)
                            gTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTCAc                                      <  1:744387/71‑1 (MQ=255)
                               gggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTCACCa                                    <  1:69128/70‑1 (MQ=255)
                               gggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTCACCa                                    <  1:2825504/70‑1 (MQ=255)
                               gggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTCACCAt                                   >  1:602803/1‑71 (MQ=255)
                                  cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTCACCAtgtg                                <  1:4069224/71‑1 (MQ=255)
                                  cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTCACCAtgtg                                <  1:3441416/71‑1 (MQ=255)
                                    aaaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCtt                                           <  1:4622680/58‑1 (MQ=255)
                                    aaaGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑G‑GCGG‑AAGTGAACTTCATTCCTTGTCACCATGTGaa                              >  1:3754836/1‑71 (MQ=255)
                                      aGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTCACCATGTGAAt                             >  1:348707/1‑70 (MQ=255)
                                      aGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTCACCATGTCAATg                            >  1:3109136/1‑71 (MQ=255)
                                       gATGAAGCGCAGGCGCTGGCAATACTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTc                                              >  1:2541057/1‑52 (MQ=39)
                                        aTGAAGCGCAGGCGCTGGCAATACTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTc                                        <  1:2396940/57‑1 (MQ=255)
                                        aTGAAGCGCAGGCGCTGGCAATACTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGcc                          <  1:433775/71‑1 (MQ=255)
                                         tGAAGCGCAGGCGCTGGCTATACTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGcc                          >  1:2707362/1‑70 (MQ=255)
                                          gAAGCGCAGGCGCTGGCAATACTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTCACCATGTGaa                              >  1:4870622/1‑65 (MQ=255)
                                           aaGCGCAGGCGCTGGCAATACTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTc                       >  1:4216309/1‑71 (MQ=255)
                                            aGCGCAGGCGCTGGCAATACTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCg                      <  1:2056170/71‑1 (MQ=255)
                                                     gcTGGCAATACTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTCa                                       >  1:5094776/1‑45 (MQ=38)
                                                         gCAATACTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg         >  1:2640908/1‑71 (MQ=255)
                                                            aTACTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt      <  1:2753406/71‑1 (MQ=255)
                                                             tACTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCa                 <  1:294237/59‑1 (MQ=255)
                                                             tACTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg         <  1:5140333/67‑1 (MQ=255)
                                                             tACTGGAGCA‑GTG‑GG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg         <  1:4294683/67‑1 (MQ=255)
                                                               cTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGcc                          >  1:2810891/1‑48 (MQ=255)
                                                               cTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGcc                          >  1:175565/1‑48 (MQ=255)
                                                               cTGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAc   <  1:4678957/71‑1 (MQ=255)
                                                                tGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg         <  1:3905491/64‑1 (MQ=255)
                                                                tGGAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACt  >  1:963000/1‑71 (MQ=255)
                                                                 ggAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCaa                >  1:3886961/1‑56 (MQ=255)
                                                                  gAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCaa                >  1:5106465/1‑55 (MQ=255)
                                                                  gAGCAGG‑G‑GG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCaa                >  1:1763327/1‑55 (MQ=255)
                                                                       |                                                                  
GCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑G‑G‑GG‑AAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACT  >  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: ‑