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

GCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGACGG  >  AP009048/3824346‑3824475
                                                              |                                                                   
gCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                                                             <  1:4479737/71‑1 (MQ=255)
gCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                               >  1:2370119/1‑69 (MQ=255)
gCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                               >  1:1520678/1‑69 (MQ=255)
 cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCagag                                                            <  1:3135069/71‑1 (MQ=255)
 cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCagag                                                            <  1:1388429/71‑1 (MQ=255)
   gTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAc                                                          <  1:1488457/71‑1 (MQ=255)
    tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                                                 <  1:5304609/63‑1 (MQ=255)
    tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc                                                         >  1:1019996/1‑71 (MQ=255)
      tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                               >  1:4177099/1‑63 (MQ=255)
      tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                                                       >  1:3415845/1‑71 (MQ=255)
         gTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACAt                                                    >  1:4176042/1‑71 (MQ=255)
          tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                               <  1:4175433/59‑1 (MQ=255)
          tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                               <  1:5133620/59‑1 (MQ=255)
          tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCTACa                                                     <  1:3334613/69‑1 (MQ=255)
           cGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg                                                   >  1:867262/1‑70 (MQ=255)
            ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGaa                                                 >  1:954892/4‑71 (MQ=255)
            ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGa                                                  <  1:3322657/67‑1 (MQ=255)
                ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                               <  1:1757383/69‑1 (MQ=255)
                ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGt                                              <  1:2453736/70‑1 (MQ=255)
                 tgtaAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGa                                                  >  1:3591983/1‑65 (MQ=255)
                 tgtaAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTa                                             <  1:4759384/70‑1 (MQ=255)
                    aaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAc                                                          >  1:3706357/2‑54 (MQ=255)
                     aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAc                                                          >  1:4059965/1‑53 (MQ=255)
                     aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                        >  1:170804/1‑71 (MQ=255)
                     aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                        >  1:3176559/1‑71 (MQ=255)
                      aCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtg                                          >  1:4758057/1‑68 (MQ=255)
                       cGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgatga                                      >  1:3872962/1‑71 (MQ=255)
                        ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTa                                             >  1:472274/4‑63 (MQ=255)
                        ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGAc                                     >  1:2319739/4‑71 (MQ=255)
                          tgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGaa                                                 >  1:912953/2‑57 (MQ=255)
                          tgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                           >  1:1121125/2‑63 (MQ=255)
                           gAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtg                                          <  1:2866328/63‑1 (MQ=255)
                                aatttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGaa                                 <  1:2339928/65‑1 (MQ=255)
                                aatttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACg                             >  1:4084591/3‑71 (MQ=255)
                                 atttGCGATCTGGAATACGTTGAAGACTCAGCCTCAGAGACCGACATGAACGTAGTGATGACCGAAGACgg                            <  1:4898941/70‑1 (MQ=255)
                                 atttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACgg                            <  1:3018851/70‑1 (MQ=255)
                                  tttGCGATCTGGAATACGTTGAAGACTCAGCCGCagag                                                            <  1:3209902/38‑1 (MQ=255)
                                  tttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGAc                                                      <  1:2388969/44‑1 (MQ=255)
                                   ttGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                                                             >  1:1522060/1‑36 (MQ=255)
                                    tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgcg                         >  1:4815722/1‑71 (MQ=255)
                                     gCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtg                                          <  1:2100922/53‑1 (MQ=255)
                                     gCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgcg                         >  1:3474697/1‑70 (MQ=255)
                                     gCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgcg                         <  1:3515450/70‑1 (MQ=255)
                                      cGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGca                       >  1:2893385/1‑71 (MQ=255)
                                      cGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGca                       <  1:1765705/71‑1 (MQ=255)
                                       gATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGt                                              <  1:372081/47‑1 (MQ=255)
                                        aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcat                      >  1:4124966/1‑70 (MQ=255)
                                         tCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTa                                             <  1:1967742/46‑1 (MQ=255)
                                         tCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatca                    >  1:3612988/1‑71 (MQ=255)
                                         tCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatca                    >  1:2691962/1‑71 (MQ=255)
                                          cTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAg                                            >  1:2967174/1‑46 (MQ=255)
                                            ggAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTa                                             >  1:722222/1‑43 (MQ=255)
                                            ggAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                        <  1:2661227/48‑1 (MQ=255)
                                             gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg                 <  1:2967960/70‑1 (MQ=255)
                                             gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg                 <  1:2524383/70‑1 (MQ=255)
                                              aaTACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGaa               <  1:5094756/71‑1 (MQ=255)
                                               aTACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGaa                                 >  1:670337/1‑52 (MQ=255)
                                               aTACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACgg                            >  1:5219455/1‑57 (MQ=255)
                                               aTACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGaa               <  1:1253419/70‑1 (MQ=255)
                                                  cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg            >  1:4715357/1‑70 (MQ=255)
                                                    ttGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGaa               >  1:5132312/1‑65 (MQ=255)
                                                    ttGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAg         <  1:978679/71‑1 (MQ=255)
                                                      gAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatca                    >  1:140825/1‑58 (MQ=255)
                                                        aGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatcat                   >  1:411754/1‑57 (MQ=255)
                                                        aGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGca                       >  1:508178/1‑53 (MQ=255)
                                                        aGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGa     >  1:3059771/1‑71 (MQ=255)
                                                         gACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAg         >  1:916932/1‑66 (MQ=255)
                                                         gACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGAc    >  1:1008338/1‑71 (MQ=255)
                                                           ctcaGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGACgg  >  1:2119850/1‑71 (MQ=255)
                                                           ctcaGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGACg   >  1:632003/1‑70 (MQ=255)
                                                              |                                                                   
GCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGACGG  >  AP009048/3824346‑3824475

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