Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,651,710 T→G 5.0% L346R (CTT→CGT)  aslA → acrylsulfatase‑like enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,651,7100TG 5.0% 324.9 / 10.9 119L346R (CTT→CGT) aslAacrylsulfatase‑like enzyme
Reads supporting (aligned to +/- strand):  ref base T (47/66);  new base G (3/3);  total (50/69)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

Only 100 of 128 total aligned reads displayed.

GCGACTGCATGGTGGAGATGAACGATGTGTTCGCTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCGAAGTACCGCCGCACGGACGCACC  >  AP009048/3651641‑3651780
                                                                     |                                                                      
gCGACTGCATGGTGGAGATGAACGATGTGTTCGCTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCtt                                                                       >  1:5258987/1‑71 (MQ=255)
gCGACTGCATGGTGGAGATGAACGATGTGTTCGCTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCt                                                                        <  1:863631/70‑1 (MQ=255)
  gACTGCATGGTGGAGATGAACGATGTGTTCGCTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCTTg                                                                      >  1:6610982/1‑70 (MQ=255)
  gACTGCATGGTGGAGATGAACGATGTGTTCGCTAATCTGTATAAAACACTGGAGAAAAACGATCAGCGTgt                                                                     >  1:6508219/1‑70 (MQ=255)
  gACTGCATGGTGGAGATGAACGATGTGTTCGCTAATCTGTATAAAACACTGGAGAAAAACGATCAGCGTgt                                                                     >  1:5962358/1‑70 (MQ=255)
   aCTGCATGGTGGAGATGAACGATGTGTTCGCTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGa                                                                     >  1:6371861/1‑70 (MQ=255)
        aTGGTGGAGATGAACGATGTGTTCGCTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAacac                                                               >  1:3828792/1‑71 (MQ=255)
         tggtggAGATGAACGATGTGTTCGCTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACg                                                              >  1:6098614/1‑71 (MQ=255)
          ggtggAGATGAACGATGTGTTCGCTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGc                                                             >  1:3420945/1‑71 (MQ=255)
           gtggAGATGAACGATGTGTTCGCTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCt                                                            >  1:1583546/1‑71 (MQ=255)
               agaTGAACGATGTGTTCGCTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATc                                                        <  1:1541950/71‑1 (MQ=255)
                 aTGAACGATGTGTTCGCTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGt                                                      <  1:7528742/71‑1 (MQ=255)
                  tGAACGATGTGTTCGCTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCtt                                                                       <  1:3288355/53‑1 (MQ=255)
                  tGAACGATGTGTTCGCTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTc                                                     <  1:1769149/71‑1 (MQ=255)
                  tGAACGATGTGTTCGCTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTc                                                     <  1:7732248/71‑1 (MQ=255)
                   gAACGATGTGTTCGCTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCt                                                    >  1:303139/1‑71 (MQ=255)
                    aaCGATGTGTTCGCTAATCTGTATAAAACACTGGAGAAAAACGATCAGCGTGTTATCAAGCTGACCGt                                                      <  1:7399474/68‑1 (MQ=21)
                     aCGATGTGTTCGCTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCg                                                       <  1:4760911/66‑1 (MQ=255)
                     aCGATGTGTTCGCTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCtt                                                   >  1:2092084/1‑70 (MQ=255)
                     aCGATGTGTTCGCTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCtt                                                   >  1:3550573/1‑70 (MQ=255)
                      cGATGTGTTCGCTAATCTGTATAAAACACTGGAGAAAGACGATCAGCGTGTTATCAAGCTGACCGTttttc                                                 <  1:4079353/71‑2 (MQ=14)
                           tgtTCGCTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTcc                                            >  1:7518265/1‑71 (MQ=255)
                            gtTCGCTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTcc                                            >  1:7003700/1‑70 (MQ=255)
                             ttCGCTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCg                                           <  1:1920022/70‑1 (MQ=255)
                             ttCGCTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCg                                           <  1:4786007/70‑1 (MQ=255)
                              tCGCTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCtt                                                                       >  1:3097622/1‑41 (MQ=255)
                               cGCTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGAt                                         <  1:3919313/70‑1 (MQ=255)
                               cGCTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGAt                                         <  1:3854636/70‑1 (MQ=255)
                                gCTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCTTg                                                                      <  1:4957072/40‑1 (MQ=255)
                                gCTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATa                                        <  1:2082236/70‑1 (MQ=255)
                                 cTAATCTGTATAAAACACTGTAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAAc                                      <  1:7900438/71‑1 (MQ=255)
                                 cTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTg                                                           <  1:7693941/50‑1 (MQ=255)
                                 cTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTg                                                           <  1:3323130/50‑1 (MQ=255)
                                 cTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTg                                                           <  1:2631101/50‑1 (MQ=255)
                                  tAATCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACg                                     >  1:6487893/1‑71 (MQ=255)
                                  aaaTCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAAc                                      >  1:3785644/2‑70 (MQ=255)
                                    aTCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACg                                                              >  1:3946687/1‑44 (MQ=255)
                                    aTCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACgg                                    <  1:7511200/70‑1 (MQ=255)
                                    aTCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACgg                                    <  1:7935701/70‑1 (MQ=255)
                                    aTCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACg                                     <  1:3373436/69‑1 (MQ=255)
                                    aTCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGa                                   >  1:6083749/1‑71 (MQ=255)
                                    aTCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGa                                   <  1:4624653/71‑1 (MQ=255)
                                     tCTTTATAAAACACTGGAGAAAAACGGTCAGCTTGATAacac                                                               <  1:7188341/42‑1 (MQ=255)
                                        gTATAAAACACTGGAGAAAAACGGTCAGCTTGATAacac                                                               <  1:1498879/39‑1 (MQ=255)
                                        gTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTg                                                           <  1:3165267/43‑1 (MQ=255)
                                         tataAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGAt                                         <  1:6840111/60‑1 (MQ=255)
                                         tataAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGAcc                                 <  1:3672755/68‑1 (MQ=255)
                                         tataAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGa                              <  1:2714119/71‑1 (MQ=255)
                                           taAAACACTGGAGAAAAACGGTCAGCTTGATAacac                                                               <  1:4393311/36‑1 (MQ=255)
                                           taAAACACTGGAGAAAAACGGTCAGCTTGATAacac                                                               <  1:4973333/36‑1 (MQ=255)
                                           taAAACACTGGAGAAAAACGGTCAGCTTGATAacac                                                               <  1:3143662/36‑1 (MQ=255)
                                           taAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAg                            >  1:7624872/1‑71 (MQ=255)
                                           taAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAg                            >  1:6687869/1‑71 (MQ=255)
                                             aaaCACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATa                                        >  1:6553007/1‑57 (MQ=255)
                                             aaaCACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATa                                        >  1:3016634/1‑57 (MQ=255)
                                               acacTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCg                         >  1:5065768/1‑70 (MQ=255)
                                                cacTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAAc                                      >  1:3957718/1‑56 (MQ=255)
                                                cacTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCg                         >  1:2239563/1‑69 (MQ=255)
                                                cacTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCg                         >  1:3210902/1‑69 (MQ=255)
                                                cacTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCg                         >  1:2617916/1‑69 (MQ=255)
                                                cacTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCGaa                       <  1:5155610/71‑1 (MQ=255)
                                                 acTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCt                                                    <  1:6897878/41‑1 (MQ=255)
                                                     gagaAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCGAAGTAcc                  <  1:6173282/71‑1 (MQ=255)
                                                     gagaAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCGAAGTAc                   <  1:1066787/70‑1 (MQ=255)
                                                       gaAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCGAAGTAccgc                >  1:7186458/1‑71 (MQ=255)
                                                        aaaaaCGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCg                                           <  1:5680289/43‑1 (MQ=255)
                                                        aaaaaCGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCg                                           <  1:5846948/43‑1 (MQ=255)
                                                        aaaaaCGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCg                                           <  1:6000524/43‑1 (MQ=255)
                                                        aaaaaCGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCg                                           <  1:4247402/43‑1 (MQ=255)
                                                        aaaaaCGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGa                                          >  1:563713/1‑44 (MQ=255)
                                                        aaaaaCGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGa                                          >  1:5704726/1‑44 (MQ=255)
                                                         aaaaCGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCGAAGTAccgccg              <  1:7693055/71‑1 (MQ=255)
                                                         aaaaCGGGCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGa                                   <  1:5844962/50‑1 (MQ=255)
                                                          aaaCGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGaa                             >  1:4940237/1‑55 (MQ=255)
                                                          aaaCGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGaa                             >  1:2747732/1‑55 (MQ=255)
                                                          aaaCGGTCAGCTTGATAACACGCTGATCGTCTTCACCTCCGATAACGGACCGGAAGCCGAAGTAccgccg              >  1:1682535/1‑70 (MQ=255)
                                                            aCGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCGAAGTACCGCCGCAc           >  1:4719627/1‑71 (MQ=255)
                                                            aCGATCAGCGTGTTATCAAGCTGACCGTTTTTACCTCCGATAACGGACCGGAAGCCGAAGTACCGCCGCAc           >  1:4769844/1‑71 (MQ=18)
                                                              ggTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCGAAGt                     <  1:6333968/59‑1 (MQ=255)
                                                              ggTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCGAAGt                     <  1:1829372/59‑1 (MQ=255)
                                                              ggTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCGAAGt                     <  1:2585041/59‑1 (MQ=255)
                                                              ggTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCGAAGTACCGCCGCACgg         <  1:1412165/71‑1 (MQ=255)
                                                              ggTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCGAAGTACCGCCGCACgg         <  1:3269563/71‑1 (MQ=255)
                                                              ggTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCGAAGTACCGCAGCACgg         <  1:2022395/71‑1 (MQ=255)
                                                               gTCAGCTTGATAACACGCTGATCGTCTTTACCTcc                                            <  1:2459520/35‑1 (MQ=255)
                                                               gTCAGCTTGATAACACGCTGATCGTCTTTACCTcc                                            <  1:3285592/35‑1 (MQ=255)
                                                               gTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCGAAGTACCGCCGCACGGa        <  1:4480103/71‑1 (MQ=255)
                                                               gTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCGAAGTACCGCCGCACGGa        >  1:5381695/1‑71 (MQ=255)
                                                               gTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCGAAGTACCGCCGCACGGa        >  1:5372847/1‑71 (MQ=255)
                                                                tCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCGAAGTACCGCCGCACg          >  1:2873539/1‑68 (MQ=255)
                                                                  aGCGTGTTATCAAGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCGAAGTACCGCCGCACGGACg      <  1:2140792/70‑1 (MQ=38)
                                                                    cTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCGAAGTAccg                 >  1:6684135/1‑57 (MQ=255)
                                                                    cTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCGAAGTAccg                 >  1:3092069/1‑57 (MQ=255)
                                                                    cTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCGAAGTAccg                 >  1:1183052/1‑57 (MQ=255)
                                                                    cTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCGAAGTACCGCCGCACGGa        >  1:1553818/1‑66 (MQ=255)
                                                                     ttGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGa                              <  1:6796697/43‑1 (MQ=255)
                                                                     ttGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGa                              <  1:2264386/43‑1 (MQ=255)
                                                                     ttGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCGAAg                      <  1:4558498/51‑1 (MQ=255)
                                                                     ttGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCGAAg                      <  1:3346516/51‑1 (MQ=255)
                                                                     ttGATAACACGCTGATCGTCTTTACATCCGATAACGGACCGGAAGCCGAAGTACCGCCGCACGGACGCAcc  >  1:1328946/1‑71 (MQ=255)
                                                                     |                                                                      
GCGACTGCATGGTGGAGATGAACGATGTGTTCGCTAATCTGTATAAAACACTGGAGAAAAACGGTCAGCTTGATAACACGCTGATCGTCTTTACCTCCGATAACGGACCGGAAGCCGAAGTACCGCCGCACGGACGCACC  >  AP009048/3651641‑3651780

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: