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.3 / NA 17pseudogene (666/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (1/16);  total (1/16)
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.

AGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGC  >  W3110S.gb/3824464‑3824593
                                                                      |                                                            
aGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                             <  1:1466172/71‑1 (MQ=255)
aGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                             <  1:1703102/71‑1 (MQ=255)
     aGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggga                                                        <  1:703582/71‑3 (MQ=255)
           cGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaa                                                   <  1:342671/71‑1 (MQ=255)
                   ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                                     <  1:67331/61‑1 (MQ=255)
                           gTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                             >  1:840105/1‑44 (MQ=38)
                           gTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                             >  1:1440321/1‑44 (MQ=38)
                           gTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                             >  1:597342/1‑44 (MQ=38)
                            ttCAGCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCa                                               <  1:1630188/58‑1 (MQ=255)
                            ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCa                                               <  1:55101/58‑1 (MQ=255)
                               aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATcc                                                <  1:1639763/54‑1 (MQ=255)
                               aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATcc                                                <  1:1725879/54‑1 (MQ=255)
                               aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATcc                                                <  1:1790946/54‑1 (MQ=255)
                                              cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGt                                           <  1:1269989/44‑1 (MQ=25)
                                              cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGt                                           <  1:219798/44‑1 (MQ=25)
                                              cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCa                  <  1:1676478/69‑1 (MQ=255)
                                                   cTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAt           >  1:1575225/1‑71 (MQ=255)
                                                       ttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAttttt       <  1:310670/71‑1 (MQ=255)
                                                            gctGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGc  <  1:1577588/70‑1 (MQ=255)
                                                            gctGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGc  <  1:921534/70‑1 (MQ=255)
                                                             ctGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGc  <  1:1574295/70‑1 (MQ=255)
                                                             ctGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGc  <  1:1463050/70‑1 (MQ=255)
                                                              tGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGc  <  1:1067719/69‑1 (MQ=255)
                                                                      |                                                            
AGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGC  >  W3110S.gb/3824464‑3824593

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: