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

CGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATC  >  W3110S.gb/3824259‑3824390
                                                                   |                                                                
cGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc                                                                <  1:2118285/70‑1 (MQ=255)
cGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc                                                                <  1:740737/70‑1 (MQ=255)
cGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCg                                                               <  1:725980/71‑1 (MQ=255)
  ggTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGAt                                                             <  1:395804/71‑1 (MQ=255)
    tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGAt                                                             <  1:784379/69‑1 (MQ=255)
    tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                           <  1:1435681/71‑1 (MQ=255)
      ccTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc                                                                >  1:1995111/1‑64 (MQ=255)
      ccTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGaa                                                          <  1:109752/70‑1 (MQ=255)
      ccTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                           >  1:3587834/1‑69 (MQ=255)
          cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAggg                                                      <  1:3551429/70‑1 (MQ=255)
              gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg                                                  <  1:487478/70‑1 (MQ=255)
              gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt                                                 <  1:215592/71‑1 (MQ=255)
               cgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg                                                  >  1:2831145/1‑69 (MQ=255)
               cgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg                                                  >  1:189376/1‑69 (MQ=255)
               cgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt                                                 <  1:509620/70‑1 (MQ=255)
                     caGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCg                                            >  1:3530331/2‑69 (MQ=255)
                     caGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCg                                            >  1:822226/2‑69 (MQ=255)
                      aGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAg                                         >  1:3989429/1‑71 (MQ=255)
                      aGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAg                                         >  1:2836336/1‑71 (MQ=255)
                        aTGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                           <  1:313897/51‑1 (MQ=255)
                          gcgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGaaa                                                         >  1:1112101/1‑51 (MQ=255)
                          gcgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATg                                                   >  1:2722420/1‑57 (MQ=255)
                          gcgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTc                                     >  1:1364936/1‑71 (MQ=255)
                           cgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGAt                                                             >  1:343662/1‑46 (MQ=255)
                           cgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCt                                    <  1:1749074/71‑1 (MQ=255)
                           cgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCt                                    <  1:223088/71‑1 (MQ=255)
                             cTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATg                                                            <  1:3946957/45‑1 (MQ=255)
                             cTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATg                                                            <  1:179968/45‑1 (MQ=255)
                              tACAGAAGCTGGTGGAAAACGGGAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTc                                 >  1:1511354/1‑71 (MQ=255)
                              tACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGc                                           <  1:49516/61‑1 (MQ=255)
                              tACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGc                                           <  1:775018/61‑1 (MQ=255)
                              tACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTc                                 >  1:3062115/1‑71 (MQ=255)
                                cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGc                                              >  1:3156292/1‑56 (MQ=255)
                                cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                               <  1:2626496/71‑1 (MQ=255)
                                cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                               >  1:3136240/1‑71 (MQ=255)
                                 aGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGAt                                                             >  1:2754358/1‑40 (MQ=255)
                                 aGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                           <  1:4137678/42‑1 (MQ=255)
                                 aGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                                      <  1:1994169/63‑1 (MQ=255)
                                 aGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                               >  1:658073/1‑70 (MQ=255)
                                  gAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCggta                             <  1:3083927/71‑3 (MQ=255)
                                  gAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCggta                             <  1:1972966/71‑3 (MQ=255)
                                  gAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCggta                             <  1:2613107/71‑3 (MQ=255)
                                    aGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGt                                        >  1:455353/1‑58 (MQ=255)
                                      cTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgt                         <  1:3434075/71‑1 (MQ=255)
                                        ggtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCggta                             >  1:968059/1‑63 (MQ=255)
                                         gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                        <  1:3480215/69‑2 (MQ=255)
                                          tggAAAACGGCAAGCTGAAAACCAACCCGATGAAAggg                                                      >  1:2774940/1‑38 (MQ=255)
                                            gAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg                   <  1:933410/71‑1 (MQ=255)
                                            gAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg                   <  1:3942924/71‑1 (MQ=255)
                                             aaaaCGGCAAGCTGAAACCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg                   >  1:1365612/1‑70 (MQ=255)
                                             aaaaCGGCAAGCTGAAAACCAACCCGATGAAAggg                                                      <  1:2097004/35‑1 (MQ=255)
                                             aaaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGc                                           <  1:538255/46‑1 (MQ=255)
                                             aaaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggt                  <  1:2504843/71‑2 (MQ=255)
                                             aaaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggt                  >  1:474976/1‑70 (MQ=255)
                                              aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                 <  1:1064469/71‑3 (MQ=255)
                                              aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                 <  1:1288315/71‑3 (MQ=255)
                                              aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                 <  1:2599340/71‑3 (MQ=255)
                                              aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggt                  >  1:1808823/1‑69 (MQ=255)
                                                aCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                 >  1:2712488/1‑67 (MQ=255)
                                                 cGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAg              >  1:2026859/1‑71 (MQ=255)
                                                    cAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                     <  1:3574529/61‑1 (MQ=255)
                                                    cAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaa           <  1:4112113/71‑3 (MQ=255)
                                                    cAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaa           <  1:2028127/71‑3 (MQ=255)
                                                     aaGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaat          <  1:3196460/71‑4 (MQ=255)
                                                     aaGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaat          <  1:659530/71‑4 (MQ=255)
                                                      aGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg                   >  1:2208659/1‑61 (MQ=255)
                                                       gCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGc             <  1:449701/66‑1 (MQ=255)
                                                        cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                      >  1:2522556/1‑56 (MQ=255)
                                                        cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGaa               >  1:433901/1‑63 (MQ=255)
                                                         tGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                 <  1:1763078/60‑3 (MQ=255)
                                                          gAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGaa               >  1:1136218/1‑61 (MQ=255)
                                                           aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGa                >  1:481760/1‑59 (MQ=255)
                                                           aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa    >  1:4048416/1‑71 (MQ=255)
                                                             aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATc  <  1:3455567/71‑1 (MQ=255)
                                                             aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATc  >  1:2731795/1‑71 (MQ=255)
                                                               ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCggta                             <  1:2746249/42‑3 (MQ=38)
                                                               ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCggta                             <  1:4059192/42‑3 (MQ=38)
                                                                   |                                                                
CGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATC  >  W3110S.gb/3824259‑3824390

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: