Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,335,995 C→T 25.0% pseudogene (3078/4605 nt)
pseudogene (3078/4104 nt)
yfaS ←
yfaS ←
ECK2220:JW2222+JW2221:b4500; hypothetical protein
ECK2220:JW2222:b2228; hypothetical protein, N‑ter fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,335,9950CT25.0% 88.2 / 28.2 48pseudogene (3078/4605 nt)
pseudogene (3078/4104 nt)
yfaS
yfaS
ECK2220:JW2222+JW2221:b4500; hypothetical protein
ECK2220:JW2222:b2228; hypothetical protein, N‑ter fragment
Reads supporting (aligned to +/- strand):  ref base C (20/16);  new base T (5/7);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-01

TCACTACTTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTG  >  W3110S.gb/2335927‑2336061
                                                                    |                                                                  
tCACTACTTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGTAt                                                                  >  1:1598101/1‑71 (MQ=255)
tCACTACTTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCAt                                                                  >  1:3752155/1‑71 (MQ=255)
tCACTACTTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCAt                                                                  >  1:229410/1‑71 (MQ=255)
  actactTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGTATGa                                                                >  1:1430687/1‑71 (MQ=255)
      ctTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACAtt                                                            >  1:972634/1‑71 (MQ=255)
               cGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTctgctg                                                   <  1:1295057/71‑1 (MQ=255)
                  atatTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTc                                                <  1:2901054/71‑1 (MQ=255)
                  atatTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTc                                                <  1:593827/71‑1 (MQ=255)
                   tatTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGTATGACAttt                                                           <  1:3833312/59‑1 (MQ=255)
                   tatTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTc                                                >  1:2063833/1‑70 (MQ=255)
                   tatTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTc                                                >  1:2430386/1‑70 (MQ=255)
                     ttGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGTATGa                                                                <  1:794466/52‑1 (MQ=255)
                     ttGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGTATGa                                                                <  1:2542273/52‑1 (MQ=255)
                     ttGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGTATGACATTTTTctg                                                      <  1:3690121/62‑1 (MQ=255)
                      tGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACg                                            >  1:791981/1‑71 (MQ=255)
                      tGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACg                                            >  1:726514/1‑71 (MQ=255)
                       gCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGTATGACATTTTTCTGCTGTTCAACgg                                           >  1:1116048/1‑71 (MQ=255)
                       gCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACg                                            <  1:23175/70‑1 (MQ=255)
                            ccTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTc                                                <  1:2130548/61‑1 (MQ=255)
                             cTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGc                                     <  1:2983489/71‑1 (MQ=255)
                                  cGGGCAACATCAGCGCGTTATCGCCACCACCGAGTATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTa                                >  1:2245693/1‑71 (MQ=255)
                                  cGGGCAACATCAGCGCATTATCGCCACCACCGAGTATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTa                                >  1:3483196/1‑71 (MQ=255)
                                         cATCAGCGCGTTATCGCCACCACCGAGTATGACAttt                                                           <  1:3229978/37‑1 (MQ=255)
                                          aTCAGCGCGTTATCGCCACCACCGAGCATGACAttttt                                                         >  1:2344433/1‑38 (MQ=255)
                                          aTCAGCGCGTTATCGCCACCACCGAGCATGACAttttt                                                         >  1:1169094/1‑38 (MQ=255)
                                          aTCAGCGCGTTATCGCCACCACCGAGCATGACAttttt                                                         >  1:3147187/1‑38 (MQ=255)
                                          aTCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACaa                         >  1:3761232/1‑70 (MQ=255)
                                           tCAGCGCGTTATCGCCACCACCGAGCATGACAttttt                                                         >  1:2876140/1‑37 (MQ=255)
                                           tCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACaaaa                       <  1:3468882/71‑1 (MQ=255)
                                                   ttATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCa                                               <  1:3790858/39‑1 (MQ=255)
                                                     ttCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAg                  <  1:1935931/65‑1 (MQ=255)
                                                     aTCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAAc                                             >  1:3149744/1‑39 (MQ=255)
                                                     aTCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGc                                     >  1:1052053/1‑47 (MQ=255)
                                                     aTCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTtg             >  1:3491358/1‑71 (MQ=255)
                                                      tCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGt                                  <  1:2451716/49‑1 (MQ=255)
                                                       cGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGt                                  <  1:2287174/48‑1 (MQ=255)
                                                         ccaccaCCGAGTATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCt         <  1:164979/71‑1 (MQ=255)
                                                         ccaccaCCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCt         <  1:3068722/71‑1 (MQ=255)
                                                         ccaccaCCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCt         <  1:3905110/71‑1 (MQ=255)
                                                             caccGAGTATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATg     <  1:613010/71‑1 (MQ=255)
                                                             caccGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTa                                >  1:1321526/1‑44 (MQ=255)
                                                             caccGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTa                                >  1:1625893/1‑44 (MQ=255)
                                                             caccGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTa                                >  1:1069084/1‑44 (MQ=255)
                                                             caccGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATg     <  1:3986540/71‑1 (MQ=255)
                                                             caccGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATg     <  1:280949/71‑1 (MQ=255)
                                                                cGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTg  <  1:3269685/71‑1 (MQ=255)
                                                                  aGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTtg             >  1:642473/1‑58 (MQ=255)
                                                                  aGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTtg             >  1:2087706/1‑58 (MQ=255)
                                                                    |                                                                  
TCACTACTTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTG  >  W3110S.gb/2335927‑2336061

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: