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

CTGCGACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGC  >  W3110S.gb/3824212‑3824346
                                                                    |                                                                  
cTGCGACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAg                                                                  <  1:1986598/71‑1 (MQ=255)
cTGCGACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAg                                                                  <  1:5130948/71‑1 (MQ=255)
  gCGACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGAt                                                                >  1:489404/1‑71 (MQ=255)
      cGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATgcgc                                                            >  1:4742837/1‑71 (MQ=255)
         gCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTa                                                          >  1:4033413/1‑70 (MQ=255)
            tCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAg                                                       >  1:2732819/1‑70 (MQ=255)
              atgCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAg                                                    <  1:149289/69‑1 (MQ=255)
              aGGCTGATGGTGGCACGCGTAGCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGa                                                                 <  1:1402241/58‑1 (MQ=39)
              aGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGaa                                                     <  1:3328833/70‑1 (MQ=255)
                gCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTCCAGAAGCt                                                  <  1:1413940/71‑1 (MQ=255)
                gCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCGGGCAGATGCGCTACAGAAGCt                                                  <  1:5065276/71‑1 (MQ=255)
                 cTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtg                                                 <  1:618147/71‑1 (MQ=255)
                  tGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtgg                                                >  1:2521894/1‑71 (MQ=255)
                    aTGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggtg                                              >  1:740968/1‑71 (MQ=255)
                    aTGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtgg                                                <  1:3948107/69‑1 (MQ=255)
                         ggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaa                                          >  1:2602929/1‑70 (MQ=255)
                          gCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaaa                                         <  1:2924412/70‑1 (MQ=255)
                          gCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaaa                                         <  1:449517/70‑1 (MQ=255)
                           cACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACg                                       <  1:3226237/71‑1 (MQ=255)
                              gcgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGc                                     >  1:390504/1‑70 (MQ=255)
                              gcgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGc                                     >  1:2533306/1‑70 (MQ=255)
                               cgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACg                                       >  1:2669063/1‑67 (MQ=255)
                                 tACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGc                                 <  1:3480152/71‑1 (MQ=255)
                                 tACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGc                                 <  1:3750792/71‑1 (MQ=255)
                                 tACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGc                                 <  1:2453361/71‑1 (MQ=255)
                                  aCCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaaa                                         >  1:3440680/1‑62 (MQ=255)
                                   ccGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCt                                >  1:2603890/1‑70 (MQ=255)
                                    cgcgTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGa                              <  1:1834262/71‑1 (MQ=255)
                                      cgTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                             <  1:2123165/70‑1 (MQ=255)
                                       gTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaaa                           >  1:3826597/1‑71 (MQ=255)
                                        tCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACgg                                      >  1:3070369/1‑59 (MQ=255)
                                        tCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAAc                          <  1:1366721/71‑1 (MQ=255)
                                        tCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAAc                          <  1:1643338/71‑1 (MQ=255)
                                         cGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAAc                          >  1:390828/1‑70 (MQ=255)
                                          gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCa                                    >  1:4357114/1‑59 (MQ=255)
                                          gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCa                                    >  1:1708013/1‑59 (MQ=255)
                                          gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCa                        <  1:971221/71‑1 (MQ=255)
                                          gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCa                        <  1:107579/71‑1 (MQ=255)
                                           aTTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCaa                       >  1:2097268/1‑71 (MQ=255)
                                             tACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtgg                                                <  1:2863478/44‑1 (MQ=255)
                                             tACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaa                                          <  1:2500321/50‑1 (MQ=255)
                                             tACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACg                                       >  1:175179/1‑53 (MQ=255)
                                             tACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaa                            >  1:143246/1‑64 (MQ=255)
                                             tACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCaacc                     >  1:1007337/1‑69 (MQ=255)
                                               cGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc                    <  1:3093114/70‑1 (MQ=255)
                                               cGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc                    <  1:4660798/70‑1 (MQ=255)
                                               cGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc                    <  1:4450720/70‑1 (MQ=255)
                                                 ggTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGAt                 <  1:4706494/71‑1 (MQ=255)
                                                 ggTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGAt                 <  1:4700747/71‑1 (MQ=255)
                                                 ggTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGAt                 <  1:4995134/71‑1 (MQ=255)
                                                  gTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGa                  <  1:1188667/69‑1 (MQ=255)
                                                   tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa               <  1:2839921/71‑1 (MQ=255)
                                                    gCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa               >  1:4331568/1‑70 (MQ=255)
                                                     ccTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGaaa             <  1:148901/71‑1 (MQ=255)
                                                     ccTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGaa              <  1:4858775/70‑1 (MQ=255)
                                                     ccTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa               <  1:2431608/69‑1 (MQ=255)
                                                      cTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAg            >  1:2152120/1‑71 (MQ=255)
                                                       tGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCt                                <  1:1862143/50‑1 (MQ=255)
                                                       tGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATg                <  1:3475845/66‑1 (MQ=255)
                                                       tGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg           >  1:256943/1‑71 (MQ=255)
                                                       tGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAg            >  1:833112/1‑70 (MQ=255)
                                                        gCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGa                              >  1:893152/1‑51 (MQ=255)
                                                        gCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa               <  1:1342032/66‑1 (MQ=255)
                                                        gCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTAGAAAACGGCAAGCTGAAAAc                          >  1:1366732/1‑55 (MQ=39)
                                                         cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaaa                           >  1:2107888/1‑53 (MQ=255)
                                                         cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaaa                           >  1:3700947/1‑53 (MQ=255)
                                                         cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGaaa             >  1:4128498/1‑67 (MQ=255)
                                                         cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa         <  1:1950969/71‑1 (MQ=255)
                                                         cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa         <  1:3308488/71‑1 (MQ=255)
                                                         cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa         >  1:2908556/1‑71 (MQ=255)
                                                           tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATg       <  1:4586568/71‑1 (MQ=255)
                                                           tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATg       <  1:841418/71‑1 (MQ=255)
                                                            ggCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGAt                 >  1:1054700/1‑60 (MQ=255)
                                                             gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt     <  1:2365600/71‑1 (MQ=255)
                                                               gcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                             >  1:4496766/1‑45 (MQ=255)
                                                               gcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa               >  1:4835398/1‑59 (MQ=255)
                                                                cTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGc  <  1:3034536/71‑1 (MQ=255)
                                                                    |                                                                  
CTGCGACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGC  >  W3110S.gb/3824212‑3824346

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: