Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,824,378 2 bp→AA 100% pseudogene (513‑514/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,3780GA100.0% 189.8 / NA 63pseudogene (513/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/32);  total (31/32)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*AP0090483,824,3790GA98.7% 182.9 / ‑5.8 63pseudogene (514/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (31/31);  minor base C (0/1);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCG  >  AP009048/3824313‑3824442
                                                                 ||                                                               
cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTg                                                             <  1:1277418/71‑1 (MQ=255)
  gAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                            >  1:1083984/1‑70 (MQ=255)
  gAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                            >  1:2099437/1‑70 (MQ=255)
  gAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                            >  1:3925391/1‑70 (MQ=255)
   aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                                           <  1:593508/70‑1 (MQ=255)
   aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                                          >  1:2164117/1‑71 (MQ=255)
    aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                                          <  1:1523096/70‑1 (MQ=255)
     aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATc                                                        >  1:4865450/1‑71 (MQ=255)
     aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATc                                                        <  1:3854110/71‑1 (MQ=255)
       ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg                                                      >  1:2753080/1‑71 (MQ=255)
           cccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                                   >  1:3876828/2‑70 (MQ=255)
             cGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTg                                                             >  1:2870340/1‑58 (MQ=255)
                  aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                            <  1:4844688/70‑1 (MQ=255)
                    aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                            >  1:4394149/1‑68 (MQ=255)
                     gggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGa                                        >  1:1353280/1‑71 (MQ=255)
                     gggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGa                                        <  1:4282583/71‑1 (MQ=255)
                      ggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAc                                       <  1:2104718/71‑1 (MQ=255)
                      ggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAAGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAc                                       <  1:4222870/71‑1 (MQ=255)
                       gATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAct                                      <  1:10972/71‑1 (MQ=255)
                         tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                            <  1:2112450/47‑1 (MQ=255)
                           gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCACTTTGCGATCTgg                                                     <  1:4880343/52‑1 (MQ=255)
                           gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                                     <  1:444889/52‑1 (MQ=255)
                           gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                                   <  1:3159820/54‑1 (MQ=255)
                           gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                                  >  1:610740/1‑71 (MQ=255)
                             aGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                          >  1:2069975/1‑61 (MQ=255)
                              gCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                            <  1:3958622/58‑1 (MQ=255)
                                cGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                >  1:156423/1‑68 (MQ=255)
                                 gCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATAc                                                <  1:1651568/51‑1 (MQ=255)
                                  cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                              >  1:3540091/1‑68 (MQ=255)
                                   aGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                >  1:1779357/1‑65 (MQ=255)
                                   aGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                >  1:1187311/1‑65 (MQ=255)
                                   aGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                >  1:1695695/1‑65 (MQ=255)
                                    gTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCt                                                       >  1:4963750/1‑41 (MQ=19)
                                     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:4957479/4‑71 (MQ=255)
                                             ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGaa                >  1:5008424/4‑71 (MQ=255)
                                               taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCagag                           <  1:4021776/57‑1 (MQ=255)
                                               taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCagag                           <  1:2779109/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:2106733/4‑71 (MQ=255)
                                                   gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt          >  1:38320/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)
                                                                 ||                                                               
CTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCG  >  AP009048/3824313‑3824442

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