Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,965,130 A→G 11.1% K204K (AAA→AAG deoA → thymidine phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,965,1300AG11.1% 170.5 / 9.9 63K204K (AAA→AAGdeoAthymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (26/30);  new base G (4/3);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGT  >  Exported/2965063‑2965198
                                                                   |                                                                    
tGATCACCGGCTCGATCCTGTCCAAGAAGCTGTCGGCGGGCCTGGATGCGCTGGTGATGGACGTGAAGGtc                                                                   >  1:3765659/1‑70 (MQ=255)
tGATCACCGGCTCGATCCTGTCCAAGAAGCTGTCGGCGGGCCTGGATGCGCTGGTGATGGACGTGAAGGtc                                                                   >  1:627910/1‑70 (MQ=255)
tGATCACCGGCTCGATCCTGTCCAAGAAGCTGTCGGCGGGCCTGGATGCGCTGGTGATGGACGTGAAGGtc                                                                   >  1:4654633/1‑70 (MQ=255)
tGATCACCGGCTCGATCCTGTCCAAGAAGCTGTCGGCGGGCCTGGATGCGCTGGTGATGGACGTGAAGGtc                                                                   >  1:5139499/1‑70 (MQ=255)
 gATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTgg                                                                  >  1:2511946/1‑71 (MQ=255)
      ccGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAgc                                                             <  1:1564036/71‑1 (MQ=255)
      ccGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAgc                                                             <  1:1572932/71‑1 (MQ=255)
       cggctcgATCCTGTCCAAGAAGCTGTCGGCGGGCCTGGATGCGCTGGTGATGGACGTGAAGGTCGGCAgcg                                                            <  1:4421560/68‑1 (MQ=255)
       cggctcgATCCTGTCCAAGAAGCTGTCGGCGGGCCTGGATGCGCTGGTGATGGACGTGAAGGTCGGCAgcg                                                            <  1:5096278/68‑1 (MQ=255)
       cggctcgATCCTGTCCAAGAAGCTGTCGGCGGGCCTGGATGCGCTGGTGATGGACGTGAAGGTCGGCAgcg                                                            <  1:5072468/68‑1 (MQ=255)
         ccTCTATTCTGGCGAAGAAACTAGCGGAAGGACTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAgcggc                                                          >  1:1677587/1‑71 (MQ=255)
           tctATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGaaa                                                                      >  1:5198074/1‑57 (MQ=255)
             tATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGt                                                      >  1:3431812/1‑71 (MQ=255)
             tATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGt                                                      >  1:3394182/1‑71 (MQ=255)
               ttCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGtt                                                     <  1:4676491/70‑1 (MQ=255)
                 cTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTa                                                   >  1:3790149/1‑70 (MQ=255)
                    gCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTAt                                                  <  1:288183/68‑1 (MQ=255)
                    gCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATg                                                 <  1:3369350/69‑1 (MQ=255)
                     cGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTa                                                   >  1:1580525/1‑66 (MQ=255)
                     cGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTa                                                   >  1:4314104/1‑66 (MQ=255)
                     cGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTa                                                   >  1:2190962/1‑66 (MQ=255)
                     cGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTa                                                   >  1:1118185/1‑66 (MQ=255)
                     cGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGAGCGTTTa                                                   >  1:1894959/1‑66 (MQ=255)
                      gaagaaACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTa                                                   <  1:3338671/65‑1 (MQ=255)
                      gaagaaACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTa                                                   <  1:3469375/65‑1 (MQ=255)
                          aaaCTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGtt                                                     <  1:1211473/59‑1 (MQ=255)
                          aaaCTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGtt                                                     <  1:414006/59‑1 (MQ=255)
                           aaCTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGaaa                                                                      >  1:1272662/1‑41 (MQ=255)
                           aaCTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGaaa                                                                      >  1:3793051/1‑41 (MQ=255)
                           aaCTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGaaa                                                                      >  1:3816621/1‑41 (MQ=255)
                           aaCTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGaaa                                                                      >  1:1756997/1‑41 (MQ=255)
                           aaCTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGaaa                                                                      >  1:2364113/1‑41 (MQ=255)
                           aaCTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGaaa                                                                      >  1:4293125/1‑41 (MQ=255)
                           aaCTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGaaa                                                                      >  1:342745/1‑41 (MQ=255)
                           aaCTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGaaa                                                                      >  1:443967/1‑41 (MQ=255)
                           aaCTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGaaa                                                                      >  1:1840212/1‑41 (MQ=255)
                           aaCTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTAc                                        >  1:2660731/1‑71 (MQ=255)
                            aCTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTAc                                        >  1:4251476/1‑70 (MQ=255)
                            aCTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTAc                                        >  1:1433811/1‑70 (MQ=255)
                             cTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGttt                                                    <  1:428748/57‑1 (MQ=255)
                             cTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGttt                                                    <  1:1124774/57‑1 (MQ=255)
                             cTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGttt                                                    <  1:2411933/57‑1 (MQ=255)
                             cTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGa                                      <  1:3479313/71‑1 (MQ=255)
                                         cTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGATCTCTCTGAAGc                          >  1:4144843/1‑71 (MQ=255)
                                         cTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGc                          >  1:4678020/1‑71 (MQ=255)
                                         cTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGc                          >  1:4575651/1‑71 (MQ=255)
                                          tGGACGCGCGGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGcc                                               <  1:2635410/49‑1 (MQ=255)
                                            gACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGttt                                                    <  1:4017494/42‑1 (MQ=255)
                                            gACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGCCCTACGAACTCTCTGAAGCCCt                       >  1:1654844/1‑71 (MQ=255)
                                            gACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAAc                                    >  1:3747509/1‑58 (MQ=255)
                                            gACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCt                       >  1:682359/1‑71 (MQ=255)
                                            gACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCt                       >  1:4248173/1‑71 (MQ=255)
                                              cgcgCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAAc                                    >  1:2877567/1‑56 (MQ=255)
                                              cgcgCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAAc                                    >  1:4830980/1‑56 (MQ=255)
                                              cgcgCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAAc                                    >  1:1480380/1‑56 (MQ=255)
                                              cgcgCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAAc                                    >  1:3056555/1‑56 (MQ=255)
                                                 gcTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCg                  <  1:4691202/71‑1 (MQ=255)
                                                   ctgTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGaa                <  1:2709613/69‑1 (MQ=255)
                                                      tGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATg                                                 <  1:3075162/35‑1 (MQ=255)
                                                      tGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCg                                              <  1:2928199/38‑1 (MQ=255)
                                                      tGATGGACGTGAAAGAGGGTAGCGGCGCGTTTATGCCg                                              <  1:1760826/38‑1 (MQ=255)
                                                        aTGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCg                                              <  1:3047542/36‑1 (MQ=255)
                                                        aTGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCg                                              <  1:2001755/36‑1 (MQ=255)
                                                        aTGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCg                                              <  1:1976977/36‑1 (MQ=255)
                                                        aTGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCg                                              <  1:4479141/36‑1 (MQ=255)
                                                        aTGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCg                                              <  1:1241880/36‑1 (MQ=255)
                                                        aTGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCg                                              <  1:1918556/36‑1 (MQ=255)
                                                        aTGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCg                                              <  1:3598476/36‑1 (MQ=255)
                                                        aTGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCg                                              <  1:8531/36‑1 (MQ=255)
                                                        aTGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGAt           <  1:4211468/71‑1 (MQ=255)
                                                          ggACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGAttg         >  1:2020882/1‑71 (MQ=255)
                                                          ggACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGAtt          >  1:1972365/1‑70 (MQ=255)
                                                               tGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAAct                                   <  1:3839671/40‑1 (MQ=255)
                                                                 aaaGTGGGTAGCGGCGCGTTTATGCCGACCTACGaa                                     >  1:3777475/1‑36 (MQ=255)
                                                                 aaaGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGt  >  1:1671103/1‑71 (MQ=255)
                                                                 aaaGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGt  <  1:5174404/71‑1 (MQ=255)
                                                                 aaaGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGt  <  1:2119485/71‑1 (MQ=255)
                                                                  aaGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGt  <  1:961750/70‑1 (MQ=255)
                                                                   |                                                                    
TGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGT  >  Exported/2965063‑2965198

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: