Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,824,324 T→C 100% pseudogene (459/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,3240TC98.3% 206.5 / ‑4.2 60pseudogene (459/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (26/33);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01
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.

ACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATC  >  AP009048/3824256‑3824388
                                                                    |                                                                
aCGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc                                                                >  1:2935435/1‑71 (MQ=255)
 cGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc                                                                <  1:2645234/70‑1 (MQ=255)
 cGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCg                                                               <  1:585790/71‑1 (MQ=255)
           cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg                                                       <  1:1415071/69‑1 (MQ=255)
           cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg                                                       <  1:2542651/69‑1 (MQ=255)
           cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAg                                                        <  1:174440/68‑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)
                   tGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGc                                              >  1:4384472/1‑70 (MQ=255)
                      caGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGcc                                             <  1:256995/67‑1 (MQ=255)
                         aTGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCTCAGtt                                       <  1:4543053/71‑1 (MQ=255)
                         aTGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGtt                                       >  1:2023255/1‑71 (MQ=255)
                          tGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                                      >  1:3721382/1‑71 (MQ=255)
                             gcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGAt                                                             <  1:651413/45‑1 (MQ=255)
                             gcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa                                                     <  1:2948026/53‑1 (MQ=255)
                             gcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGtt                                       <  1:1826880/67‑1 (MQ=255)
                              cTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGt                                  <  1:4884105/71‑1 (MQ=255)
                                aCAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCg                                >  1:477730/1‑71 (MQ=255)
                                 cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                               >  1:2431471/1‑71 (MQ=255)
                                 cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                               >  1:2369932/1‑71 (MQ=255)
                                 cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                               >  1:4910736/1‑71 (MQ=255)
                                 cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                               <  1:1248891/71‑1 (MQ=255)
                                 cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                               >  1:4058954/1‑71 (MQ=255)
                                   gAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCggta                             <  1:2569743/71‑3 (MQ=255)
                                    aaGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg                                                       <  1:933474/44‑1 (MQ=255)
                                        tggtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTa                                                <  1:596347/47‑1 (MQ=255)
                                        tggtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTa                                                <  1:3457649/47‑1 (MQ=255)
                                          gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgtaa                       <  1:1414674/70‑3 (MQ=255)
                                          gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                        <  1:2351791/69‑2 (MQ=255)
                                          gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                        <  1:760980/69‑2 (MQ=255)
                                            ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTc                                 <  1:3864843/58‑1 (MQ=255)
                                            ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACg                    >  1:5218398/1‑71 (MQ=255)
                                             gAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTg                                   >  1:2803735/1‑55 (MQ=255)
                                             gAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg                   >  1:2909608/1‑71 (MQ=255)
                                             gAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg                   <  1:4858929/71‑1 (MQ=255)
                                              aaaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                      >  1:2634112/1‑67 (MQ=255)
                                               aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggt                  <  1:4445981/70‑2 (MQ=255)
                                               aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggt                  >  1:1718140/1‑69 (MQ=255)
                                               aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACg                    <  1:3969534/68‑1 (MQ=255)
                                                aaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg                                                  <  1:1766097/37‑1 (MQ=255)
                                                aaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATg                                                   <  1:4326532/36‑1 (MQ=255)
                                                aaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCg                                            <  1:1356147/43‑1 (MQ=255)
                                                aaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg                   >  1:1505201/1‑68 (MQ=255)
                                                aaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGa                <  1:5118812/71‑1 (MQ=255)
                                                 aCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                      >  1:1355205/1‑64 (MQ=255)
                                                  cGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                 <  1:4834124/68‑3 (MQ=255)
                                                   ggCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCt                                    <  1:3229632/48‑1 (MQ=255)
                                                    gCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGc             >  1:5367070/1‑70 (MQ=255)
                                                    gCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTCTCGCTATTGTAAACGGTGAAGc             >  1:1179099/1‑70 (MQ=255)
                                                     cAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                     <  1:1067627/61‑1 (MQ=255)
                                                         cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTg       <  1:664473/71‑1 (MQ=255)
                                                         cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGAGAAGCAATTTg       >  1:2355962/1‑71 (MQ=255)
                                                           gAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc      >  1:3037079/1‑70 (MQ=255)
                                                            aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa    >  1:3451413/1‑71 (MQ=255)
                                                            aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa    >  1:180746/1‑71 (MQ=255)
                                                             aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg                   >  1:4660370/1‑55 (MQ=255)
                                                             aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt   <  1:4163764/71‑1 (MQ=255)
                                                             aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt   >  1:1900461/1‑71 (MQ=255)
                                                              aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATc  <  1:4235171/71‑1 (MQ=255)
                                                                ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCt                                    >  1:906782/1‑35 (MQ=255)
                                                                    |                                                                
ACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATC  >  AP009048/3824256‑3824388

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: