Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,824,372 C→T 100% pseudogene (507/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,3720CT99.9% 202.5 / ‑7.3 71pseudogene (507/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (33/37);  minor base A (1/0);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
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.

AACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGA  >  AP009048/3824304‑3824439
                                                                    |                                                                   
aaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGa                                                                   <  1:5118812/71‑1 (MQ=255)
    gCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGc                                                                >  1:5367070/1‑70 (MQ=255)
    gCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTCTCGCTATTGTAAACGGTGAAGc                                                                >  1:1179099/1‑70 (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)
             aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                                      <  1:4163764/71‑1 (MQ=255)
             aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                                      >  1:1900461/1‑71 (MQ=255)
              aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATc                                                     <  1:4235171/71‑1 (MQ=255)
                  aacccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                                        <  1:1276449/61‑1 (MQ=255)
                    cccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                                >  1:202895/2‑70 (MQ=255)
                     ccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAAt                                               <  1:2003655/70‑1 (MQ=255)
                     ccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa                                              <  1:1468140/71‑1 (MQ=255)
                     ccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATATGGa                                                 >  1:1635690/1‑68 (MQ=255)
                       gATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAg                                                                 >  1:494115/1‑50 (MQ=255)
                       gATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                                        <  1:2871004/59‑1 (MQ=255)
                       gATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                                        <  1:2454485/59‑1 (MQ=255)
                       gATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACg                                            <  1:1206914/71‑1 (MQ=255)
                       gATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACg                                            <  1:2633645/71‑1 (MQ=255)
                         tGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                                          <  1:4219743/71‑1 (MQ=255)
                          gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGaa                                                                  >  1:265434/1‑46 (MQ=255)
                          gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                                          >  1:272002/1‑70 (MQ=255)
                          gAAAGGGATGGTAGACGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                                          >  1:2116264/1‑70 (MQ=255)
                           aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                         <  1:968946/70‑1 (MQ=255)
                           aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGa                                        <  1:3418086/71‑1 (MQ=255)
                           aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGa                                        <  1:4534448/71‑1 (MQ=255)
                            aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                                          >  1:2561458/1‑68 (MQ=255)
                            aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                       <  1:15194/71‑1 (MQ=255)
                             aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAg                                      <  1:1766644/71‑1 (MQ=255)
                               ggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAc                                    >  1:2867265/1‑71 (MQ=255)
                                gATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAct                                   <  1:1594620/71‑1 (MQ=255)
                                 aTGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActc                                  <  1:4927624/71‑1 (MQ=255)
                                  tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGt                                           >  1:1411444/1‑61 (MQ=255)
                                    gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGACGCAATTTGCGATCt                                                    <  1:4867931/50‑1 (MQ=255)
                                    gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                       <  1:1698921/63‑1 (MQ=255)
                                        ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                                  <  1:2303058/48‑1 (MQ=25)
                                          gCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                         <  1:4479737/71‑1 (MQ=255)
                                          gCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                           >  1:2370119/1‑69 (MQ=255)
                                          gCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                           >  1:1520678/1‑69 (MQ=255)
                                           cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGt                                           <  1:2260693/52‑1 (MQ=255)
                                           cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCagag                        <  1:1388429/71‑1 (MQ=255)
                                           cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCagag                        <  1:3135069/71‑1 (MQ=255)
                                             gTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAc                      <  1:1488457/71‑1 (MQ=255)
                                              tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                                          <  1:1866857/50‑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)
                                                  tGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACg                                            >  1:2099090/1‑44 (MQ=21)
                                                   gTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATAc                                             <  1:5008720/42‑1 (MQ=19)
                                                   gTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACAt                >  1:4176042/1‑71 (MQ=255)
                                                    tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                           <  1:5133620/59‑1 (MQ=255)
                                                    tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                           <  1:4175433/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)
                                                        taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                       >  1:3462160/2‑43 (MQ=25)
                                                          ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg           <  1:1757383/69‑1 (MQ=255)
                                                          ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGt          <  1:2453736/70‑1 (MQ=255)
                                                           tgtaAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGa                                     >  1:4081326/1‑42 (MQ=25)
                                                           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)
                                                                    |                                                                   
AACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGA  >  AP009048/3824304‑3824439

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: