Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,824,280 T→C 100% pseudogene (413/687 nt) rph → ECK3633:JW3618:b3643; ribonuclease PH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,824,2800TC100.0% 178.9 / NA 52pseudogene (413/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/27);  total (25/27)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGG  >  W3110S.gb/3824215‑3824342
                                                                 |                                                              
cGACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATg                                                           >  1:4487672/1‑71 (MQ=255)
       cTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTa                                                      >  1:2374610/1‑69 (MQ=255)
       cTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTa                                                      >  1:794454/1‑69 (MQ=255)
       cTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACa                                                    <  1:2346754/71‑1 (MQ=255)
          cAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAg                                                   >  1:2566406/1‑69 (MQ=255)
          cAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAg                                                   >  1:2001185/1‑69 (MQ=255)
              cTGATGGTGGCACGCGTACCGCTTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtg                                             <  1:522889/71‑1 (MQ=255)
              cTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACa                                                    <  1:2095802/64‑1 (MQ=255)
                 aTGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggtg                                          >  1:4082151/1‑71 (MQ=255)
                  tggtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggtgg                                         <  1:3512986/71‑1 (MQ=255)
                   ggtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggtgg                                         >  1:50335/1‑70 (MQ=255)
                    gtggCACGCGTACCGCGTCTATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaa                                       <  1:2077300/71‑1 (MQ=255)
                    gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaa                                       <  1:3004090/71‑1 (MQ=255)
                    gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaa                                       <  1:3238496/71‑1 (MQ=255)
                     tggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaa                                       <  1:2497974/70‑1 (MQ=255)
                     tggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaa                                       >  1:3682488/1‑70 (MQ=255)
                      ggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaaa                                     >  1:3197207/1‑71 (MQ=255)
                       gCACGCGTACCGCGTAGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaaa                                     <  1:2267028/70‑1 (MQ=255)
                         aCGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACgg                                  >  1:170935/1‑71 (MQ=255)
                          cgcgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaa                                       <  1:1062244/65‑1 (MQ=255)
                           gcgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaaa                                     <  1:279479/66‑1 (MQ=255)
                           gcgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGc                                 >  1:1060026/1‑70 (MQ=255)
                           gcgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCa                                <  1:111486/71‑1 (MQ=255)
                            cgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAg                                                >  1:1836003/1‑54 (MQ=255)
                            cgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGc                                 >  1:1709862/1‑69 (MQ=255)
                            cgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCaa                               <  1:543172/71‑1 (MQ=255)
                             gTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtgg                                            >  1:3835250/1‑57 (MQ=255)
                              tACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaaa                                     <  1:789219/63‑1 (MQ=255)
                              tACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGc                             >  1:4518474/1‑71 (MQ=255)
                                ccGCGTCGATTACGGGTGCCTGCGCGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCt                            >  1:1512912/1‑70 (MQ=255)
                                 cgcgTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGa                          <  1:1550189/71‑1 (MQ=255)
                                  gcgTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaaa                                     >  1:2208051/1‑59 (MQ=255)
                                  gcgTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaa                                      >  1:2360245/1‑58 (MQ=255)
                                        aTTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCa                                <  1:3238938/58‑1 (MQ=255)
                                        aTTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCaa                   <  1:2653993/71‑1 (MQ=255)
                                         ttACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCaa                   >  1:2996593/1‑70 (MQ=255)
                                          tACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCaacc                 >  1:3835217/1‑69 (MQ=255)
                                                tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTg                           <  1:4106465/55‑1 (MQ=255)
                                                 gCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGa                                        >  1:2920593/1‑41 (MQ=255)
                                                  ccTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACg                                   >  1:999131/1‑45 (MQ=255)
                                                  ccTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGaa          <  1:1909233/70‑1 (MQ=255)
                                                   cTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATg            <  1:4110753/67‑1 (MQ=255)
                                                    tGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCa                                <  1:1983623/46‑1 (MQ=255)
                                                     gCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGAt             <  1:2950411/64‑1 (MQ=255)
                                                      cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg       <  1:3349747/69‑1 (MQ=255)
                                                      cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg       <  1:3004055/69‑1 (MQ=255)
                                                      cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg       <  1:2271140/69‑1 (MQ=255)
                                                       ggggCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGAt             <  1:4482454/60‑1 (MQ=255)
                                                       gTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGAt    >  1:2745132/1‑71 (MQ=255)
                                                          gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                         >  1:1434950/1‑47 (MQ=255)
                                                          gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                         >  1:954815/1‑47 (MQ=255)
                                                          gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg  <  1:2485265/70‑1 (MQ=255)
                                                                 |                                                              
CGACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGG  >  W3110S.gb/3824215‑3824342

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: