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

CGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAG  >  W3110S.gb/3824308‑3824419
                                                      |                                                         
cGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAg                                           >  1:1567859/1‑71 (MQ=255)
cGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAg                                           >  1:759081/1‑71 (MQ=255)
   cAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaa                                        <  1:2081560/71‑3 (MQ=255)
           aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                <  1:2394451/71‑1 (MQ=255)
               cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                            <  1:739129/71‑1 (MQ=255)
                      aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                         <  1:1374735/51‑2 (MQ=255)
                      aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                         <  1:887256/51‑2 (MQ=255)
                      aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                         <  1:424417/51‑2 (MQ=255)
                         aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                         <  1:2493548/48‑2 (MQ=255)
                         aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                   <  1:98981/70‑1 (MQ=255)
                         aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                   <  1:2329738/70‑1 (MQ=255)
                           aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATATGCGATCTGGAATACGTTGAAg                <  1:1683518/71‑1 (MQ=255)
                            gggATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                                                   <  1:165997/35‑1 (MQ=37)
                            gggATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                                                   <  1:630153/35‑1 (MQ=37)
                            gggATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                                                   <  1:622924/35‑1 (MQ=37)
                                tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                   <  1:2628388/63‑1 (MQ=255)
                                tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                   <  1:1885686/63‑1 (MQ=255)
                                tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                   <  1:1758700/63‑1 (MQ=255)
                                tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActca           <  1:2489363/71‑2 (MQ=255)
                                 ggTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActca           <  1:887380/70‑2 (MQ=255)
                                  gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                 <  1:1245086/63‑1 (MQ=255)
                                  gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                 <  1:1756343/63‑1 (MQ=255)
                                  gTAGCCCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                 <  1:1088448/62‑1 (MQ=255)
                                   tAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACg                      >  1:1718350/1‑57 (MQ=255)
                                   tAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACg                      >  1:1611925/1‑57 (MQ=255)
                                   tAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACg                      >  1:1329422/1‑57 (MQ=255)
                                     gCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg       <  1:2147498/70‑1 (MQ=255)
                                            tttCTGTCTGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCag    <  1:2460481/66‑1 (MQ=255)
                                              tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCagag  >  1:1864132/1‑66 (MQ=255)
                                                      |                                                         
CGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAG  >  W3110S.gb/3824308‑3824419

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: