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

ACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGG  >  W3110S.gb/3824217‑3824342
                                                               |                                                              
aCGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATgcg                                                         >  1:2748913/1‑71 (MQ=255)
aCGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGACTGCGTGGCGCTGGCAGATgcg                                                         >  1:221252/1‑71 (MQ=255)
aCGTGCTTCAGGCTGATGGTGGAACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATgcg                                                         >  1:3908445/1‑71 (MQ=255)
     cTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTa                                                      >  1:3329952/1‑69 (MQ=255)
     cTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTa                                                      >  1:1646785/1‑69 (MQ=255)
     cTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCCCTGGCAGATGCGCTa                                                      >  1:4933024/1‑69 (MQ=255)
       tCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGa                                                             >  1:1470512/1‑60 (MQ=255)
        cAGGCTGAAGGTGGCACGCGTACCGCGACGAATACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAg                                                   >  1:2174787/1‑69 (MQ=255)
         aGGCTGATGGTGGCACGCGTACCGCGTCCATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGaa                                                 <  1:1034464/70‑1 (MQ=255)
           gCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCt                                              <  1:184335/71‑1 (MQ=255)
           gCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCt                                              <  1:2109933/71‑1 (MQ=255)
           gCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCt                                              >  1:923798/1‑71 (MQ=255)
           gCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCt                                              >  1:1202865/1‑71 (MQ=255)
           gCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCt                                              <  1:5141750/71‑1 (MQ=255)
           gCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCt                                              <  1:4099768/71‑1 (MQ=255)
            cTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtg                                             >  1:3986966/1‑71 (MQ=255)
            cTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtg                                             >  1:291786/1‑71 (MQ=255)
            cTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtg                                             >  1:1676969/1‑71 (MQ=255)
             tGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGaa                                                 <  1:1914477/66‑1 (MQ=255)
             tGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGaa                                                 <  1:401021/66‑1 (MQ=255)
               aTGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggtg                                          >  1:252537/1‑71 (MQ=255)
               aTGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggtg                                          >  1:845658/1‑71 (MQ=255)
               aTGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggtg                                          >  1:3459863/1‑71 (MQ=255)
                 ggtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggtgg                                         >  1:2290317/1‑70 (MQ=255)
                 ggtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggtgg                                         >  1:3464006/1‑70 (MQ=255)
                   tggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATgc                                                          <  1:939383/51‑1 (MQ=255)
                            tACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggtg                                          <  1:4942075/58‑1 (MQ=255)
                                 cgTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                         >  1:2278963/1‑70 (MQ=255)
                                   tCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                         <  1:1417619/68‑1 (MQ=255)
                                     gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaaa                       >  1:2532353/1‑68 (MQ=255)
                                     gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaaa                       >  1:301203/1‑68 (MQ=255)
                                     gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaaa                       >  1:3385562/1‑68 (MQ=255)
                                     gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAAcc                     >  1:1068651/1‑70 (MQ=255)
                                        tACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCaacc                 >  1:59054/1‑69 (MQ=255)
                                         aCGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaaa                       >  1:4376093/1‑64 (MQ=255)
                                          cGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc                <  1:3254817/70‑1 (MQ=255)
                                          cGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc                <  1:1811421/70‑1 (MQ=255)
                                              tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa           <  1:3921587/71‑1 (MQ=255)
                                              tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa           <  1:2313830/71‑1 (MQ=255)
                                                ccTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCt                            >  1:3520482/1‑52 (MQ=255)
                                                ccTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTg                           >  1:902144/1‑53 (MQ=255)
                                                 cTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCaa                   >  1:2736004/1‑60 (MQ=255)
                                                   gCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg       >  1:2893815/1‑70 (MQ=255)
                                                    cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAggg      <  1:399752/70‑1 (MQ=255)
                                                    cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa     >  1:2343706/1‑71 (MQ=255)
                                                       ggCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGa                          <  1:5082570/47‑1 (MQ=255)
                                                       ggCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGa                          <  1:4312257/47‑1 (MQ=255)
                                                        gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg  >  1:3842272/1‑70 (MQ=255)
                                                          gcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                         <  1:2759748/45‑1 (MQ=255)
                                                          gcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                         <  1:1746509/45‑1 (MQ=255)
                                                           cTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg       <  1:543844/62‑1 (MQ=255)
                                                           cTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg       <  1:765549/62‑1 (MQ=255)
                                                           cTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg       <  1:3047220/62‑1 (MQ=255)
                                                               |                                                              
ACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGG  >  W3110S.gb/3824217‑3824342

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: