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

GTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACG  >  REL606/1329446‑1329584                                                                       |                                                                     gTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATc                                                                 >  2:671429/1‑76 (MQ=255)  tCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATcc                                                                >  1:495132/1‑76 (MQ=255)   ccccGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCcg                                                               <  2:800286/76‑1 (MQ=255)    cccGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCcgc                                                              >  2:973211/1‑76 (MQ=255)       ggCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAt                                                           <  2:138453/76‑1 (MQ=255)         cAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAttt                                                         <  1:764439/76‑1 (MQ=255)         cAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAttt                                                         >  2:77181/1‑76 (MQ=255)         cAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAttt                                                         >  1:584734/1‑76 (MQ=255)         cAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTAAATCCGCGAttt                                                         <  1:72715/76‑1 (MQ=255)         cAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGATTAGATCCGCGAttt                                                         <  1:72702/76‑1 (MQ=255)           aaaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGc                                                       <  1:219786/76‑1 (MQ=255)             aGCCAAAACGATCACCAAGTATCTCCGTAGTGACTACGTGGTCAAATCCAGCGTCGGTTACATCTGCGATTTGCCg                                                     <  1:31791/76‑1 (MQ=255)                cAAAACGATCAACAAGTATCTGGGTAGTGACGACNTGGTGAANTCNAGNGTCGCTTAGATCCGCGATTTGCCGAcc                                                  >  1:1156253/1‑76 (MQ=255)                 aaaaCGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCa                                                 <  1:979915/76‑1 (MQ=255)                  aaaCGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAg                                                <  1:317194/76‑1 (MQ=255)                    aCGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTg                                              <  2:237552/76‑1 (MQ=255)                     cGATCAAGAAGTATGTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTgg                                             <  1:31681/76‑1 (MQ=255)                     cGATCAACAAGTATCTGGTTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTgg                                             <  2:726591/76‑1 (MQ=255)                     cGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTgg                                             >  2:606097/1‑76 (MQ=255)                      gATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGc                                            <  2:485834/76‑1 (MQ=255)                      gATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGc                                            <  2:678609/76‑1 (MQ=255)                        tCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTc                                          >  1:1017669/1‑76 (MQ=255)                        tCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTc                                          <  2:125399/76‑1 (MQ=255)                          aacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAg                                        <  2:681325/76‑1 (MQ=255)                                tATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCa                                  >  2:405258/1‑76 (MQ=255)                                   cTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCCGCCaaaa                               >  1:417002/1‑76 (MQ=255)                                      ggTAGTGACTACGTGGTGAACTCCAGCGTCGGTTACATCCGCGATTTGCCGACTAGTGGCTCAGCTGCCAAAAaga                            >  2:109337/1‑76 (MQ=255)                                        tAGTGACTACTTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGt                          <  1:320737/76‑1 (MQ=255)                                           tGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGcc                       >  1:410483/1‑76 (MQ=255)                                            gACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCg                      <  1:895438/76‑1 (MQ=255)                                               tACGTGGTGAAATCCAGCCTCTGTTAAAGCCGCGTATTGTCGATCAGTAGCTCAGCTGCCGATAAGGGTGCCGAct                   >  1:1339491/1‑76 (MQ=255)                                                 cGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGActct                 <  2:1267969/76‑1 (MQ=255)                                                 cGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGActct                 >  1:394961/1‑76 (MQ=255)                                                   tggtgAAATCCAGCGTCGGTTCCATCCNCGCTTTGNNGACCACTGGCTCAGCTGCCATAAAGAGTACCTACTCTAc               >  1:1322838/1‑76 (MQ=255)                                                     gtgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCt             <  2:772402/76‑1 (MQ=255)                                                      tgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTc            <  1:688932/76‑1 (MQ=255)                                                        aaaTCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTcca          <  2:1233642/76‑1 (MQ=255)                                                        aaaTCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTcca          >  1:1029638/1‑76 (MQ=255)                                                        aaaCCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTcca          <  2:653488/76‑1 (MQ=255)                                                         aaTCCAGCGTCGGTTACATCCGCGATTNGCCGACCACTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTccac         >  1:1169051/1‑76 (MQ=255)                                                           tCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTccacca       >  1:308629/1‑76 (MQ=255)                                                             cAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAg     >  1:752799/1‑76 (MQ=255)                                                                cGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACg  >  1:813972/1‑76 (MQ=255)                                                                       |                                                                     GTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACG  >  REL606/1329446‑1329584

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