Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,824,410 T→A 100% pseudogene (543/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,4100TA100.0% 96.0 / NA 32pseudogene (543/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (16/16);  total (16/16)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGC  >  W3110S.gb/3824342‑3824468
                                                                    |                                                          
gTAGCCGCAGTTTCTGTCTGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                                                          <  1:1784946/71‑1 (MQ=255)
 tAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGcc                                                         >  1:117338/1‑71 (MQ=255)
   gCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                                        >  1:2733032/1‑70 (MQ=255)
      gCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                      >  1:1970621/1‑69 (MQ=255)
       cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCagag                                                   <  1:142135/71‑1 (MQ=255)
                    taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                      >  1:585810/2‑71 (MQ=255)
                    taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                      <  1:3103781/70‑1 (MQ=255)
                    taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                      >  1:1000291/2‑71 (MQ=255)
                     aTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                      >  1:1994275/1‑70 (MQ=255)
                     aTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                      >  1:1488061/1‑70 (MQ=255)
                      ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                      <  1:2925350/53‑1 (MQ=255)
                           aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                               <  1:2355211/71‑1 (MQ=255)
                            aCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                               >  1:1783843/1‑70 (MQ=255)
                            aCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                               >  1:460502/1‑70 (MQ=255)
                               gtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATAAc                            >  1:1356933/3‑70 (MQ=255)
                                tgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCg                          <  1:1102035/70‑1 (MQ=255)
                                tgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCg                          <  1:2280246/70‑1 (MQ=255)
                                       atttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGaa                        <  1:1894713/65‑1 (MQ=255)
                                          tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGaa                                        >  1:2714771/1‑47 (MQ=255)
                                          tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGaa                                        >  1:1725330/1‑47 (MQ=255)
                                          tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgcg                >  1:2685013/1‑71 (MQ=255)
                                            cGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgcgc               <  1:1106917/70‑1 (MQ=255)
                                                cTGGAATACGTTGAAGACTCAGCCGCAGAGACAGACATGAACGTAGTGATGACCGaa                        <  1:1929000/57‑1 (MQ=255)
                                                   gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg        <  1:2853955/70‑1 (MQ=255)
                                                     aTACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACTGGCGCATCATTGAAg     <  1:634648/71‑1 (MQ=255)
                                                      tACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgatga                             >  1:1979884/1‑46 (MQ=255)
                                                        cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCAtt         >  1:2346702/1‑64 (MQ=255)
                                                        cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg   >  1:931813/1‑70 (MQ=255)
                                                        cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGc  <  1:207546/71‑1 (MQ=255)
                                                         gTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                               <  1:2213934/41‑1 (MQ=255)
                                                          ttGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCg                          <  1:506859/45‑1 (MQ=255)
                                                          ttGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGaa                        <  1:638766/47‑1 (MQ=255)
                                                                    |                                                          
GTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGC  >  W3110S.gb/3824342‑3824468

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: