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%	163.3 / NA	41	H33Y (CAC→TAC)	topA	DNA topoisomerase I
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/24);  total (17/24)

GAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTA  >  REL606/1329444‑1329588                                                                         |                                                                         gAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTtaca                                                                       <  2:947803/76‑1 (MQ=255)     ccccGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCcg                                                                   <  1:981246/76‑1 (MQ=255)     ccccGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCcg                                                                   >  2:827855/1‑76 (MQ=255)     ccccGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCcg                                                                   <  2:1438616/76‑1 (MQ=255)           cAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAttt                                                             <  2:1259156/76‑1 (MQ=255)           cAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAttt                                                             <  2:1056797/76‑1 (MQ=255)              aaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAACCCAGCGTCGGTTACATCCGCGATTTGcc                                                          <  2:212308/76‑1 (MQ=255)                   aaaaCGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCa                                                     >  2:853517/1‑76 (MQ=255)                   aaaaCGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCa                                                     <  1:745436/76‑1 (MQ=255)                           caacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa                                             <  1:901253/76‑1 (MQ=255)                           caacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa                                             <  1:404812/76‑1 (MQ=255)                                 gTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGcc                                       <  1:159499/76‑1 (MQ=255)                                  tATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTGAGCTGCCa                                      <  1:183062/76‑1 (MQ=255)                                  tATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCa                                      >  1:751004/1‑76 (MQ=255)                                   aTCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaa                                     >  1:1546791/1‑76 (MQ=255)                                   aTCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaa                                     >  2:1175619/1‑76 (MQ=255)                                    tCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaa                                    <  2:1399343/76‑1 (MQ=255)                                      tGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaaaa                                  >  1:306081/1‑76 (MQ=255)                                         gTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAagag                               >  2:1565019/1‑76 (MQ=255)                                           aGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTg                             >  2:1518216/1‑76 (MQ=255)                                            gTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCTAAAAAGAGTGc                            >  2:1447069/1‑76 (MQ=255)                                             tGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGcc                           >  1:410391/1‑76 (MQ=255)                                              gACTACGTGGTGAAATCCAGCGTCGGTTACATCTGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCg                          >  1:1096227/1‑76 (MQ=255)                                              gACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCg                          >  1:470163/1‑76 (MQ=255)                                               aCTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGa                         <  1:719388/76‑1 (MQ=255)                                                cTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGAc                        <  1:1272633/76‑1 (MQ=255)                                                 tACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGAct                       <  2:581897/76‑1 (MQ=255)                                                   cGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGActct                     <  1:1194745/76‑1 (MQ=255)                                                   cGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGActct                     <  1:1366309/76‑1 (MQ=255)                                                   cGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGActct                     <  1:295233/76‑1 (MQ=255)                                                    gtggtgAAATCCAGCGTCTGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTa                    <  2:29806/76‑1 (MQ=255)                                                     tggtgAACTCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTAc                   <  2:505268/76‑1 (MQ=255)                                                         gAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTcc               <  1:309843/76‑1 (MQ=255)                                                            aTCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTccacc            <  1:166154/76‑1 (MQ=255)                                                             tCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTccacca           <  2:391314/76‑1 (MQ=255)                                                              ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCCGCTGCCAAAAAGAGCGCCGACTCTACCTccaccac          >  1:300254/1‑75 (MQ=255)                                                               cAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAg         <  2:623857/76‑1 (MQ=255)                                                                 gCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGAc       <  2:93607/76‑1 (MQ=255)                                                                 gCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGAc       <  1:170442/76‑1 (MQ=255)                                                                   gTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACgg     <  1:1550138/76‑1 (MQ=255)                                                                    tCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGc    >  2:201194/1‑76 (MQ=255)                                                                    tCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGc    >  2:157708/1‑76 (MQ=255)                                                                     cGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCt   >  2:1414874/1‑76 (MQ=255)                                                                      ggTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTa  >  2:278499/1‑76 (MQ=255)                                                                      ggTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTa  <  1:343824/76‑1 (MQ=255)                                                                         |                                                                         GAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTA  >  REL606/1329444‑1329588

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