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

GAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTA  >  REL606/1329444‑1329588                                                                         |                                                                         gAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTTGTGACTACGTGGTGAAATCCAGCGTCGGTtaca                                                                       >  1:949288/1‑76 (MQ=255)   gTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATc                                                                     >  2:870736/1‑76 (MQ=255)   gTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATc                                                                     <  2:489482/76‑1 (MQ=255)     ccccGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCcg                                                                   <  1:264721/76‑1 (MQ=255)       ccGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCcgcg                                                                 <  2:911336/76‑1 (MQ=255)             aaaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCNGTTACATCCGCGATTTGc                                                           >  1:941449/1‑76 (MQ=255)                gCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAACTCAGCGTCGGTTACATCCGCGATTTGCCGa                                                        <  1:1169979/76‑1 (MQ=255)                  cAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGAcc                                                      <  1:360590/76‑1 (MQ=255)                    aaaCGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAg                                                    >  2:616520/1‑76 (MQ=255)                           naacaaGTATCTGGGTAGTGANTACGTGNTGAAAACCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa                                             <  1:1005556/75‑1 (MQ=255)                           caacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa                                             <  2:253192/76‑1 (MQ=255)                           caacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa                                             >  1:680818/1‑76 (MQ=255)                           caacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa                                             <  2:486315/76‑1 (MQ=255)                                aGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGc                                        >  2:423136/1‑76 (MQ=255)                                    tCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaa                                    >  1:29808/1‑76 (MQ=255)                                     cTGGGTAGTGCCTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaaa                                   <  1:65569/76‑1 (MQ=255)                                     cTGGGTAGTGCCTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaaa                                   >  1:577101/1‑76 (MQ=255)                                      tGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaaaa                                  >  1:733531/1‑76 (MQ=255)                                      tGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaaaa                                  >  2:1185625/1‑76 (MQ=255)                                        ggcaGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAaga                                <  1:563066/73‑1 (MQ=255)                                          tAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGt                              <  1:421878/76‑1 (MQ=255)                                              gACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGGTTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCg                          <  1:833553/76‑1 (MQ=255)                                              gACTACGTAGTGAAATCCAGCGTCGGTTACATCCGCGATTTGNCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCg                          >  1:957846/1‑76 (MQ=255)                                                ctccGTGGTGAAATCCAGCGTCGGTTNCATCCGTGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCAGAc                        >  2:1026787/4‑76 (MQ=255)                                                    gtggtgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTa                    <  1:1158221/76‑1 (MQ=255)                                                          aaaTCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTcca              >  1:189052/1‑76 (MQ=255)                                                          aaaTCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGAGCCGACTCTACCTcca              >  2:793534/1‑76 (MQ=255)                                                                     cGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCt   <  1:905806/76‑1 (MQ=255)                                                                      ggTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTa  >  1:825853/1‑76 (MQ=255)                                                                      ggTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTa  <  1:146665/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: ‑