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

CGACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCC  >  W3110S.gb/3824215‑3824347
                                                                 |                                                                   
cGACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGa                                                                  <  1:3065413/69‑1 (MQ=255)
          cAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGa                                                       <  1:1643744/70‑1 (MQ=255)
           aGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATgc                                                               <  1:1113323/61‑1 (MQ=255)
                gATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggt                                                >  1:2485621/1‑71 (MQ=255)
                gATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggt                                                >  1:1775167/1‑71 (MQ=255)
                gATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggt                                                >  1:1344937/1‑71 (MQ=255)
                 aTGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggtg                                               >  1:948298/1‑71 (MQ=255)
                 aTGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggtg                                               >  1:708813/1‑71 (MQ=255)
                    gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaa                                            <  1:1236952/71‑1 (MQ=255)
                    gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaa                                            <  1:3280016/71‑1 (MQ=255)
                    gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGa                                             <  1:1550942/70‑1 (MQ=255)
                           gcgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaaa                                          >  1:3038427/1‑66 (MQ=255)
                               aCCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCt                                 >  1:179802/1‑71 (MQ=255)
                               aCCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCt                                 >  1:1012212/1‑71 (MQ=255)
                                ccGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCt                                 >  1:1298585/1‑70 (MQ=255)
                                     tCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggtgg                                              <  1:754343/52‑1 (MQ=255)
                                     tCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggtgg                                              <  1:3364923/52‑1 (MQ=255)
                                     tCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggtgg                                              <  1:1097633/52‑1 (MQ=255)
                                       gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAAcc                          <  1:45016/70‑1 (MQ=255)
                                       gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAAcc                          <  1:3642515/70‑1 (MQ=255)
                                       gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAAcc                          >  1:3282497/1‑70 (MQ=255)
                                       gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAAcc                          <  1:1275078/70‑1 (MQ=255)
                                         ttACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCaa                                    <  1:2454565/58‑1 (MQ=255)
                                          tACGGGTGCCTGCGTGGCGCTGGCATATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCaacc                      >  1:917324/1‑69 (MQ=255)
                                          tACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCaacc                      >  1:1007236/1‑69 (MQ=255)
                                           aCGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCaacc                      <  1:3357242/70‑3 (MQ=255)
                                                tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                <  1:2159072/71‑1 (MQ=255)
                                                tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                <  1:1560826/71‑1 (MQ=255)
                                                tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                <  1:1543915/71‑1 (MQ=255)
                                                    tGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAAcc                          <  1:1429375/57‑1 (MQ=255)
                                                    tGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATg                 <  1:2975730/66‑1 (MQ=255)
                                                     gCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGaaa              >  1:197076/1‑68 (MQ=255)
                                                       gTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGTAAAACGGCAATCTGAAAACCAACCCGATGAAAgg            <  1:953326/68‑1 (MQ=255)
                                                       gTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATg                 >  1:2477765/1‑63 (MQ=255)
                                                          gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt      <  1:2863484/71‑1 (MQ=255)
                                                            gcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAg    <  1:2041328/71‑1 (MQ=255)
                                                             cTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                <  1:1762999/58‑1 (MQ=255)
                                                             cTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                <  1:850618/58‑1 (MQ=255)
                                                              tGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATg                 <  1:2896570/56‑1 (MQ=255)
                                                              tGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGAt         >  1:617483/1‑64 (MQ=255)
                                                              tGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGcc  >  1:1470786/1‑71 (MQ=255)
                                                                 |                                                                   
CGACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCC  >  W3110S.gb/3824215‑3824347

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: