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

TCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTA  >  W3110S.gb/3824457‑3824589
                                                                   |                                                                  
tcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTg                                                                 <  1:29291/71‑1 (MQ=255)
  atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                               >  1:1301098/1‑71 (MQ=255)
  atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                               >  1:2658195/1‑71 (MQ=255)
  atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                               <  1:215542/71‑1 (MQ=255)
     gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                             <  1:4919231/70‑1 (MQ=255)
      aaGTGCAGGGGACGGCATAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGa                                                           <  1:3869615/71‑1 (MQ=255)
      aaGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                            <  1:155805/70‑1 (MQ=255)
      aaGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                            <  1:4507173/70‑1 (MQ=255)
          gCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTg                                                                 <  1:1870265/61‑1 (MQ=255)
             ggggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggggaa                                                   >  1:1219573/1‑68 (MQ=255)
                gACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                                 <  1:2702811/71‑1 (MQ=255)
                  cGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaa                                               >  1:2179771/1‑71 (MQ=255)
                   ggCAGAAGGCGAGACGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                               >  1:1423527/1‑54 (MQ=255)
                     cAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATcc                                            <  1:336040/71‑1 (MQ=255)
                       gAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                                          >  1:1498362/1‑71 (MQ=255)
                       gAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                                          >  1:1947035/1‑71 (MQ=255)
                       gAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                                          <  1:39532/71‑1 (MQ=255)
                       gAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                                          >  1:2681350/1‑71 (MQ=255)
                       gAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                                          >  1:1709149/1‑71 (MQ=255)
                         aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCTAGGGGGAATCGAATCCATTg                                        <  1:3007926/71‑1 (MQ=255)
                         aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATcc                                            >  1:1318397/1‑67 (MQ=255)
                         aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                        >  1:4716825/1‑71 (MQ=255)
                         aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                        >  1:810483/1‑71 (MQ=255)
                         aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                        <  1:5035234/71‑1 (MQ=255)
                         aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                        >  1:4259546/1‑71 (MQ=255)
                           gCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                         <  1:852341/51‑1 (MQ=255)
                             gAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc                                    <  1:1341509/71‑1 (MQ=255)
                             gAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc                                    <  1:2957866/71‑1 (MQ=255)
                              aGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGTCCCg                                                            <  1:3361969/46‑1 (MQ=255)
                              aGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                            <  1:4967712/46‑1 (MQ=255)
                               gCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                                  >  1:2741922/1‑71 (MQ=255)
                               gCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                                  >  1:5182989/1‑71 (MQ=255)
                                ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                        <  1:43815/64‑1 (MQ=255)
                                ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                        <  1:367303/64‑1 (MQ=255)
                                   ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                             <  1:4033860/40‑1 (MQ=38)
                                   ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATcc                                            <  1:3240957/57‑1 (MQ=255)
                                    tCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACg                                >  1:3132652/1‑68 (MQ=255)
                                     cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                             <  1:406229/70‑1 (MQ=255)
                                        ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                             <  1:4977382/51‑1 (MQ=255)
                                        ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                             <  1:4058662/51‑1 (MQ=255)
                                        ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                             <  1:2884469/51‑1 (MQ=255)
                                        ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAgg                         <  1:3437704/71‑1 (MQ=255)
                                         cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                            <  1:230744/35‑1 (MQ=37)
                                         cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggc                        >  1:231081/1‑71 (MQ=255)
                                          aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcg                       >  1:1632755/1‑71 (MQ=255)
                                           tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc                                    <  1:2434941/57‑1 (MQ=255)
                                           tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcg                       >  1:3581657/1‑70 (MQ=255)
                                              agagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGgcgc                   <  1:2117491/71‑1 (MQ=255)
                                               gagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCt                  <  1:3176865/71‑1 (MQ=255)
                                                   tACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                             <  1:1756081/40‑1 (MQ=25)
                                                   tACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                             <  1:4921390/40‑1 (MQ=25)
                                                     cTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                             <  1:3024893/38‑1 (MQ=21)
                                                     cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCa              <  1:2116972/69‑1 (MQ=255)
                                                     cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCa              <  1:1633471/69‑1 (MQ=255)
                                                        accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTg         <  1:2510922/69‑1 (MQ=255)
                                                         ccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc                                    <  1:2154355/42‑1 (MQ=37)
                                                         ccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc                                    <  1:4558368/42‑1 (MQ=37)
                                                         ccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAgg                         <  1:1234504/53‑1 (MQ=255)
                                                          cTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAt       >  1:927948/1‑71 (MQ=255)
                                                           ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                                  >  1:2806385/1‑43 (MQ=37)
                                                           ttgttgGCGCTGGCCCGAGGGGGAATAGAATCCAt                                          >  1:711374/1‑35 (MQ=19)
                                                            tgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                                  >  1:3896057/1‑42 (MQ=37)
                                                            tgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa                           >  1:3252714/1‑49 (MQ=255)
                                                              ttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAttttt   <  1:2404871/71‑1 (MQ=255)
                                                              ttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAttttt   <  1:944066/71‑1 (MQ=255)
                                                                ggCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTa  <  1:4514327/70‑1 (MQ=255)
                                                                   |                                                                  
TCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTA  >  W3110S.gb/3824457‑3824589

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: