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% 151.7 4097topADNA topoisomerase I
Reads supporting (aligned to +/- strand):  new base (13/27):  ref base (0/0):  total (13/27)

GGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCT  >  REL606/1329483‑1329548
                                 |                                
ggTAGTGACTACGTGGTGAAATCCAGCGTCGGTtac                                <  1:5006799/36‑1 (MQ=255)
ggTAGTGACTACGTGGTGAAATCCAGCGTCGGTtac                                >  1:3050817/1‑36 (MQ=255)
ggTAGTGACTACGTGGTGAAATCCAGCGTCGGTtac                                <  1:2495211/36‑1 (MQ=255)
 gTAGTGACTACGTGGTGAAATCCAGCGTCGGTtaca                               <  1:6425374/36‑1 (MQ=255)
  tAGTGACTACGTGGTGAAATCCAGCGTCGGTTACAt                              >  1:935619/1‑36 (MQ=255)
   aGTGACTACGTGGTGAAATCCAGCGTCGGTTACATc                             >  1:2869315/1‑36 (MQ=255)
    gTGACTACGTGGTGAAATCCAGCGTCGGTTACATcc                            <  1:4686100/36‑1 (MQ=255)
     tGACTACGTGGTGAAATCCAGCGTCGGTTACATCcg                           >  1:1006097/1‑36 (MQ=255)
      gACTACGTGGTGAAACCCAGCGTCGGTTACATCcgc                          <  1:2892509/36‑1 (MQ=37)
       aCTACGTGGTGAAATCCAGCGTCGGTTACATCcgcg                         <  1:4291220/36‑1 (MQ=255)
       aCTACGTGGTGAAATCCAGCGTCGGTTACATCcgcg                         <  1:3641561/36‑1 (MQ=255)
        cTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGa                        <  1:4053838/36‑1 (MQ=255)
         tACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAt                       <  1:1199001/36‑1 (MQ=255)
          aCGTGGTGAAATCCAGCGTCGGTTACATCCTCGAtt                      >  1:4409329/1‑36 (MQ=37)
          aCGTGGTGAAATCCAGCGTCGGTTACATCCGCGAtt                      >  1:4836721/1‑36 (MQ=255)
           cGTGGTGAAATCCAGCGTCGGTTACATCCGCGAttt                     <  1:7795978/36‑1 (MQ=255)
           cGTGGTGAAATCCAGCGTCGGTTACATCCGCGAttt                     >  1:7390734/1‑36 (MQ=255)
           cGTGGTGAAATCCAGCGTCGGTTACATCCGCGAttt                     <  1:6936130/36‑1 (MQ=255)
            gtggtgAAATCCAGCGTCGGTTACATCCGCGATTTg                    <  1:2660296/36‑1 (MQ=255)
             tggtgAAATCCAGCGTCGGTTACATCCGCGAttttc                   >  1:2506700/1‑34 (MQ=37)
             tggtgAAATCCAGCGTCGGTTACATCCGCGATTTGc                   <  1:7742789/36‑1 (MQ=255)
              ggtgAAATCTAGCGTCGGTTACATCCGCGATTTGcc                  <  1:3365775/36‑1 (MQ=37)
              ggtgAAATCCAGCGTCGGTTACATCCGCGATTTGcc                  <  1:619650/36‑1 (MQ=255)
               gtgAAATCCAGCGTCGGTTACATCCGCGATTTTCCg                 >  1:7500639/1‑36 (MQ=37)
               gtgAAATCCAGCGTCGGTTACATCCGCGATTTGCCg                 <  1:6686867/36‑1 (MQ=255)
               gtgAAATCCAGCGTCGGTTACATCCGCGATTTGCCg                 >  1:3732521/1‑36 (MQ=255)
                tgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGa                <  1:7750184/36‑1 (MQ=255)
                tgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGa                >  1:6098661/1‑36 (MQ=255)
                tgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGa                <  1:4700295/36‑1 (MQ=255)
                 gAAATCCAGCGTCGGTTACATCCGCGATTTGCCGAc               >  1:2554995/1‑36 (MQ=255)
                  aaaTCCAGCGTCGGTTACATCCGCGATTTGCCGAcc              <  1:6164691/36‑1 (MQ=255)
                   aaTCCAGCGTCGGTTTCATCCGCGATTTGCCGACCa             >  1:5330545/1‑36 (MQ=37)
                     tCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGt           <  1:2623903/36‑1 (MQ=255)
                      ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTg          >  1:6039866/1‑36 (MQ=255)
                      ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTg          <  1:1421929/36‑1 (MQ=255)
                      ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTg          <  1:4223566/36‑1 (MQ=255)
                       cAGCGTCGGTTACATCCGCGATTTGCCGACCAGTgg         <  1:3004074/36‑1 (MQ=255)
                         gCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCt       <  1:4326642/36‑1 (MQ=255)
                          cGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTc      <  1:3282079/36‑1 (MQ=255)
                           gTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa     <  1:1943290/36‑1 (MQ=255)
                            tCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAg    <  1:1253430/36‑1 (MQ=255)
                             cGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGc   <  1:1775359/36‑1 (MQ=255)
                              ggTTACATCCGCGATTTGCCGACCAGTGGCTCAGCt  <  1:6384120/36‑1 (MQ=255)
                              ggTTACATCCGCGATTTGCCGACCAGTGGCTCAGCt  <  1:7409082/36‑1 (MQ=255)
                                 |                                
GGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCT  >  REL606/1329483‑1329548

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 26 ≤  ATCG  < 34 ≤  ATCG  < 37 ≤  ATCG