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

TGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCT  >  REL606/1329443‑1329587                                                                          |                                                                        tGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTtac                                                                       >  2:1584152/1‑76 (MQ=255)        ccGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCcgcg                                                                <  2:970910/76‑1 (MQ=255)              aaaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGc                                                          >  2:1257343/1‑76 (MQ=255)              aaaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGc                                                          <  1:1176218/76‑1 (MQ=255)              aaaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGGGGTGAAATCCAGCGTCGGTTACATCCGCAATTTGc                                                          >  2:374902/1‑76 (MQ=255)                aGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGCTTTGCCg                                                        >  1:1152062/1‑76 (MQ=255)                  ccAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGAc                                                      >  1:1411272/1‑76 (MQ=255)                   cAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGAcc                                                     <  1:1629403/76‑1 (MQ=255)                   cAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGAcc                                                     >  2:953867/1‑76 (MQ=255)                     aaaCGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAg                                                   >  1:381957/1‑76 (MQ=255)                        cGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTgg                                                >  1:739593/1‑76 (MQ=255)                         gATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGc                                               <  2:94903/76‑1 (MQ=255)                         gATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGc                                               >  1:322144/1‑76 (MQ=255)                         gATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGc                                               <  2:537098/76‑1 (MQ=255)                          aTCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCt                                              >  1:192371/1‑76 (MQ=255)                          aTCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCt                                              >  2:655344/1‑76 (MQ=255)                            caacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa                                            >  2:353835/1‑76 (MQ=255)                            caacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa                                            <  2:655054/76‑1 (MQ=255)                            caacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa                                            <  2:1473112/76‑1 (MQ=255)                            caacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa                                            <  2:1415860/76‑1 (MQ=255)                            caacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa                                            <  1:1080452/76‑1 (MQ=255)                            caacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa                                            <  1:1560563/76‑1 (MQ=255)                                aaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCg                                                        <  2:533949/60‑1 (MQ=255)                                aaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCg                                                        >  1:533949/1‑60 (MQ=255)                                aaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTg                                        <  2:1168027/76‑1 (MQ=255)                                    aTCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaa                                    >  1:991255/1‑76 (MQ=255)                                     tCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaa                                   >  1:212604/1‑76 (MQ=255)                                      cTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaaa                                  >  1:423811/1‑76 (MQ=255)                                       tGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaaaa                                 >  2:808258/1‑76 (MQ=255)                                        gggTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTTCCGACCAGTTGCTCAGCTGCCaaaaca                                >  2:321755/1‑74 (MQ=255)                                        gggTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAag                                <  1:98693/76‑1 (MQ=255)                                           tAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGt                             >  1:499149/1‑76 (MQ=255)                                                aCTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGa                        >  2:1238401/1‑76 (MQ=255)                                                 cTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGAc                       <  1:468175/76‑1 (MQ=255)                                                 cTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGAc                       >  2:647235/1‑76 (MQ=255)                                                 cTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGAc                       <  2:396361/76‑1 (MQ=255)                                                  tACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGAct                      <  2:84655/76‑1 (MQ=255)                                                    cGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGActct                    <  2:23943/76‑1 (MQ=255)                                                    cGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGActct                    <  2:428150/76‑1 (MQ=255)                                                     gtggtgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTa                   >  1:1025991/1‑76 (MQ=255)                                                     gtggtgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTa                   <  2:1144763/76‑1 (MQ=255)                                                        gtgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCt                >  2:103002/1‑76 (MQ=255)                                                         tgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTc               <  2:1241147/76‑1 (MQ=255)                                                            aaTCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTccac            >  2:785309/1‑76 (MQ=255)                                                               ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCaa         <  2:1698287/76‑1 (MQ=255)                                                                cAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAg        <  1:1262226/76‑1 (MQ=255)                                                                  gCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGAc      <  1:253826/76‑1 (MQ=255)                                                                   cGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACg     <  1:1699426/76‑1 (MQ=255)                                                                    ggcGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACgg    <  2:1484615/74‑1 (MQ=255)                                                                    gTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACgg    <  2:205746/76‑1 (MQ=255)                                                                      cGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCt  <  1:1141363/76‑1 (MQ=255)                                                                          |                                                                        TGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCT  >  REL606/1329443‑1329587

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