Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,824,512 T→C 100% pseudogene (647/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,5120TC100.0% 263.9 / NA 76pseudogene (647/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/43);  total (33/43)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GACGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAA  >  AP009048/3824445‑3824576
                                                                   |                                                                   
gACGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCtt                                                                  <  1:332203/71‑1 (MQ=255)
  cGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttg                                                                 <  1:2257791/70‑1 (MQ=255)
   gggCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttgtt                                                               <  1:4953375/71‑1 (MQ=255)
     gcgcATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTgg                                                             <  1:2872687/71‑1 (MQ=255)
        catcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGcgc                                                          >  1:4893853/1‑69 (MQ=255)
        catcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGcgc                                                          <  1:14027/71‑3 (MQ=255)
        catcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGcgc                                                          >  1:3574104/1‑69 (MQ=255)
         atcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGcgct                                                         >  1:1649746/1‑71 (MQ=255)
         atcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGcgct                                                         >  1:1766467/1‑71 (MQ=255)
          tcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTg                                                        >  1:1686129/1‑71 (MQ=255)
          tcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTg                                                        >  1:3638616/1‑71 (MQ=255)
           catTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTgg                                                       >  1:2127348/1‑71 (MQ=255)
             ttGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGcc                                                     >  1:808876/1‑71 (MQ=255)
               gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                    >  1:1543959/1‑70 (MQ=255)
               gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                    <  1:4011340/70‑1 (MQ=255)
                 aGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                <  1:3091346/71‑1 (MQ=255)
                 aGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                <  1:3998775/71‑1 (MQ=255)
                 aGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                <  1:1450519/71‑1 (MQ=255)
                  gTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgg                                              <  1:1215560/71‑1 (MQ=255)
                  gTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgg                                              <  1:2937648/71‑1 (MQ=255)
                   tGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTg                                                        <  1:4100916/62‑1 (MQ=255)
                   tGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggg                                             <  1:3200147/71‑1 (MQ=255)
                      aGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggga                                          <  1:1625377/71‑3 (MQ=255)
                      aGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgg                                              >  1:401086/1‑67 (MQ=255)
                       ggggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggggaa                                         >  1:3472279/1‑68 (MQ=255)
                        gggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggg                                            <  1:388939/67‑1 (MQ=255)
                         ggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGTGGgaatc                                       <  1:2204984/71‑1 (MQ=255)
                          gACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                      <  1:2059803/71‑1 (MQ=255)
                          gACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                      <  1:4456662/71‑1 (MQ=255)
                           aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttgt                                                                <  1:4272010/46‑1 (MQ=255)
                             ggCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                <  1:2738779/59‑1 (MQ=255)
                              gCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                  >  1:4274375/1‑71 (MQ=255)
                               cAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                  >  1:3612343/1‑70 (MQ=255)
                                 gAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                   <  1:2156009/53‑1 (MQ=255)
                                  aaGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAtt                              <  1:3388543/71‑1 (MQ=255)
                                  aaGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAtt                              >  1:2757799/1‑71 (MQ=255)
                                   aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                             <  1:2909006/71‑1 (MQ=255)
                                   aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggga                                          <  1:4119877/58‑3 (MQ=255)
                                   aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggga                                          <  1:633878/58‑3 (MQ=255)
                                       gAGCCGTTCACCCATGAAGAGCTACTCACCttgtt                                                               <  1:1767320/35‑1 (MQ=255)
                                       gAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTg                                                        >  1:2572589/1‑42 (MQ=38)
                                       gAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc                         >  1:1702098/1‑71 (MQ=255)
                                        aGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCg                        <  1:1448893/71‑1 (MQ=255)
                                         gCCGTTCACGCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCg                        >  1:1712121/1‑70 (MQ=255)
                                         gCCGTTCACCCATGAAGAGCTACTCACCTTGTTgg                                                             >  1:2977950/1‑35 (MQ=255)
                                         gCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggg                                             >  1:726700/1‑49 (MQ=255)
                                         gCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggg                                             >  1:1947211/1‑49 (MQ=255)
                                             ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                   <  1:349125/41‑1 (MQ=38)
                                             ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                   <  1:814858/41‑1 (MQ=38)
                                             ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                   <  1:2498737/41‑1 (MQ=38)
                                              tCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACg                     >  1:2140511/1‑68 (MQ=255)
                                               cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAatccaatcca                                >  1:1904658/1‑56 (MQ=255)
                                               cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCa                                >  1:3302379/1‑56 (MQ=255)
                                               cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                  <  1:856093/70‑1 (MQ=255)
                                               cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                  <  1:5323230/70‑1 (MQ=255)
                                               cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                  <  1:5293998/70‑1 (MQ=255)
                                               cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGAGACGCAg                  <  1:2262496/70‑1 (MQ=255)
                                                aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                  <  1:5300224/69‑1 (MQ=255)
                                                aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa                <  1:1001932/71‑1 (MQ=255)
                                                aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGa                 <  1:950119/70‑1 (MQ=255)
                                                  ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAgg              <  1:4097046/71‑1 (MQ=255)
                                                  ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAgg              >  1:2778095/1‑71 (MQ=255)
                                                   cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                       >  1:2043063/1‑61 (MQ=255)
                                                    aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCa                                >  1:132913/1‑51 (MQ=255)
                                                    aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACg                     <  1:85282/62‑1 (MQ=255)
                                                    aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcg            >  1:3040977/1‑71 (MQ=255)
                                                    aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcg            <  1:2774330/71‑1 (MQ=255)
                                                    aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACACAGAAggcg            >  1:3598421/1‑71 (MQ=255)
                                                      gAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACg                     <  1:2852550/60‑1 (MQ=255)
                                                       aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa                >  1:1324126/1‑64 (MQ=255)
                                                       aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCAACGCAGAAGgcggcg         >  1:2717052/1‑71 (MQ=255)
                                                         gagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGAAATCGAATCCATTGTAGc                         >  1:2390187/1‑53 (MQ=255)
                                                            cTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTa                           >  1:4981970/1‑48 (MQ=255)
                                                            cTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGc                    <  1:2534944/55‑1 (MQ=255)
                                                              aCTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa                >  1:4265161/1‑57 (MQ=255)
                                                               cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaa  <  1:4369559/70‑1 (MQ=255)
                                                                   |                                                                   
GACGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAA  >  AP009048/3824445‑3824576

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: