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

ACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGAT  >  W3110S.gb/3824217‑3824340
                                                               |                                                            
aCGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATgct                                                       <  1:5482053/71‑2 (MQ=255)
aCGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATgcg                                                       <  1:4725189/71‑1 (MQ=255)
aCGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATgcg                                                       <  1:3567545/71‑1 (MQ=255)
aCGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATgcg                                                       <  1:3172986/71‑1 (MQ=255)
aCGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATgcg                                                       <  1:2811189/71‑1 (MQ=255)
        cAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAg                                                 >  1:4653513/1‑69 (MQ=255)
        cAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAg                                                 >  1:4763796/1‑69 (MQ=255)
        cAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAg                                                 >  1:5242201/1‑69 (MQ=255)
        cAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAg                                                 >  1:2534678/1‑69 (MQ=255)
        cAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAg                                                 >  1:2126340/1‑69 (MQ=255)
        cAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAg                                                 >  1:1791246/1‑69 (MQ=255)
                     gCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaaa                                   <  1:1012473/70‑1 (MQ=255)
                     gCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaaa                                   <  1:2972653/70‑1 (MQ=255)
                     gCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaaa                                   <  1:5411697/70‑1 (MQ=255)
                     gCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaaa                                   <  1:2653793/70‑1 (MQ=255)
                     gCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaaa                                   <  1:1660662/70‑1 (MQ=255)
                     gCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaaa                                   <  1:1579000/70‑1 (MQ=255)
                     gCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaaa                                   <  1:1517591/70‑1 (MQ=255)
                     gCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaaa                                   <  1:795790/70‑1 (MQ=255)
                     gCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaaa                                   <  1:1207603/70‑1 (MQ=255)
                        cgcgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCTCTACAGAAGCTGGTGGaaaa                                   >  1:1672811/1‑67 (MQ=255)
                        cgcgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaaa                                   >  1:1743153/1‑67 (MQ=255)
                        cgcgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaaa                                   >  1:1027817/1‑67 (MQ=255)
                        cgcgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaaa                                   >  1:4578746/1‑67 (MQ=255)
                        cgcgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaa                                    >  1:1006857/1‑66 (MQ=255)
                         gcgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGc                               <  1:5557870/70‑1 (MQ=255)
                           gTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAg                            >  1:3895797/1‑71 (MQ=255)
                                gcgTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                       <  1:356158/71‑1 (MQ=255)
                                gcgTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                       <  1:811598/71‑1 (MQ=255)
                                         aCGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTg                         >  1:974389/1‑60 (MQ=255)
                                         aCGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaaa                     >  1:4107890/1‑64 (MQ=255)
                                         aCGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaaa                     >  1:4027052/1‑64 (MQ=255)
                                         aCGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaaa                     >  1:3774787/1‑64 (MQ=255)
                                         aCGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaaa                     >  1:5330761/1‑64 (MQ=255)
                                         aCGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaaa                     >  1:5600043/1‑64 (MQ=255)
                                         aCGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaaa                     >  1:1901052/1‑64 (MQ=255)
                                         aCGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc              >  1:4475720/1‑71 (MQ=255)
                                         aCGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc              >  1:1457664/1‑71 (MQ=255)
                                         aCGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc              >  1:4316506/1‑71 (MQ=255)
                                         aCGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGATGaaaa                     >  1:1390589/1‑64 (MQ=255)
                                              tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa         <  1:5390511/71‑1 (MQ=255)
                                              tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa         <  1:1322077/71‑1 (MQ=255)
                                              tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa         <  1:5544971/71‑1 (MQ=255)
                                              tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa         <  1:3633971/71‑1 (MQ=255)
                                               gCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggtgg                                       >  1:1473942/1‑40 (MQ=255)
                                               gCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa         <  1:2946350/70‑1 (MQ=255)
                                               gCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa         >  1:4423688/1‑70 (MQ=255)
                                                     gTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGAt  <  1:3780255/71‑1 (MQ=255)
                                                               |                                                            
ACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGAT  >  W3110S.gb/3824217‑3824340

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: