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

CTGCGACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGC  >  AP009048/3824210‑3824344
                                                                    |                                                                  
cTGCGACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAg                                                                  <  1:4808915/71‑1 (MQ=255)
cTGCGACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAg                                                                  <  1:4802779/71‑1 (MQ=255)
      cGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATgcgc                                                            >  1:3687788/1‑71 (MQ=255)
       gTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATgcgc                                                            >  1:1357188/1‑70 (MQ=255)
       gTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATgcgc                                                            >  1:5359556/1‑70 (MQ=255)
       gTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATgcgc                                                            >  1:343839/1‑70 (MQ=255)
       gTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCt                                                           <  1:651262/71‑1 (MQ=255)
             cAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAg                                                       >  1:3661084/1‑69 (MQ=255)
             cAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAg                                                       >  1:4955341/1‑69 (MQ=255)
              aGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAg                                                    <  1:2352072/71‑1 (MQ=255)
               ggCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGc                                                   >  1:4178408/1‑71 (MQ=255)
                gCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCt                                                  <  1:1513036/71‑1 (MQ=255)
                 cTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtg                                                 <  1:737038/71‑1 (MQ=255)
                    aTGGTGGCACGCGTACCGCGTCGATTACTGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtgg                                                <  1:992892/69‑1 (MQ=255)
                    aTGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggt                                               <  1:4094455/70‑1 (MQ=255)
                    aTGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggt                                               <  1:707778/70‑1 (MQ=255)
                     tggtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTa                                                          >  1:3297656/1‑58 (MQ=255)
                       gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtgg                                                <  1:2382886/66‑1 (MQ=255)
                       gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaa                                           <  1:1366962/71‑1 (MQ=255)
                       gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaa                                           <  1:3770437/71‑1 (MQ=255)
                        tggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaa                                          <  1:3700154/71‑1 (MQ=255)
                        tggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaa                                          <  1:3032936/71‑1 (MQ=255)
                          gCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaaa                                         <  1:4703329/70‑1 (MQ=255)
                          gCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaaa                                         <  1:1779742/70‑1 (MQ=255)
                           cACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtgg                                                >  1:84044/1‑62 (MQ=255)
                             cgcgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAAc                                        <  1:2509239/68‑1 (MQ=255)
                              gcgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAg                                                    >  1:2535033/1‑55 (MQ=255)
                              gcgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGc                                     <  1:636739/70‑1 (MQ=255)
                               cgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGc                                     >  1:1850040/1‑69 (MQ=255)
                                gTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAAAAGCTGGTGGAAAACGGCAAg                                  >  1:303945/1‑71 (MQ=255)
                                 tACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtgg                                                <  1:5258668/56‑1 (MQ=255)
                                 tACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGc                                 >  1:2692979/1‑71 (MQ=255)
                                  aCCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtgg                                                >  1:4474663/1‑55 (MQ=255)
                                  aCCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCt                                >  1:3105445/1‑71 (MQ=255)
                                  aCCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCt                                >  1:820706/1‑71 (MQ=255)
                                   ccGCGTCGATTACGTGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTg                               <  1:3784493/71‑1 (MQ=255)
                                   ccGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtgg                                                >  1:5161509/1‑54 (MQ=255)
                                   ccGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCt                                >  1:1061013/1‑70 (MQ=255)
                                   ccGCGTCGATTACGGGTGCCTGCGTGGCGCAGGCAGATGCGCTACAGAAGCTGGTGGTAAACTGCAAGCTg                               <  1:1900554/71‑1 (MQ=255)
                                    cgcgTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGa                              <  1:2958790/71‑1 (MQ=255)
                                     gcgTCGATTACGGGTGCCTGCGTGGCGTTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                             <  1:4348726/71‑1 (MQ=255)
                                     gcgTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                             <  1:3211478/71‑1 (MQ=255)
                                     gcgTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                             <  1:5316409/71‑1 (MQ=255)
                                     gcgTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                             <  1:4033246/71‑1 (MQ=255)
                                     gcgTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                             <  1:830990/71‑1 (MQ=255)
                                     gcgTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                             <  1:1466695/71‑1 (MQ=255)
                                        tCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCa                                    <  1:2413359/61‑1 (MQ=255)
                                         cGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaaa                           >  1:2577537/1‑69 (MQ=255)
                                         cGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAAc                          >  1:2701498/1‑70 (MQ=255)
                                          gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAAcc                         >  1:1205140/1‑70 (MQ=255)
                                          gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAAcc                         <  1:2998931/70‑1 (MQ=255)
                                           aTTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTa                                                          <  1:421376/36‑1 (MQ=255)
                                           aTTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAg                                                       <  1:1451815/39‑1 (MQ=255)
                                           aTTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtg                                                 <  1:3506205/45‑1 (MQ=255)
                                           aTTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGc                                     >  1:2357937/1‑57 (MQ=255)
                                           aTTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCaa                       >  1:711131/1‑71 (MQ=255)
                                            ttACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaaa                                         >  1:2815599/1‑52 (MQ=255)
                                            ttACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCaac                      >  1:4582197/1‑70 (MQ=255)
                                              aCGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc                    >  1:2935435/1‑71 (MQ=255)
                                               cGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc                    <  1:2645234/70‑1 (MQ=255)
                                               cGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCg                   <  1:585790/71‑1 (MQ=255)
                                                      cTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCt                                <  1:1730589/51‑1 (MQ=255)
                                                       tGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaa                                           <  1:908195/39‑1 (MQ=255)
                                                       tGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCa                                    <  1:4981164/46‑1 (MQ=255)
                                                        gCGTGGCGCTGGCAGATGCGCTACAGTAGCTGGTGGAAAACGGCAAGCTGaa                             <  1:237809/52‑1 (MQ=39)
                                                        gCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                             >  1:2868051/1‑52 (MQ=255)
                                                         cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg           <  1:1415071/69‑1 (MQ=255)
                                                         cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg           <  1:2542651/69‑1 (MQ=255)
                                                         cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAg            <  1:174440/68‑1 (MQ=255)
                                                             gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGa                              <  1:2501400/46‑1 (MQ=255)
                                                             gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg      <  1:3370479/70‑1 (MQ=255)
                                                             gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt     <  1:5098955/71‑1 (MQ=255)
                                                              cgcTGGCAGATGCGCTACATAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt     >  1:3706329/1‑70 (MQ=255)
                                                              cgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg      >  1:981643/1‑69 (MQ=255)
                                                              cgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTa    <  1:822376/71‑1 (MQ=255)
                                                               gcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                             <  1:280291/45‑1 (MQ=255)
                                                                 tGGCAGATGCTCTACAGAAGCTGGTGGAAAACGGCAAGCTGa                              <  1:4712254/42‑1 (MQ=38)
                                                                 tGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGc  >  1:4384472/1‑70 (MQ=255)
                                                                    |                                                                  
CTGCGACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGC  >  AP009048/3824210‑3824344

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: