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

GGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGAC  >  REL606/1329452‑1329583                                                                 |                                                                    ggCAAAATCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAt                                                          <  2:635690/76‑1 (MQ=255)      aaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGcg                                                     >  2:659127/1‑75 (MQ=255)       aGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCg                                                    >  1:513309/1‑76 (MQ=255)              aCGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTg                                             >  1:573558/1‑76 (MQ=255)                      caaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCt                                     <  1:498647/76‑1 (MQ=255)                       aaGTATCTGGGTAGTGACTACGTGGCGAAATACAGCGTCGNNTACATCCGCGNNTGGCCGACCAGTGGATCAGCTg                                    >  1:1053741/1‑76 (MQ=255)                         gTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGcc                                  <  2:440416/76‑1 (MQ=255)                         gTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGcc                                  <  2:240841/76‑1 (MQ=255)                           aTCTGGGAAGTGACTACGTGCTNNAATCCAGCGTNNCTTACATCTGCGATGTGCCGACCAGTGGCTCAGCTGCCaa                                <  1:1045034/76‑1 (MQ=255)                            tCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaa                               <  2:787193/76‑1 (MQ=255)                              tGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaaaa                             >  2:3902/1‑76 (MQ=255)                                 gTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAagag                          <  2:201531/76‑1 (MQ=255)                                  tAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGt                         >  2:442417/1‑76 (MQ=255)                                        cTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGAc                   >  1:370490/1‑76 (MQ=255)                                         tACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGAct                  >  2:255084/1‑76 (MQ=255)                                         tACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGAct                  >  1:1073714/1‑76 (MQ=255)                                           cGTGGTGAAAGCCAGGGTCGTTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGActct                <  1:343458/76‑1 (MQ=255)                                                 gAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTcc          <  2:195568/76‑1 (MQ=255)                                                 gAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTcc          >  2:323839/1‑76 (MQ=255)                                                         gCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGAc  <  1:242590/76‑1 (MQ=255)                                                                 |                                                                    GGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGAC  >  REL606/1329452‑1329583

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