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

CTGCGACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAG  >  W3110S.gb/3824212‑3824345
                                                                    |                                                                 
cTGCGACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAg                                                                 <  1:1246287/71‑1 (MQ=255)
   cGACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGa                                                                <  1:4878867/69‑1 (MQ=255)
      cGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATgcgc                                                           <  1:3064544/71‑1 (MQ=255)
      cGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATgcgc                                                           >  1:4270890/1‑71 (MQ=255)
        tGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTa                                                         <  1:638871/71‑1 (MQ=255)
          cTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTAc                                                        <  1:2518781/70‑1 (MQ=255)
           ttCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATg                                                              <  1:4788309/63‑1 (MQ=255)
           ttCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATg                                                              <  1:2689315/63‑1 (MQ=255)
               ggCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGc                                                  <  1:3380461/71‑1 (MQ=255)
                      ggtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggtgg                                            >  1:3164198/1‑70 (MQ=255)
                       gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaa                                          <  1:4851484/71‑1 (MQ=255)
                          gCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAAc                                       <  1:3104521/71‑1 (MQ=255)
                            aCGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACgg                                     <  1:1836978/71‑1 (MQ=255)
                             cgcgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACg                                      <  1:3288753/69‑1 (MQ=255)
                             cgcgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGc                                    <  1:34196/71‑1 (MQ=255)
                              gcgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCa                                   <  1:1577826/71‑1 (MQ=255)
                               cgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGc                                    >  1:3243451/1‑69 (MQ=255)
                               cgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGc                                    >  1:2722966/1‑69 (MQ=255)
                                   ccGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCt                               >  1:999028/1‑70 (MQ=255)
                                   ccGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCt                               >  1:1182047/1‑70 (MQ=255)
                                     gcgTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                            <  1:4180513/71‑1 (MQ=255)
                                     gcgTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                            <  1:4857348/71‑1 (MQ=255)
                                     gcgTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                            <  1:1727339/71‑1 (MQ=255)
                                      cgTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaa                           <  1:4429795/71‑1 (MQ=255)
                                      cgTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaa                           <  1:2048060/71‑1 (MQ=255)
                                        tCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAg                                                   <  1:4931915/45‑1 (MQ=255)
                                        tCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAAc                         >  1:4651473/1‑71 (MQ=255)
                                        tCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGAGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAAc                         >  1:3413757/1‑71 (MQ=255)
                                          gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGc                                    >  1:1457683/1‑58 (MQ=255)
                                          gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaa                           >  1:2821527/1‑67 (MQ=255)
                                             tACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCaacc                    >  1:2511212/1‑69 (MQ=255)
                                             tACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCaacc                    >  1:4324177/1‑69 (MQ=255)
                                             tACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCaacc                    >  1:4418787/1‑69 (MQ=255)
                                             tACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCaacc                    >  1:67200/1‑69 (MQ=255)
                                              aCGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaaa                                        >  1:4571947/1‑50 (MQ=255)
                                               cGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc                   <  1:3232627/70‑1 (MQ=255)
                                               cGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc                   <  1:1661034/70‑1 (MQ=255)
                                                   tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCaa                                  <  1:3096392/51‑1 (MQ=255)
                                                   tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa              <  1:4711698/71‑1 (MQ=255)
                                                   tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa              <  1:4598021/71‑1 (MQ=255)
                                                       tGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaa                                          <  1:1433340/39‑1 (MQ=255)
                                                       tGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaa                                          <  1:1436348/39‑1 (MQ=255)
                                                        gCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCaa                      >  1:2841970/1‑58 (MQ=255)
                                                        gCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg          >  1:1071725/1‑70 (MQ=255)
                                                         cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa        >  1:3368245/1‑71 (MQ=255)
                                                         cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa        >  1:920586/1‑71 (MQ=255)
                                                         cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa        >  1:731573/1‑71 (MQ=255)
                                                          gTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                            <  1:4344282/50‑1 (MQ=255)
                                                          gTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGAt       >  1:4030951/1‑71 (MQ=255)
                                                           tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTg                              <  1:3483057/47‑1 (MQ=255)
                                                           tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCACCCCGATGAAAGGGATg      <  1:407112/71‑1 (MQ=255)
                                                           tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGAt                <  1:4444260/61‑1 (MQ=255)
                                                           tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATg      <  1:4596870/71‑1 (MQ=255)
                                                            ggCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                            >  1:1492533/1‑48 (MQ=255)
                                                               gcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa              <  1:3120995/59‑1 (MQ=255)
                                                               gcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAg  <  1:3572879/71‑1 (MQ=255)
                                                                 tGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                            <  1:5205875/43‑1 (MQ=255)
                                                                 tGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                            <  1:1090842/43‑1 (MQ=255)
                                                                 tGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg          <  1:3261338/61‑1 (MQ=255)
                                                                 tGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg          <  1:2895742/61‑1 (MQ=255)
                                                                    |                                                                 
CTGCGACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAG  >  W3110S.gb/3824212‑3824345

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: