Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,824,533:1 +G 100% pseudogene (666/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,5331.G100.0% 182.0 / NA 54pseudogene (666/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (21/33);  total (21/33)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

AAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACT  >  W3110S.gb/3824463‑3824597
                                                                       |                                                                
aaGTGCAGGGGACGGCATAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGa                                                                   <  1:3869615/71‑1 (MQ=255)
       ggggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggggaa                                                           >  1:1219573/1‑68 (MQ=255)
          gACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                                         <  1:2702811/71‑1 (MQ=255)
            cGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaa                                                       >  1:2179771/1‑71 (MQ=255)
               cAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATcc                                                    <  1:336040/71‑1 (MQ=255)
                 gAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                                                  >  1:2681350/1‑71 (MQ=255)
                 gAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                                                  >  1:1947035/1‑71 (MQ=255)
                 gAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                                                  >  1:1709149/1‑71 (MQ=255)
                 gAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                                                  <  1:39532/71‑1 (MQ=255)
                 gAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                                                  >  1:1498362/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:4259546/1‑71 (MQ=255)
                   aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                                >  1:810483/1‑71 (MQ=255)
                   aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                                <  1:5035234/71‑1 (MQ=255)
                     gCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                                 <  1:852341/51‑1 (MQ=255)
                       gAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc                                            <  1:2957866/71‑1 (MQ=255)
                       gAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc                                            <  1:1341509/71‑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)
                             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:2884469/51‑1 (MQ=255)
                                  ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                                     <  1:4058662/51‑1 (MQ=255)
                                  ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAgg                                 <  1:3437704/71‑1 (MQ=255)
                                   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:4921390/40‑1 (MQ=25)
                                             tACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                                     <  1:1756081/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:4558368/42‑1 (MQ=37)
                                                   ccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc                                            <  1:2154355/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)
                                                             gctGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGc      <  1:3922488/70‑1 (MQ=255)
                                                             gctGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGc      <  1:222305/70‑1 (MQ=255)
                                                               tGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGa    >  1:2265170/1‑71 (MQ=255)
                                                                 gCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt  <  1:2789960/71‑1 (MQ=255)
                                                                       |                                                                
AAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACT  >  W3110S.gb/3824463‑3824597

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: