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

GACGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAA  >  W3110S.gb/3824447‑3824564
                                                                   |                                                    
gACGGGCGAATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCtt                                                   >  1:113323/1‑71 (MQ=255)
   gggCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttttt                                                >  1:1740989/1‑71 (MQ=255)
      cgcATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGc                                             <  1:502912/71‑1 (MQ=255)
          tcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTg                                         >  1:726896/1‑71 (MQ=255)
           catTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCATGAAGAGCTACTCACCttg                                                  >  1:1827341/1‑60 (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)
                                      cGAGCCGTTCACCCATGAAGAGCTACTCACCttgtt                                                >  1:2663854/1‑36 (MQ=255)
                                      cGAGCCGTTCACCCATGAAGAGCTACTCACCttgtt                                                >  1:1173645/1‑36 (MQ=255)
                                      cGAGCCGTTCACCCATGAAGAGCTACTCACCttgtt                                                >  1:369863/1‑36 (MQ=255)
                                      cGAGCCGTTCACCCATGAAGAGCTACTCACCttgtt                                                >  1:1546481/1‑36 (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:2356098/40‑1 (MQ=38)
                                             ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                     <  1:1486067/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:1095182/1‑64 (MQ=255)
                                                       aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa  >  1:1070846/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)
                                                                   |                                                    
GACGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAA  >  W3110S.gb/3824447‑3824564

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: