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

GCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCT  >  REL606/1329453‑1329587                                                                |                                                                        gCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAtt                                                             <  2:1072638/76‑1 (MQ=255)    aaaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGc                                                          >  2:1006260/1‑76 (MQ=255)    aaaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGc                                                          >  2:718602/1‑76 (MQ=255)                aTCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCt                                              >  1:970885/1‑76 (MQ=255)                 tCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTc                                             <  2:1392995/76‑1 (MQ=255)                   aacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAg                                           <  1:922482/76‑1 (MQ=255)                        gTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGcc                                      <  2:72163/76‑1 (MQ=255)                          aTCTGGGTCGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTGAGCTGCCaa                                    >  1:45581/1‑76 (MQ=255)                              gggTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAag                                >  1:877194/1‑76 (MQ=255)                                 tAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGt                             >  1:558482/1‑76 (MQ=255)                                            tggtgAAATCCAGTGTCGGTTACATCCGCGATTNGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTAc                  >  2:1208261/1‑76 (MQ=255)                                                gAAATCCANCATCGGTTACNTCCGCGATTTGCCGACCAGTGGNTCAGCTGCCAAAAAGAGTGCCGACTCTACCTcc              <  2:1185570/76‑1 (MQ=255)                                                    tCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTccacca          >  2:1002149/1‑76 (MQ=255)                                                    tCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTccacca          <  2:637091/76‑1 (MQ=255)                                                      cAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAg        <  2:605786/76‑1 (MQ=255)                                                       aGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGa       <  1:62138/76‑1 (MQ=255)                                                       aGCATCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGa       <  2:1365028/76‑1 (MQ=255)                                                          gTCGGTTACATCCGCGATTTGCNGNCCAGAGGCTCNGCTGCCAAAAAGAGNGCCGACTCTGCCTCCGCCAAGACgg    <  2:23729/76‑1 (MQ=255)                                                           tCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGc   >  1:289050/1‑76 (MQ=255)                                                            cGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCt  >  1:264392/1‑76 (MQ=255)                                                            cGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCt  <  2:753158/76‑1 (MQ=255)                                                                |                                                                        GCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCT  >  REL606/1329453‑1329587

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