Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,329,516 C→T 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% 170.3 5797topADNA topoisomerase I
Reads supporting (aligned to +/- strand):  new base (27/29):  ref base (0/0):  total (27/30)

GTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCT  >  REL606/1329484‑1329548
                                |                                
gTAGTGACTACGTGGTGAAATCCAGCGTCGGTtaca                               <  2:536308/36‑1 (MQ=255)
 tAGTGACTTCGTGGTGAAATCCAGCGTCGGTTACAt                              <  1:3001701/36‑1 (MQ=37)
  aGTGACTACGTGGTGAAATCCAGCGTCGGTTACATc                             <  1:1289749/36‑1 (MQ=255)
  aGTGACTACGTGGTGAAATCCAGCGTCGGTTACATc                             >  2:2372947/1‑36 (MQ=255)
   gTGACTACGTGGTGAAATCCAGCGTCGGTTACATcc                            <  2:1826867/36‑1 (MQ=255)
   gTGACTACGTGGTGAAATCCAGCGTCGGGTACATcc                            >  2:638003/1‑36 (MQ=37)
    tGACTACGTGGTGAAATCCAGCGTCGGTTACATCcg                           <  2:917386/36‑1 (MQ=255)
     gACTACGTGGTGAAATCCAGCGTCGGTTATATcccc                          >  1:2551201/1‑34 (MQ=25)
      acgaCGTGGTGATATCCGGCGTCGGTTACATCcgcg                         <  1:3459140/33‑1 (MQ=19)
      aCTACGTGGTGAAATCCAGCGTCGGTTACATCcgcg                         <  2:2925665/36‑1 (MQ=255)
        tACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAt                       <  2:1019401/36‑1 (MQ=255)
        tACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAt                       >  1:961844/1‑36 (MQ=255)
        tACGTGGTGAAATCCAGCGTCGGGTACATTCGCGAt                       >  1:770289/1‑36 (MQ=21)
          cGTGGTGTAATCCAGCGTCGGTTACATCCGCGAttt                     <  2:1784468/36‑1 (MQ=37)
          cGTGGTGAAATCCAGCGTCGGTTACATCCGCGAttt                     <  2:741418/36‑1 (MQ=255)
          cGTGGTGAAATCCAGCGTCGGTTACATCCGCGAttt                     <  1:3299145/36‑1 (MQ=255)
          cGTGGTGAAATCCAGCGTCGGTTACATCCGCGAttt                     <  2:3725306/36‑1 (MQ=255)
          cGTGGTGAAATCCAGCGTCGGGTACATCCGCGAttt                     >  1:3021222/1‑36 (MQ=37)
           gtggtgATATCCAGCGTTGGTTACATCCGCGATTTg                    <  2:3511299/36‑1 (MQ=25)
           gtggtgAAATCCAGCGTCGGTTACATCCGTGAtttc                    >  2:3416557/1‑35 (MQ=25)
           gtggtgAAATCCAGCGTCGGTTACATCCGCGCtttt                    >  2:2021684/1‑35 (MQ=25)
            tggtgAAAGCCAGCGTCGGTTACATCCGCGATTTGc                   <  1:2899738/36‑1 (MQ=37)
             ggtgAAATCCAGCGTCGGTTTCATCCGGGATTTGcc                  >  1:3597597/1‑36 (MQ=21)
             ggtgAAATCCAGCGTCGGTTTCATCCGCGGTTTGcc                  >  2:196985/1‑36 (MQ=21)
             ggtgAAATCCAGCGTCGGTTACATCCGCGATTTGcc                  <  1:442916/36‑1 (MQ=255)
             ggtgAAATCCAGCGTCGGTTACATCCGCGATTTGcc                  >  1:2350665/1‑36 (MQ=255)
              gtgAAATCCAACGTCGGTTACATCCGCGATTTGCCg                 <  1:56434/36‑1 (MQ=37)
               tgAAAGCCGGCGGCGGTTACATCCGCGATTTGCCGa                <  1:3021935/36‑1 (MQ=19)
                gAAATCCAGCGTCGGTTACATCCGCGATTTGCCGAc               <  1:2512046/36‑1 (MQ=255)
                 aaaTCCAGCGTCGGTTACATCCGCGACTTGCCGAcg              >  2:1821645/1‑35 (MQ=37)
                 aaaTCCAGCGTCGGTTACAACCGCGATTTGCCGAcc              >  1:904040/1‑36 (MQ=37)
                  aaTCCAGCGTCGGTTACATTCGCGATTTGTCGACCa             >  1:1676873/1‑36 (MQ=25)
                  aaTCCAGCGTCGGTTACATCCGCGATTTGCCGACCa             <  1:2885107/36‑1 (MQ=255)
                   aTCCAGCTTCGGTTACATCCGCGATTTGCCGACCAg            <  2:859986/36‑1 (MQ=37)
                   aTCCAGCGTCGGTTACATCCGCGATTTTCCGACCAg            >  2:3211556/1‑36 (MQ=37)
                   aTCCACCGTCGGTTACATCCGCGATTTGCCGACCAg            <  2:523755/36‑1 (MQ=37)
                    tCCAGCGTCGGTTACATTCGTGATTTTCCTAGCAGt           >  1:614753/1‑36 (MQ=12)
                    tCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGt           <  2:3220258/36‑1 (MQ=255)
                     ccAGCGTCGGTTACATCCGCGGTTTTCCGACCTGTg          >  2:2201847/1‑36 (MQ=16)
                     ccAGCGTCGGTTACATCCGCGATTTTTCGACCAGTg          >  2:466689/1‑36 (MQ=25)
                     ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTg          <  2:3379580/36‑1 (MQ=255)
                     ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTg          <  1:589067/36‑1 (MQ=255)
                     ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTg          <  1:502295/36‑1 (MQ=255)
                     ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTg          <  1:3614789/36‑1 (MQ=255)
                     ccAGCGTCGGGTACATCCGCGATTTGCCGACCAGTg          >  2:1252775/1‑36 (MQ=37)
                       aGCGTTGGTTACATCCGTGATTTGCCGACCAGTGGc        <  2:2077838/36‑1 (MQ=25)
                        gCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCt       <  2:583720/36‑1 (MQ=255)
                         cGTCGGTTACATCCGCGATTTGCCGACCCGTGGccc      >  2:3459958/1‑34 (MQ=37)
                         cGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTc      <  2:39012/36‑1 (MQ=255)
                         cGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTc      >  1:636767/1‑36 (MQ=255)
                         cGTCGGTTACATCCGCGAATTTCCGACCAGCGACTc      >  2:129392/1‑36 (MQ=16)
                          gTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa     >  1:1446545/1‑36 (MQ=255)
                          gTCGGTTACATCCGCGATTTGCCGACCACAGGCTCa     >  2:2821785/1‑36 (MQ=25)
                          gTCGGTAACATCCGCGATTTACCGACCAGTGGCTCa     <  2:2857948/36‑1 (MQ=37)
                           tCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAg    <  1:576345/36‑1 (MQ=255)
                            cGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGc   >  2:886622/1‑36 (MQ=255)
                            cGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGc   >  2:1824174/1‑36 (MQ=255)
                             ggTTACATCCGCGATTTGCCGACCAGTGGCTCAGCt  >  1:2072150/1‑36 (MQ=255)
                                |                                
GTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCT  >  REL606/1329484‑1329548

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 12 ≤  ATCG  < 25 ≤  ATCG  < 32 ≤  ATCG