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% 146.3 4197topADNA topoisomerase I
Reads supporting (aligned to +/- strand):  new base (17/24):  ref base (0/0):  total (17/24)

GGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCT  >  REL606/1329482‑1329548
                                  |                                
gggTAGTGACTACGTGGTGAAATCCAGCGTCGGTca                                 >  1:27014/1‑36 (MQ=255)
 ggTAGTGACTACGTGGTGAACTCCAGCGTCGGTtac                                >  1:4303437/1‑36 (MQ=37)
 ggTAGTGACTACGTGGTGAAATCCAGCGTCGGTtac                                <  1:1860564/36‑1 (MQ=255)
 ggTAGTGACTACGTGGTGAAATCCAGCGTCGGTtac                                >  1:2550592/1‑36 (MQ=255)
  gTAGTGCCTACGTGGTGAAATCCAGCGTCGGTtaca                               <  1:6252245/36‑1 (MQ=37)
  gTAGTGACTACGTGGTGAAATCCAGCGTCGGTtaca                               <  1:2981406/36‑1 (MQ=255)
   tAGTGACTACGTGGTGAAATCCAGCGTCGGTTACAt                              >  1:5825873/1‑36 (MQ=255)
   tAGTGACTACGTGGTGAAATCCAGCGTCGGTTACAt                              <  1:3387784/36‑1 (MQ=255)
    aGTGGCTCCGTGGTGAAATCCAGCGTCGGTTACATc                             <  1:2892711/36‑1 (MQ=21)
    aGTGACTACGTGGTGAAATCCAGCGTCGGTTACATc                             <  1:2032215/36‑1 (MQ=255)
     gTGATTACGTGGTGAAATCCAGCGTCGGTTACATcc                            <  1:7418947/36‑1 (MQ=37)
     gTGACTACGTGGTGAAATCCAGCGTCGGTTACATcc                            <  1:6747444/36‑1 (MQ=255)
     gTGACTACGTGGTGAAATCCAGCGTCGGTTACATcc                            <  1:2135302/36‑1 (MQ=255)
       gACTACGTGGTGAAATCCAGCGTCGGTTACATCcgc                          <  1:2478742/36‑1 (MQ=255)
         cTACGTGGTGAAATCCAGCGTCGGTTACATCcgcgg                        >  1:4031995/1‑35 (MQ=255)
         cTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGa                        >  1:5150676/1‑36 (MQ=255)
         cTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGa                        <  1:7620701/36‑1 (MQ=255)
          tACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAt                       <  1:3375349/36‑1 (MQ=255)
          tACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAt                       <  1:666115/36‑1 (MQ=255)
          tACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAt                       >  1:729790/1‑36 (MQ=255)
           aCGTGGTGAAATCCAGCGTCGGTTACATCCGCGAtt                      <  1:6461233/36‑1 (MQ=255)
             gtggtgAAATCCAGCGTCGGTTACATCCGCGATTTg                    >  1:6315922/1‑36 (MQ=255)
             gtggtgAAATCCAGCGTCGGTTACATCCGCGATTTg                    >  1:6495344/1‑36 (MQ=255)
              tggtgAAATCCAGCGTCGGTTACATCCGCGATTTGc                   >  1:3266559/1‑36 (MQ=255)
              tggtgAAATCCAGCGTCGGTTACATCCGCGATTTGc                   >  1:4707116/1‑36 (MQ=255)
               ggtgAAATCCAGCGTCGGTTACATCCGCGATTTGcc                  <  1:7184143/36‑1 (MQ=255)
                 tgAAATCCAGCGTCGGTTACATCCGCGCTTTTCCGa                >  1:4119959/1‑36 (MQ=25)
                 tgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGa                <  1:1047572/36‑1 (MQ=255)
                 tgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGa                <  1:2548212/36‑1 (MQ=255)
                    aaTCCAGCGTCGGTTACATCCGCGATTTGCCGACCa             <  1:7452875/36‑1 (MQ=255)
                    aaTCCAGCGTCGGTTACATCCGCGATTTGCCGACCa             >  1:1067880/1‑36 (MQ=255)
                    aaTCCAGCGTCGGTTACATCCGCGATTTGCCGACCa             <  1:6354400/36‑1 (MQ=255)
                       ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTg          >  1:5336081/1‑36 (MQ=255)
                       ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTg          >  1:398651/1‑36 (MQ=255)
                       ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTg          <  1:1645610/36‑1 (MQ=255)
                        cAGCGTCGGTTACATCCGCGATTTGCCGACCAGTgg         >  1:3781760/1‑36 (MQ=255)
                         aGCGTCGGTTACATACGCGATTTTCCGACCAGGggg        >  1:838978/1‑35 (MQ=19)
                           cGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTc      <  1:1195698/36‑1 (MQ=255)
                            gTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa     >  1:6288696/1‑36 (MQ=255)
                            gTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa     <  1:7456129/36‑1 (MQ=255)
                            gTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa     <  1:4935979/36‑1 (MQ=255)
                             tCGGTTACATCCGCGGTTTGCCGCCCAGTGGCTCAg    <  1:1913015/36‑1 (MQ=25)
                             tCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAg    <  1:1056885/36‑1 (MQ=255)
                             tCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAg    <  1:6026340/36‑1 (MQ=255)
                              cGGTTACATCCTCGATTTGCCGACCAGTGGCTCAGc   <  1:7607907/36‑1 (MQ=38)
                              cGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGc   >  1:5839566/1‑36 (MQ=255)
                               ggTTACATCCGCGATTTGCCGACCAGTGGCTCAGCt  >  1:3322622/1‑36 (MQ=255)
                                  |                                
GGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCT  >  REL606/1329482‑1329548

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG