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

GACGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAA  >  W3110S.gb/3824447‑3824579
                                                                   |                                                                   
gACGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCtt                                                                  <  1:2313456/71‑1 (MQ=255)
  cGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttgt                                                                <  1:3025357/71‑1 (MQ=255)
  cGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttgt                                                                <  1:3996992/71‑1 (MQ=255)
   gggCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttgt                                                                <  1:3647412/70‑1 (MQ=255)
      cgcATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGc                                                            <  1:1428808/71‑1 (MQ=255)
      cgcATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGc                                                            <  1:4023225/71‑1 (MQ=255)
            atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCCCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                      >  1:3974297/1‑71 (MQ=255)
            atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                      >  1:3966724/1‑71 (MQ=255)
            atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                      >  1:1946438/1‑71 (MQ=255)
            atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                      >  1:303039/1‑71 (MQ=255)
               gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                    <  1:1032396/70‑1 (MQ=255)
               gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                    <  1:3529524/70‑1 (MQ=255)
               gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                    <  1:3479879/70‑1 (MQ=255)
               gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                    <  1:3382178/70‑1 (MQ=255)
               gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                    <  1:2927978/70‑1 (MQ=255)
               gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                    <  1:2705091/70‑1 (MQ=255)
               gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                    <  1:2137022/70‑1 (MQ=255)
               gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                    <  1:1801507/70‑1 (MQ=255)
               gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                    <  1:1319736/70‑1 (MQ=255)
                 aGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                <  1:2542525/71‑1 (MQ=255)
                 aGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                <  1:4534074/71‑1 (MQ=255)
                       ggggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggggaa                                         <  1:4070384/71‑4 (MQ=255)
                       ggggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggggaa                                         <  1:3374437/71‑4 (MQ=255)
                       ggggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggggaa                                         <  1:1728560/71‑4 (MQ=255)
                       ggggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggggaa                                         <  1:765373/71‑4 (MQ=255)
                       ggggACCGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggggaa                                         <  1:2743841/71‑4 (MQ=255)
                        gggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggggaa                                         <  1:2553953/70‑4 (MQ=255)
                              gCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                   >  1:5599192/1‑71 (MQ=255)
                              gCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                   >  1:928857/1‑71 (MQ=255)
                                               cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                   <  1:460388/70‑1 (MQ=255)
                                               cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGTGGGAATCGAATCCATTGTAGCGACGCAg                   <  1:4644712/70‑1 (MQ=255)
                                                   cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggc              <  1:81577/71‑1 (MQ=255)
                                                   cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggc              <  1:193377/71‑1 (MQ=255)
                                                   cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggc              <  1:2613572/71‑1 (MQ=255)
                                                    aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCa                                 >  1:3291855/1‑51 (MQ=255)
                                                      gAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcgg            >  1:826736/1‑70 (MQ=255)
                                                       aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcggc           >  1:4065948/1‑70 (MQ=255)
                                                       aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGgcggcg          >  1:4750198/1‑71 (MQ=255)
                                                       aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGgcggcg          >  1:4756516/1‑71 (MQ=255)
                                                       aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGgcggcg          >  1:675291/1‑71 (MQ=255)
                                                       aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGgcggcg          >  1:1481673/1‑71 (MQ=255)
                                                           gCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTgg      <  1:5503657/71‑1 (MQ=255)
                                                           gCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTgg      <  1:3230419/71‑1 (MQ=255)
                                                              aCTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCa    >  1:2565575/1‑70 (MQ=255)
                                                              aCTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCa    >  1:2970863/1‑70 (MQ=255)
                                                              aCTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCa    >  1:3451124/1‑70 (MQ=255)
                                                               cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaaa  >  1:1117958/1‑71 (MQ=255)
                                                               cTCACCTTGTTGGCGCGGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaaa  <  1:2360676/71‑1 (MQ=255)
                                                                   |                                                                   
GACGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAA  >  W3110S.gb/3824447‑3824579

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: