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

GCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCA  >  W3110S.gb/3824454‑3824577
                                                            |                                                                
gcATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGcg                                                        >  1:3679166/1‑70 (MQ=255)
gcATCATTGAAGTGCAGGGGACGGCAGAAGGAGAGCCGTTCACCCATGAAGAGCTAATCACCttgt                                                             >  1:2105601/1‑66 (MQ=255)
   tcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCCCCttgt                                                             <  1:476055/63‑1 (MQ=255)
     atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTCTTGGCGCTGGc                                                   >  1:1643638/1‑71 (MQ=255)
        gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttg                                                              >  1:2726636/1‑57 (MQ=255)
        gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttg                                                              >  1:1593705/1‑57 (MQ=255)
        gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                 >  1:1725287/1‑70 (MQ=255)
                           aaGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCa                               <  1:1198614/69‑1 (MQ=255)
                           aaGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCa                               <  1:3304901/69‑1 (MQ=255)
                                      ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                              <  1:2149201/59‑1 (MQ=255)
                                      ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                              <  1:3465505/59‑1 (MQ=255)
                                      ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                              <  1:358155/59‑1 (MQ=255)
                                      ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                              <  1:550770/59‑1 (MQ=255)
                                      ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                              <  1:172658/59‑1 (MQ=255)
                                       tCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                              <  1:1885036/58‑1 (MQ=255)
                                         aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa               <  1:2376607/71‑1 (MQ=255)
                                         aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa               <  1:1994305/71‑1 (MQ=255)
                                         aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa               <  1:1491893/71‑1 (MQ=255)
                                         aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa               <  1:2911268/71‑1 (MQ=255)
                                         aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa               <  1:2960663/71‑1 (MQ=255)
                                         aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa               <  1:1476005/71‑1 (MQ=255)
                                         aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa               <  1:698876/71‑1 (MQ=255)
                                               gAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaat                                  >  1:2499300/1‑46 (MQ=255)
                                               gAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaat                                  >  1:3618079/1‑46 (MQ=255)
                                                aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGAGACGCAGAAggcgg          >  1:1361894/1‑69 (MQ=255)
                                                aaGAGCTACTCACCTTCTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGgcggcg        >  1:2307142/1‑71 (MQ=255)
                                                 agagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCAATTGTAGCGACGCAGAAGGCGgcgc       >  1:1163830/1‑71 (MQ=255)
                                                    ggtACTCACCTTGTTGGCGCTGGCCCGAGGGGGAAGCGAATCCAt                              <  1:371261/43‑1 (MQ=255)
                                                        cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATcc                                <  1:2971979/39‑1 (MQ=25)
                                                        cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATcc                                <  1:2855204/39‑1 (MQ=25)
                                                        cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATcc                                <  1:2565044/39‑1 (MQ=25)
                                                        cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATcc                                <  1:1040105/39‑1 (MQ=25)
                                                        cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCTCTGGCa  <  1:2426/69‑1 (MQ=255)
                                                        cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCa  <  1:2574147/69‑1 (MQ=255)
                                                        cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCa  <  1:838604/69‑1 (MQ=255)
                                                            |                                                                
GCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCA  >  W3110S.gb/3824454‑3824577

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: