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

TCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTA  >  W3110S.gb/3824457‑3824589
                                                                   |                                                                   
tcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTg                                                                  >  1:726896/1‑71 (MQ=255)
    tGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                              >  1:3253443/1‑71 (MQ=255)
      aaGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                              >  1:425592/1‑69 (MQ=255)
       aGTTCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                          <  1:2574685/71‑1 (MQ=255)
       aGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                          <  1:3077997/71‑1 (MQ=255)
       aGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                          <  1:3073790/71‑1 (MQ=255)
        gTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgg                                                         >  1:3015130/1‑71 (MQ=255)
                gACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGTgaatc                                                  >  1:906219/1‑70 (MQ=255)
                     cAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                              <  1:3678080/54‑1 (MQ=255)
                     cAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                             <  1:1584609/55‑1 (MQ=255)
                     cAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                             >  1:2589909/1‑70 (MQ=255)
                     cAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGACCGAGGGGgaatcgaatc                                             <  1:1088597/70‑1 (MQ=255)
                        aaGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCa                                           >  1:1135389/1‑69 (MQ=255)
                          ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                        >  1:2955416/1‑70 (MQ=255)
                           gCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggga                                                    <  1:603375/56‑3 (MQ=255)
                                ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                              <  1:2443089/43‑1 (MQ=38)
                                ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                              <  1:3303550/43‑1 (MQ=38)
                                ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGAc                                 <  1:3128187/71‑1 (MQ=255)
                                  gTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                        <  1:3146124/62‑1 (MQ=255)
                                  gTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                                  >  1:2925693/1‑68 (MQ=255)
                                   ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                              <  1:1486067/40‑1 (MQ=38)
                                   ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                              <  1:2356098/40‑1 (MQ=38)
                                   ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                              <  1:2917101/40‑1 (MQ=38)
                                     cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGAATCGAATCCAt                                          >  1:1964553/1‑56 (MQ=255)
                                        ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                                 <  1:2815257/47‑1 (MQ=255)
                                           tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                          <  1:3672158/35‑1 (MQ=37)
                                            gAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaat                                              >  1:704675/1‑46 (MQ=255)
                                             aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa                           >  1:1070846/1‑64 (MQ=255)
                                             aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa                           >  1:1095182/1‑64 (MQ=255)
                                              agagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                                  >  1:1263837/1‑56 (MQ=255)
                                                     cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc                                    <  1:131486/47‑1 (MQ=255)
                                                     cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCa                              <  1:2816232/53‑1 (MQ=255)
                                                      tCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGa                            >  1:168848/1‑54 (MQ=255)
                                                       caccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACt          >  1:2908440/4‑71 (MQ=255)
                                                        accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaaa            >  1:1913107/3‑68 (MQ=255)
                                                        accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaaa            >  1:1220052/3‑68 (MQ=255)
                                                          cTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGAc                                 <  1:3096261/45‑1 (MQ=255)
                                                           ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa                           <  1:2587878/50‑1 (MQ=255)
                                                                ggCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTa  <  1:3049618/70‑1 (MQ=255)
                                                                   |                                                                   
TCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTA  >  W3110S.gb/3824457‑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: