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

CTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGA  >  W3110S.gb/3824315‑3824424
                                               |                                                              
cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTg                                         >  1:1128407/1‑71 (MQ=255)
cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTg                                         >  1:1829888/1‑71 (MQ=255)
cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTg                                         >  1:267940/1‑71 (MQ=255)
        cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                 <  1:692401/71‑1 (MQ=255)
        cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                 <  1:690142/71‑1 (MQ=255)
                 gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGt                          >  1:764295/1‑69 (MQ=255)
                  aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                        >  1:681826/1‑70 (MQ=255)
                  aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                        >  1:1830059/1‑70 (MQ=255)
                  aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                        >  1:248900/1‑70 (MQ=255)
                    aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAg                     <  1:1134052/71‑1 (MQ=255)
                     gggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                              <  1:200206/45‑2 (MQ=255)
                     gggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                              <  1:99291/45‑2 (MQ=255)
                     gggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                              <  1:708658/45‑2 (MQ=255)
                     gggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                              <  1:34742/45‑2 (MQ=255)
                     gggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                              <  1:1491618/45‑2 (MQ=255)
                      ggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAc                   <  1:61938/71‑1 (MQ=255)
                           gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc              <  1:1386288/71‑1 (MQ=255)
                                 gCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActca                >  1:1395535/1‑62 (MQ=255)
                                      ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGttaa                       >  1:1387703/1‑49 (MQ=255)
                                      ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                      >  1:1114141/1‑52 (MQ=255)
                                      ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                      >  1:1802188/1‑52 (MQ=255)
                                      ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                      >  1:46528/1‑52 (MQ=255)
                                      ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                      >  1:1485878/1‑52 (MQ=255)
                                      ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                      >  1:112359/1‑52 (MQ=255)
                                        cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa  <  1:373701/70‑1 (MQ=255)
                                               |                                                              
CTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGA  >  W3110S.gb/3824315‑3824424

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: