Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,608,514 T→C V351A (GTC→GCC)  deoA → thymidine phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,608,5140TC99.3% 613.6 / ‑5.2 152V351A (GTC→GCC) deoAthymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (79/72);  total (80/72)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 157 total aligned reads displayed.

GTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGA  >  REL606/4608441‑4608586
                                                                         |                                                                         
gTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCa                                                                         <  2:2691652/76‑1 (MQ=255)
gTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCa                                                                         >  1:5407051/1‑76 (MQ=255)
gTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCa                                                                         <  1:5304710/76‑1 (MQ=255)
  tGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGt                                                                       <  2:3912317/76‑1 (MQ=255)
  tGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGt                                                                       >  2:2046321/1‑76 (MQ=255)
  tGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGt                                                                       <  1:307154/76‑1 (MQ=255)
   gagaACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTg                                                                      >  2:3771405/1‑76 (MQ=255)
   gagaACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTg                                                                      <  1:32331/76‑1 (MQ=255)
    agaACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGa                                                                     >  2:5030573/1‑76 (MQ=255)
    agaACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGa                                                                     >  2:1457233/1‑76 (MQ=255)
     gaACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGaa                                                                    <  1:1071129/76‑1 (MQ=255)
     gaACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGaa                                                                    <  1:1563552/76‑1 (MQ=255)
      aaCTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGaaa                                                                   <  1:1038139/76‑1 (MQ=255)
      aaCTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGaaa                                                                   <  2:3794302/76‑1 (MQ=255)
      aaCTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTGTGCCCGTTaaa                                                                   >  2:754751/1‑76 (MQ=255)
       aCTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAAt                                                                  <  1:6065820/76‑1 (MQ=255)
        cTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGCAATg                                                                 >  2:1361203/1‑76 (MQ=255)
        cTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATg                                                                 >  2:5624948/1‑76 (MQ=255)
           cgcgAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGAt                                                              <  1:6302549/76‑1 (MQ=255)
            gcgAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATa                                                             <  1:3280223/76‑1 (MQ=255)
             cgAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATAc                                                            <  2:2125097/76‑1 (MQ=255)
             cgAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATAc                                                            >  1:5443245/1‑76 (MQ=255)
              gAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATAcc                                                           <  1:2907979/76‑1 (MQ=255)
              gAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATAcc                                                           <  1:818916/76‑1 (MQ=255)
              gAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATAcc                                                           >  2:2674854/1‑76 (MQ=255)
               aaGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATAccc                                                          >  2:3194176/1‑76 (MQ=255)
               aaGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATAacc                                                          >  2:4214797/1‑76 (MQ=255)
                aGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCcg                                                         <  1:5177632/76‑1 (MQ=255)
                  tATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCcgcg                                                       <  1:4181755/76‑1 (MQ=255)
                  tATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATAACcgcg                                                       >  2:3143974/1‑76 (MQ=255)
                  tATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATAACcgcg                                                       >  2:5028437/1‑76 (MQ=255)
                  tATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATAACcccg                                                       >  1:208448/1‑76 (MQ=255)
                    tCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCcgcgcg                                                     <  2:4441325/76‑1 (MQ=255)
                    tCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCcgcgcg                                                     <  2:5958970/76‑1 (MQ=255)
                     cTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATAACCccgcgc                                                    >  2:5828126/1‑76 (MQ=255)
                       gCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATAACCCCGCGCTg                                                  >  1:3953958/1‑76 (MQ=255)
                         cGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTggg                                                >  1:1617/1‑76 (MQ=255)
                         cGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTggg                                                <  1:3926567/76‑1 (MQ=255)
                          gACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTgggg                                               >  1:3503489/1‑76 (MQ=255)
                          gACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTgggg                                               >  2:1141267/1‑76 (MQ=255)
                          gACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATAACCGCGCGCTgggg                                               >  1:5184665/1‑76 (MQ=255)
                           aCAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGa                                              <  2:5251220/76‑1 (MQ=255)
                           aCAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGa                                              <  2:255177/76‑1 (MQ=255)
                           aCAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGa                                              <  2:4744453/76‑1 (MQ=255)
                            cAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGAt                                             <  1:2809413/76‑1 (MQ=255)
                            cAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGAt                                             <  2:6127246/76‑1 (MQ=255)
                            cAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATAACCGCGCGCTggggga                                             >  2:2731107/1‑74 (MQ=255)
                             aGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATg                                            <  1:1998989/76‑1 (MQ=255)
                               cGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGAgggg                                          >  2:5236961/1‑75 (MQ=255)
                               cGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATggt                                          <  1:5328987/76‑2 (MQ=255)
                               cGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATAACCGCGCGCTGGGGATggg                                          >  2:4093329/1‑75 (MQ=255)
                                gATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGGTGGCa                                         >  2:1133366/1‑76 (MQ=255)
                                gATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCa                                         >  1:4957412/1‑76 (MQ=255)
                                 aTGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAg                                        <  1:558310/76‑1 (MQ=255)
                                     tGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGAGGGCTGCGGt                                    >  2:686149/1‑76 (MQ=255)
                                      gACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGGTGGCAgggggt                                   >  2:1427321/1‑74 (MQ=255)
                                      gACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGtt                                   <  1:2715096/76‑1 (MQ=255)
                                      gACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATAACCGCGCGCTTGGGATGGGAGTgggt                                   >  2:1977874/1‑74 (MQ=255)
                                       aCGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTg                                  <  2:6011137/76‑1 (MQ=255)
                                       aCGAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTg                                  <  2:5418353/76‑1 (MQ=255)
                                         gAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCa                               <  1:3782233/76‑1 (MQ=255)
                                         gAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCa                               <  1:3974152/76‑1 (MQ=255)
                                         gAAAGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCa                               >  1:1369563/1‑76 (MQ=255)
                                          aaaGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCCCGCGCTGGGGATGGTAGCTGTTTCaa                              >  2:1872885/1‑76 (MQ=255)
                                           aaGCAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGGAAt                             >  1:2769538/1‑76 (MQ=255)
                                              cAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATggg                          >  1:4024116/1‑76 (MQ=255)
                                              cAGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATAACCGCGCGCTGGGGATGGCAGTGGGTGGAATtgg                          >  2:3372166/1‑76 (MQ=255)
                                               aGTCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCCCGCGCTGGGGATGGCAGTGGTTGGAATgggg                         >  2:3873543/1‑75 (MQ=255)
                                                 tCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGGTGGCAGTGGTTGCAATgggggg                       >  2:4347150/1‑76 (MQ=255)
                                                 tCTATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGggcgg                       >  1:2803074/1‑76 (MQ=255)
                                                   tATGCTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGGAATgggggggg                     >  2:5553178/1‑74 (MQ=255)
                                                       cTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGGTGGCAGTGGTTGCCAATGgggggggggg                  >  2:5443090/1‑74 (MQ=255)
                                                       cTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGGGGGGggcg                 >  2:5192841/1‑76 (MQ=255)
                                                       cTGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGGTGCAATGGGGggcgggcg                 >  1:3865962/1‑76 (MQ=255)
                                                       cTGATACCGAAGGTTTTGCCAGTGAAATGGATAACCGCGCGCTGGGGATGGCAGTGGTTGGCATGggggggggggc                 >  1:2384212/1‑73 (MQ=255)
                                                       cTGATACCGAAGGTTTTGCCAGTGAAATGGATAACCGCGCGCTGGGGATGGCAGTGGTTGCAATTGGGGGGGggcg                 >  2:1091224/1‑76 (MQ=255)
                                                        tGATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGgtcgt                <  2:2971312/76‑1 (MQ=255)
                                                         gATACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGggcgtc               >  1:4717729/1‑76 (MQ=255)
                                                          aTACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGggcggcgggcggcg              >  1:2336866/1‑76 (MQ=255)
                                                          aTACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGgtcgtcg              <  2:4797738/76‑1 (MQ=255)
                                                          aTACCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGAGGGCAGTGGTTGGAATGGGGGGCGggggtcg              >  2:2836778/1‑76 (MQ=255)
                                                            aCCGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGgtcgtcgtc            <  2:5960662/76‑1 (MQ=255)
                                                             ccGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGGTGCAATGGGCGGGGggcgtcggcc           >  2:183517/1‑73 (MQ=255)
                                                             ccGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGGGGTTGCAATGGGCGGCGGGCGTCGTCa           >  1:910674/1‑76 (MQ=255)
                                                              cGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAg          <  2:4735388/76‑1 (MQ=255)
                                                              cGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAg          <  1:6287136/76‑1 (MQ=255)
                                                              cGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGGCGTCGTCAg          >  1:720057/1‑76 (MQ=255)
                                                              cGAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGGGGTTGCCATGGGCGGGGggcgtcgtcca          >  2:3682340/1‑74 (MQ=255)
                                                               gAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGTGCGGCGGTCGTCGTCAgg         <  2:4807949/76‑1 (MQ=255)
                                                               gAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAgg         <  1:572375/76‑1 (MQ=255)
                                                               gAAGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGGGGTTGCAATGGGCGGCGGCGGGCGTCAgg         >  1:1387109/1‑76 (MQ=255)
                                                                aaGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGGCGTCGCCAGGc        >  2:4534626/1‑76 (MQ=255)
                                                                 aGGTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCa       <  2:6307449/76‑1 (MQ=255)
                                                                  ggTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCAt      <  1:996582/76‑1 (MQ=255)
                                                                  ggTTTTGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCAt      <  1:6158602/76‑1 (MQ=255)
                                                                    ttttGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCt    <  2:2300022/76‑1 (MQ=255)
                                                                    ttttGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGCCGGCCGGCctct    >  1:6147638/1‑76 (MQ=255)
                                                                    ttttGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGGCGTCGTCAGGCATCt    >  1:4743674/1‑76 (MQ=255)
                                                                      ttGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGa  <  2:3326126/76‑1 (MQ=255)
                                                                      ttGCCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGGTGCAATGGGCGGCGGGCGTCCGCAGGCCTCTGa  >  1:1755801/1‑76 (MQ=255)
                                                                         |                                                                         
GTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGA  >  REL606/4608441‑4608586

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: