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

GTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGA  >  W3110S.gb/3824300‑3824424
                                                              |                                                              
gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                                                          <  1:2560864/69‑2 (MQ=255)
gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                                                          <  1:2850429/69‑2 (MQ=255)
gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                                                        >  1:1640399/1‑71 (MQ=255)
 tggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                                                        >  1:2743344/1‑70 (MQ=255)
 tggAAAACGGCAAGCCGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                                                       <  1:2246619/71‑1 (MQ=255)
  ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                                                       <  1:1785498/70‑1 (MQ=255)
   gAAAACGGCAAGCTGAAAACCCACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAt                                                              >  1:619153/1‑62 (MQ=255)
   gAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAt                                                              >  1:2619663/1‑62 (MQ=255)
        cGGCAAGCTGAAAACCAACCCGATTAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                                                          <  1:2790087/61‑2 (MQ=255)
        cGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                                                          <  1:1324876/61‑2 (MQ=255)
            aaGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACg                                                      >  1:835850/1‑61 (MQ=255)
              gCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                   <  1:713490/62‑3 (MQ=255)
                 gAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                                                        <  1:233126/54‑1 (MQ=255)
                 gAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                        >  1:2322201/1‑70 (MQ=255)
                  aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                              >  1:1680602/1‑62 (MQ=255)
                  aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                       <  1:2978456/70‑1 (MQ=255)
                  aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                       <  1:2938087/70‑1 (MQ=255)
                  aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                      >  1:952803/1‑71 (MQ=255)
                   aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                                                       >  1:1595356/1‑53 (MQ=255)
                   aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                      <  1:2198851/70‑1 (MQ=255)
                    aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATc                                    >  1:1453574/1‑71 (MQ=255)
                     aCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAg                                                >  1:1311104/1‑58 (MQ=255)
                          cccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                                                          <  1:858432/42‑2 (MQ=38)
                          cccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                                                          <  1:2773498/42‑2 (MQ=38)
                          cccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                                                          <  1:3464266/42‑2 (MQ=38)
                           ccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa                             <  1:2311671/71‑1 (MQ=255)
                              aTGAAAGGGATGGTAGCCTCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaa                                             >  1:3100767/1‑50 (MQ=255)
                              aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                 <  1:1344436/64‑1 (MQ=255)
                                gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACg                                                      >  1:1485738/1‑41 (MQ=38)
                                gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                         >  1:124476/1‑70 (MQ=255)
                                 aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                        <  1:3238901/70‑1 (MQ=255)
                                 aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                        <  1:3311843/70‑1 (MQ=255)
                                  aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                        >  1:2753446/1‑69 (MQ=255)
                                  aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                      <  1:2148535/71‑1 (MQ=255)
                                   aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAg                     <  1:359915/71‑1 (MQ=255)
                                    gggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGAGATCTg                                  >  1:80718/1‑57 (MQ=255)
                                     ggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAc                   <  1:2153997/71‑1 (MQ=255)
                                             gCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg            <  1:1637724/70‑1 (MQ=255)
                                             gCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg            >  1:3338537/1‑70 (MQ=255)
                                              ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg            <  1:3126838/69‑1 (MQ=255)
                                              ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa          <  1:268550/71‑1 (MQ=255)
                                                    tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                 <  1:324375/42‑1 (MQ=19)
                                                     ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc    >  1:3429196/1‑70 (MQ=255)
                                                      tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                               >  1:2728424/1‑42 (MQ=19)
                                                      tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa  >  1:3337848/1‑71 (MQ=255)
                                                      tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa  >  1:297211/1‑71 (MQ=255)
                                                      tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa  >  1:769137/1‑71 (MQ=255)
                                                      tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa  >  1:2055310/1‑71 (MQ=255)
                                                       cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                        <  1:3000520/48‑1 (MQ=25)
                                                       cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa          >  1:2608084/1‑62 (MQ=255)
                                                          tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGc           <  1:3194150/58‑1 (MQ=255)
                                                              |                                                              
GTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGA  >  W3110S.gb/3824300‑3824424

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: