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

GTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGACG  >  W3110S.gb/3824342‑3824476
                                                                    |                                                                  
gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                                                                  >  1:967319/1‑71 (MQ=255)
         gTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAc                                                         >  1:1306627/1‑71 (MQ=255)
            tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                                                      >  1:782049/1‑71 (MQ=255)
                tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACAt                                                   >  1:1376429/1‑70 (MQ=255)
                       tgtaAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAg                                           >  1:1626945/1‑71 (MQ=255)
                        gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACAt                                                   >  1:80629/4‑62 (MQ=255)
                        gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                          >  1:1643029/4‑71 (MQ=255)
                        gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                          >  1:1669191/4‑71 (MQ=255)
                        gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                          >  1:518804/4‑71 (MQ=255)
                                       atttGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                              <  1:2021183/35‑1 (MQ=38)
                                       atttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACgg                           <  1:1277048/70‑1 (MQ=255)
                                       atttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACgg                           <  1:2534500/70‑1 (MQ=255)
                                          tGCGATGTGGAATACGTTGAAGACTCAGCCGCAGAGAcc                                                        >  1:6212/1‑39 (MQ=38)
                                          tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgatga                                     >  1:500806/1‑58 (MQ=255)
                                          tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgcg                        >  1:1743395/1‑71 (MQ=255)
                                            cGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgcgc                       <  1:967115/70‑1 (MQ=255)
                                                  ggAATACGTTGAAGACTCAGCCGCAGAGACCGACAt                                                   <  1:681280/36‑1 (MQ=255)
                                                   gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatcat                  <  1:934713/68‑1 (MQ=255)
                                                   gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg                <  1:1450015/70‑1 (MQ=255)
                                                   gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg                <  1:585212/70‑1 (MQ=255)
                                                   gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg                >  1:283442/1‑70 (MQ=255)
                                                    aaTACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGaa              <  1:1627447/71‑1 (MQ=255)
                                                      tACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgatg                                      <  1:2315738/45‑1 (MQ=255)
                                                      tACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCg                                  <  1:285062/49‑1 (MQ=255)
                                                      tACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGaa                                <  1:2207939/51‑1 (MQ=255)
                                                      tACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCAtt                 <  1:1695032/66‑1 (MQ=255)
                                                      tACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGa               <  1:221443/68‑1 (MQ=255)
                                                        cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGa                              >  1:2360122/1‑51 (MQ=255)
                                                        cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGttc          <  1:295217/71‑3 (MQ=255)
                                                          ttGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGaa                                <  1:2611733/47‑1 (MQ=255)
                                                          ttGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGa                              >  1:2269440/1‑49 (MQ=255)
                                                          ttGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGa                              >  1:2266406/1‑49 (MQ=255)
                                                           tGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGc          <  1:1012708/68‑1 (MQ=255)
                                                            gAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAggg      >  1:2452836/1‑71 (MQ=255)
                                                            gAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAggg      >  1:2414950/1‑71 (MQ=255)
                                                            gAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAggg      >  1:2393492/1‑71 (MQ=255)
                                                            gAAGACTCAGCCGCAGAGACCGACATGAACGTACTGATGACCGAAGACGGGCGCATCATTGAAGTGCAgg       >  1:2434719/1‑70 (MQ=255)
                                                               gACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGaa                                >  1:579822/1‑42 (MQ=255)
                                                               gACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGaa                                >  1:2127389/1‑42 (MQ=255)
                                                                 ctcaGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGACg  <  1:900491/70‑1 (MQ=255)
                                                                    |                                                                  
GTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGACG  >  W3110S.gb/3824342‑3824476

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: