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

CGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGG  >  W3110S.gb/3824269‑3824393
                                                         |                                                                   
cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg                                                          <  1:1085217/69‑1 (MQ=255)
cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg                                                          <  1:2123490/69‑1 (MQ=255)
cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg                                                          <  1:1568277/69‑1 (MQ=255)
cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg                                                          <  1:1430768/69‑1 (MQ=255)
 gTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGAt                                                       >  1:2507836/1‑71 (MQ=255)
 gTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGAt                                                       >  1:2253892/1‑71 (MQ=255)
           caGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCg                                               >  1:1740802/2‑69 (MQ=255)
              aTGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGtt                                          <  1:3460312/71‑1 (MQ=255)
                gcgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTc                                        >  1:3331424/1‑71 (MQ=255)
                gcgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTc                                        >  1:2119532/1‑71 (MQ=255)
                     aCAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa                                                        >  1:3649345/1‑50 (MQ=255)
                     aCAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa                                                        >  1:1184751/1‑50 (MQ=255)
                      cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                                  <  1:1352030/71‑1 (MQ=255)
                      cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                                  <  1:3635472/71‑1 (MQ=255)
                      cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                                  <  1:1836671/71‑1 (MQ=255)
                      cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                                  <  1:836630/71‑1 (MQ=255)
                      cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                                  <  1:3173894/71‑1 (MQ=255)
                      cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                                  <  1:3094908/71‑1 (MQ=255)
                      cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                                  <  1:3049483/71‑1 (MQ=255)
                      cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                                  <  1:2850928/71‑1 (MQ=255)
                          aGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGaaa                                                            >  1:1559669/1‑41 (MQ=255)
                                tggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                        >  1:973517/1‑71 (MQ=255)
                                tggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                        >  1:1059108/1‑71 (MQ=255)
                                tggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                        >  1:2577962/1‑71 (MQ=255)
                                tggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                        >  1:1925007/1‑71 (MQ=255)
                                          cAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaa              <  1:1228974/71‑3 (MQ=255)
                                          cAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaa              <  1:1195975/71‑3 (MQ=255)
                                                 aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa       <  1:2498491/71‑1 (MQ=255)
                                                 aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa       <  1:2082027/71‑1 (MQ=255)
                                                     ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg   >  1:1497214/1‑71 (MQ=255)
                                                      cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                    >  1:1195597/1‑51 (MQ=255)
                                                      cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg  <  1:2539083/71‑1 (MQ=255)
                                                      cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg  <  1:2534936/71‑1 (MQ=255)
                                                      cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg  <  1:2305458/71‑1 (MQ=255)
                                                      cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg  <  1:3434670/71‑1 (MQ=255)
                                                      cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg  <  1:889236/71‑1 (MQ=255)
                                                      cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg   >  1:739639/1‑70 (MQ=255)
                                                         |                                                                   
CGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGG  >  W3110S.gb/3824269‑3824393

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: