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% 194.4 / NA 56pseudogene (459/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/31);  total (25/31)
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:3232627/70‑1 (MQ=255)
cGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc                                                                   <  1:1661034/70‑1 (MQ=255)
    tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                              <  1:4711698/71‑1 (MQ=255)
    tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                              <  1:4598021/71‑1 (MQ=255)
         gCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg                                                          >  1:1071725/1‑70 (MQ=255)
          cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa                                                        >  1:920586/1‑71 (MQ=255)
          cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa                                                        >  1:3368245/1‑71 (MQ=255)
          cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa                                                        >  1:731573/1‑71 (MQ=255)
           gTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGAt                                                       >  1:4030951/1‑71 (MQ=255)
            tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCACCCCGATGAAAGGGATg                                                      <  1:407112/71‑1 (MQ=255)
            tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGAt                                                                <  1:4444260/61‑1 (MQ=255)
            tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATg                                                      <  1:4596870/71‑1 (MQ=255)
                gcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                              <  1:3120995/59‑1 (MQ=255)
                gcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAg                                                  <  1:3572879/71‑1 (MQ=255)
                  tGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg                                                          <  1:3261338/61‑1 (MQ=255)
                  tGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg                                                          <  1:2895742/61‑1 (MQ=255)
                       gATTCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                              <  1:4223515/52‑1 (MQ=255)
                       gATGCGATACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGAAGCCGCAGt                                           >  1:2258886/1‑71 (MQ=255)
                        aTGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg                                                          <  1:3761759/55‑1 (MQ=255)
                        aTGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGtt                                          >  1:4829269/1‑71 (MQ=255)
                          gcgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTc                                        >  1:4485927/1‑71 (MQ=255)
                           cgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCt                                       <  1:758694/71‑1 (MQ=255)
                              tACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAg                                            >  1:1019948/1‑63 (MQ=255)
                               aCAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTc                                    <  1:1051670/70‑1 (MQ=255)
                               aCAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCg                                   >  1:1332392/1‑71 (MQ=255)
                               aCAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCg                                   >  1:52434/1‑71 (MQ=255)
                                cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                                  <  1:796280/71‑1 (MQ=255)
                                cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                                  <  1:3405025/71‑1 (MQ=255)
                                     gCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtg                             >  1:2484455/1‑71 (MQ=255)
                                       tggtggAAAACGGCAAGCTGAAAACCAACCCGATg                                                               <  1:5133761/35‑1 (MQ=255)
                                       tggtggAAAACGGCAAGCTGAAAACCAACCCGATg                                                               <  1:2967513/35‑1 (MQ=255)
                                         gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                         >  1:4168232/1‑71 (MQ=255)
                                          tggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGt                                           <  1:599720/52‑1 (MQ=255)
                                               aaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                         >  1:1327000/1‑65 (MQ=255)
                                                aCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                                  >  1:4451338/1‑55 (MQ=255)
                                                 cGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAg                                                  >  1:526101/1‑38 (MQ=255)
                                                 cGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAg                 <  1:3436888/71‑1 (MQ=255)
                                                 cGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAg                 <  1:2976049/71‑1 (MQ=255)
                                                 cGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAg                 >  1:2827780/1‑71 (MQ=255)
                                                 cGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAg                 >  1:1997023/1‑71 (MQ=255)
                                                  ggCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                         >  1:2273313/1‑62 (MQ=255)
                                                  ggCAAGATGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                         >  1:5165/1‑62 (MQ=255)
                                                   gCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                                  <  1:4902488/52‑1 (MQ=255)
                                                   gCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                                  <  1:104312/52‑1 (MQ=255)
                                                   gCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                                  <  1:473267/52‑1 (MQ=255)
                                                   gCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca               <  1:4080151/71‑2 (MQ=255)
                                                     aaGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGtt                                          <  1:4753023/42‑1 (MQ=255)
                                                      aGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgt                            <  1:2675823/55‑1 (MQ=255)
                                                      aGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgt                            <  1:2435440/55‑1 (MQ=255)
                                                      aGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgt                            <  1:1251802/55‑1 (MQ=255)
                                                       gCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                    <  1:875016/62‑3 (MQ=255)
                                                        cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgtaa                          >  1:1632003/1‑53 (MQ=255)
                                                         tGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGaa                  >  1:3535515/1‑62 (MQ=255)
                                                          gAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGtt                                          <  1:3534936/37‑1 (MQ=255)
                                                          gAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGc                >  1:1105785/1‑63 (MQ=255)
                                                           aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                         >  1:4783576/1‑53 (MQ=255)
                                                            aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt      >  1:3678523/1‑71 (MQ=255)
                                                             aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATATGCGATc     >  1:3541660/1‑71 (MQ=255)
                                                                cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg  <  1:1837864/71‑1 (MQ=255)
                                                                   |                                                                   
CGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGG  >  W3110S.gb/3824259‑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: