Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,824,326 T→C 100% pseudogene (459/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,3260TC100.0% 154.1 / NA 47pseudogene (459/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/29);  total (18/29)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGG  >  W3110S.gb/3824261‑3824393
                                                                 |                                                                   
ggTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGAt                                                                <  1:2424002/71‑1 (MQ=255)
  tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                              <  1:1337312/71‑1 (MQ=255)
     cTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCg                                                                  <  1:1786158/64‑1 (MQ=255)
        cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa                                                        >  1:19330/1‑71 (MQ=255)
          tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGAt                                                       <  1:2817274/70‑1 (MQ=255)
          tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATg                                                      <  1:3427517/71‑1 (MQ=255)
            gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg                                                     >  1:1650181/1‑70 (MQ=255)
            gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg                                                     <  1:4051569/70‑1 (MQ=255)
            gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt                                                    <  1:1716395/71‑1 (MQ=255)
                   caGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGa                                                                 >  1:2654422/2‑51 (MQ=255)
                     gATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGt                                           <  1:3108271/71‑1 (MQ=255)
                        gcgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGaaa                                                            >  1:1769496/1‑51 (MQ=255)
                          gcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGt                                           <  1:1370318/66‑1 (MQ=255)
                              cAGAAGCTGGTGGAGAACGGCAAGCTGAAAACCAACCCGATGAAAgg                                                          <  1:2922694/47‑1 (MQ=39)
                              cAGAAGCTGGTGGAAACCGGCAAGCTGAAAACCAACCCGATGAAAgg                                                          <  1:2124760/47‑1 (MQ=39)
                              cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                                  <  1:2211256/71‑1 (MQ=255)
                              cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                                  <  1:4050499/71‑1 (MQ=255)
                                aaaGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCg                                                                  <  1:2885222/36‑1 (MQ=38)
                                 aaGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAt                               >  1:1163782/1‑71 (MQ=255)
                                  aGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTTTCGGTAtt                              >  1:647909/1‑71 (MQ=255)
                                    cTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                              <  1:3080577/37‑1 (MQ=255)
                                     tggtggAAAACGGCAAGCTGAAAACCAACCCGATGaaa                                                            >  1:3281293/1‑38 (MQ=255)
                                     tggtggAAAACGGCAAGCTGAAAACCAACCCGATGa                                                              >  1:1755939/1‑36 (MQ=255)
                                      ggtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                           >  1:1047394/1‑69 (MQ=255)
                                       gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                           <  1:1029863/69‑2 (MQ=255)
                                          gAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACg                       <  1:3073902/70‑1 (MQ=255)
                                            aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGtt                                          <  1:2132721/49‑1 (MQ=255)
                                            aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                           >  1:907218/1‑63 (MQ=255)
                                            aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                           >  1:1910478/1‑63 (MQ=255)
                                            aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                    <  1:3928350/71‑3 (MQ=255)
                                            aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                    <  1:2024184/71‑3 (MQ=255)
                                             aaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                    >  1:274404/1‑68 (MQ=255)
                                              aCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                    >  1:4196372/1‑67 (MQ=255)
                                                ggCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAg                 >  1:1330687/1‑70 (MQ=255)
                                                   aaGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                                  <  1:191196/50‑1 (MQ=255)
                                                   aaGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                           <  1:3124732/57‑2 (MQ=255)
                                                   aaGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaat             <  1:4306097/71‑4 (MQ=255)
                                                        gAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc         >  1:232767/1‑70 (MQ=255)
                                                        gAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc         >  1:3694425/1‑70 (MQ=255)
                                                         aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa       >  1:2329814/1‑71 (MQ=255)
                                                         aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa       <  1:1813908/71‑1 (MQ=255)
                                                         aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa       <  1:2884486/71‑1 (MQ=255)
                                                             ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg   <  1:202790/71‑1 (MQ=255)
                                                              cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg  <  1:3809293/71‑1 (MQ=255)
                                                              cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg  <  1:1846119/71‑1 (MQ=255)
                                                              cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg  <  1:2392574/71‑1 (MQ=255)
                                                              cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg   >  1:1456726/1‑70 (MQ=255)
                                                                 |                                                                   
GGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGG  >  W3110S.gb/3824261‑3824393

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: