Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 64,346 C→G 10.7% A575P (GCC→CCC)  hepA ← ATP‑dependent helicase HepA

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL60664,3460C→G10.7% 2.2 ‑3.0 ‑0.0451723hepAATP‑dependent helicase HepA
Reads supporting (aligned to +/- strand):  new base (4/1):  ref base (12/27):  total (16/29)
Fisher's exact test strand distribution p-value = 5.09e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.11e-01

GTTGAATGGCAGGTCAAACATCACCATGTGGCTGGCGAACTGGAAGTTACGTCCTTCAGAACCGATTTCTG  >  REL606/64311‑64381
                                   |                                   
gTTGAATGGCAGGTCAAACATCACCATGtggctggc                                     <  1:1609569/36‑1 (MQ=255)
gTTGAATGGCAGGTCAAACATCACCATGtggctggc                                     <  1:7110130/36‑1 (MQ=255)
  tGAATGGCAGGTCAAACATCACCATGTGGCTGGCGa                                   <  1:4968527/36‑1 (MQ=255)
   gAATGGCAGGTCAAACATCACCATGTGGCTGGCGaa                                  <  1:1553176/36‑1 (MQ=255)
   gAATGGCAGGTCAAACATCACCATGTGGCTGGCGaa                                  <  1:4098812/36‑1 (MQ=255)
    aaTGGCAGGTCAAACATCACCATGTGGCTGGGGAAc                                 <  1:5999892/36‑1 (MQ=255)
    aaTGGCAGGTCAAACATCACCATGTGGCTGGCGAAc                                 <  1:271436/36‑1 (MQ=255)
      tGGCAGGTCAAACATCACCATGTGGCTGGGGGACTg                               >  1:4889685/1‑36 (MQ=37)
       ggCAGGTCAAACATCACCATGTGGCTGGCGAACTgg                              <  1:6001603/36‑1 (MQ=255)
       ggCAGGTCAAACATCACCATGTGGCTGGCGAACTgg                              <  1:5833448/36‑1 (MQ=255)
        gCAGGTCAAACATCACCATGTGGCTGGGGGACTggg                             >  1:7686456/1‑35 (MQ=37)
        gCAGGTCAAACATCACCATGTGGCTGGCGAACTGGa                             <  1:6112884/36‑1 (MQ=255)
          aGGTCAAACATCACCATGTGGCTGGCGGACTGGGAg                           >  1:6560756/1‑36 (MQ=37)
          aGGTCAAACATCACCATGTGGCTGGCGAACCGGGAg                           >  1:6547997/1‑36 (MQ=37)
           ggTCAAACATCACCATGTGGCTGGCGGACTGGAAgg                          >  1:201237/1‑35 (MQ=255)
           ggTCAAACATCACCATGTGGCTGGCGAACTGGAAGt                          <  1:2346282/36‑1 (MQ=255)
             tCAAACATCACCATGTGGCTGGCGAACTGGGAGttt                        >  1:647881/1‑35 (MQ=255)
             tCAAACATCACCATGTGGCTGGCGAACTGGAAGTTa                        <  1:2733820/36‑1 (MQ=255)
              cAAACATCACCATGTGGCTGGCGGACTCGAAGTTaa                       >  1:3527037/1‑35 (MQ=37)
               aaaCATCTCCAGGTGGCTGGAGAGCTGGAAGTTACg                      <  1:6783525/36‑1 (MQ=19)
               aaaCATCACCATGTGGCTGGCGAACTGGAAGTTACg                      <  1:3732241/36‑1 (MQ=255)
                aaCATCACCATGTGGCTGGCGAACTTGAAGTTACGt                     >  1:6811242/1‑36 (MQ=255)
                aaCATCACCATGTGGCTGGCGAACTGGAAGTTACGt                     <  1:5207694/36‑1 (MQ=255)
                aaCATCACCATGTGGCTGGCGAACTGGAAGTTACGt                     <  1:22366/36‑1 (MQ=255)
                aaCATCACCATGTGGCTGGCGAACTGGAAGTTACGt                     <  1:7736959/36‑1 (MQ=255)
                 aCATCACCATGTGGCTGGCGAACTGGAAGTTACGTc                    <  1:3745730/36‑1 (MQ=255)
                  cATCACCATGTGGCTGGCGGACTGGAAGTTACGTcc                   >  1:2368169/1‑36 (MQ=255)
                  cATCACCATGTGGCTGGCGAACTGGGAGTTACGTcc                   >  1:6181848/1‑36 (MQ=255)
                   aTCACCATGTGGCTGGGGAACTGGAAGTTACGTTCt                  >  1:2697829/1‑36 (MQ=37)
                   aTCACCATGTGGCTGGCGGACTTGAAGTTTCTTCCt                  >  1:5040044/1‑36 (MQ=19)
                    tCACCATGTGGCTGGCGGGCTTGGAGTTTCGTCCtt                 >  1:2840601/1‑36 (MQ=16)
                    tCACCATGTGGCTGGCGAACTGGAAGTTACGTCCtt                 <  1:7822684/36‑1 (MQ=255)
                     ctccATGTGGCTGGCGAACTGGAAGTTACGTCCTTc                <  1:534287/34‑1 (MQ=255)
                     cACCATGTGGCTGGGCGACTGGAAGTTTCCTCCTTc                >  1:5941767/1‑36 (MQ=16)
                      aCCATGTGGCTGGCGAACTGGAAGTTACGTCCTTCa               <  1:3677741/36‑1 (MQ=255)
                      aCCATGTGGCTGGCGAACTGGAAGTTACGTCCTTCa               <  1:6163292/36‑1 (MQ=255)
                      aCCATGTGGCTGGCGAACTGGAAGTTACGTCCTTCa               <  1:799755/36‑1 (MQ=255)
                      aCCATGTGGCTGGCCAACTGGAAGGTACGTCCTTcc               >  1:8320948/1‑35 (MQ=37)
                       ccATGTGGCTGGCGAACTGGAAGTTACGTCCTTCAg              <  1:6268166/36‑1 (MQ=255)
                        cATGTGGCTGGCGAACTGGAAGTTACGTCCTTCAGa             <  1:6103647/36‑1 (MQ=255)
                        cATGTGGCTGGCGAACTGGAAGTTACGTCCTTCAGa             <  1:4438713/36‑1 (MQ=255)
                        cATGTGGCTGGCGAACTGGAAGTTACGTCCTTCAGa             <  1:6531303/36‑1 (MQ=255)
                        cATGTGGCTGGCGAACTGGAAGTTACGTCCTTCAGa             <  1:7047233/36‑1 (MQ=255)
                        cATGTGGCTGGCGAACTGGAAGTTACGTCCTTCAGa             >  1:5665776/1‑36 (MQ=255)
                         aTGTGGCTGGCGAACTGGAAGTTACGTCCTTCAGaa            <  1:7301085/36‑1 (MQ=255)
                           gtgGCTGGCGAACTGGAAGTTACGTCCTTCAGAAcc          <  1:1186028/36‑1 (MQ=255)
                           gtgGCTGGCGAACTGGAAGTTACGTCCTTCAGAAcc          <  1:2633364/36‑1 (MQ=255)
                            tggctggcGAACTGGAAGTTACGTCCTTCAGAACCg         <  1:956809/36‑1 (MQ=255)
                            tggctggcGAACTGGAAGTTACGTCCTTCAGAACCg         >  1:2726856/1‑36 (MQ=255)
                            tggctggcGAACTGGAAGTTACGTCCTTCAGAACCg         >  1:5329739/1‑36 (MQ=255)
                             ggctggcGAACTGGAAGTTACGTCCTTCAGAACCGa        <  1:8459627/36‑1 (MQ=255)
                             ggctggcGAACTGGAAGTTACGTCCTTCAGAACCGa        <  1:216924/36‑1 (MQ=255)
                             ggctggcGAACTGGAAGTTACGTCCTTCAGAACCGa        <  1:295633/36‑1 (MQ=255)
                             ggctggcGAACTGGAAGTTACGTCCTTCAGAACCGa        >  1:614676/1‑36 (MQ=255)
                             ggctggcGAACTGGAAGTTACGTCCTTCAGAACCGa        >  1:6168808/1‑36 (MQ=255)
                              gctggtGAACTGGAAGTTACGTCCTTCAGAACCGAt       <  1:7161316/36‑1 (MQ=255)
                              gctggcGAACTGGAAGTTACGTCCTTCAGAACCGAt       >  1:3338744/1‑36 (MQ=255)
                              gctggcGAACTGGAAGTTACGTCCTTCAGAACCGAt       <  1:4924990/36‑1 (MQ=255)
                                tggcGAACTGGAAGTTACGTCCTTCAGAACCGAttt     <  1:596707/36‑1 (MQ=255)
                                tggcGAACTGGAAGTTACGTCCTTCAGAACCGAttt     >  1:1722994/1‑36 (MQ=255)
                                tggcGAACTGGAAGTTACGTCCTTCAGAACCGAttt     >  1:4989923/1‑36 (MQ=255)
                                tggcGAACTGGAAGTTACGTCCTTCAGAACCGAttt     >  1:7218104/1‑36 (MQ=255)
                                tggcGAACTGGAAGTTACGTCCTTCAGAACCGAttt     <  1:5270875/36‑1 (MQ=255)
                                tggcGAACTGGAAGTTACGTCCTTCAGAACCGAttt     <  1:3127794/36‑1 (MQ=255)
                                tggcGAACTGGAAGTTACGTCCTTCAGAACCGAttt     <  1:4412236/36‑1 (MQ=255)
                                tggcGAACTGGAAGTTACGTCCTTCAGAACCGAttt     >  1:3769590/1‑36 (MQ=255)
                                tggcGAACTGGAAGTTACGTCCTTCAGAACCGAttt     <  1:6387483/36‑1 (MQ=255)
                                 ggcGAACTGGAAGTTACGTCCTTCAGAACCGATTTc    >  1:3229622/1‑36 (MQ=255)
                                 ggcGAACTGGAAGTTACGTCCTTCAGAACCGATTTc    <  1:4873767/36‑1 (MQ=255)
                                 ggcGAACTGGAAGTTACGTCCTTCAGAACCGATTTc    >  1:3653565/1‑36 (MQ=255)
                                 ggcGAACTGGAAGTTACGTCCTTCAGAACCGATTTc    <  1:5893111/36‑1 (MQ=255)
                                 ggcGAACTGGAAGTTACGTCCTTCAGAACCGATTTc    >  1:1884200/1‑36 (MQ=255)
                                  gcGAACTGGAAGTTACGTCCTTCAGAACCGATTTCt   <  1:3547070/36‑1 (MQ=255)
                                  gcGAACTGGAAGTTACGTCCTTCAGAACCGATTTCt   <  1:2005342/36‑1 (MQ=255)
                                  gcGAACTGGAAGTTACGTCCTTCAGAACCGATTTCt   <  1:3823782/36‑1 (MQ=255)
                                  gcGAACTGGAAGTTACGTCCTTCAGAACCGATTTCt   <  1:586477/36‑1 (MQ=255)
                                  gcGAACTGGAAGTTACGTCCTTCAGAACCGATTTCt   >  1:5258329/1‑36 (MQ=255)
                                   cGAACTGGAAGTTACGTCCTTCAGAACCGATTTCTg  >  1:6268957/1‑36 (MQ=255)
                                   cGAACTGGAAGTTACGTCCTTCAGAACCGATTTCTg  <  1:4035636/36‑1 (MQ=255)
                                   |                                   
GTTGAATGGCAGGTCAAACATCACCATGTGGCTGGCGAACTGGAAGTTACGTCCTTCAGAACCGATTTCTG  >  REL606/64311‑64381

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG