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% 123.4 / NA 37pseudogene (459/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/19);  total (18/19)
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:127195/70‑1 (MQ=255)
    tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                              <  1:1062620/71‑1 (MQ=255)
      ccTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                              >  1:873491/1‑69 (MQ=255)
      ccTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                              >  1:2697341/1‑69 (MQ=255)
       cTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                              <  1:1725850/68‑1 (MQ=255)
          cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa                                                        >  1:966303/1‑71 (MQ=255)
          cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa                                                        <  1:2605963/71‑1 (MQ=255)
              gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg                                                     >  1:1702763/1‑70 (MQ=255)
              gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg                                                     <  1:2603283/70‑1 (MQ=255)
              gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt                                                    <  1:682440/71‑1 (MQ=255)
                      aGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAg                                            <  1:1871363/71‑1 (MQ=255)
                       gATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                              <  1:848177/52‑1 (MQ=255)
                         tGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa                                                        >  1:326378/1‑56 (MQ=255)
                         tGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa                                                        >  1:1683089/1‑56 (MQ=255)
                         tGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGAt                                                       >  1:1645809/1‑57 (MQ=255)
                         tGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGAt                                                       >  1:218698/1‑57 (MQ=255)
                         tGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                                         >  1:2543200/1‑71 (MQ=255)
                          gcgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGt                                           >  1:300590/1‑68 (MQ=255)
                          gcgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTc                                        >  1:1804217/1‑71 (MQ=255)
                                cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTa                                                   <  1:2040098/54‑1 (MQ=255)
                                cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                                  <  1:1454158/71‑1 (MQ=255)
                                cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGAGGAAAGGGATGGTa                                                   <  1:1373620/54‑1 (MQ=255)
                                   aaGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAt                               <  1:1753082/71‑1 (MQ=255)
                                         gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg                                                          <  1:1739180/38‑1 (MQ=255)
                                            gAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa                                                        >  1:2060947/1‑37 (MQ=255)
                                            gAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa                                                        >  1:2361285/1‑37 (MQ=255)
                                            gAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg                      >  1:794802/1‑71 (MQ=255)
                                                 cGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                         >  1:838214/1‑63 (MQ=255)
                                                 cGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAg                 >  1:1538054/1‑71 (MQ=255)
                                                        cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTg          >  1:2184014/1‑71 (MQ=255)
                                                           aaaaCCAACCCGATGAAAGGGATGGTAGCCGCATTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa       <  1:1100532/71‑1 (MQ=255)
                                                           aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg        <  1:188839/70‑1 (MQ=255)
                                                           aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa       <  1:196637/71‑1 (MQ=255)
                                                            aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa       <  1:1160993/70‑1 (MQ=255)
                                                               ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg   >  1:2304483/1‑71 (MQ=255)
                                                                cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg  <  1:1520686/71‑1 (MQ=255)
                                                                cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg  <  1:1455672/71‑1 (MQ=255)
                                                                cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg  <  1:2576568/71‑1 (MQ=255)
                                                                   |                                                                   
CGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGG  >  W3110S.gb/3824259‑3824393

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