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

TCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTT  >  W3110S.gb/3824457‑3824588
                                                                   |                                                                  
tcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTg                                                                 <  1:1071970/71‑1 (MQ=255)
 catTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTgg                                                                >  1:844147/1‑71 (MQ=255)
  atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                               >  1:1472200/1‑71 (MQ=255)
  atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                               >  1:433815/1‑71 (MQ=255)
  atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                               >  1:4294499/1‑71 (MQ=255)
  atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTCGc                                                               >  1:2884109/1‑71 (MQ=255)
      aaGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                            <  1:1925989/70‑1 (MQ=255)
          gCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggg                                                     >  1:2090835/1‑71 (MQ=255)
            aGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgg                                                       <  1:320288/67‑1 (MQ=255)
              gggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                             <  1:2464491/61‑1 (MQ=255)
                  cGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaa                                              <  1:3976126/71‑1 (MQ=255)
                    gCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                            <  1:1468296/71‑1 (MQ=255)
                      aGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCa                                          <  1:1019521/71‑1 (MQ=255)
                       gAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTg                                                                 <  1:2461655/48‑1 (MQ=255)
                        aaGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                         <  1:1823944/54‑1 (MQ=255)
                         aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                       <  1:2992418/71‑1 (MQ=255)
                          ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGt                                      <  1:459663/71‑1 (MQ=255)
                          ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGt                                      <  1:952348/71‑1 (MQ=255)
                               gCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                                 >  1:3968067/1‑71 (MQ=255)
                               gCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                                 >  1:3596064/1‑71 (MQ=255)
                                  gTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATcc                                           >  1:4191603/1‑58 (MQ=255)
                                  gTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACg                               >  1:741107/1‑70 (MQ=255)
                                   ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCa                             <  1:1419014/71‑1 (MQ=255)
                                   ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCa                             <  1:4434611/71‑1 (MQ=255)
                                     cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc                                   >  1:393889/1‑63 (MQ=255)
                                     cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                            <  1:2680082/70‑1 (MQ=255)
                                     cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggg                                                      >  1:15250/1‑43 (MQ=38)
                                     cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggg                                                      >  1:2206608/1‑43 (MQ=38)
                                      aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                            <  1:3327371/38‑1 (MQ=38)
                                       cccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAg                         <  1:480963/71‑1 (MQ=255)
                                        ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa                          >  1:3925294/1‑69 (MQ=255)
                                        ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAgg                        <  1:4197056/71‑1 (MQ=255)
                                         cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggc                       <  1:2837832/71‑1 (MQ=255)
                                           tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcgg                     <  1:15012/71‑1 (MQ=255)
                                            gAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcggc                    >  1:471214/1‑71 (MQ=255)
                                             aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGgcggcg                   >  1:2356664/1‑71 (MQ=255)
                                               gagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCt                 <  1:4356356/71‑1 (MQ=255)
                                                agCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaat                                             >  1:1979395/1‑42 (MQ=25)
                                                agCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaat                                             >  1:1870504/1‑42 (MQ=25)
                                                agCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                       <  1:3520350/48‑1 (MQ=255)
                                                 gCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTgg               >  1:871215/1‑71 (MQ=255)
                                                  cTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCa                             <  1:3813840/56‑1 (MQ=255)
                                                  cTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCa                             <  1:1349279/56‑1 (MQ=255)
                                                     cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaaa           <  1:3481175/71‑1 (MQ=255)
                                                     cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaaa           >  1:3997264/1‑71 (MQ=255)
                                                      tCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                                         >  1:864592/1‑40 (MQ=25)
                                                      tCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                                 >  1:4263843/1‑48 (MQ=255)
                                                       caccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACt         >  1:458988/4‑71 (MQ=255)
                                                        accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACt         >  1:4326285/3‑70 (MQ=255)
                                                        accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACt         >  1:125311/3‑70 (MQ=255)
                                                        accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCAATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACt         >  1:3976496/3‑70 (MQ=255)
                                                         ccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                       <  1:2365548/38‑1 (MQ=25)
                                                         ccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                            <  1:1558109/49‑1 (MQ=255)
                                                          cTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAt      <  1:2176766/71‑1 (MQ=255)
                                                           ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa                          >  1:1159887/1‑50 (MQ=255)
                                                           ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa                          >  1:1164930/1‑50 (MQ=255)
                                                           ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggc                       <  1:2430966/53‑1 (MQ=255)
                                                           ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTg        <  1:1179984/68‑1 (MQ=255)
                                                           ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAtt     >  1:1862383/1‑71 (MQ=255)
                                                            tgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcgg                     >  1:3640854/1‑54 (MQ=255)
                                                            tgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcgg                     >  1:1352416/1‑54 (MQ=255)
                                                             gttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGACAAACTGAttt    >  1:3947978/1‑70 (MQ=255)
                                                              ttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                            <  1:2468947/45‑1 (MQ=255)
                                                              ttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                            <  1:2833094/45‑1 (MQ=255)
                                                              ttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAttttt  <  1:4259466/71‑1 (MQ=255)
                                                              ttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAttttt  <  1:3103152/71‑1 (MQ=255)
                                                              ttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCCGCGCTGGCAAACTGAttttt  <  1:490716/71‑1 (MQ=255)
                                                              ttgGCGCTGGCCCGAGTGGGAATCGAATCCATTGTAGCGACGCAg                            <  1:1692655/45‑1 (MQ=255)
                                                                   |                                                                  
TCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTT  >  W3110S.gb/3824457‑3824588

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: