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

GGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGC  >  REL606/1329482‑1329550                                   |                                   gggTAGTGGCTACGTGGTGAAATCCAGCGTCGGTtac                                  <  2:927980/37‑1 (MQ=37) gggTAGTGACTACGTGGTGAAATCCAGCGTCGGTtac                                  <  2:9221855/37‑1 (MQ=255) gggTAGTGACTACGTGGTGAAATCCAGCGTCGGTtac                                  >  2:9114058/1‑37 (MQ=255)  ggTAGTGACTACGTGGTGAAATCCAGCGTCGGTtaca                                 <  2:3195010/37‑1 (MQ=255)   gTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACAt                                <  1:7529178/37‑1 (MQ=255)   gTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACAt                                <  1:2404965/37‑1 (MQ=255)   gTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACAt                                <  2:6354869/37‑1 (MQ=255)    tAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATc                               <  1:10813350/37‑1 (MQ=255)    tAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATc                               >  2:5301234/1‑37 (MQ=255)     aGTGACTACGTGGTGAAATCCAGCGTCGGTTACATcc                              <  2:3856923/37‑1 (MQ=255)     aGTGACTACGTGGTGAAATCCAGCGTCGGTTACATcc                              <  2:1421693/37‑1 (MQ=255)      gTGACTAGGTGGTGAAATCCAGCGTCGGTTCCATCcg                             <  1:4077439/37‑1 (MQ=21)      gTGACTACGTGGTGAAATCCAGCGTCGGTTACATCcg                             <  2:551667/37‑1 (MQ=255)      gTGACTACGTGGTGAAATCCAGCGTCGGTTACATCcg                             <  2:4528981/37‑1 (MQ=255)      gTGACTACGTGGTGAAATCCAGCGTCGGTTACATCcg                             <  2:2651331/37‑1 (MQ=255)       tGACTACGTGGTGAAATCCAGCGTCGGTTACATCcgc                            >  2:10332620/1‑37 (MQ=255)       tGACTACGTGGTGAAATCCAGCGTCGGTTACATCcgc                            <  1:7330732/37‑1 (MQ=255)       tGACTACGTGGTGAAATCCAGCGTCGGTTACATCcgc                            >  1:3625738/1‑37 (MQ=255)       tGACTACGTGGTGAAATCCAGCGTCGGTTACATCcgc                            <  1:2619715/37‑1 (MQ=255)       tGACTACGTGGTGAAATCCAGCGTCGGTTACATCcgc                            <  2:9659318/37‑1 (MQ=255)        gACTACGTGGTGAAATCCGGCGTCGGTTACCNCcgcg                           <  2:10592265/37‑1 (MQ=21)        gACTACGTGGTGAAATCCAGCGTCGGTTACATCcgcg                           <  2:1576572/37‑1 (MQ=255)        gACTACGTGGTGAAATCCAGCGTCGGTTACATCcgcg                           <  2:10121320/37‑1 (MQ=255)         aCTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGa                          <  2:9869076/37‑1 (MQ=255)          cTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAt                         <  2:5715930/37‑1 (MQ=255)           tACGTGGTGAAATCCAGCGTCGGTTACATCCGCGCtt                        >  2:10490409/1‑37 (MQ=37)           tACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAtt                        <  1:10237402/37‑1 (MQ=255)           tACGTGGCGAAATCCAGCGTCGGTTACATCCGCGAtt                        <  2:172733/37‑1 (MQ=38)            aCGTGGTGAAATCCAGCGTCGGTTACATCCGCGAttt                       <  1:2439587/37‑1 (MQ=255)             cGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTg                      <  1:6872087/37‑1 (MQ=255)              gtggtgAAATCCAGCGTCGGTTACATCCGCGATTTGc                     >  2:10253294/1‑37 (MQ=255)              gtggtgAAATCCAGCGTCGGTTACATCCGCGATTTGc                     <  2:8541099/37‑1 (MQ=255)              gtggtgAAATCCAGCGTCGGTTACATCCGCGATTTGc                     >  2:3064759/1‑37 (MQ=255)               tggtgAAATCCAGCGTCGGTTACATCCGCGATTTGcc                    <  1:4852670/37‑1 (MQ=255)               tggtgAAATCCAGCGTCGGTTACATCCGCGATTTGcc                    <  2:1729946/37‑1 (MQ=255)               tggtgAAATCCAGCGTCGGTTACATCCGCCACTTGcc                    >  1:9145227/1‑37 (MQ=25)                ggtgAAATCCAGCGTCGGTTACATCCGCGATTTGcca                   <  2:2956813/37‑2 (MQ=38)                ggtgAAATCCAGCGTCGGTTACATCCGCGATTTGCCg                   >  2:10881595/1‑37 (MQ=255)                 gtgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGa                  <  2:8735471/37‑1 (MQ=255)                 gtgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGa                  <  1:8328933/37‑1 (MQ=255)                 gtgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGa                  <  2:3319473/37‑1 (MQ=255)                  tgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGAc                 <  1:7488329/37‑1 (MQ=255)                  tgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGAc                 <  2:8608784/37‑1 (MQ=255)                  tgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGAc                 <  2:9384193/37‑1 (MQ=255)                  tgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGAc                 >  2:1112069/1‑37 (MQ=255)                   gAAGTCCAGCGTCGGTTACATCCGCGATTTGCCGAcc                <  2:10333900/37‑1 (MQ=38)                   gAACTCCCGCGTCGGTTACATCCGCGATTTGCCGAcc                <  2:7412881/37‑1 (MQ=25)                   gAAATCCAGCGTCGGTTACATCCGCGATTTGCCGAcc                <  2:5055143/37‑1 (MQ=255)                   gAAATCCAGCGTCGGTTACATCCGCGATTTGCCGAcc                <  2:1370679/37‑1 (MQ=255)                   gAAATCCAGCGTCGGTTACATCCGCGATTTGCCGAcc                <  1:2721101/37‑1 (MQ=255)                   gAAATCCAGCGTCGGTTACATCCGCGATTTGCCGAcc                <  1:1583189/37‑1 (MQ=255)                   gAAAGCCAGCGTCGGTTACATCCGCGATTTGCCGAcc                <  2:3156647/37‑1 (MQ=38)                   gAAACCCAGCGTCGGTTACATCCGCGATTTGCCGAcc                <  1:4396825/37‑1 (MQ=37)                    aaaTCCAGCGTCGGTTACATCCGCGATTTGCCGATCa               <  2:5611127/37‑1 (MQ=38)                    aaaTCCAGCGTCGGTTACATCCGCGATTTGCCGACCa               <  2:1073324/37‑1 (MQ=255)                    aaaTCCAGCGTCGGTTACATCCGCGATTTGCCGACCa               <  1:8430733/37‑1 (MQ=255)                    aaaTCCAGCGTCGGTTACATCCGCGATTTGCCGACCa               <  2:4158390/37‑1 (MQ=255)                    aaaCCCAGCGTCGGTTACATCCGCGATTTGCCGACCa               <  1:6654914/37‑1 (MQ=37)                      aTCCATCGTCGGTTACATCCGCGATTTGCCGACCAGt             <  1:10687245/37‑1 (MQ=38)                      aTCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGt             <  2:8573798/37‑1 (MQ=255)                       tCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTg            <  2:7011209/37‑1 (MQ=255)                       gccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTg            <  2:3732681/36‑1 (MQ=38)                        ccAGCGTCGGTTACATCCGCGATTTGCCGCCCAGTgg           <  2:10835550/37‑1 (MQ=38)                        ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTgg           <  2:2700108/37‑1 (MQ=255)                        ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTgg           >  2:320512/1‑37 (MQ=255)                          aGCGTCGGTTACATCCGCGATTTGCCGACCATTGGCt         <  2:452316/37‑1 (MQ=37)                          aGCGGCGGTTACATCCGCGATTTGCCGACCAGTGGCt         <  1:5997308/37‑1 (MQ=38)                           gCGTCGGTTACATCCGCGCTTTGCCGACCAGTGGCTc        <  1:1313306/37‑1 (MQ=38)                           gCGTCGGTTACATCCGCGATTTTCCGACCCGTGGCTc        >  2:5738183/1‑37 (MQ=37)                           gCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTc        <  2:2554034/37‑1 (MQ=255)                            cGTCGGTTACATCCGCGATTTGCCGGCCAGTGGCTCa       <  1:4129102/37‑1 (MQ=38)                            cGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa       <  2:5899799/37‑1 (MQ=255)                            cGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa       <  2:6219487/37‑1 (MQ=255)                            cGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa       <  1:6427533/37‑1 (MQ=255)                             gTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAg      >  2:925156/1‑37 (MQ=255)                             gTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAg      <  2:2504108/37‑1 (MQ=255)                             gTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAg      <  1:10245692/37‑1 (MQ=255)                              tCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGc     <  1:4864278/37‑1 (MQ=255)                              tCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGc     >  1:7414739/1‑37 (MQ=255)                              tCGGTTACATCCGCGATTTGCCGACCAGTGGCGCAGc     >  1:2385390/1‑37 (MQ=38)                               cGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCt    >  2:7389596/1‑37 (MQ=255)                               cGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCt    <  1:936493/37‑1 (MQ=255)                               cGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCt    <  2:5030417/37‑1 (MQ=255)                                ggTTACATCCGCGATTTGCCGAGCAGTGGCTCAGCTg   <  2:154977/37‑1 (MQ=38)                                ggTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTg   <  2:8820670/37‑1 (MQ=255)                                ggTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTg   <  1:1422374/37‑1 (MQ=255)                                ggTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTg   >  1:2887190/1‑37 (MQ=255)                                ggTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTg   <  2:18969/37‑1 (MQ=255)                                ggTTACATCCGCGATTTGCCGA‑CAGTGGCTCAGCTGc  <  2:6007611/37‑1 (MQ=25)                                   |                                   GGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGC  >  REL606/1329482‑1329550

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