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% 106.2 3097topADNA topoisomerase I
Reads supporting (aligned to +/- strand):  new base (12/18):  ref base (0/0):  total (12/18)

GGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCT  >  REL606/1329483‑1329548
                                 |                                
ggTAGTGACTACGTGGTGAAATCCAGCGTCGGTtac                                <  1:120455/36‑1 (MQ=255)
   aGTGACTACGTGGTGAAATCCAGCGTCGGTTACATc                             <  1:2250567/36‑1 (MQ=255)
    gTGACTACGTGGTGAAATCCAGCGTCGGTTACATcc                            <  1:2270150/36‑1 (MQ=255)
    gTGACTACGTGGTGAAATCCAGCGTCGGTTACATcc                            >  1:2834441/1‑36 (MQ=255)
     tGACTACGTGGTGAAATCCAGCGTCAGTTACATCcg                           <  1:2874732/36‑1 (MQ=37)
      gACTACGTGGTGAAATCCAGCGTCGGTTACATCcgc                          >  1:2963316/1‑36 (MQ=255)
       aCTACGTGGTGAAATCCAGCGTCGGTTACATCcgcg                         >  1:5126374/1‑36 (MQ=255)
       aCTACGTGGTGAAATCCAGCGTCGGTTACATCcgcg                         <  1:4445907/36‑1 (MQ=255)
         tACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAt                       >  1:1662517/1‑36 (MQ=255)
          aCGTGGTGAAATCCAGCGTCGGTTACATCCGCGAtt                      <  1:904450/36‑1 (MQ=255)
           cGTGGTGAAATCCAGGGTCTTTTACATCCGCGAttt                     <  1:4717911/36‑1 (MQ=19)
           cGTGGTGAAATCCAGCGTCGGTTACATCCGCGAttt                     <  1:1317438/36‑1 (MQ=255)
           cGTGGTGAAATCCAGCGTCGGTTACATCCGCGAttt                     <  1:3499397/36‑1 (MQ=255)
              ggtgAAATCCATCGTCGGTTACATCCGCGATTTGcc                  <  1:3305966/36‑1 (MQ=37)
              ggtgAAATCCAGCGTCGGTTACATCCGCGATTTGcc                  <  1:3675053/36‑1 (MQ=255)
                 gAAATCCAGCGTCGGTTAAATCCGCGATTTTTCGAc               >  1:2438457/1‑36 (MQ=16)
                 gAAAACCAGCGTCGGTTACATCCGCGATTTGCCGAc               <  1:2966458/36‑1 (MQ=37)
                  aaaTCCAGCGTCGGTTACATCCGCGACTTGCCGAcc              >  1:5402750/1‑36 (MQ=37)
                     tCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGt           >  1:5293056/1‑36 (MQ=255)
                     cccAGCGTCGGTTACATCCCCGCTTTGCCGACCAGt           <  1:3950125/35‑1 (MQ=25)
                      ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTg          >  1:1003773/1‑36 (MQ=255)
                       cAGCGTCGGTTACATCCGCGATTTGCCGACCAGTgg         <  1:3901073/36‑1 (MQ=255)
                        aGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGc        >  1:4066705/1‑36 (MQ=255)
                         gCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCt       >  1:3374604/1‑36 (MQ=255)
                          catCGGTTACATCCTCGATTTGCCGCCCAGTGGCTc      <  1:5036225/34‑1 (MQ=21)
                          cGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTc      >  1:5012644/1‑36 (MQ=255)
                           gTCGGTTACATCCGCGATTTTCCGACCAGTGGCTCa     >  1:4036047/1‑36 (MQ=37)
                           gTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa     <  1:422495/36‑1 (MQ=255)
                            tCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAg    <  1:3590963/36‑1 (MQ=255)
                              ggTTACATCCGCGATTTGCCGACCAGTGGCTCAGCt  <  1:442742/36‑1 (MQ=255)
                              ggTTACATCCGCGATTTGCCGACCAGTGGCTCAGCt  <  1:5275584/36‑1 (MQ=255)
                                 |                                
GGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCT  >  REL606/1329483‑1329548

Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 18 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG