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.

TCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGA‑GGG‑GAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTA  >  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)
       aGTTCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑g                                                          <  1:2574685/71‑1 (MQ=255)
       aGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑g                                                          <  1:3077997/71‑1 (MQ=255)
       aGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑g                                                          <  1:3073790/71‑1 (MQ=255)
        gTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑gg                                                         >  1:3015130/1‑71 (MQ=255)
                gACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑GGGTgaatc                                                  >  1:906219/1‑70 (MQ=255)
                     cAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                              <  1:3678080/54‑1 (MQ=255)
                     cAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                             <  1:1584609/55‑1 (MQ=255)
                     cAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGG‑gaatcgaatc                                             >  1:2589909/1‑70 (MQ=255)
                     cAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGACCGAGGGG‑gaatcgaatc                                             <  1:1088597/70‑1 (MQ=255)
                        aaGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCa                                           >  1:1135389/1‑69 (MQ=255)
                          ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTg                                        >  1:2955416/1‑70 (MQ=255)
                           gCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑ggg‑gga                                                    <  1:603375/56‑3 (MQ=255)
                                ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                              <  1:2443089/43‑1 (MQ=38)
                                ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                              <  1:3303550/43‑1 (MQ=38)
                                ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGAc                                 <  1:3128187/71‑1 (MQ=255)
                                  gTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTg                                        <  1:3146124/62‑1 (MQ=255)
                                  gTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGa                                  >  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)
                                     cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑GGG‑GAATCGAATCCAt                                          >  1:1964553/1‑56 (MQ=255)
                                        ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGG‑gaatcg                                                 <  1:2815257/47‑1 (MQ=255)
                                           tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑g                                                          <  1:3672158/35‑1 (MQ=37)
                                            gAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGG‑gaatcgaat                                              >  1:704675/1‑46 (MQ=255)
                                             aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAGaa                           >  1:1070846/1‑64 (MQ=255)
                                             aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAGaa                           >  1:1095182/1‑64 (MQ=255)
                                              agagCTACTCACCTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGa                                  >  1:1263837/1‑56 (MQ=255)
                                                     cTCACCTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGc                                    <  1:131486/47‑1 (MQ=255)
                                                     cTCACCTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCa                              <  1:2816232/53‑1 (MQ=255)
                                                      tCACCTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAGa                            >  1:168848/1‑54 (MQ=255)
                                                       caccTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACt          >  1:2908440/4‑71 (MQ=255)
                                                        accTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaaa            >  1:1913107/3‑68 (MQ=255)
                                                        accTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaaa            >  1:1220052/3‑68 (MQ=255)
                                                          cTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGAc                                 <  1:3096261/45‑1 (MQ=255)
                                                           ttgttgGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAGaa                           <  1:2587878/50‑1 (MQ=255)
                                                                ggCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTa  <  1:3049618/70‑1 (MQ=255)
                                                                   |                                                                   
TCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGA‑GGG‑GAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTA  >  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: ‑