Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 547,835 +G pseudogene (255/261 nt)
pseudogene (255/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.2% 128.2 / ‑3.1 36pseudogene (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/0);  new base G (24/11);  total (25/11)
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.86e-01

GCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGG  >  W3110S.gb/547773‑547887
                                                           |                                                         
gcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggg                                                        <  1:3658952/62‑1 (MQ=255)
gcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggg                                                        <  1:2460653/62‑1 (MQ=255)
gcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggg                                                        <  1:4321430/62‑1 (MQ=255)
gcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggg                                                        <  1:2172386/62‑1 (MQ=255)
gcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggg                                                        <  1:3158120/62‑1 (MQ=255)
gcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggg                                                        <  1:1079776/62‑1 (MQ=255)
gcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggg                                                        <  1:1636962/62‑1 (MQ=255)
gcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggg                                                        <  1:1380509/62‑1 (MQ=255)
gcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggg                                                        <  1:841724/62‑1 (MQ=255)
gcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggg                                                        <  1:3789608/62‑1 (MQ=255)
 cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgggg                                                       <  1:4576502/62‑1 (MQ=255)
 cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgggg                                                       <  1:4535208/62‑1 (MQ=255)
 cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgggg                                                       <  1:226640/62‑1 (MQ=255)
 cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgggg                                                       <  1:2376275/62‑1 (MQ=255)
 cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgggg                                                       <  1:449364/62‑1 (MQ=255)
       tGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTg                                                 >  1:2543428/1‑62 (MQ=255)
            gAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                                            >  1:3931158/1‑62 (MQ=255)
            gAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                                            >  1:727754/1‑62 (MQ=255)
            gAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                                            >  1:4031169/1‑62 (MQ=255)
            gAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACtt                                            >  1:1205238/1‑62 (MQ=255)
            gAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACt                                             >  1:3446511/1‑61 (MQ=255)
            gAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACt                                             >  1:844200/1‑61 (MQ=255)
                        aaaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCa                                <  1:2045817/62‑1 (MQ=255)
                         aaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGaa                                               >  1:2046019/1‑46 (MQ=38)
                         aaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCa                                >  1:1261892/1‑61 (MQ=255)
                         aaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAc                               >  1:4151559/1‑62 (MQ=255)
                         aaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAc                               >  1:2240112/1‑62 (MQ=255)
                         aaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAc                               >  1:1560790/1‑62 (MQ=255)
                         aaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAc                               >  1:3895682/1‑62 (MQ=255)
                         aaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAc                               >  1:1320215/1‑62 (MQ=255)
                              gAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtg                           >  1:3328202/1‑61 (MQ=255)
                              gAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtg                           >  1:3437347/1‑61 (MQ=255)
                              gAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtg                           >  1:1835032/1‑61 (MQ=255)
                              gAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtg                           >  1:2153426/1‑61 (MQ=255)
                              gAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtg                           >  1:2504137/1‑61 (MQ=255)
                              gAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtg                           >  1:2409042/1‑61 (MQ=255)
                                      ggCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCg                  <  1:443585/62‑1 (MQ=255)
                                          cTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCgg              >  1:421106/1‑62 (MQ=255)
                                          cTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCgg              >  1:2412805/1‑62 (MQ=255)
                                               aaTACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCaa         <  1:3574527/62‑1 (MQ=255)
                                               aaTACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCaa         <  1:1342631/62‑1 (MQ=255)
                                                  aCTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCg                  <  1:3364215/50‑1 (MQ=255)
                                                  aCTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATgg      <  1:1080931/62‑1 (MQ=255)
                                                  aCTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATgg      <  1:4023656/62‑1 (MQ=255)
                                                  aCTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATgg      <  1:396140/62‑1 (MQ=255)
                                                   cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATggg     >  1:670669/1‑62 (MQ=255)
                                                   cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATggg     >  1:171950/1‑62 (MQ=255)
                                                     ggAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCg   <  1:3036903/62‑1 (MQ=255)
                                                     ggAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCg   <  1:739427/62‑1 (MQ=255)
                                                     ggAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCg   <  1:3792904/62‑1 (MQ=255)
                                                       aGCAGGGGTAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg  >  1:752618/1‑61 (MQ=255)
                                                           |                                                         
GCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGG  >  W3110S.gb/547773‑547887

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: