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% 144.7 / NA 43pseudogene (413/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/16);  total (27/16)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCGACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGT  >  W3110S.gb/3824214‑3824343
                                                                  |                                                               
gCGACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAg                                                               >  1:702835/1‑69 (MQ=255)
gCGACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAg                                                               >  1:1490693/1‑69 (MQ=255)
      tGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTATGGGTGCCTGCGTGGCGCTGGCAGATgcgc                                                         >  1:1296586/1‑69 (MQ=255)
      tGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATgcgc                                                         >  1:1039049/1‑69 (MQ=255)
      tGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATgcgc                                                         >  1:494708/1‑69 (MQ=255)
        cTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATgcgc                                                         >  1:2460750/1‑67 (MQ=255)
        cTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATgcgc                                                         >  1:93232/1‑67 (MQ=255)
        cTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATgcgc                                                         >  1:315209/1‑67 (MQ=255)
             ggCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGaa                                                  >  1:1602179/1‑69 (MQ=255)
                     gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggtg                                           >  1:1831971/1‑68 (MQ=255)
                     gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggtg                                           >  1:3239280/1‑68 (MQ=255)
                     gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggtg                                           >  1:496077/1‑68 (MQ=255)
                     gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggtg                                           >  1:605905/1‑68 (MQ=255)
                     gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggtg                                           >  1:1060064/1‑68 (MQ=255)
                             cgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACg                                    >  1:2969350/1‑67 (MQ=255)
                             cgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACg                                    >  1:1637114/1‑67 (MQ=255)
                             cgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACg                                    >  1:1764112/1‑67 (MQ=255)
                                 ccGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAg                               >  1:2971662/1‑68 (MQ=255)
                                 ccGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAg                               >  1:955062/1‑68 (MQ=255)
                                 ccGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAg                               >  1:1846183/1‑68 (MQ=255)
                                     gTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                          <  1:2485771/69‑1 (MQ=255)
                                        gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaaa                        >  1:1201104/1‑68 (MQ=255)
                                        gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaaa                        >  1:1322688/1‑68 (MQ=255)
                                        gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaaa                        >  1:1349259/1‑68 (MQ=255)
                                        gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaaa                        >  1:3519141/1‑68 (MQ=255)
                                        gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaaa                        >  1:3225336/1‑68 (MQ=255)
                                          ttACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAAcc                      <  1:1764927/68‑1 (MQ=255)
                                          ttACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAAcc                      <  1:2772312/68‑1 (MQ=255)
                                          ttACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAAcc                      <  1:161352/68‑1 (MQ=255)
                                          ttACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAAcc                      <  1:1604124/68‑1 (MQ=255)
                                           tACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGGTGGTGGAAAACGGCAAGCTGAAAACCaa                    >  1:530537/1‑69 (MQ=255)
                                               ggTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaa                         <  1:2969318/60‑1 (MQ=255)
                                               ggTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCg                <  1:2041035/69‑1 (MQ=255)
                                                     tGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAAcc                      >  1:1541582/1‑57 (MQ=255)
                                                       cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAAcc                      <  1:1337443/55‑1 (MQ=255)
                                                       cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCaac                   <  1:3393801/58‑2 (MQ=255)
                                                       cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCaac                   <  1:1510400/58‑2 (MQ=255)
                                                       cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCaac                   <  1:1615425/58‑2 (MQ=255)
                                                         tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                          <  1:1156158/49‑1 (MQ=255)
                                                         tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg        <  1:2332350/67‑1 (MQ=255)
                                                             gcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt  <  1:3050358/69‑1 (MQ=255)
                                                             gcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt  <  1:2917845/69‑1 (MQ=255)
                                                             gcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt  <  1:2837151/69‑1 (MQ=255)
                                                                  |                                                               
GCGACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGT  >  W3110S.gb/3824214‑3824343

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