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

GCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAG  >  W3110S.gb/3824296‑3824419
                                                                  |                                                         
gCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtg                                                       <  1:2205601/71‑1 (MQ=255)
  tggtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgt                                                      >  1:1220407/1‑70 (MQ=255)
      ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgt                                                      >  1:2420233/1‑66 (MQ=255)
      ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACg                                                 >  1:1225382/1‑71 (MQ=255)
       gAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg                                                >  1:1809538/1‑71 (MQ=255)
        aaaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggt                                               <  1:328515/71‑2 (MQ=255)
        aaaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggt                                               >  1:2010030/1‑70 (MQ=255)
          aaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                                                     <  1:723956/63‑2 (MQ=255)
                aaGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaat                                       <  1:967616/71‑4 (MQ=255)
                  gCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaa                                        <  1:1847149/68‑3 (MQ=255)
                   cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACg                                                 >  1:108782/1‑58 (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)
                             acccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAt                                                         <  1:1679583/38‑1 (MQ=38)
                             acccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                          <  1:1997451/69‑1 (MQ=255)
                             acccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                          <  1:2972358/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:222107/1‑42 (MQ=19)
                                                          tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                          >  1:1550379/1‑42 (MQ=19)
                                                          tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                   <  1:1837339/49‑1 (MQ=25)
                                                                  |                                                         
GCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAG  >  W3110S.gb/3824296‑3824419

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: