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

GACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCG  >  W3110S.gb/3824473‑3824594
                                                             |                                                             
gACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggggaa                                                        >  1:1224168/1‑65 (MQ=255)
gACGGCAGAAGGCGAACCGTTCACCCACGAAGAGTTACTCACCTTGTTGGCGCTGGCCCGAgggggaa                                                        >  1:3574819/1‑65 (MQ=255)
gACGGCAGAAGGCGAACCGTTCACCCACGAAGAGTTACTCACCTTGTTGGCGCTGGCCCGAgggggaa                                                        >  1:2494651/1‑65 (MQ=255)
      aGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                                  <  1:1441789/69‑1 (MQ=255)
      aGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                                  <  1:1799044/69‑1 (MQ=255)
      aGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                                  <  1:1820027/69‑1 (MQ=255)
      aGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                                  <  1:3359025/69‑1 (MQ=255)
         aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                                               <  1:690926/69‑1 (MQ=255)
            cGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGt                                            >  1:855324/1‑69 (MQ=255)
            cGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGt                                            >  1:1885121/1‑69 (MQ=255)
            cGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGt                                            >  1:1082857/1‑69 (MQ=255)
            cGAGCCGTTCACCCATGAAGAGCTACTCACATTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGt                                            >  1:3597858/1‑69 (MQ=255)
                        ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGa                                 <  1:2912290/68‑1 (MQ=255)
                             aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcgg                           >  1:3725584/1‑69 (MQ=255)
                                    aCTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTgg                     >  1:30672/1‑68 (MQ=255)
                                       caccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaaa                 >  1:3077547/4‑69 (MQ=255)
                                       caccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaaa                 >  1:2491559/4‑69 (MQ=255)
                                       caccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaaa                 >  1:1395943/4‑69 (MQ=255)
                                       caccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaa                  >  1:2660993/4‑68 (MQ=255)
                                       caccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCa                   <  1:3100367/64‑1 (MQ=255)
                                        accTTGTTGTCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaaa                 >  1:2771683/3‑68 (MQ=255)
                                        accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaaa                 >  1:3595513/3‑68 (MQ=255)
                                          cTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTg              <  1:2917880/69‑1 (MQ=255)
                                           ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                  <  1:819935/48‑1 (MQ=255)
                                           ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                  <  1:924944/48‑1 (MQ=255)
                                             gttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAAc                >  1:2384017/1‑64 (MQ=255)
                                                      ggCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCg  >  1:1914403/1‑69 (MQ=255)
                                                             |                                                             
GACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCG  >  W3110S.gb/3824473‑3824594

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: