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

GTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGG  >  W3110S.gb/3824342‑3824472
                                                                    |                                                              
gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                                                              >  1:3980273/1‑71 (MQ=255)
gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                                                              >  1:302832/1‑71 (MQ=255)
 tAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGcc                                                             >  1:264769/1‑71 (MQ=255)
     cGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                                            >  1:724100/1‑68 (MQ=255)
     cGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                                            >  1:4445335/1‑68 (MQ=255)
       cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                          >  1:3491147/1‑68 (MQ=255)
         gTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAc                                                     >  1:3924004/1‑71 (MQ=255)
          tttCTGTCTGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAca                                                    <  1:4383232/71‑2 (MQ=255)
          tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                                            <  1:1541307/63‑1 (MQ=255)
            tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                                                              >  1:973091/1‑59 (MQ=255)
            tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCagag                                                       >  1:3623071/1‑66 (MQ=255)
            tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                                                  >  1:2337851/1‑71 (MQ=255)
             cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                                                  >  1:1676582/1‑70 (MQ=255)
             cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                                                  >  1:1968961/1‑70 (MQ=255)
               gTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGAc                                                 >  1:1331496/1‑69 (MQ=255)
                tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg                                              <  1:3662472/71‑1 (MQ=255)
                 cGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg                                              >  1:3417052/1‑70 (MQ=255)
                  ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGaa                                            >  1:3695592/4‑71 (MQ=255)
                  ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGaa                                            >  1:562478/4‑71 (MQ=255)
                      ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                                            <  1:357173/51‑1 (MQ=255)
                      ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGAc                                                 >  1:675131/1‑62 (MQ=255)
                        gtaaaCGGTGAAGCAGTTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                      >  1:3019551/4‑71 (MQ=255)
                        gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                      >  1:3517319/4‑71 (MQ=255)
                         taaaCGGTGAAGCAATTTGCGATCTTGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                      <  1:2654906/68‑1 (MQ=255)
                         taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                      <  1:3041565/68‑1 (MQ=255)
                           aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                   >  1:3138434/1‑71 (MQ=255)
                           aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                   >  1:4074333/1‑71 (MQ=255)
                                tgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCg                              <  1:2761148/70‑1 (MQ=255)
                                tgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCg                              <  1:487381/70‑1 (MQ=255)
                                        tttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc                                                    <  1:2913977/41‑1 (MQ=255)
                                         ttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtga                                    >  1:2019722/1‑56 (MQ=255)
                                          tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                                                  >  1:4310624/1‑41 (MQ=255)
                                          tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgcg                    >  1:2677615/1‑71 (MQ=255)
                                          tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAAAGTAGTGATGACCGAAGACGGgc                     >  1:2638393/1‑70 (MQ=255)
                                           gCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                   <  1:3419110/55‑1 (MQ=255)
                                           gCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                   <  1:3330079/55‑1 (MQ=255)
                                           gCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGaa                            <  1:4214771/62‑1 (MQ=255)
                                            cGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgcgc                   <  1:2675511/70‑1 (MQ=255)
                                             gATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcat                 >  1:3523330/1‑71 (MQ=255)
                                             gATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcat                 >  1:3528587/1‑71 (MQ=255)
                                               tCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGAc                         <  1:1583700/61‑1 (MQ=255)
                                                   gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatc                <  1:2923011/66‑1 (MQ=255)
                                                   gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg            <  1:1534669/70‑1 (MQ=255)
                                                   gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGa           <  1:715560/71‑1 (MQ=255)
                                                   gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGa           <  1:2756431/71‑1 (MQ=255)
                                                      tACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg            <  1:3311912/67‑1 (MQ=255)
                                                        cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGaa          >  1:804026/1‑67 (MQ=255)
                                                        cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg       <  1:2599545/70‑1 (MQ=255)
                                                        cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg       >  1:1877168/1‑70 (MQ=255)
                                                        cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg       <  1:4102335/70‑1 (MQ=255)
                                                         gTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCa     <  1:4278859/71‑1 (MQ=255)
                                                         gTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCa     <  1:324978/71‑1 (MQ=255)
                                                         gTTGAAGACTCAGCCGAAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcat                 <  1:790907/59‑1 (MQ=255)
                                                           tGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAgg   >  1:35674/1‑71 (MQ=255)
                                                            gAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcat                 <  1:3291220/56‑1 (MQ=255)
                                                            gAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAggg  >  1:976323/1‑71 (MQ=255)
                                                            gAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAgg   <  1:1594658/70‑1 (MQ=255)
                                                            gAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAgg   <  1:1449355/70‑1 (MQ=255)
                                                                aCTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCg                              <  1:2729323/39‑1 (MQ=255)
                                                                aCTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg       >  1:1826692/1‑62 (MQ=255)
                                                                    |                                                              
GTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGG  >  W3110S.gb/3824342‑3824472

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: