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% 177.5 / NA 53pseudogene (647/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/29);  total (24/29)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAA  >  W3110S.gb/3824451‑3824579
                                                               |                                                                  
ggCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttg                                                                <  1:704476/68‑1 (MQ=255)
ggCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttg                                                                <  1:2260616/68‑1 (MQ=255)
     atcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGcgct                                                        <  1:2464369/71‑1 (MQ=255)
     atcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGcgct                                                        >  1:3802376/1‑71 (MQ=255)
      tcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttgtt                                                              >  1:3320752/1‑64 (MQ=255)
      tcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGc                                                           <  1:538418/67‑1 (MQ=255)
      tcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTg                                                       <  1:2921874/71‑1 (MQ=255)
       catTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTgg                                                      >  1:3623908/1‑71 (MQ=255)
        atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttgt                                                               >  1:3235649/1‑61 (MQ=255)
        atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                     >  1:4256091/1‑71 (MQ=255)
        atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                     >  1:1764122/1‑71 (MQ=255)
        atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                     >  1:3927245/1‑71 (MQ=255)
              gTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgg                                              <  1:2054071/71‑1 (MQ=255)
                  aGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggg                                           >  1:4317960/1‑69 (MQ=255)
                      gACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttgtt                                                              >  1:3645350/1‑48 (MQ=255)
                      gACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGa                                                 >  1:1686194/1‑61 (MQ=255)
                      gACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                       <  1:1946085/71‑1 (MQ=255)
                      gACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                       <  1:162854/71‑1 (MQ=255)
                            aGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGcgc                                                         >  1:478905/1‑45 (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:3306939/1‑70 (MQ=255)
                                ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                              >  1:1423235/1‑70 (MQ=255)
                                   gAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGcg                                                          <  1:1500056/39‑2 (MQ=39)
                                   gAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGcg                                                          <  1:412078/39‑2 (MQ=39)
                                      ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                                <  1:1579121/62‑1 (MQ=255)
                                         ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCa                                 >  1:2352885/1‑58 (MQ=255)
                                           cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                   <  1:1530146/70‑1 (MQ=255)
                                           cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                   >  1:2173365/1‑70 (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)
                                                               |                                                                  
GGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAA  >  W3110S.gb/3824451‑3824579

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: