Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 548,801 G→A 54.5% pseudogene (1221/1259 nt)
pseudogene (964/1002 nt)
ylbE →
ylbE →
ECK0512:JW0508+JW0507:b4572; hypothetical protein
ECK0512:JW0508:b0519; hypothetical protein, C‑ter fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb548,8010GA54.5% ‑2.5 / 27.5 22pseudogene (1221/1259 nt)
pseudogene (964/1002 nt)
ylbE
ylbE
ECK0512:JW0508+JW0507:b4572; hypothetical protein
ECK0512:JW0508:b0519; hypothetical protein, C‑ter fragment
Reads supporting (aligned to +/- strand):  ref base G (6/4);  new base A (7/5);  total (13/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

GGGGATCGGGCAGATTGGCGCAGGCACCGTGCGGGCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCT  >  W3110S.gb/548743‑548867
                                                          |                                                                  
ggggATCGGGCAGATTGGCGCAGGCACCGTGCGGGCACCGCTGGCGTGCTTTGAACAGACGCTGGAAGCa                                                         >  1:1369960/1‑70 (MQ=255)
  ggATCGGGCAGATTGGCGCAGGCACCGTGCGGGCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTg                                                      >  1:1650049/1‑71 (MQ=255)
              ttGGCGCAGGCACCGTGCGGGCACCGCTGGCGTGCTTTGAACAGACGCTGGAAGCACTGGCTGAAAGCATg                                          >  1:569715/1‑71 (MQ=255)
                            gTGCGGGCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAg                             >  1:781645/1‑70 (MQ=255)
                              gCGGGCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGaa                          >  1:1293453/1‑71 (MQ=255)
                               cGGGCACCGCTGGCGTGCTTTGAACAGACGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAAc                         <  1:1486892/71‑1 (MQ=255)
                               cGGGCACCGCTGGCGTGCTTTGAACAGACGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAAc                         <  1:1654485/71‑1 (MQ=255)
                                   cACCGCTGGCGTGCTTTGAACAGACGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAAcgcg                      >  1:1308450/1‑70 (MQ=255)
                                    aCCGCTGGCGTGCTTTGAACAGACGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACgcgc                     >  1:1720488/1‑70 (MQ=255)
                                    aCCGCTGGCGTGCTTTGAACAGACGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACgcgc                     >  1:214583/1‑70 (MQ=255)
                                       gCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAAt                  >  1:975598/1‑70 (MQ=255)
                                        cTGGCGTGCTTTGAACAGACGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGa                <  1:662504/71‑1 (MQ=255)
                                             gTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGatc           <  1:1309246/71‑1 (MQ=255)
                                             gTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGatc           <  1:834546/71‑1 (MQ=255)
                                             gTGCTTTGAACAGACGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGat            <  1:1091761/70‑1 (MQ=255)
                                              tGCTTTGAACAGACGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGa                           <  1:584621/54‑1 (MQ=255)
                                               gCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGa                >  1:581505/1‑64 (MQ=255)
                                                cTTTGAACAGACGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGatcatc        >  1:336712/1‑71 (MQ=255)
                                                cTTTGAACAGACGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGatcatc        >  1:499625/1‑71 (MQ=255)
                                                   tGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATcatcc       <  1:247058/69‑1 (MQ=255)
                                                   tGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATcatcc       <  1:1050933/69‑1 (MQ=255)
                                                        aGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGAACATGCATCct  >  1:199350/1‑69 (MQ=255)
                                                          |                                                                  
GGGGATCGGGCAGATTGGCGCAGGCACCGTGCGGGCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCT  >  W3110S.gb/548743‑548867

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: