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

CCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTA  >  REL606/1329448‑1329588                                                                     |                                                                         ccccGGCAAAAGCNAAAACGATCAACAAGTATCTGGGTAGTGACTGCGTGGTGAAATCCAGCGTCGGTTACATCcg                                                                   <  2:1273435/76‑1 (MQ=255)      gCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAtt                                                              <  1:1056139/76‑1 (MQ=255)        aaaaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTg                                                            <  1:141619/76‑1 (MQ=255)             ccAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGAc                                                       <  2:540555/76‑1 (MQ=255)               aaaaCGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCa                                                     <  1:253756/76‑1 (MQ=255)                   cGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTgg                                                 <  1:1134895/76‑1 (MQ=255)                       caacaaGTATCTGGGTAGTGACTACGTGGTGAANTCCAGCGTCGGTTACATCCGCGNTTTGCCGACCNGTGGCTCa                                             >  2:1168571/1‑76 (MQ=255)                       caacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa                                             <  1:651492/76‑1 (MQ=255)                          caaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCt                                          >  1:636972/1‑76 (MQ=255)                            aGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGc                                        <  1:898925/76‑1 (MQ=255)                              tATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCa                                      <  1:253109/76‑1 (MQ=255)                                tCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaa                                    >  2:704375/1‑76 (MQ=255)                                  tGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaaaa                                  >  2:894179/1‑76 (MQ=255)                                  tGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaaaa                                  >  1:530390/1‑76 (MQ=255)                                  tGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGCGGCTCAGCTGCCcaaaa                                  >  2:552843/1‑76 (MQ=255)                                     gTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAagag                               <  1:254091/76‑1 (MQ=255)                                      tAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGt                              <  1:846603/76‑1 (MQ=255)                                          gACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCg                          >  1:475711/1‑76 (MQ=255)                                           aCTACGTGNTGAAATCCAGCGTCGGTTACATCCGCGATTGGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGa                         <  1:1225625/76‑1 (MQ=255)                                           aCTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGa                         >  2:446296/1‑76 (MQ=255)                                             tACGTGGTGAAATCCAGGGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGAct                       <  1:149045/76‑1 (MQ=255)                                             tACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGAct                       >  1:196430/1‑76 (MQ=255)                                             tACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGAct                       <  2:347822/76‑1 (MQ=255)                                              aCGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGActc                      >  1:527464/1‑76 (MQ=255)                                                      aaaTCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTcca              <  1:143957/76‑1 (MQ=255)                                                      aaaTCCAGCGTCGGTTACATCCGCGATATGCCGACCAGTGGCTCNGCNGCCAAAAAGAGTGCNGACTCTACCTcca              >  2:1261111/1‑76 (MQ=255)                                                       aaTCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTccac             <  2:1216430/76‑1 (MQ=255)                                                       aaTCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTccac             >  2:257713/1‑76 (MQ=255)                                                          ccAGCGTCGGTTACATCGGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCACCACCaa          >  1:73876/1‑76 (MQ=255)                                                                tCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGc    <  1:577941/76‑1 (MQ=255)                                                                  ggTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTa  <  2:1115452/76‑1 (MQ=255)                                                                  ggTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTa  >  1:367114/1‑76 (MQ=255)                                                                     |                                                                         CCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTA  >  REL606/1329448‑1329588

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