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

GTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACTG  >  W3110S.gb/3824465‑3824598
                                                                     |                                                                 
gTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgg                                                                 <  1:2054071/71‑1 (MQ=255)
    aGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggg                                                              >  1:4317960/1‑69 (MQ=255)
        gACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGa                                                                    >  1:1686194/1‑61 (MQ=255)
        gACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                                          <  1:162854/71‑1 (MQ=255)
        gACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                                          <  1:1946085/71‑1 (MQ=255)
              aGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCa                                                    >  1:3977566/1‑71 (MQ=255)
               gAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                                                   <  1:2096616/71‑1 (MQ=255)
                aaGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAtt                                                  <  1:4050467/71‑1 (MQ=255)
                aaGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                                                   >  1:4177249/1‑70 (MQ=255)
                aaGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggg                                                                <  1:789549/56‑1 (MQ=255)
                  ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGa                                                                    >  1:687218/1‑51 (MQ=255)
                  ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                                 >  1:1423235/1‑70 (MQ=255)
                  ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                                 >  1:3306939/1‑70 (MQ=255)
                        ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                                                   <  1:1579121/62‑1 (MQ=255)
                           ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCa                                                    >  1:2352885/1‑58 (MQ=255)
                             cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                      >  1:2173365/1‑70 (MQ=255)
                             cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                      <  1:1530146/70‑1 (MQ=255)
                              aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                      >  1:2146382/1‑69 (MQ=255)
                                ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAgg                                  <  1:4428092/71‑1 (MQ=255)
                                 cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggc                                 >  1:4248508/1‑71 (MQ=255)
                                  aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaat                                                       >  1:2830920/1‑48 (MQ=255)
                                  aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTa                                               >  1:60448/1‑56 (MQ=255)
                                   tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                                                   <  1:1451569/51‑1 (MQ=255)
                                   tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCg                                            <  1:2971390/58‑1 (MQ=255)
                                     aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAg                                              >  1:3209600/1‑54 (MQ=255)
                                     aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcggc                              <  1:3558123/70‑1 (MQ=255)
                                          cTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                      <  1:750796/57‑1 (MQ=255)
                                           tACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaa                                                        <  1:4511194/38‑1 (MQ=21)
                                            aCTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                      <  1:120805/55‑1 (MQ=255)
                                            aCTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCa                       >  1:4397794/1‑70 (MQ=255)
                                             cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAgg                                  <  1:731188/58‑1 (MQ=255)
                                             cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaaa                     <  1:2920199/71‑1 (MQ=255)
                                             cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaaa                     <  1:2638683/71‑1 (MQ=255)
                                             cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCa                       <  1:1583088/69‑1 (MQ=255)
                                             cTCACCTTGTTGGCGCTGGACCGAGGGGGAATCGAATcc                                                     <  1:2329575/39‑1 (MQ=19)
                                                 ccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCa                                       >  1:4030821/2‑49 (MQ=255)
                                                  cTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAt                >  1:3830321/1‑71 (MQ=255)
                                                  cTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAt                <  1:3670890/71‑1 (MQ=255)
                                                   ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAt                >  1:440652/1‑70 (MQ=255)
                                                      ttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAttttt            <  1:2753349/71‑1 (MQ=255)
                                                      ttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAttttt            <  1:2561841/71‑1 (MQ=255)
                                                      ttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAtttt             <  1:3558198/70‑1 (MQ=255)
                                                       tgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAtttct            <  1:2567070/70‑3 (MQ=255)
                                                        ggCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAttt              >  1:1422153/1‑67 (MQ=255)
                                                                cccGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACtg  <  1:1731375/71‑1 (MQ=255)
                                                                     |                                                                 
GTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACTG  >  W3110S.gb/3824465‑3824598

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: