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		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	1,329,516	0	C	T	100.0%	114.4 / NA	32	H33Y (CAC→TAC)	topA	DNA topoisomerase I
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/18);  total (14/18)

GAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCT  >  REL606/1329444‑1329587                                                                         |                                                                        gAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTtaca                                                                      >  1:238807/1‑76 (MQ=255)  aGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACAt                                                                     <  2:757396/76‑1 (MQ=255)   gTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTTCGTGGTGAAATCCAGCGTCGGTTACATc                                                                    <  2:1004143/76‑1 (MQ=255)    tCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATcc                                                                   >  1:440718/1‑76 (MQ=255)     ccccGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCcg                                                                  <  1:57524/76‑1 (MQ=255)         ggCAAGAGCCAAGACGGTCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAt                                                              <  2:1221313/76‑1 (MQ=255)         ggCAAAAGCCAAAACGATCAACACGTATCTGGGTAGTGACTACGTGGGGAAATCCAGCGTCGGTTACATCCGCGAt                                                              >  1:402158/1‑76 (MQ=255)         ggCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAt                                                              <  2:702958/76‑1 (MQ=255)          gCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCGAGCGTCGGTTACATCCGCGAtt                                                             <  2:982102/76‑1 (MQ=255)          gCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAtt                                                             <  1:241527/76‑1 (MQ=255)            aaaaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTg                                                           <  2:910112/76‑1 (MQ=255)                gCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGa                                                       <  2:517116/76‑1 (MQ=255)                 ccAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGAc                                                      >  1:1274420/1‑76 (MQ=255)                      aCGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTg                                                 >  2:479809/1‑76 (MQ=255)                        gATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGc                                               <  2:946128/76‑1 (MQ=255)                           caacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa                                            <  1:1477416/76‑1 (MQ=255)                             acaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTAAAACCGCGATTTGGCGGCCAGTGGCTCAGc                                          >  1:947964/1‑76 (MQ=255)                               aaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTg                                        >  2:1438567/1‑76 (MQ=255)                                aGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGc                                       <  2:188408/76‑1 (MQ=255)                                aGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACAACCGCGATTTGCCGACCCGTGGCTCCGCTGc                                       >  2:927644/1‑76 (MQ=255)                                 gTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGcc                                      >  2:1579333/1‑76 (MQ=255)                                  tATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCa                                     >  1:347401/1‑76 (MQ=255)                                   aTCTGGGTAGTGACAACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaa                                    >  1:853806/1‑76 (MQ=255)                                          tAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGt                             >  1:754934/1‑76 (MQ=255)                                            gTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGc                           <  2:982256/76‑1 (MQ=255)                                            gTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGc                           >  1:1172392/1‑76 (MQ=255)                                               aCTGCGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGa                        <  1:1049840/76‑1 (MQ=255)                                                  ncGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCCGCCAAAAAGAGTGCCGActc                     >  1:1222994/2‑76 (MQ=255)                                                   cGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGActct                    >  1:1400874/1‑76 (MQ=255)                                                    gtggtgAAATCCAGCGTCGGTTACGTCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTa                   <  2:1072287/76‑1 (MQ=255)                                                    gtggtgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTa                   <  2:402158/76‑1 (MQ=255)                                                          aaaTCCAGCGTCGTTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTcca             <  1:607574/76‑1 (MQ=255)                                                          aaaTCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGGGCCGACTCTACCTcca             >  1:185897/1‑76 (MQ=255)                                                             tCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTccacca          >  2:574733/1‑76 (MQ=255)                                                              ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCaa         <  1:1587710/76‑1 (MQ=255)                                                               cAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAg        <  1:1047374/76‑1 (MQ=255)                                                                     cGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCt  >  2:1513744/1‑76 (MQ=255)                                                                     cGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCt  <  2:419161/76‑1 (MQ=255)                                                                         |                                                                        GAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCT  >  REL606/1329444‑1329587

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑