Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,824,408 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
*AP0090483,824,4080TA100.0% 250.4 / NA 82pseudogene (543/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (39/43);  total (39/43)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGACG  >  AP009048/3824340‑3824474
                                                                    |                                                                  
gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                                                                  >  1:610740/1‑71 (MQ=255)
     cGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                                                >  1:156423/1‑68 (MQ=255)
       cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                              >  1:3540091/1‑68 (MQ=255)
        aGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                                                >  1:1187311/1‑65 (MQ=255)
        aGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                                                >  1:1779357/1‑65 (MQ=255)
        aGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                                                >  1:1695695/1‑65 (MQ=255)
          tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc                                                        <  1:1355882/71‑1 (MQ=255)
            tCTGTCGGTATTGTAAACGGTTAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                                                      >  1:2261763/1‑71 (MQ=255)
             cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                                                      <  1:3272012/70‑1 (MQ=255)
              tGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                                                      >  1:2403161/1‑69 (MQ=255)
                tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACAt                                                   <  1:3395292/70‑1 (MQ=255)
                tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg                                                  <  1:4912068/71‑1 (MQ=255)
                tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg                                                  <  1:490448/71‑1 (MQ=255)
                 cGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg                                                  >  1:2850926/1‑70 (MQ=255)
                  ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGaa                                                >  1:5008424/4‑71 (MQ=255)
                  ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGaa                                                >  1:4957479/4‑71 (MQ=255)
                    taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCagag                                                           <  1:2779109/57‑1 (MQ=255)
                    taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCagag                                                           <  1:4021776/57‑1 (MQ=255)
                    taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                              <  1:2827096/70‑1 (MQ=255)
                     aTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                              >  1:700162/1‑70 (MQ=255)
                       tgtaAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                                                <  1:4340959/50‑1 (MQ=255)
                        gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAg                                           >  1:1104829/4‑70 (MQ=255)
                        gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                          >  1:38320/4‑71 (MQ=255)
                        gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                          >  1:2106733/4‑71 (MQ=255)
                         taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                          <  1:2454168/68‑1 (MQ=255)
                           aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGa                                                 <  1:2146643/61‑1 (MQ=255)
                           aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                       <  1:4475537/71‑1 (MQ=255)
                           aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                       >  1:1351523/1‑71 (MQ=255)
                           aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                       <  1:1614929/71‑1 (MQ=255)
                            aCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                       >  1:2913282/1‑70 (MQ=255)
                              ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGATTCAGCCGCAGAGACCGACATGAACGTAGTGATGAc                                    >  1:3562501/4‑71 (MQ=255)
                                tgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGAcc                                   >  1:3150832/2‑70 (MQ=255)
                                tgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCg                                  <  1:1648802/70‑1 (MQ=255)
                                tgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCg                                  <  1:1956093/70‑1 (MQ=255)
                                tgAAGCAATTTGCGATCTGGAATACGTTCAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCg                                  <  1:4546954/70‑1 (MQ=255)
                                 gAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCg                                                       <  1:1579488/49‑1 (MQ=38)
                                     caatttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGAc                             <  1:5001230/68‑1 (MQ=255)
                                     caatttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGAc                             <  1:4858328/68‑1 (MQ=255)
                                      aatttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg                                                  <  1:1693283/47‑1 (MQ=39)
                                      aatttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACg                            <  1:29258/69‑1 (MQ=255)
                                       atttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGa                              >  1:1467635/2‑68 (MQ=255)
                                       atttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACgg                           <  1:1993508/70‑1 (MQ=255)
                                       atttGCGATCTGGAATACGTTGAAGACTCAGACGCAGAGACCGACATGAACGTAGTGATGACCGAAGa                              >  1:66397/2‑68 (MQ=255)
                                        tttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGt                                             <  1:1725559/52‑1 (MQ=255)
                                        tttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCg                                  <  1:495362/63‑1 (MQ=255)
                                        tttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg                          <  1:2623132/71‑1 (MQ=255)
                                          tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                       >  1:2873629/1‑56 (MQ=255)
                                          tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGAc                                    >  1:4605309/1‑59 (MQ=255)
                                             gATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcat                     >  1:2890906/1‑71 (MQ=255)
                                              aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                          <  1:1558789/49‑1 (MQ=255)
                                              aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgc                         <  1:1762211/66‑1 (MQ=255)
                                              aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgc                         <  1:2805482/66‑1 (MQ=255)
                                               tCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGaa                                                >  1:4729248/1‑42 (MQ=255)
                                               tCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                              <  1:3872732/44‑1 (MQ=255)
                                               tCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGa                              >  1:1173655/1‑60 (MQ=255)
                                               tCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatca                   >  1:2759572/1‑71 (MQ=255)
                                               tCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatca                   >  1:4591060/1‑71 (MQ=255)
                                                 tGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACgg                           >  1:3308947/1‑61 (MQ=255)
                                                   taaTACGTTGAAGACTCAGCCGCAGAGACCGACATg                                                  <  1:1572428/35‑1 (MQ=39)
                                                   gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatca                   <  1:4867748/67‑1 (MQ=255)
                                                   gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg                <  1:3744537/70‑1 (MQ=255)
                                                   gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg                <  1:2806874/70‑1 (MQ=255)
                                                   gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg                >  1:2760594/1‑70 (MQ=255)
                                                   gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGa               <  1:3971981/71‑1 (MQ=255)
                                                     aTACGTTGAAGACTCAGGCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcat                     >  1:1289500/1‑63 (MQ=255)
                                                     aTACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCg                                  <  1:1591902/50‑1 (MQ=255)
                                                      tACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGt            >  1:4263991/1‑71 (MQ=255)
                                                        cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtga                                        >  1:4881456/1‑41 (MQ=255)
                                                        cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg           <  1:536904/70‑1 (MQ=255)
                                                         gTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCa         <  1:851115/71‑1 (MQ=255)
                                                         gTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCa         <  1:3948023/71‑1 (MQ=255)
                                                         gTTGAAGACTCAGCCGCAGAGACCGACATCAACGTAGTGATGACCGa                                 <  1:1918347/47‑1 (MQ=38)
                                                          ttGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCTAAGACGGGCGCATCATTGaa              >  1:600124/1‑65 (MQ=255)
                                                          ttGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACTGgcgc                       <  1:2888038/56‑1 (MQ=255)
                                                          ttGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg                <  1:1046276/63‑1 (MQ=255)
                                                          ttGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAg        <  1:5047291/71‑1 (MQ=255)
                                                           tGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGa                                 >  1:4686169/1‑45 (MQ=255)
                                                            gAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGa               >  1:542119/1‑62 (MQ=255)
                                                            gAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAggg      >  1:952205/1‑71 (MQ=255)
                                                             aaGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGa                              >  1:2517676/1‑46 (MQ=255)
                                                             aaGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACg                            >  1:3654096/1‑48 (MQ=255)
                                                               gACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGAc   <  1:1085181/71‑1 (MQ=255)
                                                                 ctcaGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGACg  <  1:4345700/70‑1 (MQ=255)
                                                                    |                                                                  
GTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGACG  >  AP009048/3824340‑3824474

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: