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

AGACGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGA‑G‑GGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGC  >  W3110S.gb/3824446‑3824576
                                                                    |                                                                
aGACGGGCGCATCATTGAAGTGCAGGGGACGGAAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCt                                                                <  1:204601/71‑1 (MQ=255)
  aCGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttg                                                              <  1:1159534/71‑1 (MQ=255)
  aCGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttg                                                              <  1:4335651/71‑1 (MQ=255)
  aCGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttg                                                              <  1:4223322/71‑1 (MQ=255)
   cGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttgt                                                             <  1:2739852/71‑1 (MQ=255)
      gcgcATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTgg                                                          <  1:1299143/71‑1 (MQ=255)
       cgcATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGc                                                         <  1:1493717/71‑1 (MQ=255)
             atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGcgc                                                       >  1:3301282/1‑65 (MQ=255)
             atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                   >  1:8617/1‑71 (MQ=255)
             atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                   >  1:4699056/1‑71 (MQ=255)
             atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                   >  1:364170/1‑71 (MQ=255)
             atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                   >  1:3176598/1‑71 (MQ=255)
             atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                   >  1:2977582/1‑71 (MQ=255)
             atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                   >  1:2108754/1‑71 (MQ=255)
             atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                   >  1:2343339/1‑71 (MQ=255)
             atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACACACCTTGTAGGCGCTGGc                                                   >  1:1043744/1‑71 (MQ=255)
              ttGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttg                                                              <  1:4575865/59‑1 (MQ=255)
              ttGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttg                                                              <  1:4433060/59‑1 (MQ=255)
              ttGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTg                                                     <  1:4091672/68‑1 (MQ=255)
                gAAGTGCAGGGGCCGGCAGACGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                 <  1:5037967/70‑1 (MQ=255)
                gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                 <  1:2408548/70‑1 (MQ=255)
                gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                 <  1:3816372/70‑1 (MQ=255)
                 aaGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                <  1:2738937/70‑1 (MQ=255)
                      cAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑g‑ggg                                         >  1:3721371/1‑70 (MQ=255)
                        ggggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑g‑ggggaa                                      >  1:1168623/1‑68 (MQ=255)
                          ggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑g‑ggg                                         >  1:3405667/1‑66 (MQ=255)
                            aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                <  1:4809717/59‑1 (MQ=255)
                                 aGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑g‑gggg                                        >  1:2085607/1‑59 (MQ=255)
                                  gAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCAt                             <  1:879039/71‑1 (MQ=255)
                                     ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGt                          <  1:3656726/71‑1 (MQ=255)
                                      gCGAGCCGTTCACCCATGAAGAGCTACTCACCttgt                                                             >  1:4137953/1‑36 (MQ=255)
                                          gCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCGGGCCCGA‑GTGGGAATCGAATCCATTGTAGCg                      <  1:1530798/70‑1 (MQ=255)
                                           ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCg                      <  1:857031/69‑1 (MQ=255)
                                           ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCg                      <  1:1945019/69‑1 (MQ=255)
                                           ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCg                      <  1:567631/69‑1 (MQ=255)
                                           ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGAc                    >  1:3987079/1‑71 (MQ=255)
                                             gTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTg                                                     >  1:2486123/1‑37 (MQ=38)
                                              ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTg                                                     <  1:2457734/36‑1 (MQ=37)
                                              ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCa                 <  1:5073320/71‑1 (MQ=255)
                                                cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGgaatcga                                   >  1:1425469/1‑51 (MQ=255)
                                                cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGc                  >  1:1024424/1‑68 (MQ=255)
                                                cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAg                <  1:344011/70‑1 (MQ=255)
                                                cACCAATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGc                  >  1:3730094/1‑68 (MQ=255)
                                                     aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAggcg          >  1:2811059/1‑71 (MQ=255)
                                                      tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑g‑gggga                                       <  1:4557815/40‑3 (MQ=37)
                                                        aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCTAATCCATTGTAGCGACGCAGAAGgcggcg       <  1:4102588/71‑1 (MQ=255)
                                                        aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAGgcggcg       >  1:2437565/1‑71 (MQ=255)
                                                          gagCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCt     <  1:1234289/71‑1 (MQ=255)
                                                           agCTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTg    >  1:1049439/1‑71 (MQ=255)
                                                             cTACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAggcgg         >  1:3666631/1‑64 (MQ=255)
                                                              tACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAggc           <  1:3083935/61‑1 (MQ=255)
                                                              tACTCACCTTGTTGGCGCTGGCCCGAGG‑GGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGc  <  1:377164/70‑1 (MQ=255)
                                                                 tCACCTTGTTGGCGCTGGCCCGAGG‑GGgaatcgaat                                 >  1:2029082/1‑36 (MQ=19)
                                                                    |                                                                
AGACGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGA‑G‑GGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGC  >  W3110S.gb/3824446‑3824576

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: ‑