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

TCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTT  >  W3110S.gb/3824457‑3824587
                                                                   |                                                                
tcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTg                                                               <  1:2915491/71‑1 (MQ=255)
tcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTg                                                               >  1:1114912/1‑71 (MQ=255)
  atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                             >  1:1116224/1‑71 (MQ=255)
  atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                             >  1:1987528/1‑71 (MQ=255)
   ttGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGcc                                                            >  1:1014091/1‑71 (MQ=255)
   ttGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGcc                                                            >  1:410209/1‑71 (MQ=255)
   ttGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                             >  1:2731189/1‑70 (MQ=255)
     gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                           <  1:257820/70‑1 (MQ=255)
      aaGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                          <  1:920468/70‑1 (MQ=255)
         tGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggg                                                     >  1:3163094/1‑71 (MQ=255)
         tGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggg                                                     >  1:1691050/1‑71 (MQ=255)
         tGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggg                                                     <  1:2412182/71‑1 (MQ=255)
              gggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggggaa                                                 <  1:1274175/70‑4 (MQ=255)
                gACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                               <  1:2998399/71‑1 (MQ=255)
                 aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                           <  1:558798/58‑1 (MQ=255)
                     cAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                          <  1:610146/55‑1 (MQ=255)
                             gAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAg                                   <  1:2626193/70‑1 (MQ=255)
                             gAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc                                  <  1:1552483/71‑1 (MQ=255)
                               gCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                           <  1:2284359/44‑1 (MQ=38)
                                 cGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGa                                                         >  1:610520/1‑44 (MQ=38)
                                  gTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                       <  1:1662968/44‑1 (MQ=38)
                                   ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCa                            <  1:1110814/71‑1 (MQ=255)
                                     cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                           >  1:2558227/1‑70 (MQ=255)
                                       cccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAg                        <  1:2240134/71‑1 (MQ=255)
                                        ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                                <  1:1002641/62‑1 (MQ=255)
                                        ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggg                                                     <  1:31796/40‑1 (MQ=38)
                                          aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcg                     >  1:2926341/1‑71 (MQ=255)
                                          aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCAGa                                                         >  1:2262146/1‑35 (MQ=25)
                                              agagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                                <  1:1343858/56‑1 (MQ=255)
                                                agCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCa                                         <  1:1930223/45‑1 (MQ=255)
                                                     cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCa            <  1:2877691/69‑1 (MQ=255)
                                                        accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACt        >  1:1248484/3‑70 (MQ=255)
                                                            tgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaa           >  1:1150065/1‑63 (MQ=255)
                                                            tgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAttt   >  1:1075948/1‑71 (MQ=255)
                                                              ttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAtttt  <  1:1578450/70‑1 (MQ=255)
                                                                   |                                                                
TCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTT  >  W3110S.gb/3824457‑3824587

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: