Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,365,218 A→G 63.4% pseudogene (179/480 nt)
pseudogene (179/1745 nt)
yhcE →
yhcE →
ECK3207:JW3184:b3217; hypothetical protein, N‑ter fragment
ECK3207:JW3184+JW3187:b4569; hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,365,2180AG63.4% 39.4 / 38.4 41pseudogene (179/480 nt)
pseudogene (179/1745 nt)
yhcE
yhcE
ECK3207:JW3184:b3217; hypothetical protein, N‑ter fragment
ECK3207:JW3184+JW3187:b4569; hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (4/11);  new base G (13/13);  total (17/24)
Fisher's exact test for biased strand distribution p-value = 1.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAACAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGA  >  W3110S.gb/3365152‑3365285
                                                                  |                                                                   
ccATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAACAATCACGCTGGTaatgcgatgc                                                                 <  1:3454508/71‑1 (MQ=255)
ccATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAACAATCACGCTGGTaatgcgatgc                                                                 <  1:1135548/71‑1 (MQ=255)
ccATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAACAATCACGCTGGTaatgcgatgc                                                                 <  1:1481663/71‑1 (MQ=255)
      gcCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAACAATCACGCTGGTAATGCGATGCGATCAg                                                           <  1:42745/71‑1 (MQ=255)
      gcCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAACAATCACGCTGGTAATGCGATGCGATCAg                                                           <  1:1624172/71‑1 (MQ=255)
      gcCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAACAATCACGCTGGTAATGCGATGCGATCAg                                                           <  1:1326758/71‑1 (MQ=255)
      gcCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAACAATCACGCTGGTAATGCAATGCGATCAg                                                           <  1:3440186/71‑1 (MQ=255)
                 cGGGGTAAAACAACACAATTAGCAGACAGAACAATCACGCTGGTaatgcgatgc                                                                 <  1:3620511/54‑1 (MQ=255)
                        aaaCAACACAATTAGCAGACAGAACAATCACGCTGGTAATGCAATGCGATCAGGACGCTCata                                                 >  1:1062502/1‑63 (MQ=255)
                        aaaCAACACAATTAGCAGACAGAACAATCACGCTGGTAATGCAATGCGATAGGACGCTCata                                                 >  1:3037198/1‑62 (MQ=255)
                         aacaacACAATTAGCAGACAGAACAATCACGCTGGTAATGCGATGCGATCAGGACGCTCATATTCGTGTTc                                        <  1:1431467/71‑1 (MQ=255)
                         aacaacACAATTAGCAGACAGAACAATCACGCTGGTAATGCGATGCGATCAGGACGCTCATATTCGTGTTc                                        <  1:3437907/71‑1 (MQ=255)
                                    tAGCAGACAGAACAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATAcc                             >  1:1249614/1‑71 (MQ=255)
                                    tAGCAGACAGAACAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATAcc                             >  1:2717870/1‑71 (MQ=255)
                                     aGCAGACAGAACAATCACGCTGGTAATGCGATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCg                            >  1:1261226/1‑71 (MQ=255)
                                     aGCAGACAGAACAATCACGCTGGTAATGCGATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCg                            >  1:1652123/1‑71 (MQ=255)
                                      gCAGACAGAACAATCACGCTGGTAATGCGATGCGATCAg                                                           <  1:1831199/39‑1 (MQ=255)
                                       cagacagaACAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCt                          <  1:3404452/71‑1 (MQ=255)
                                       cagacagaACAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCt                          <  1:2063457/71‑1 (MQ=255)
                                       cagacagaACAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCt                          <  1:3008040/71‑1 (MQ=255)
                                       cagacagaACAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCt                          <  1:1115081/71‑1 (MQ=255)
                                           cagaACAATCACGCTGGTAATGCGATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAata                      >  1:3103056/1‑71 (MQ=255)
                                           cagaACAATCACGCTGGTAATGCGATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAata                      >  1:2253380/1‑71 (MQ=255)
                                            agaACAATCACGCTGGTAATGCGATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAatat                     <  1:2908974/71‑1 (MQ=255)
                                            agaACAATCACGCTGGTAATGCGATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAatat                     <  1:3563694/71‑1 (MQ=255)
                                            agaACAATCACGCTGGTAATGCGATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAatat                     <  1:1301333/71‑1 (MQ=255)
                                            agaACAATCACGCTGGTAATGCGATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAata                      >  1:3574158/1‑70 (MQ=255)
                                            agaACAATCACGCTGGTAATGCGATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAata                      >  1:3283124/1‑70 (MQ=255)
                                            agaACAATCACGCTGGTAATGCGATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAat                       >  1:3089118/1‑69 (MQ=255)
                                              aaCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTg                                   <  1:2306734/55‑1 (MQ=255)
                                              aaCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTg                                   <  1:2841082/55‑1 (MQ=255)
                                              aaCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTg                                   <  1:914185/55‑1 (MQ=255)
                                              aaCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTg                                   <  1:1718991/55‑1 (MQ=255)
                                              aaCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTg                                   <  1:1737206/55‑1 (MQ=255)
                                              aaCAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTg                                   <  1:622800/55‑1 (MQ=255)
                                                 aaTCACGCTGGTAATGCGATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTcc                 >  1:84294/1‑70 (MQ=255)
                                                 aaTCACGCTGGTAATGCGATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTcc                 >  1:3575454/1‑70 (MQ=255)
                                                 aaTCACGCTGGTAATGCGATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTcc                 >  1:2950855/1‑70 (MQ=255)
                                                     aCGCTGGTAATGCGATGCGATCAGGACGCTCATATTCGTGTTCAg                                      >  1:727751/1‑45 (MQ=255)
                                                     aCGCTGGTAATGCGATGCGATCAGGACGCTCATATTCGTGTTCAg                                      >  1:2128187/1‑45 (MQ=255)
                                                     aCGCTGGTAATGCGATGCGATCAGGACGCTCATATTCGTGTTCAg                                      >  1:2176970/1‑45 (MQ=255)
                                                       gCTGGTAATGCGATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACg          <  1:2336797/71‑1 (MQ=255)
                                                       gCTGGTAATGCGATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACg          <  1:2738942/71‑1 (MQ=255)
                                                       gCTGGTAATGCGATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACg          <  1:1936745/71‑1 (MQ=255)
                                                       gCTGGTAATGCGATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACg          <  1:1549409/71‑1 (MQ=255)
                                                              atgcgatgcgATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTaa                        <  1:1343862/50‑1 (MQ=255)
                                                              atgcgatgcgATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTaa                        <  1:1063587/50‑1 (MQ=255)
                                                              atgcgatgcgATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACggatttg    >  1:2610822/1‑70 (MQ=255)
                                                              atgcgatgcgATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACggatttg    >  1:2173988/1‑70 (MQ=255)
                                                               tgcgatgcgaTCAGGACGCTCATATTCGTGTTCAGTTGaa                                 <  1:3230108/40‑1 (MQ=255)
                                                               tgcgatgcgaTCAGGACGCTCATATTCGTGTTCAGTTGaa                                 <  1:2477007/40‑1 (MQ=255)
                                                               tgcgatgcgaTCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCca                <  1:3623678/57‑1 (MQ=255)
                                                               tgcgatgcgaTCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCca                <  1:2556605/57‑1 (MQ=255)
                                                               tgcgatgcgaTCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACggatttgga  <  1:929100/71‑1 (MQ=255)
                                                                  |                                                                   
CCATGCGCCAAGTTGACCGGGGTAAAACAACACAATTAGCAGACAGAACAATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGA  >  W3110S.gb/3365152‑3365285

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: