Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,658,871 G→A 12.7% P357P (CCG→CCA ppx → exopolyphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,658,8710GA12.7% 167.0 / 12.6 63P357P (CCG→CCAppxexopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base G (6/49);  new base A (3/5);  total (9/54)
Fisher's exact test for biased strand distribution p-value = 7.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

ACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGTTGAATATCAAC  >  Exported/1658802‑1658931
                                                                     |                                                            
aCGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCTGGAACAGCAACCGAAGCTGGCGCATCCGc                                                             <  1:3420278/71‑1 (MQ=255)
     tGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTg                                                        <  1:1571405/71‑1 (MQ=255)
          gATACCACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGgc                                                   <  1:3824303/71‑1 (MQ=255)
          gATACCACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGgc                                                   <  1:2192599/71‑1 (MQ=255)
            tACCACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGgcgc                                                 >  1:5181152/1‑71 (MQ=255)
               cACTATGCAAATGTACGAACATTGGCGGGAACAGCAACCGAAGCTGGCGCATCCACAACTgg                                                       <  1:1998725/62‑1 (MQ=255)
               cACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCACAACTgg                                                       <  1:4721775/62‑1 (MQ=255)
               cACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCACAACTgg                                                       <  1:3641685/62‑1 (MQ=255)
               cACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCACAACTgg                                                       <  1:4440943/62‑1 (MQ=255)
                aCTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCACAACTgg                                                       <  1:2112406/61‑1 (MQ=255)
                  tATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCACAACTGGAGGCGCTACTGc                                           >  1:475909/1‑71 (MQ=255)
                  tATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCACAACTGGAGGCGCTACTGc                                           >  1:2633753/1‑71 (MQ=255)
                  tATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCACAACTGGAGGCGCTACTGc                                           >  1:3755705/1‑71 (MQ=255)
                    tGCAAATGTACGAACAGTGGCGGGAACAGCAACTGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGa                                         >  1:2379050/1‑71 (MQ=255)
                    tGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTg                                            <  1:3480246/68‑1 (MQ=255)
                      cAAATGTACGAACAGTTGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGAt                                        <  1:518715/70‑1 (MQ=255)
                      cAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGAt                                        <  1:4203805/70‑1 (MQ=255)
                      cAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATg                                       <  1:1900269/71‑1 (MQ=255)
                      cAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATg                                       <  1:4420442/71‑1 (MQ=255)
                       aaaTGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGAt                                        <  1:251663/69‑1 (MQ=255)
                       aaaTGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATgg                                      <  1:4786114/71‑1 (MQ=255)
                       aaaTGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATg                                       <  1:1403103/70‑1 (MQ=255)
                       aaaTGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATg                                       <  1:2884352/70‑1 (MQ=255)
                        aaTGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATggg                                     <  1:1588765/71‑1 (MQ=255)
                        aaTGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATggc                                     <  1:2485885/71‑2 (MQ=255)
                        aaTGTACGAACAGTGGCGGGAACAGCAACCGAAGATGGCGCATCCGCAACTGGAGGCGCTACTGCGATggg                                     <  1:959902/71‑1 (MQ=255)
                         aTGTACGAACAGTGGCGGGAACAGCAACCTAAGCTGGCGCATCCGc                                                             <  1:327614/46‑1 (MQ=255)
                         aTGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGc                                                             <  1:3133695/46‑1 (MQ=255)
                         aTGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGc                                                             <  1:1833610/46‑1 (MQ=255)
                         aTGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGc                                                             <  1:544182/46‑1 (MQ=255)
                                   aGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCt                                                <  1:4892124/49‑1 (MQ=255)
                                   aGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCt                                                <  1:5007386/49‑1 (MQ=255)
                                   aGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCt                                                <  1:565121/49‑1 (MQ=255)
                                   aGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCt                                                <  1:1355999/49‑1 (MQ=255)
                                   aGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCt                                                <  1:844200/49‑1 (MQ=255)
                                    gTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCAtgctg                          >  1:3739563/1‑70 (MQ=255)
                                    gTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCAtgctg                          >  1:3863871/1‑70 (MQ=255)
                                       gCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGgc                                                   <  1:4149296/42‑1 (MQ=255)
                                         gggAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTATTGCGAt                                        <  1:9836/51‑1 (MQ=255)
                                         gggAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATg                                       <  1:135283/52‑1 (MQ=255)
                                         gggAACAGCAACCGAAGCGGGCGCATCCGCAACTGGAGGCGCTACTGCGATg                                       <  1:3409147/52‑1 (MQ=255)
                                           gAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGt                  <  1:796742/71‑1 (MQ=255)
                                            aaCAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATg                                       <  1:5013665/49‑1 (MQ=255)
                                             aCAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCg                <  1:4610877/71‑1 (MQ=255)
                                             aCAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCg                <  1:753733/71‑1 (MQ=255)
                                             aCAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCg                <  1:1521391/71‑1 (MQ=255)
                                             aCAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCg                <  1:1652537/71‑1 (MQ=255)
                                              cAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTg                                            <  1:4785590/42‑1 (MQ=255)
                                                gCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTg                                            <  1:1262017/40‑1 (MQ=255)
                                                 cAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGt                  <  1:3886468/65‑1 (MQ=255)
                                                 cAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCg                <  1:4744732/67‑1 (MQ=255)
                                                  aaCCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGTTg           <  1:1251660/71‑1 (MQ=255)
                                                  aaCCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGTTg           <  1:1515770/71‑1 (MQ=255)
                                                  aaCCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGTTg           <  1:4691287/71‑1 (MQ=255)
                                                  aaCCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGTTg           <  1:3589504/71‑1 (MQ=255)
                                                   aCCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGttag          <  1:2150800/71‑3 (MQ=255)
                                                   aCCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGttag          <  1:3098259/71‑3 (MQ=255)
                                                   aCCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGttag          <  1:4432055/71‑3 (MQ=255)
                                                   aCCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGTTGa          <  1:4571143/71‑1 (MQ=255)
                                                       taGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGt                  <  1:3279671/58‑1 (MQ=255)
                                                       aaGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGTTg           <  1:3921552/66‑1 (MQ=255)
                                                                    cGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGTTGAATATCAAc  >  1:2924265/1‑62 (MQ=255)
                                                                    cGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGTTGAATATCAAc  >  1:1417490/1‑62 (MQ=255)
                                                                     |                                                            
ACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGTTGAATATCAAC  >  Exported/1658802‑1658931

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: