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

GTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGA  >  W3110S.gb/3824342‑3824474
                                                                    |                                                                
gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                                                                >  1:1692238/1‑71 (MQ=255)
gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                                                                >  1:4437080/1‑71 (MQ=255)
  aGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                                              >  1:969216/1‑71 (MQ=255)
  aGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                                              >  1:2969710/1‑71 (MQ=255)
  aGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                                              >  1:1403140/1‑71 (MQ=255)
  aGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                                              >  1:4756674/1‑71 (MQ=255)
    ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                            <  1:2343926/71‑1 (MQ=255)
      gCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                            >  1:1084617/1‑69 (MQ=255)
          tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                                                                <  1:88627/61‑1 (MQ=255)
          tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                            <  1:2808130/65‑1 (MQ=255)
           ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAgaga                                                        >  1:4257332/1‑68 (MQ=255)
           ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAgaga                                                        >  1:4251805/1‑68 (MQ=255)
           ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc                                                      >  1:4250726/1‑70 (MQ=255)
            tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTTAAGACTCAGCCGCAGAGACCGa                                                    >  1:4375180/1‑71 (MQ=255)
            tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                            >  1:61881/1‑63 (MQ=255)
            tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                            >  1:3595046/1‑63 (MQ=255)
            tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCg                                                     <  1:3101014/70‑1 (MQ=255)
            tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCg                                                     <  1:3095358/70‑1 (MQ=255)
            tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCg                                                     <  1:3410738/70‑1 (MQ=255)
            tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                                                    >  1:2294059/1‑71 (MQ=255)
            tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                                                    >  1:3160188/1‑71 (MQ=255)
              tGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACa                                                  >  1:3662788/1‑71 (MQ=255)
                tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACAt                                                 <  1:3682119/70‑1 (MQ=255)
                tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg                                                <  1:2158607/71‑1 (MQ=255)
                tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg                                                <  1:2876150/71‑1 (MQ=255)
                 cGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg                                                >  1:4333780/1‑70 (MQ=255)
                  ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGa                                               <  1:4503438/67‑1 (MQ=255)
                      ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                            >  1:2619281/1‑69 (MQ=255)
                      ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                            >  1:2963666/1‑69 (MQ=255)
                      ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                            >  1:3350290/1‑69 (MQ=255)
                        gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                        >  1:4469711/4‑71 (MQ=255)
                        gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                        >  1:1365419/4‑71 (MQ=255)
                         taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtg                                       <  1:1075244/69‑1 (MQ=255)
                         taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                        <  1:1887434/68‑1 (MQ=255)
                         taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                        <  1:679653/68‑1 (MQ=255)
                         taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                        <  1:1027049/68‑1 (MQ=255)
                           aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                     <  1:3792754/71‑1 (MQ=255)
                              ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                        <  1:4449869/62‑1 (MQ=255)
                                tgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGc                                                             <  1:1629147/41‑1 (MQ=37)
                                tgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                                                          >  1:3259276/2‑45 (MQ=37)
                                 gAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCg                                >  1:1341514/1‑70 (MQ=255)
                                     caatttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGa                            >  1:3218615/4‑70 (MQ=255)
                                     caatttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGAc                           <  1:1803846/68‑1 (MQ=255)
                                     caatttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGAc                           >  1:799258/4‑71 (MQ=255)
                                       atttGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                            >  1:3216119/2‑36 (MQ=38)
                                       atttGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                            >  1:641903/2‑36 (MQ=38)
                                       atttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAg                                         >  1:3654416/2‑55 (MQ=255)
                                       atttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACgg                         >  1:2427180/2‑71 (MQ=255)
                                       atttGCGATCTGGAATACGTTGAAGACTCAGACGCAGAGACCGACATGAACGTAGTGATGACCGAAGACgg                         >  1:4373505/2‑71 (MQ=255)
                                         ttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGAc                                  >  1:3781644/1‑60 (MQ=255)
                                         ttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGAc                                  >  1:135864/1‑60 (MQ=255)
                                         ttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGAc                                  >  1:4141442/1‑60 (MQ=255)
                                         ttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgc                       >  1:1813757/1‑71 (MQ=255)
                                              aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                        <  1:2561243/49‑1 (MQ=255)
                                              aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                        <  1:1488909/49‑1 (MQ=255)
                                               tCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                     <  1:3131371/51‑1 (MQ=255)
                                               tCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                     <  1:3597137/51‑1 (MQ=255)
                                               tCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGa                            >  1:547781/1‑60 (MQ=255)
                                                cTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatcat                >  1:80936/1‑71 (MQ=255)
                                                 tGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                        <  1:1310492/46‑1 (MQ=255)
                                                 tGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGaa                              <  1:2786080/56‑1 (MQ=255)
                                                 tGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGaa                              <  1:695860/56‑1 (MQ=255)
                                                   gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGAGCGAAGACGGGCGCATCATTg              >  1:1341068/1‑70 (MQ=255)
                                                   gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg              <  1:2719608/70‑1 (MQ=255)
                                                   gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg              <  1:4376014/70‑1 (MQ=255)
                                                   gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg              <  1:83120/70‑1 (MQ=255)
                                                   gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg              >  1:677128/1‑70 (MQ=255)
                                                   gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg              <  1:424270/70‑1 (MQ=255)
                                                    aaTACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg              >  1:1879794/1‑69 (MQ=255)
                                                    aaTACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg              >  1:4097619/1‑69 (MQ=255)
                                                     aTACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAg           >  1:4216887/1‑71 (MQ=255)
                                                      tACGTTGAAGACTCAGCCGGAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGt          >  1:4227544/1‑71 (MQ=255)
                                                      tACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGAc                                  <  1:4247831/47‑1 (MQ=255)
                                                        cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGaa            >  1:3225710/1‑67 (MQ=255)
                                                         gTTGAAGACTCAGCCGCAGAGACCGACATGAACGt                                           <  1:925655/35‑1 (MQ=255)
                                                         gTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCa       <  1:3415804/71‑1 (MQ=255)
                                                          ttGAAGACTCAGCCGCAGAGACCGACATGAACGTa                                          <  1:1139875/35‑1 (MQ=255)
                                                            gAAGACTCAGCCGCAGAGACCGACATGAACGTAGtga                                      >  1:316346/1‑37 (MQ=255)
                                                            gAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                        <  1:3839685/35‑1 (MQ=255)
                                                            gAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAggg    >  1:4734057/1‑71 (MQ=255)
                                                            gAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAggg    >  1:1231702/1‑71 (MQ=255)
                                                            gAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAggg    >  1:1575777/1‑71 (MQ=255)
                                                            gAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAggg    >  1:741090/1‑71 (MQ=255)
                                                             aaGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAggg    <  1:548352/70‑1 (MQ=255)
                                                              aGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGa  >  1:2754795/1‑71 (MQ=255)
                                                                aCTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg              <  1:3311709/57‑1 (MQ=255)
                                                                    |                                                                
GTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGA  >  W3110S.gb/3824342‑3824474

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: