Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,824,410 T→A 100% pseudogene (543/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,4100TA100.0% 192.4 / NA 66pseudogene (543/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (34/32);  total (34/32)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGACGG  >  W3110S.gb/3824343‑3824477
                                                                   |                                                                   
tAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                                                                   <  1:2969009/70‑1 (MQ=255)
   ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                               <  1:2591528/71‑1 (MQ=255)
   ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                               <  1:3312934/71‑1 (MQ=255)
   ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                               <  1:4256024/71‑1 (MQ=255)
     gCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                               >  1:795706/1‑69 (MQ=255)
      cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                               >  1:4011926/1‑68 (MQ=255)
        gTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                               <  1:1526810/66‑1 (MQ=255)
           tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                                                             >  1:1660460/1‑65 (MQ=255)
           tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAccta                                                       >  1:3942372/1‑69 (MQ=255)
           tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                                                       >  1:3277817/1‑71 (MQ=255)
            cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                                                       <  1:4522083/70‑1 (MQ=255)
               tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg                                                   <  1:773061/71‑1 (MQ=255)
               tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg                                                   <  1:3782366/71‑1 (MQ=255)
                  gtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAc                                                          <  1:538908/59‑1 (MQ=255)
                   taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                                                             >  1:1362981/2‑57 (MQ=255)
                   taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                               <  1:829932/70‑1 (MQ=255)
                    aTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACAt                                                    >  1:2758869/1‑65 (MQ=255)
                    aTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                               >  1:2579504/1‑70 (MQ=255)
                     ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg                                                   <  1:3587321/65‑1 (MQ=255)
                     ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                               >  1:1788059/1‑69 (MQ=255)
                           aCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgatg                                       <  1:1483511/71‑1 (MQ=255)
                           aCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                        >  1:1049470/1‑70 (MQ=255)
                              gtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                               <  1:1668691/42‑1 (MQ=25)
                              gtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                               <  1:3152273/58‑1 (MQ=255)
                               tgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGAcc                                    >  1:1161518/2‑70 (MQ=255)
                                gAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCg                                   >  1:4245443/1‑70 (MQ=255)
                                     aatttGCGATCTGGAATACGTTGAAGACTGAGCCGCAGAGACCGACATGAACGTAGtgatg                                       <  1:150451/59‑1 (MQ=255)
                                     aatttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGAc                              <  1:1481654/68‑1 (MQ=255)
                                      atttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGAc                              <  1:3495254/68‑1 (MQ=255)
                                      atttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACg                             >  1:2107723/2‑70 (MQ=255)
                                       tttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGa                                                  <  1:1233308/48‑1 (MQ=255)
                                       tttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACTGAAGAc                              <  1:584848/68‑1 (MQ=255)
                                        ttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTa                                             >  1:2206059/1‑52 (MQ=255)
                                         tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGa                               >  1:2883728/1‑65 (MQ=255)
                                         tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGAc                              >  1:2933888/1‑66 (MQ=255)
                                          gCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgcg                         >  1:2329622/1‑70 (MQ=255)
                                          gCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgcg                         >  1:3413700/1‑70 (MQ=255)
                                           cGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg                           >  1:4397491/1‑67 (MQ=255)
                                           cGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgcgc                        <  1:1905389/70‑1 (MQ=255)
                                           cGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGAAGACCGAAGACggg                           >  1:1049845/1‑67 (MQ=255)
                                             aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                           <  1:4186445/49‑1 (MQ=255)
                                              tCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                        >  1:1749905/1‑51 (MQ=255)
                                               cTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACg                             >  1:3237454/1‑61 (MQ=255)
                                               cTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatcat                   >  1:3963532/1‑71 (MQ=255)
                                                tGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGa                                                  <  1:286740/39‑1 (MQ=255)
                                                 ggAATACGTTGAAGACTCAGCCGCAGAGACCGACATGa                                                  <  1:922894/38‑1 (MQ=255)
                                                  gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatca                    <  1:561108/67‑1 (MQ=255)
                                                  gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatca                    <  1:3127039/67‑1 (MQ=255)
                                                  gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCAtt                  >  1:3661891/1‑69 (MQ=255)
                                                  gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg                 >  1:2092264/1‑70 (MQ=255)
                                                   aaTACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCAtt                  <  1:1699866/68‑1 (MQ=255)
                                                     tACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCg                                   >  1:92574/1‑49 (MQ=255)
                                                      aCGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg            >  1:429880/1‑71 (MQ=255)
                                                       cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGAc                                     >  1:4233352/1‑45 (MQ=255)
                                                        gTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg            <  1:775328/69‑1 (MQ=255)
                                                         ttGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                        <  1:3269263/40‑1 (MQ=255)
                                                         ttGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                        <  1:4199087/40‑1 (MQ=255)
                                                         ttGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                        <  1:3826804/40‑1 (MQ=255)
                                                         ttGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAg         <  1:486531/71‑1 (MQ=255)
                                                          tGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAg              >  1:4283872/1‑65 (MQ=255)
                                                           gAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatc                     >  1:3635515/1‑57 (MQ=255)
                                                           gAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAggg       >  1:3597218/1‑71 (MQ=255)
                                                            aaGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGAc                                     >  1:171072/1‑40 (MQ=255)
                                                            aaGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAggg       <  1:559336/70‑1 (MQ=255)
                                                             aGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGa     >  1:652963/1‑71 (MQ=255)
                                                             aGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGa     >  1:833106/1‑71 (MQ=255)
                                                                ctcaGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGACgg  >  1:3484374/1‑71 (MQ=255)
                                                                   |                                                                   
TAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGACGG  >  W3110S.gb/3824343‑3824477

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: