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.G98.6% 234.0 / ‑3.9 69pseudogene (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 (40/28);  total (40/29)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-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.

CGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTC  >  W3110S.gb/547762‑547896
                                                                      |                                                                  
cgtgcgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAg                                                                   <  1:3404683/71‑1 (MQ=255)
    cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggggg                                                              <  1:7518171/71‑2 (MQ=255)
    cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgggg                                                               >  1:4404801/1‑70 (MQ=255)
            cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAAc                                                       >  1:3184785/1‑71 (MQ=255)
               gTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                                                     >  1:6898249/1‑70 (MQ=255)
                    ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCa                                                   >  1:583065/1‑67 (MQ=255)
                      cGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCa                                                                     <  1:769673/48‑1 (MQ=255)
                      cGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCtt                                             >  1:4270000/1‑71 (MQ=255)
                       gAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTg                                            >  1:2876220/1‑71 (MQ=255)
                        aaGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGt                                           >  1:504549/1‑71 (MQ=255)
                         aGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgg                                                                  >  1:7238333/1‑47 (MQ=255)
                          ggTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCa                                         <  1:39252/71‑1 (MQ=255)
                          ggTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCa                                         >  1:4605997/1‑71 (MQ=255)
                          ggTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCa                                         <  1:1724504/71‑1 (MQ=255)
                             tGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCGCCa                                      >  1:7933981/1‑71 (MQ=255)
                             tGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCa                                      >  1:3537693/1‑71 (MQ=255)
                                gCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGa                                                         <  1:582097/49‑1 (MQ=38)
                                gCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGa                                                         <  1:3617416/49‑1 (MQ=38)
                                 cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                                  <  1:797001/71‑1 (MQ=255)
                                 cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                                  <  1:3361019/71‑1 (MQ=255)
                                 cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgt                                   >  1:5921742/1‑70 (MQ=255)
                                 cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGAGAAGTGAACTTCATTCCTTGTCACCAtgtg                                  <  1:6063386/71‑1 (MQ=255)
                                  gAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGaa                                                        <  1:4142000/48‑1 (MQ=38)
                                  gAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGaa                                                        <  1:5077274/48‑1 (MQ=38)
                                  gAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGaa                                                        <  1:1236655/48‑1 (MQ=38)
                                   aaaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgt                                   >  1:101206/1‑68 (MQ=255)
                                    aaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAAt                               >  1:5536742/1‑71 (MQ=255)
                                    aaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAAt                               <  1:3703023/71‑1 (MQ=255)
                                     aGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATg                              >  1:6186211/1‑71 (MQ=255)
                                       aTGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGaa                                <  1:7344852/67‑1 (MQ=255)
                                       aTGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAAt                               >  1:1441996/1‑68 (MQ=255)
                                        tGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCg                           >  1:3885399/1‑71 (MQ=255)
                                         gAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATg                              >  1:3855521/1‑67 (MQ=255)
                                           aGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCAtt                                                 <  1:5388332/46‑1 (MQ=38)
                                           aGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCg                        >  1:2064920/1‑71 (MQ=255)
                                            gcgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAAc                                                       <  1:7225166/39‑1 (MQ=37)
                                             cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCtt                                             <  1:1977406/48‑1 (MQ=38)
                                             cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCgg                       <  1:5237695/70‑1 (MQ=255)
                                              gcAGGAGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGc                     >  1:2414919/1‑71 (MQ=255)
                                               cAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTc                         <  1:2875068/66‑1 (MQ=255)
                                               cAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGc                     >  1:2725857/1‑70 (MQ=255)
                                                 ggCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCaa                  >  1:2430764/1‑71 (MQ=255)
                                                 ggCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCaa                  <  1:3619932/71‑1 (MQ=255)
                                                  gcgcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAAt                 >  1:7809436/1‑71 (MQ=255)
                                                   cgcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATg                >  1:5681196/1‑71 (MQ=255)
                                                    gcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATgg               >  1:5128617/1‑71 (MQ=255)
                                                      tGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATggg              <  1:1722262/70‑1 (MQ=255)
                                                       ggCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCg            >  1:6179588/1‑71 (MQ=255)
                                                        gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGaa                                >  1:3377476/1‑50 (MQ=255)
                                                        gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGaa                                >  1:6988622/1‑50 (MQ=255)
                                                        gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGaa                                >  1:4521346/1‑50 (MQ=255)
                                                        gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGaa                                >  1:3208213/1‑50 (MQ=255)
                                                        gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGc                             >  1:7367344/1‑53 (MQ=255)
                                                         cAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg           >  1:3196243/1‑70 (MQ=255)
                                                         cAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGt          <  1:5314996/71‑1 (MQ=255)
                                                          aaTACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTa         >  1:4304963/1‑71 (MQ=255)
                                                            tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAcc                                       <  1:1662089/39‑1 (MQ=38)
                                                            tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAcc                                       <  1:6608481/39‑1 (MQ=38)
                                                            tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTc                         >  1:7667677/1‑53 (MQ=255)
                                                            tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTc                         >  1:7483499/1‑53 (MQ=255)
                                                            tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAtt       <  1:1670846/71‑1 (MQ=255)
                                                              cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAt                                     >  1:5309585/1‑39 (MQ=38)
                                                              cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTa      <  1:6399169/70‑1 (MQ=255)
                                                              cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAc     <  1:3944217/71‑1 (MQ=255)
                                                              cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAc     >  1:5357327/1‑71 (MQ=255)
                                                               tGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGaa                                >  1:5180797/1‑43 (MQ=38)
                                                               tGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATgg               <  1:5336765/60‑1 (MQ=255)
                                                               tGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATgg               <  1:6266936/60‑1 (MQ=255)
                                                               tGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATgg               <  1:2085814/60‑1 (MQ=255)
                                                                 gAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATg                >  1:7551715/1‑57 (MQ=255)
                                                                 gAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACt    <  1:3498588/69‑1 (MQ=255)
                                                                 gAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTc  >  1:423441/1‑71 (MQ=255)
                                                                 gAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTc  >  1:1977199/1‑71 (MQ=255)
                                                                 gAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTc  >  1:3657858/1‑71 (MQ=255)
                                                                      |                                                                  
CGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTC  >  W3110S.gb/547762‑547896

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: