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

ACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTG  >  W3110S.gb/3824258‑3824392
                                                                    |                                                                  
aCGGGTGCCTGCGTGGCGCTGGCAGATGAGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc                                                                  >  1:126460/1‑71 (MQ=255)
   ggTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGa                                                                >  1:448786/1‑70 (MQ=255)
       ccTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGaa                                                            <  1:4202291/70‑1 (MQ=255)
         tGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGa                                                                >  1:1314635/1‑64 (MQ=255)
         tGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGa                                                                >  1:2117336/1‑64 (MQ=255)
           cGTGGCGCTGGCAGATGCTCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAggg                                                        <  1:1404347/70‑1 (MQ=255)
            gTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGAt                                                      <  1:1797476/71‑1 (MQ=255)
             tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                             <  1:4791122/63‑1 (MQ=255)
             tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATg                                                     <  1:1593009/71‑1 (MQ=255)
               gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg                                                    <  1:1433368/70‑1 (MQ=255)
               gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt                                                   <  1:4318730/71‑1 (MQ=255)
                cgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg                                                    >  1:2205033/1‑69 (MQ=255)
                 gcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAg                                                 <  1:4242978/71‑1 (MQ=255)
                  cTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGc                                                >  1:634445/1‑71 (MQ=255)
                     gcaGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGc                                             <  1:2539217/69‑1 (MQ=255)
                     gcaGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCAGATGAAAGGGATGGTAGCCg                                              >  1:2672794/3‑70 (MQ=255)
                      caGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCg                                              >  1:1846939/2‑69 (MQ=255)
                      caGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCa                                            <  1:2399793/70‑1 (MQ=255)
                         aTGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg                                                         <  1:814563/55‑1 (MQ=255)
                         aTGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg                                                         <  1:979145/55‑1 (MQ=255)
                         aTGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGtt                                         >  1:4146270/1‑71 (MQ=255)
                          tGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                                        >  1:2005765/1‑71 (MQ=255)
                          tGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                                        >  1:2322245/1‑71 (MQ=255)
                           gcgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATg                                                              <  1:356109/48‑1 (MQ=255)
                           gcgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTc                                       >  1:1591470/1‑71 (MQ=255)
                           gcgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTc                                       >  1:2586992/1‑71 (MQ=255)
                           gcgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCAGATGAAAGGGATGGTAGCCGCAGtt                                         >  1:3735038/1‑69 (MQ=255)
                            cgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCt                                      <  1:4186392/71‑1 (MQ=255)
                             gcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                             <  1:4149957/47‑1 (MQ=255)
                             gcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                             <  1:1351307/47‑1 (MQ=255)
                              cTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                             <  1:2292316/46‑1 (MQ=255)
                              cTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                             <  1:873513/46‑1 (MQ=255)
                              cTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                             <  1:1069469/46‑1 (MQ=255)
                              cTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTg                                     >  1:1114626/1‑70 (MQ=255)
                                aCAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGaaa                                                           >  1:1408632/1‑46 (MQ=255)
                                aCAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAggg                                                        <  1:2895395/49‑1 (MQ=255)
                                aCAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAggg                                                        <  1:742973/49‑1 (MQ=255)
                                 cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGAt                                                      <  1:1577910/50‑1 (MQ=255)
                                  aGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGAt                                                      >  1:4741268/1‑49 (MQ=255)
                                   gAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGt                                          <  1:2744711/60‑1 (MQ=255)
                                    aaGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                             <  1:3497084/40‑1 (MQ=255)
                                    aaGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTa                                                  >  1:4503567/1‑51 (MQ=255)
                                    aaGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTa                                                  >  1:25148/1‑51 (MQ=255)
                                     aGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAtt                             >  1:3680275/1‑71 (MQ=255)
                                           tggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                          <  1:2144706/68‑2 (MQ=255)
                                             gAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg                     >  1:3952797/1‑71 (MQ=255)
                                               aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                   <  1:4236526/71‑3 (MQ=255)
                                                aaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt                                                   >  1:512133/1‑38 (MQ=255)
                                                aaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGa                  <  1:1135845/71‑1 (MQ=255)
                                                 aCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt                                                   <  1:2983898/37‑1 (MQ=255)
                                                   ggCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGaa                 >  1:756269/1‑69 (MQ=255)
                                                   ggCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGc               <  1:79219/71‑1 (MQ=255)
                                                   ggCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGc               <  1:2292187/71‑1 (MQ=255)
                                                     cAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCggt                                <  1:1288781/52‑2 (MQ=255)
                                                     cAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                          <  1:2525647/58‑2 (MQ=255)
                                                      aaGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaat            >  1:4226816/1‑68 (MQ=255)
                                                         cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                   <  1:1375978/61‑3 (MQ=255)
                                                         cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAg                <  1:494591/64‑1 (MQ=255)
                                                         cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTg         <  1:2035815/71‑1 (MQ=255)
                                                          tGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc        >  1:2328291/1‑71 (MQ=255)
                                                          tGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc        >  1:329634/1‑71 (MQ=255)
                                                           gAAAACCAACCCGATGAAAGGGATGGTAGCCGAAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc        >  1:3915910/1‑70 (MQ=255)
                                                            aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                   <  1:1856311/58‑3 (MQ=255)
                                                            aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg       <  1:383498/70‑1 (MQ=255)
                                                              aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt     <  1:620597/70‑1 (MQ=255)
                                                              aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATc    <  1:1702977/71‑1 (MQ=255)
                                                              aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATc    >  1:4108682/1‑71 (MQ=255)
                                                               aCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCt   <  1:319233/71‑1 (MQ=255)
                                                               aCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCt   >  1:4502242/1‑71 (MQ=255)
                                                                ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg  >  1:4411400/1‑71 (MQ=255)
                                                                ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg  <  1:4296311/71‑1 (MQ=255)
                                                                 cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCt   >  1:2810242/1‑69 (MQ=255)
                                                                    |                                                                  
ACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTG  >  W3110S.gb/3824258‑3824392

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: