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

TGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCT  >  REL606/1329443‑1329587                                                                          |                                                                        tGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTtac                                                                       <  1:797495/76‑1 (MQ=255) tGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTaac                                                                       >  1:546008/1‑76 (MQ=255)         cGGCAAAAGCCAAAACGATCATCAAGTATCTNGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGa                                                               >  1:102147/1‑76 (MQ=255)           gCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGAGTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAtt                                                             <  2:350775/76‑1 (MQ=255)            cAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAttt                                                            <  2:1137944/76‑1 (MQ=255)              aaaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGc                                                          >  2:148429/1‑76 (MQ=255)                 gCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGa                                                       >  2:344620/1‑76 (MQ=255)                    aaaaCGATCAACAAGTATCTGGGTNGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCa                                                    >  2:201925/1‑76 (MQ=255)                    aaaaCGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCTGACCa                                                    >  2:939373/1‑76 (MQ=255)                     aaaCGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAg                                                   <  2:739389/76‑1 (MQ=255)                       aCGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTg                                                 <  2:426044/76‑1 (MQ=255)                          aTCAACAAGTATCTGGGTCGTGACTACGTGGTGAAATCCAGCGTCGGATACATCCGCGATTTGCCGACCAGTGGCt                                              >  2:1051459/1‑76 (MQ=255)                          aTCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCt                                              <  1:158270/76‑1 (MQ=255)                          aTCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCt                                              <  2:641934/76‑1 (MQ=255)                                  gTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGcc                                      >  1:605532/1‑76 (MQ=255)                                   tATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCa                                     >  1:1144451/1‑76 (MQ=255)                                    aTCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaa                                    >  2:366755/1‑76 (MQ=255)                                    aTCTGGGTAGTGACTACGTGGTGAAAGCCAGCGNCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaa                                    <  1:929296/76‑1 (MQ=255)                                      cTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaaa                                  >  1:330624/1‑76 (MQ=255)                                          gTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAagag                              >  1:284810/1‑76 (MQ=255)                                                         tgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTc               <  1:68893/76‑1 (MQ=255)                                                                 aGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAGAGAGTGCCGACTCTACCTCCACCACGa       >  1:550852/1‑76 (MQ=255)                                                                  gCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGAc      <  2:471195/76‑1 (MQ=255)                                                                  gCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGAc      <  2:631843/76‑1 (MQ=255)                                                                      cGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGCGTGCCGACTCTACCTCCACCAAGACGGCt  <  1:1135284/76‑1 (MQ=255)                                                                          |                                                                        TGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCT  >  REL606/1329443‑1329587

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