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% 141.5 4097topADNA topoisomerase I
Reads supporting (aligned to +/- strand):  new base (20/19):  ref base (1/0):  total (21/19)

GGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCT  >  REL606/1329483‑1329548
                                 |                                
ggTAGTGACTACGTGGTGAAATCCAGCGTCGGTtac                                <  1:1838617/36‑1 (MQ=255)
ggTAGTGACTACGTGGTGAAATCCAGCGTCGGTtac                                <  1:6479982/36‑1 (MQ=255)
 gTAGTGACTACGTGGTGAAATCCAGCGTCGGTtaca                               >  1:7351859/1‑36 (MQ=255)
  tAGTGGCTACGTGGTGAAATCCAGCGTCGGTTACAt                              <  1:5369843/36‑1 (MQ=37)
   aGTGACTACGTGGTGAAATCCAGCGTCGGTTACATc                             <  1:4399252/36‑1 (MQ=255)
   aGTGACTACGTGGTGAAATCCAGCGTCGGTTACATc                             <  1:7274882/36‑1 (MQ=255)
    gTGACTACGTGGTGAAATCCAGCGTCGGTTACATcc                            >  1:543160/1‑36 (MQ=255)
    gTGACTACGTGGTGAAATCCAGCGTCGGTTACATcc                            <  1:7447043/36‑1 (MQ=255)
    gTGACTACGTGGTGAAATCCAGCGTCGGTTACATcc                            <  1:2639790/36‑1 (MQ=255)
    gTGACTACGTGGTGAAATCCAGCGTCGGTTAAATcc                            >  1:1528204/1‑36 (MQ=37)
     tGACTACGTGGTGAAATCCAGCGTCGGTCACATCcg                           >  1:3510440/1‑36 (MQ=255)
      gACTACGTGGTGAAATCCAGCGTCGGTTACATCcgg                          >  1:7588676/1‑35 (MQ=255)
      gACTACGTGGTGAAATCCAGCGTCGGTTACATCcgc                          >  1:7916453/1‑36 (MQ=255)
      gACTACGTGGTGAAATCCAGCGTCGGTTACATCcgc                          >  1:843407/1‑36 (MQ=255)
       aCTACGTGGTGAAATCCAGCGTCGGTTACATCcgcg                         >  1:2685863/1‑36 (MQ=255)
       aCTACGTGGTGAAATCCAGCGTCGGCTACCTCcgcg                         >  1:2640538/1‑36 (MQ=21)
         tACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAt                       <  1:2647539/36‑1 (MQ=255)
         tACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAt                       >  1:1161896/1‑36 (MQ=255)
           cGTGGTGAAATCCAGCGTCGGTTACATCCGCGAttt                     >  1:4265186/1‑36 (MQ=255)
             tggtgAAATCCAGCGTCGGTTACATCCGCGATTTGc                   <  1:5220381/36‑1 (MQ=255)
             tggtgAAATCCAGCGTCGGTTACATCCGCGATTTGc                   <  1:305683/36‑1 (MQ=255)
              ggtgAAATCCAGCGTCGGTTACATCCGCGATTTTcc                  <  1:4063918/36‑1 (MQ=37)
              ggtgAAACCCAGCGTCGGTTACATCCGCGATTTGcc                  <  1:7919173/36‑1 (MQ=37)
               gtgAAATCCAGCGTCGGTTACATCCGCGATTTGCCg                 <  1:7248614/36‑1 (MQ=255)
                tgAAATCCAGCGTCGGTTACATCCGCGATTTTCCGa                >  1:6986361/1‑36 (MQ=37)
                tgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGa                >  1:5909900/1‑36 (MQ=255)
                  aaaTCCAGCGTCGGTTACATCCGCGATTTGCCGAcc              >  1:2845317/1‑36 (MQ=255)
                   aaTCCAGCGTCGGTTACATCCGCGATTTGCCGACCa             >  1:6943511/1‑36 (MQ=255)
                    aTCCAGCGTCGGTTACATCCGCGATTTTCCGACCAg            >  1:230380/1‑36 (MQ=37)
                    aTCCAGCGTCGGTTACATCCCCGATTTGCCGACCAg            >  1:7406749/1‑36 (MQ=38)
                     tCCAGCGTCGGTTAAATCCGCGATTTGCCGACCAGt           >  1:5675857/1‑36 (MQ=25)
                      ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTg          <  1:5341366/36‑1 (MQ=255)
                       cAGCGTCGGTTACATCCGCGATTTGCCGACCAGTgg         <  1:2117340/36‑1 (MQ=255)
                       cAGCGTCGGTTACATCCGCGATTTCCCGACCAGTgg         >  1:1026888/1‑36 (MQ=38)
                         tcctCGGTTACATCCGCGATTTGCCGCCCAGTGGCt       <  1:1583289/33‑1 (MQ=25)
                         gCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCt       <  1:1353234/36‑1 (MQ=255)
                          cGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTc      <  1:6798214/36‑1 (MQ=255)
                          cGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTc      <  1:6557229/36‑1 (MQ=255)
                             cGGTTACATCCGCGATTTGCCGACCAGGGGCTCAGc   >  1:352676/1‑36 (MQ=37)
                             cGGTTACATCCGAGATTTGCCGACCAGTTGCTCAGc   >  1:4523098/1‑36 (MQ=25)
                              ggTTACATCCGCGATTTGCCGACCAGTGGCTCAGCt  <  1:6608914/36‑1 (MQ=255)
                              ggTTACATCCGCGATTTGCCGACCAGTGGCTCAGCt  <  1:6127808/36‑1 (MQ=255)
                              ggTTACATCCGCGATTTGCCGACCAGTGGCTCAGCt  >  1:3807387/1‑36 (MQ=255)
                                 |                                
GGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCT  >  REL606/1329483‑1329548

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