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

ATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTT  >  W3110S.gb/3824456‑3824588
                                                                    |                                                                 
atcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGcgct                                                                 >  1:3802376/1‑71 (MQ=255)
atcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGcgct                                                                 <  1:2464369/71‑1 (MQ=255)
 tcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTg                                                                <  1:2921874/71‑1 (MQ=255)
  catTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTgg                                                               >  1:3623908/1‑71 (MQ=255)
   atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                              >  1:3927245/1‑71 (MQ=255)
   atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                              >  1:1764122/1‑71 (MQ=255)
   atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                              >  1:4256091/1‑71 (MQ=255)
         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)
                        gAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGcgct                                                                 <  1:716322/47‑1 (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)
                                       aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                            <  1:1315759/37‑1 (MQ=38)
                                       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:2638683/71‑1 (MQ=255)
                                                      cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaaa           <  1:2920199/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:2561841/71‑1 (MQ=255)
                                                               ttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAttttt  <  1:2753349/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)
                                                                    |                                                                 
ATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTT  >  W3110S.gb/3824456‑3824588

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: