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

TGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTA  >  REL606/1329443‑1329588                                                                          |                                                                         tGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTtac                                                                        <  2:176070/76‑1 (MQ=255) tGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTtac                                                                        >  1:1149290/1‑76 (MQ=255)  gAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTtaca                                                                       <  1:358713/76‑1 (MQ=255)     tCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGAGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATcc                                                                    >  1:38389/1‑76 (MQ=255)     tCCCCGGCAAAAGCCAAAACGATCAACAAGTATCCGGGTAGNGACTACGTGGTGAAATCCNGCNTCGGTTACATcc                                                                    >  1:1498381/1‑76 (MQ=255)            cAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAttt                                                             <  2:1044892/76‑1 (MQ=255)            cAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAttt                                                             >  1:618721/1‑76 (MQ=255)            cAAAAGCCAAAACGATCAACAAGTATCTGGGTAGAGANTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAttt                                                             >  1:1780648/1‑76 (MQ=255)              aaaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGc                                                           >  2:913726/1‑76 (MQ=255)              aaaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGc                                                           >  2:747979/1‑76 (MQ=255)                   cAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGATATCCAGCGTCGGNTACATCCGAGNTTTGCCGAcc                                                      <  2:1804831/76‑1 (MQ=255)                    aaaaCGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCa                                                     <  1:499826/76‑1 (MQ=255)                      aaCGATCAACAAGTNTCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGTGAGTTGCCGACCAGt                                                   <  1:1562059/76‑1 (MQ=255)                          aTCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCt                                               <  2:735349/76‑1 (MQ=255)                           tCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTc                                              <  2:123584/76‑1 (MQ=255)                           tCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTc                                              <  2:1149289/76‑1 (MQ=255)                           tCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTc                                              <  2:450321/76‑1 (MQ=255)                            caacaaGTATCTGGGTAGTGACTACGTGGTGACATCCAGCGTCGGTTACATCCGCGATTTACCGACCAGTGGCTcc                                             >  1:1520157/1‑75 (MQ=255)                            caacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa                                             >  1:257330/1‑76 (MQ=255)                            caacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa                                             <  1:1388200/76‑1 (MQ=255)                               caaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCt                                          <  2:1180708/76‑1 (MQ=255)                               caaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGCTACATCCGCGATTTGCCGACCAGTGGCTCAGCt                                          >  1:1337270/1‑76 (MQ=255)                                  gTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGcc                                       >  1:1508529/1‑76 (MQ=255)                                  gTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGcc                                       >  1:914168/1‑76 (MQ=255)                                    aTCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaa                                     >  2:1115026/1‑76 (MQ=255)                                     tCTGGGTAGTGAGTACGGAGTGAAATCCAACGTCGGTTACGTCAGCGATTTGCCGACCAGTGGCTCAGCTGCCaaa                                    <  2:30037/76‑1 (MQ=255)                                     tCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaa                                    <  1:1377105/76‑1 (MQ=255)                                     tCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaa                                    <  1:141589/76‑1 (MQ=255)                                     tCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaa                                    <  2:1621543/76‑1 (MQ=255)                                     tCTGGGTAGTAACTACGAGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTCGCTCAGCTGCCaaa                                    >  1:144646/1‑76 (MQ=255)                                      cTGGGTAGTGACTACGTGGTGAAATCCTGCGACGNTTACCTCNGCNATNTGCCTACCAGTGGCTCAGCTGCCaaaa                                   >  1:1591324/1‑76 (MQ=255)                                      cTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaaa                                   <  2:1678091/76‑1 (MQ=255)                                        gggTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAag                                 <  2:417879/76‑1 (MQ=255)                                           tAGTGACTACGTAGTCAAATCCAGCGTCGGTTAAATCCGCGATTTGCCGACCAGTGGCTCANCTGCCAAAAAGAGt                              >  1:1567732/1‑76 (MQ=255)                                            aGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTg                             >  1:659440/1‑76 (MQ=255)                                               gACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCg                          >  2:1199270/1‑76 (MQ=255)                                               gACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCg                          <  2:1439695/76‑1 (MQ=255)                                                  tACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGAct                       >  1:455412/1‑76 (MQ=255)                                                    cGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGActct                     <  2:734769/76‑1 (MQ=255)                                                    cGTGGTAAAATCCAGCGTCGTTTACATCCGCGATTTGCCGACCAGTGGATCAGCTGCCAAAAAGAGTGCCGActct                     <  1:1625599/76‑1 (MQ=255)                                                     gtggtgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTa                    >  1:1102675/1‑76 (MQ=255)                                                     gtggtgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTa                    <  1:176932/76‑1 (MQ=255)                                                         tgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTc                <  2:1708899/76‑1 (MQ=255)                                                         tgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTc                <  2:1688816/76‑1 (MQ=255)                                                         tgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTc                <  2:618720/76‑1 (MQ=255)                                                          gAAATCCAGCGTCGTTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTcc               <  2:455411/76‑1 (MQ=255)                                                           aaaTCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTcca              >  1:1099810/1‑76 (MQ=255)                                                            aaTCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTccac             <  1:476272/76‑1 (MQ=255)                                                               ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCNCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCaa          <  1:1778693/76‑1 (MQ=255)                                                                    gTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACgg     >  1:1169746/1‑76 (MQ=255)                                                                       ggTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTa  <  2:1440800/76‑1 (MQ=255)                                                                       ggTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTa  >  2:963041/1‑76 (MQ=255)                                                                       ggTTACATCCGCGATTTGCCAACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTa  <  1:1811460/76‑1 (MQ=255)                                                                          |                                                                         TGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTA  >  REL606/1329443‑1329588

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