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% 214.2 / NA 59pseudogene (657/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (24/35);  total (24/35)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTA  >  W3110S.gb/3824459‑3824589
                                                                 |                                                                   
atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCCCCCATGAAGAGCTACTCCCCTTGTTGGCGCTGGc                                                                >  1:1465924/1‑71 (MQ=255)
atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                                >  1:1709945/1‑71 (MQ=255)
 ttGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGcc                                                               <  1:3131606/71‑1 (MQ=255)
  tGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGcc                                                               <  1:2354812/70‑1 (MQ=255)
   gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                              >  1:395542/1‑70 (MQ=255)
   gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                              >  1:4123170/1‑70 (MQ=255)
   gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                              >  1:4417783/1‑70 (MQ=255)
   gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                              >  1:567389/1‑70 (MQ=255)
   gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                              <  1:806098/70‑1 (MQ=255)
   gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                              >  1:1029565/1‑70 (MQ=255)
    aaGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                             <  1:4523435/70‑1 (MQ=255)
         cAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggg                                                       <  1:2243050/70‑1 (MQ=255)
          aGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggga                                                     <  1:325776/71‑3 (MQ=255)
              gACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                                 <  1:2437694/71‑1 (MQ=255)
                  gCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGACCCGAGGGGgaatcgaatc                                             <  1:4499115/71‑1 (MQ=255)
                   cAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                             >  1:2640708/1‑70 (MQ=255)
                    aGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCa                                           >  1:1454572/1‑71 (MQ=255)
                    aGAAGGCGAGCAGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCa                                           >  1:3344862/1‑71 (MQ=255)
                    aGAAGGCGAGCAGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCa                                           >  1:3339346/1‑71 (MQ=255)
                     gAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                              <  1:296282/52‑1 (MQ=255)
                      aaGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCTAATCCAtt                                         >  1:1623665/1‑71 (MQ=255)
                       aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                        <  1:2896200/71‑1 (MQ=255)
                            aGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCg                                   <  1:833556/71‑1 (MQ=255)
                             gCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                                  >  1:2021230/1‑71 (MQ=255)
                                 ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                              <  1:479488/40‑1 (MQ=38)
                                 ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGt                                       <  1:3903044/62‑1 (MQ=255)
                                 ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc                                    <  1:3579384/65‑1 (MQ=255)
                                 ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggga                                                     <  1:2668799/48‑3 (MQ=255)
                                  tCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACg                                >  1:1014895/1‑68 (MQ=255)
                                   cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                             <  1:1253262/70‑1 (MQ=255)
                                   cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                             >  1:2363465/1‑70 (MQ=255)
                                   cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                             <  1:3662907/70‑1 (MQ=255)
                                    aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa                           <  1:1395267/71‑1 (MQ=255)
                                       cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggc                        >  1:3968211/1‑71 (MQ=255)
                                         tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCa                                           >  1:712740/1‑50 (MQ=255)
                                         tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCa                                           >  1:4114454/1‑50 (MQ=255)
                                         tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                        <  1:3738464/53‑1 (MQ=255)
                                         tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc                                    <  1:382667/57‑1 (MQ=255)
                                         tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggc                        <  1:1861839/69‑1 (MQ=255)
                                          gAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaat                                              <  1:2683285/46‑1 (MQ=255)
                                           aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCATAAGgcggcg                    >  1:3084796/1‑71 (MQ=255)
                                           aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcggc                     <  1:1469392/70‑1 (MQ=255)
                                           aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcggc                     <  1:3730553/70‑1 (MQ=255)
                                           aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcg                       <  1:4199753/68‑1 (MQ=255)
                                            agagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGgcgc                   <  1:2303595/71‑1 (MQ=255)
                                            agagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGgcgc                   <  1:1183897/71‑1 (MQ=255)
                                             gagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCt                  <  1:1566880/71‑1 (MQ=255)
                                                  aCTCACCTTGTTGGCGCTGGCCCGAGGGGTAATCGAATCCATTGTAGCGACg                                >  1:3708060/1‑52 (MQ=255)
                                                   cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCa              <  1:1455440/69‑1 (MQ=255)
                                                   cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCa              <  1:642195/69‑1 (MQ=255)
                                                        cTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                        <  1:2291068/38‑1 (MQ=25)
                                                          tgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAttt     >  1:2357567/1‑71 (MQ=255)
                                                          tgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAtt      >  1:1521075/1‑70 (MQ=255)
                                                          tgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAtt      >  1:2991052/1‑69 (MQ=255)
                                                           gttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGa                            >  1:2467183/1‑47 (MQ=255)
                                                            ttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAACCTGAttttt   <  1:1262096/71‑1 (MQ=255)
                                                            ttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAttttt   <  1:3787248/71‑1 (MQ=255)
                                                            ttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAttttt   <  1:3770800/71‑1 (MQ=255)
                                                              ggCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTa  <  1:118110/70‑1 (MQ=255)
                                                                 |                                                                   
ATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTA  >  W3110S.gb/3824459‑3824589

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