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.

AGACGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGC  >  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)
                      cAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggg                                         >  1:3721371/1‑70 (MQ=255)
                        ggggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggggaa                                      >  1:1168623/1‑68 (MQ=255)
                          ggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggg                                         >  1:3405667/1‑66 (MQ=255)
                            aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                <  1:4809717/59‑1 (MQ=255)
                                 aGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggg                                        >  1:2085607/1‑59 (MQ=255)
                                  gAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                             <  1:879039/71‑1 (MQ=255)
                                     ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGt                          <  1:3656726/71‑1 (MQ=255)
                                      gCGAGCCGTTCACCCATGAAGAGCTACTCACCttgt                                                             >  1:4137953/1‑36 (MQ=255)
                                          gCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCGGGCCCGAGTGGGAATCGAATCCATTGTAGCg                      <  1:1530798/70‑1 (MQ=255)
                                           ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCg                      <  1:857031/69‑1 (MQ=255)
                                           ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCg                      <  1:1945019/69‑1 (MQ=255)
                                           ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCg                      <  1:567631/69‑1 (MQ=255)
                                           ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGAc                    >  1:3987079/1‑71 (MQ=255)
                                             gTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTg                                                     >  1:2486123/1‑37 (MQ=38)
                                              ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTg                                                     <  1:2457734/36‑1 (MQ=37)
                                              ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCa                 <  1:5073320/71‑1 (MQ=255)
                                                cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcga                                   >  1:1425469/1‑51 (MQ=255)
                                                cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGc                  >  1:1024424/1‑68 (MQ=255)
                                                cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                <  1:344011/70‑1 (MQ=255)
                                                cACCAATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGc                  >  1:3730094/1‑68 (MQ=255)
                                                     aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcg          >  1:2811059/1‑71 (MQ=255)
                                                      tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggga                                       <  1:4557815/40‑3 (MQ=37)
                                                        aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCTAATCCATTGTAGCGACGCAGAAGgcggcg       <  1:4102588/71‑1 (MQ=255)
                                                        aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGgcggcg       >  1:2437565/1‑71 (MQ=255)
                                                          gagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCt     <  1:1234289/71‑1 (MQ=255)
                                                           agCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTg    >  1:1049439/1‑71 (MQ=255)
                                                             cTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcgg         >  1:3666631/1‑64 (MQ=255)
                                                              tACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggc           <  1:3083935/61‑1 (MQ=255)
                                                              tACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGc  <  1:377164/70‑1 (MQ=255)
                                                                 tCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaat                                 >  1:2029082/1‑36 (MQ=19)
                                                                    |                                                                
AGACGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGC  >  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: