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

GCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAG  >  W3110S.gb/3824345‑3824465
                                                                 |                                                       
gCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                                     <  1:697092/70‑1 (MQ=255)
         tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                                     <  1:1802931/61‑1 (MQ=255)
                  aTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                                                 >  1:2165406/1‑56 (MQ=255)
                  aTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                                                 >  1:1700654/1‑56 (MQ=255)
                  aTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                                                 >  1:2562644/1‑56 (MQ=255)
                  aTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                                                 >  1:1655156/1‑56 (MQ=255)
                   ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGt                                  <  1:2985217/70‑1 (MQ=255)
                   ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGt                                  <  1:1213667/70‑1 (MQ=255)
                   ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGt                                  <  1:1574387/70‑1 (MQ=255)
                   ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGt                                  <  1:1541783/70‑1 (MQ=255)
                    tgtaAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGt                                  <  1:52733/69‑1 (MQ=255)
                     gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGaa                                     >  1:3269202/4‑65 (MQ=255)
                     gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAg                                >  1:858409/4‑70 (MQ=255)
                     gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                               >  1:2628082/4‑71 (MQ=255)
                     gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                               >  1:2583715/4‑71 (MQ=255)
                     gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                               >  1:690408/4‑71 (MQ=255)
                     gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                               >  1:2935466/4‑71 (MQ=255)
                     gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                               >  1:2008813/4‑71 (MQ=255)
                      taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                               <  1:1649458/68‑1 (MQ=255)
                      taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                               <  1:3177597/68‑1 (MQ=255)
                      taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                               <  1:3191354/68‑1 (MQ=255)
                      taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                               <  1:762755/68‑1 (MQ=255)
                              gAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTa                                 >  1:2786641/1‑60 (MQ=255)
                              gAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGAc                         >  1:475880/1‑68 (MQ=255)
                              gAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGAc                         >  1:1808667/1‑68 (MQ=255)
                              gAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGAc                         >  1:565637/1‑68 (MQ=255)
                              gAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGAc                         >  1:679057/1‑68 (MQ=255)
                                      ttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGAc                  >  1:3351789/1‑67 (MQ=255)
                                      ttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGAc                  >  1:2470346/1‑67 (MQ=255)
                                      ttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGAc                  >  1:2753659/1‑67 (MQ=255)
                                      ttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGAc                  >  1:1158868/1‑67 (MQ=255)
                                       tGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                                                 >  1:557093/1‑35 (MQ=255)
                                       tGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                                                 >  1:2309136/1‑35 (MQ=255)
                                       tGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                                                 >  1:561663/1‑35 (MQ=255)
                                       tGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                                                 >  1:1774427/1‑35 (MQ=255)
                                       tGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                                                 >  1:3078003/1‑35 (MQ=255)
                                       tGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                                                 >  1:1697492/1‑35 (MQ=255)
                                       tGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                                                 >  1:1664414/1‑35 (MQ=255)
                                       tGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                                                 >  1:1388840/1‑35 (MQ=255)
                                       tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACg                 >  1:2471769/1‑67 (MQ=255)
                                       tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACg                 >  1:3043613/1‑67 (MQ=255)
                                       tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACg                 >  1:2822716/1‑67 (MQ=255)
                                       tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACg                 >  1:2202799/1‑67 (MQ=255)
                                       tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACg                 >  1:1445325/1‑67 (MQ=255)
                                             cTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatcat       >  1:744863/1‑71 (MQ=255)
                                             cTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatcat       >  1:2286470/1‑71 (MQ=255)
                                             cTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatcat       >  1:2984748/1‑71 (MQ=255)
                                             cTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatcaa       >  1:1241765/1‑70 (MQ=255)
                                              tGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGaa                     <  1:1848505/56‑1 (MQ=255)
                                                  aTACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGaa                     <  1:1799448/52‑1 (MQ=255)
                                                  aTACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAg  <  1:3353531/71‑1 (MQ=255)
                                                   tACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatca        >  1:2175660/1‑64 (MQ=255)
                                                   tACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatca        >  1:886473/1‑64 (MQ=255)
                                                   tACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatca        >  1:1269200/1‑64 (MQ=255)
                                                   tACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatca        >  1:1386710/1‑64 (MQ=255)
                                                   tACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatca        >  1:1457633/1‑64 (MQ=255)
                                                   tACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatca        >  1:1521178/1‑64 (MQ=255)
                                                   tACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatca        >  1:2401703/1‑64 (MQ=255)
                                                   tACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatca        >  1:2415276/1‑64 (MQ=255)
                                                   tACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatca        >  1:249812/1‑64 (MQ=255)
                                                       ttGAAGACTCAGCCGCAGAGACCGACATGAAGGTAGTGATGACCGaa                     <  1:1094414/47‑1 (MQ=255)
                                                       ttGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                               <  1:3041961/37‑1 (MQ=255)
                                                                 |                                                       
GCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAG  >  W3110S.gb/3824345‑3824465

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: