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

CGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGG  >  W3110S.gb/3824259‑3824393
                                                                   |                                                                   
cGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc                                                                   <  1:3556117/70‑1 (MQ=255)
    tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                              <  1:3171568/71‑1 (MQ=255)
    tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                              <  1:1327274/71‑1 (MQ=255)
    tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                              <  1:3997744/71‑1 (MQ=255)
          cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg                                                          <  1:1061818/69‑1 (MQ=255)
           gTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGAt                                                       >  1:315723/1‑71 (MQ=255)
           gTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGAt                                                       >  1:2796944/1‑71 (MQ=255)
           gTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGAt                                                       >  1:163422/1‑71 (MQ=255)
              gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt                                                    <  1:2385991/71‑1 (MQ=255)
              gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt                                                    <  1:2986312/71‑1 (MQ=255)
              gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt                                                    <  1:2400660/71‑1 (MQ=255)
              gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt                                                    <  1:1150996/71‑1 (MQ=255)
              gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt                                                    <  1:1433399/71‑1 (MQ=255)
                    gcaGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGc                                              <  1:2035453/69‑1 (MQ=255)
                    gcaGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGc                                              <  1:2779873/69‑1 (MQ=255)
                      aGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt                                                    >  1:3476561/1‑63 (MQ=255)
                      aGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt                                                    >  1:385805/1‑63 (MQ=255)
                      aGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt                                                    >  1:1587433/1‑63 (MQ=255)
                      aGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAg                                            >  1:193352/1‑71 (MQ=255)
                        aTGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                              <  1:2013033/51‑1 (MQ=255)
                         tGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGaaa                                                            >  1:1585176/1‑52 (MQ=255)
                         tGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg                                                          <  1:3475892/54‑1 (MQ=255)
                         tGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg                                                          <  1:3477082/54‑1 (MQ=255)
                         tGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg                                                          <  1:605175/54‑1 (MQ=255)
                         tGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                                         >  1:3102867/1‑71 (MQ=255)
                          gcgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGaaa                                                            >  1:3832479/1‑51 (MQ=255)
                          gcgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGaaa                                                            >  1:3959554/1‑51 (MQ=255)
                           cgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGaaa                                                            >  1:107490/1‑50 (MQ=255)
                           cgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTc                                        >  1:2842011/1‑70 (MQ=255)
                            gcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCt                                       <  1:3648395/70‑1 (MQ=255)
                               aCAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCg                                   >  1:760031/1‑71 (MQ=255)
                                cAGAAGCTGGTGGAAAACGGGAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                                  >  1:313391/1‑71 (MQ=255)
                                cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAg                                            >  1:404821/1‑61 (MQ=255)
                                cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAg                                            >  1:2327700/1‑61 (MQ=255)
                                cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAg                                            >  1:2581874/1‑61 (MQ=255)
                                cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                                  <  1:469022/71‑1 (MQ=255)
                                cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                                  <  1:2680115/71‑1 (MQ=255)
                                    aGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAtt                              >  1:3641157/1‑71 (MQ=255)
                                    aGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAtt                              >  1:2333468/1‑71 (MQ=255)
                                    aGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAt                               <  1:3849773/70‑1 (MQ=255)
                                     gCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                              <  1:3529589/38‑1 (MQ=255)
                                          tggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                        <  1:3459786/71‑1 (MQ=255)
                                             aaaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGAt                                                       <  1:4038334/37‑1 (MQ=255)
                                             aaaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGAt                                                       <  1:2609751/37‑1 (MQ=255)
                                                   gCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGc                >  1:3943012/1‑70 (MQ=255)
                                                      aGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCt                                       <  1:852665/44‑1 (MQ=255)
                                                      aGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCt                                       <  1:265199/44‑1 (MQ=255)
                                                         tGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc         >  1:610660/1‑71 (MQ=255)
                                                           aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg        <  1:350950/70‑1 (MQ=255)
                                                           aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg        <  1:3496848/70‑1 (MQ=255)
                                                              aCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg        <  1:3361217/67‑1 (MQ=255)
                                                                cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg  <  1:1021298/71‑1 (MQ=255)
                                                                cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg  <  1:2033232/71‑1 (MQ=255)
                                                                cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGGCGGTATTGTAACCGGTGAAGCAATTTGCGATCTgg  <  1:3730322/71‑1 (MQ=255)
                                                                cAACCCGATGAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg  <  1:1766045/70‑1 (MQ=255)
                                                                   |                                                                   
CGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGG  >  W3110S.gb/3824259‑3824393

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: