Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,824,368 G→A 100% pseudogene (501/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,3680GA100.0% 119.3 / NA 40pseudogene (501/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (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.

GGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACG  >  W3110S.gb/3824302‑3824432
                                                                  |                                                                
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)
          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:138151/1‑71 (MQ=255)
                aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                              >  1:47592/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:222107/1‑42 (MQ=19)
                                                    tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                       >  1:1550379/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:1994275/1‑70 (MQ=255)
                                                             aTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg  >  1:1488061/1‑70 (MQ=255)
                                                              ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                  <  1:2925350/53‑1 (MQ=255)
                                                                  |                                                                
GGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACG  >  W3110S.gb/3824302‑3824432

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: