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

TGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATC  >  W3110S.gb/3824263‑3824390
                                                               |                                                                
tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                           <  1:2159072/71‑1 (MQ=255)
tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                           <  1:1560826/71‑1 (MQ=255)
tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                           <  1:1543915/71‑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:850618/58‑1 (MQ=255)
             cTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                           <  1:1762999/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)
                   gATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa                                                     <  1:1946538/58‑1 (MQ=255)
                   gATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGt                                        <  1:2618996/71‑1 (MQ=255)
                    aTGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAACCGCAGtt                                       <  1:2747403/71‑1 (MQ=255)
                     tGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACACGATGAAAGGGATg                                                   >  1:3561484/1‑58 (MQ=255)
                         cTACTGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGt                                  <  1:1886949/71‑1 (MQ=255)
                         cTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                           <  1:3192599/46‑1 (MQ=255)
                         cTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGt                                  >  1:3431467/1‑71 (MQ=255)
                            cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                           <  1:1931399/43‑1 (MQ=255)
                            cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                               <  1:311461/71‑1 (MQ=255)
                            cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                               <  1:2353987/71‑1 (MQ=255)
                            cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                               >  1:358197/1‑71 (MQ=255)
                             aGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                               >  1:1533798/1‑70 (MQ=255)
                              gAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCggta                             >  1:1367151/1‑69 (MQ=255)
                                 gCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                           <  1:3610015/38‑1 (MQ=255)
                                   tggtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGc                                           <  1:2792680/52‑1 (MQ=255)
                                     gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                        <  1:2560864/69‑2 (MQ=255)
                                     gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                        <  1:2850429/69‑2 (MQ=255)
                                     gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                      >  1:1640399/1‑71 (MQ=255)
                                      tggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTc                                     <  1:3474123/55‑1 (MQ=255)
                                      tggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTc                                     <  1:1235114/55‑1 (MQ=255)
                                      tggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                      >  1:2743344/1‑70 (MQ=255)
                                      tggAAAACGGCAAGCCGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                     <  1:2246619/71‑1 (MQ=255)
                                       ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTc                                     <  1:3652198/54‑1 (MQ=255)
                                       ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                     <  1:1785498/70‑1 (MQ=255)
                                        gAAAACGGCAAGCTGAAAACCCACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAt                            >  1:619153/1‑62 (MQ=255)
                                        gAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGt                                  <  1:2635259/56‑1 (MQ=255)
                                        gAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAt                            >  1:2619663/1‑62 (MQ=255)
                                          aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGa                                                     >  1:859371/1‑35 (MQ=255)
                                           aaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                                      <  1:153073/49‑1 (MQ=255)
                                             cGGCAAGCTGAAAACCAACCCGATTAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                        <  1:2790087/61‑2 (MQ=255)
                                             cGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                        <  1:1324876/61‑2 (MQ=255)
                                                cAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                               <  1:387336/51‑1 (MQ=255)
                                                 aaGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACg                    >  1:835850/1‑61 (MQ=255)
                                                   gCTGAAAACCAACCCGATGAAAGGGATTGTAGCCGCAGTTTCTGTCgg                               <  1:306996/48‑1 (MQ=255)
                                                   gCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                 <  1:713490/62‑3 (MQ=255)
                                                      gAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                      <  1:233126/54‑1 (MQ=255)
                                                      gAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc      >  1:2322201/1‑70 (MQ=255)
                                                       aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca            >  1:1680602/1‑62 (MQ=255)
                                                       aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg     <  1:2938087/70‑1 (MQ=255)
                                                       aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg     <  1:2978456/70‑1 (MQ=255)
                                                       aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa    >  1:952803/1‑71 (MQ=255)
                                                        aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                     >  1:1595356/1‑53 (MQ=255)
                                                        aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa    <  1:2198851/70‑1 (MQ=255)
                                                         aaccaaccCGATTAAAGGGTTGGTAGCCGCAGTTTCTGTCggt                              <  1:1327648/43‑2 (MQ=37)
                                                         aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATc  >  1:1453574/1‑71 (MQ=255)
                                                          aCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAg              >  1:1311104/1‑58 (MQ=255)
                                                               |                                                                
TGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATC  >  W3110S.gb/3824263‑3824390

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: