Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,329,516 C→T 100% H33Y (CAC→TAC)  topA → DNA topoisomerase I

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,329,5160C→T100.0% 181.1 4997topADNA topoisomerase I
Reads supporting (aligned to +/- strand):  new base (18/31):  ref base (0/0):  total (18/31)

GGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCT  >  REL606/1329483‑1329548
                                 |                                
ggTAGTTACTACGTGGTGAAAACCCGCTTCGGTtac                                <  1:437577/36‑1 (MQ=12)
 gTAGTGACTACGTGGTGAAATCCAGCGTCGGTtaca                               >  1:2547732/1‑36 (MQ=255)
  tAGTGACTACGTGGTGAAATCCAGCGTCGGTTACAt                              <  1:6992968/36‑1 (MQ=255)
  tAGTGACTACGTGGTGAAATCCAGCGTCGGTTACAt                              <  1:2645915/36‑1 (MQ=255)
  tAGTGACTACGTGGTGAAATCCAGCGTCGGTTACAt                              <  1:2061012/36‑1 (MQ=255)
    gTGACTACGTGGTGAAATCCAGCGTCGGTTACATcc                            >  1:504603/1‑36 (MQ=255)
    gTGACTACGTGGTGAAATCCAGCGTCGGTTACATcc                            <  1:1852079/36‑1 (MQ=255)
    gTGACTACGTGGTGAAATCCAGCGTCGGTTACATcc                            >  1:1289883/1‑36 (MQ=255)
      gACTACGTGGTGAAATCCAGCGTCGGTTACATCcgc                          <  1:1501207/36‑1 (MQ=255)
       aCTACGTGGTGAAATCCAGCGTCGGTTACATCcgcg                         >  1:3258546/1‑36 (MQ=255)
        cTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGa                        <  1:6422146/36‑1 (MQ=255)
        cTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGa                        <  1:1587913/36‑1 (MQ=255)
        cTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGa                        <  1:2944211/36‑1 (MQ=255)
         tACGTGGTGAAATCCAGCGTCGGTTACATCCTCGAt                       <  1:4285566/36‑1 (MQ=37)
         tACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAt                       <  1:6852632/36‑1 (MQ=255)
         tACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAt                       >  1:4318076/1‑36 (MQ=255)
         tACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAt                       <  1:3257905/36‑1 (MQ=255)
          aCGTGGTGAAATCCAGCGTCGGTTACATCCGCGAtt                      >  1:5740563/1‑36 (MQ=255)
          aCGTGGTGAAATCCAGCGTCGGTTACATCCGCGAtt                      >  1:4763019/1‑36 (MQ=255)
          aCGTGGTAAAATCCAGCGTCGGTTACATCCGCGAtt                      >  1:3968323/1‑36 (MQ=37)
           cGTGTTGAAATCCAGCGTCGGTTACATCCGCGAttt                     <  1:7444143/36‑1 (MQ=37)
           cGTGGTGAAATCCAGCGTCGGTTACATCCGCGAttt                     <  1:7234702/36‑1 (MQ=255)
           cGTGGTGAAATCCAGCGTCGGTTACATCCGCGAttt                     <  1:699251/36‑1 (MQ=255)
             tggtgAAATCCAGCGTCGGTTACATCCGCGATTTGc                   <  1:5270537/36‑1 (MQ=255)
             tggtgAAATCCAGCGTCGATTACATCCGCGATTTGc                   <  1:4621165/36‑1 (MQ=37)
              ggtgAAATCCAGCGTCGGTTACATCCGCGATTTGcc                  <  1:4306499/36‑1 (MQ=255)
              ggtgAAATCCAGCGTCGGTTACATCCGCGATTTGcc                  <  1:1350335/36‑1 (MQ=255)
              ggtgAAATCCAGCGTCGGTTACATCCGCGATTTGcc                  >  1:2088478/1‑36 (MQ=255)
               gtgAAATCCAGCGTCGGTTACATCCGCGATTTTCCg                 >  1:6862770/1‑36 (MQ=37)
                  aaaTCCAGCGTCGGTTACATCCGCGATTTGCCGAcc              >  1:4035341/1‑36 (MQ=255)
                  aaaTCCAGCGTCGGTTACATCCGCGATTTGCCGAcc              <  1:7393825/36‑1 (MQ=255)
                  aaaTCCAGCGTCGGTTACATCCGCGATTTGCCGAcc              >  1:2422873/1‑36 (MQ=255)
                   aaTCCAGCGTCGGTTACATCCGCGATTTGCCGACCa             <  1:953933/36‑1 (MQ=255)
                    aTCCAGCGTCGGTTACATCCGTGATTTGCCGACCAg            >  1:772287/1‑36 (MQ=37)
                     tCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGt           >  1:1742216/1‑36 (MQ=255)
                     tCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGt           <  1:2356089/36‑1 (MQ=255)
                     tCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGt           <  1:1539793/36‑1 (MQ=255)
                       cAGCGTCGGTTACATCCGCGATTTGCCGACCAGTgg         >  1:753103/1‑36 (MQ=255)
                        aGCGTCGGTTACATCCGCGATTTTCCGACCAGTGGc        >  1:784254/1‑36 (MQ=38)
                        aGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGc        <  1:4797591/36‑1 (MQ=255)
                           gTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa     <  1:3919500/36‑1 (MQ=255)
                            tCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAg    <  1:5380252/36‑1 (MQ=255)
                            tCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAg    >  1:1951653/1‑36 (MQ=255)
                            tCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAg    <  1:3031025/36‑1 (MQ=255)
                            tCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAg    <  1:891021/36‑1 (MQ=255)
                             cGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGc   <  1:6226473/36‑1 (MQ=255)
                             cGGTTACATCCGCGATTTGCCGACCAGGGGCTCAGc   >  1:5068161/1‑36 (MQ=37)
                              ggTTACATCCGCGATTTGCCGACCAGTGGCTCAGCt  <  1:6519930/36‑1 (MQ=255)
                              ggTTACATCCGCGATTTGCCGACCAGTGGCTCAGCt  >  1:4784908/1‑36 (MQ=255)
                              ggTTACATCCGCGATTTGCCGACCAGTGGCTCAGCt  <  1:4685911/36‑1 (MQ=255)
                              ggTTACATCCGCGATTTGCCGACCAGTGGCTCAGCt  <  1:2444566/36‑1 (MQ=255)
                                 |                                
GGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCT  >  REL606/1329483‑1329548

Base quality scores:  ATCG  < 3 ≤  ATCG  < 12 ≤  ATCG  < 23 ≤  ATCG  < 31 ≤  ATCG  < 37 ≤  ATCG