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
*	REL606	1,329,516	0	C→T	100.0%	154.0	53	H33Y (CAC→TAC)	topA	DNA topoisomerase I
Reads supporting (aligned to +/- strand):  new base (24/28):  ref base (0/0):  total (24/29)

GTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCT  >  REL606/1329484‑1329548                                 |                                 gTAGTGACTACGTGGTGAAATCCAGCGTCGGTtaca                               <  1:4275518/36‑1 (MQ=255)  tAGTGACTTCGTGGTGAAATCCAGCGTCGGTTACAt                              <  1:3001701/36‑1 (MQ=21)   aGTGACTACGTGGTGAAATCCAGCGTCGGTTACATc                             <  1:1289749/36‑1 (MQ=255)   aGTGACTACGTGGTGAAATCCAGCGTCGGTTACATc                             >  1:6112157/1‑36 (MQ=255)    gTGACTACGTGGTGAAATCCAGCGTCGGTTACATcc                            <  1:5566077/36‑1 (MQ=255)    gTGACTACGTGGTGAAATCCAGCGTCGGGTACATcc                            >  1:4377213/1‑36 (MQ=21)     tGACTACGTGGTGAAATCCAGCGTCGGTTACATCcg                           <  1:4656596/36‑1 (MQ=255)      gACTACGTGGTGAAATCCAGCGTCGGTTATATcccc                          >  1:2551201/1‑34 (MQ=21)       aCTACGTGGTGAAATCCAGCGTCGGTTACATCcgcg                         <  1:6664875/36‑1 (MQ=255)         tACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAt                       >  1:961844/1‑36 (MQ=255)         tACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAt                       <  1:4758611/36‑1 (MQ=255)         tACGTGGTGAAATCCAGCGTCGGGTACATTCGCGAt                       >  1:770289/1‑36 (MQ=18)           cGTGGTGTAATCCAGCGTCGGTTACATCCGCGAttt                     <  1:5523678/36‑1 (MQ=25)           cGTGGTGAAATCCAGCGTCGGTTACATCCGCGAttt                     <  1:3299145/36‑1 (MQ=255)           cGTGGTGAAATCCAGCGTCGGTTACATCCGCGAttt                     <  1:7464516/36‑1 (MQ=255)           cGTGGTGAAATCCAGCGTCGGTTACATCCGCGAttt                     <  1:4480628/36‑1 (MQ=255)           cGTGGTGAAATCCAGCGTCGGGTACATCCGCGAttt                     >  1:3021222/1‑36 (MQ=21)            gtggtgATATCCAGCGTTGGTTACATCCGCGATTTg                    <  1:7250509/36‑1 (MQ=21)            gtggtgAAATCCAGCGTCGGTTACATCCGTGAtttc                    >  1:7155767/1‑35 (MQ=21)            gtggtgAAATCCAGCGTCGGTTACATCCGCGCtttt                    >  1:5760894/1‑35 (MQ=21)             tggtgAAAGCCAGCGTCGGTTACATCCGCGATTTGc                   <  1:2899738/36‑1 (MQ=25)              ggtgAAATCCAGCGTCGGTTTCATCCGGGATTTGcc                  >  1:3597597/1‑36 (MQ=18)              ggtgAAATCCAGCGTCGGTTTCATCCGCGGTTTGcc                  >  1:3936195/1‑36 (MQ=18)              ggtgAAATCCAGCGTCGGTTACATCCGCGATTTGcc                  <  1:442916/36‑1 (MQ=255)              ggtgAAATCCAGCGTCGGTTACATCCGCGATTTGcc                  >  1:2350665/1‑36 (MQ=255)               gtgAAATCCAACGTCGGTTACATCCGCGATTTGCCg                 <  1:56434/36‑1 (MQ=25)                tgAAAGCCGGCGGCGGTTACATCCGCGATTTGCCGa                <  1:3021935/36‑1 (MQ=18)                 gAAATCCAGCGTCGGTTACATCCGCGATTTGCCGAc               <  1:2512046/36‑1 (MQ=255)                  aaaTCCAGCGTCGGTTACATCCGCGACTTGCCGAcg              >  1:5560855/1‑35 (MQ=25)                  aaaTCCAGCGTCGGTTACAACCGCGATTTGCCGAcc              >  1:904040/1‑36 (MQ=21)                   aaTCCAGCGTCGGTTACATTCGCGATTTGTCGACCa             >  1:1676873/1‑36 (MQ=21)                   aaTCCAGCGTCGGTTACATCCGCGATTTGCCGACCa             <  1:2885107/36‑1 (MQ=255)                    aTCCAGCTTCGGTTACATCCGCGATTTGCCGACCAg            <  1:4599196/36‑1 (MQ=25)                    aTCCAGCGTCGGTTACATCCGCGATTTTCCGACCAg            >  1:6950766/1‑36 (MQ=25)                    aTCCACCGTCGGTTACATCCGCGATTTGCCGACCAg            <  1:4262965/36‑1 (MQ=25)                     tCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGt           <  1:6959468/36‑1 (MQ=255)                      ccAGCGTCGGTTACATCCGCGATTTTTCGACCAGTg          >  1:4205899/1‑36 (MQ=21)                      ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTg          <  1:502295/36‑1 (MQ=255)                      ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTg          <  1:589067/36‑1 (MQ=255)                      ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTg          <  1:7118790/36‑1 (MQ=255)                      ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTg          <  1:3614789/36‑1 (MQ=255)                      ccAGCGTCGGGTACATCCGCGATTTGCCGACCAGTg          >  1:4991985/1‑36 (MQ=25)                        aGCGTTGGTTACATCCGTGATTTGCCGACCAGTGGc        <  1:5817048/36‑1 (MQ=21)                         gCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCt       <  1:4322930/36‑1 (MQ=255)                          cGTCGGTTACATCCGCGATTTGCCGACCCGTGGccc      >  1:7199168/1‑34 (MQ=25)                          cGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTc      >  1:636767/1‑36 (MQ=255)                          cGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTc      <  1:3778222/36‑1 (MQ=255)                           gTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa     >  1:1446545/1‑36 (MQ=255)                           gTCGGTTACATCCGCGATTTGCCGACCACAGGCTCa     >  1:6560995/1‑36 (MQ=21)                           gTCGGTAACATCCGCGATTTACCGACCAGTGGCTCa     <  1:6597158/36‑1 (MQ=21)                            tCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAg    <  1:576345/36‑1 (MQ=255)                             cGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGc   >  1:5563384/1‑36 (MQ=255)                             cGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGc   >  1:4625832/1‑36 (MQ=255)                              ggTTACATCCGCGATTTGCCGACCAGTGGCTCAGCt  >  1:2072150/1‑36 (MQ=255)                                 |                                 GTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCT  >  REL606/1329484‑1329548 Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 13 ≤  ATCG  < 26 ≤  ATCG  < 32 ≤  ATCG