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

AAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGAT  >  W3110S.gb/3824312‑3824439
                                                              |                                                                 
aaGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaat                                                           <  1:967616/71‑4 (MQ=255)
  gCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaa                                                            <  1:1847149/68‑3 (MQ=255)
      aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                                      >  1:2939689/1‑70 (MQ=255)
      aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                                     >  1:47592/1‑71 (MQ=255)
      aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                                     >  1:138151/1‑71 (MQ=255)
       aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                       >  1:1231733/1‑68 (MQ=255)
       aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                       >  1:3084135/1‑68 (MQ=255)
        aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGaa                                                                >  1:177619/1‑58 (MQ=255)
          ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg                                                 <  1:869668/71‑1 (MQ=255)
           cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaa                                                            <  1:846170/59‑3 (MQ=255)
           cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                                <  1:404031/71‑1 (MQ=255)
             acccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                              <  1:2972358/69‑1 (MQ=255)
             acccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                              <  1:1997451/69‑1 (MQ=255)
             acccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGa                                               >  1:1615969/3‑70 (MQ=255)
              cccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                                             <  1:1136491/54‑2 (MQ=255)
               ccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGTTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa                                            <  1:3000435/71‑1 (MQ=255)
                     aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                       <  1:2778503/54‑1 (MQ=255)
                     aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                       <  1:2404834/70‑1 (MQ=255)
                       aGGGATGGTAGCCGCAGATTCTGTCGGTATTGTAAACGGTGAAg                                                               >  1:2293483/1‑44 (MQ=25)
                         ggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAc                                  >  1:1137035/1‑71 (MQ=255)
                              gTAGCCGCAGTTTCTGTCTGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                             <  1:1784946/71‑1 (MQ=255)
                               tAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGcc                            >  1:117338/1‑71 (MQ=255)
                                 gCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                           >  1:2733032/1‑70 (MQ=255)
                                   cGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                       >  1:536245/1‑40 (MQ=19)
                                    gCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                         >  1:1970621/1‑69 (MQ=255)
                                     cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCagag                      <  1:142135/71‑1 (MQ=255)
                                          tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                              >  1:1550379/1‑42 (MQ=19)
                                          tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                              >  1:222107/1‑42 (MQ=19)
                                          tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                       <  1:1837339/49‑1 (MQ=25)
                                                  taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg         <  1:3103781/70‑1 (MQ=255)
                                                  taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg         >  1:585810/2‑71 (MQ=255)
                                                  taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg         >  1:1000291/2‑71 (MQ=255)
                                                   aTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg         >  1:1488061/1‑70 (MQ=255)
                                                   aTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg         >  1:1994275/1‑70 (MQ=255)
                                                    ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                         <  1:2925350/53‑1 (MQ=255)
                                                         aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAc                                  >  1:1203539/1‑39 (MQ=25)
                                                         aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat  <  1:2355211/71‑1 (MQ=255)
                                                          aCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAct                                 >  1:3122730/1‑39 (MQ=25)
                                                          aCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat  >  1:1783843/1‑70 (MQ=255)
                                                          aCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat  >  1:460502/1‑70 (MQ=255)
                                                              |                                                                 
AAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGAT  >  W3110S.gb/3824312‑3824439

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: