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

TGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCT  >  REL606/1329443‑1329587                                                                          |                                                                        tGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTtac                                                                       >  1:672149/1‑76 (MQ=255)  gAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTtaca                                                                      >  1:625235/1‑76 (MQ=255)     tCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATcc                                                                   >  2:354182/1‑76 (MQ=255)      ccccGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCcg                                                                  <  1:139533/76‑1 (MQ=255)         cGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGa                                                               >  1:793712/1‑76 (MQ=255)              aaaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGc                                                          >  2:1005615/1‑76 (MQ=255)              aaaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGc                                                          <  2:565391/76‑1 (MQ=255)              aaaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGc                                                          >  2:719795/1‑76 (MQ=255)               aaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGcc                                                         >  1:744885/1‑76 (MQ=255)                  ccAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGNTGCAATCCAGCGTCGGTTACATCCGCGATTTGCCGAc                                                      <  1:1180420/76‑1 (MQ=255)                   cAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGAcc                                                     >  2:408611/1‑76 (MQ=255)                     aaaCGATCAACNAGNATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATGCGCGATTAGCCGACCAg                                                   <  2:193505/76‑1 (MQ=255)                          aTCGACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCt                                              >  1:904526/1‑76 (MQ=255)                          aTCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCt                                              >  1:904545/1‑76 (MQ=255)                           tCAACAAGTATCTGGGTAGTTACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTc                                             <  2:498737/76‑1 (MQ=255)                             aacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAg                                           >  2:859299/1‑76 (MQ=255)                             aacaaGTACCTGGGNAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAg                                           >  1:82913/1‑76 (MQ=255)                                  gTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGAGCAGTGGCTCAGCTGcc                                      <  1:841468/76‑1 (MQ=255)                                    aTCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaa                                    >  1:412107/1‑76 (MQ=255)                                    aTCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaa                                    <  1:240139/76‑1 (MQ=255)                                      cTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaaa                                  <  2:1165999/76‑1 (MQ=255)                                        gggTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAag                                >  2:822433/1‑76 (MQ=255)                                          gTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTCGCTCAGCTGCCAAAAagag                              <  1:1023874/76‑1 (MQ=255)                                            aGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTg                            >  2:1026761/1‑76 (MQ=255)                                                  tACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGAct                      >  1:467552/1‑76 (MQ=255)                                                   aCGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGActc                     <  1:973773/76‑1 (MQ=255)                                                     gtggtgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGGCTACTCTa                   <  1:1121194/76‑1 (MQ=255)                                                        gtgAAATCCAGCGTCGGTCGCATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCt                <  2:541585/76‑1 (MQ=255)                                                         tgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGGGGGTCAGCTGCCCAAAAGAGTGCCGAGTCTACCTc               >  1:467835/1‑76 (MQ=255)                                                              tCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTccacca          <  2:937333/76‑1 (MQ=255)                                                               ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCaa         <  2:744485/76‑1 (MQ=255)                                                                cAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAg        <  1:454264/76‑1 (MQ=255)                                                                 aGCGTCGTTTACACCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGa       <  1:609360/76‑1 (MQ=255)                                                                 aGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGActct                    <  1:857943/63‑1 (MQ=255)                                                                 aGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGActct                    >  2:857943/1‑63 (MQ=255)                                                                   cGTCGGTTACATCCGCGCTTTGCGGACCAGTGGCTCAGCNGCCAACAAGAGTGCCNACTCTATCTCCACCAAGACg     >  2:1199479/1‑76 (MQ=255)                                                                      cGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCt  <  1:303478/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 < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑