Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,824,533:1 +G 100% pseudogene (666/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,5331.G100.0% 55.6 / NA 16pseudogene (666/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (5/11);  total (5/11)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

TGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTT  >  W3110S.gb/3824466‑3824587
                                                                    |                                                      
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)
                    gAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAg                                   <  1:2626193/70‑1 (MQ=255)
                    gAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc                                  <  1:1552483/71‑1 (MQ=255)
                        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)
                                                           ctGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCa                            <  1:3020590/38‑1 (MQ=38)
                                                                    |                                                      
TGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTT  >  W3110S.gb/3824466‑3824587

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: