BRESEQ :: Evidence

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

GAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGC  >  REL606/1329444‑1329586

gAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTtac                                                                      <  1:1248574/75‑1 (MQ=255)

gAGTCCCCGGCAAAAGCCAAAACGATCAACAAATATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTtac                                                                      >  1:720216/1‑75 (MQ=255)

 aGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTtaca                                                                     <  2:711012/75‑1 (MQ=255)

  gTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACAt                                                                    <  1:29857/75‑1 (MQ=255)

     cccGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCcg                                                                 <  1:163165/75‑1 (MQ=255)

     cccGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTAAATCcg                                                                 >  2:638293/1‑75 (MQ=255)

        ggCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGa                                                              <  1:455417/75‑1 (MQ=255)

         gCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAt                                                             <  2:280797/75‑1 (MQ=255)

         gCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAt                                                             <  1:1179458/75‑1 (MQ=255)

           aaaaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAttt                                                           <  1:111929/75‑1 (MQ=255)

            aaaGCCAAAACGATCAACAAGTATCTGGGTAGTGACCACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTg                                                          <  2:1806220/75‑1 (MQ=255)

             aaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGc                                                         >  2:556178/1‑75 (MQ=255)

              aGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGcc                                                        <  1:1591336/75‑1 (MQ=255)

              aGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGcc                                                        <  2:1790661/75‑1 (MQ=255)

                 cAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGAc                                                     >  1:1627792/1‑75 (MQ=255)

                    aaCGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAg                                                  <  2:829313/75‑1 (MQ=255)

                     aCGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGTCCAGt                                                 >  1:658178/1‑75 (MQ=255)

                     aCGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGt                                                 <  2:313061/75‑1 (MQ=255)

                     aCGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGt                                                 <  1:846827/75‑1 (MQ=255)

                        aTCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGCCCAGTGGc                                              <  2:217101/75‑1 (MQ=255)

                        aTCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGc                                              >  1:1197690/1‑75 (MQ=255)

                         tCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCt                                             >  1:1243446/1‑75 (MQ=255)

                           aacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa                                           <  2:971501/75‑1 (MQ=255)

                           aacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa                                           <  1:1667315/75‑1 (MQ=255)

                            acaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAg                                          <  2:609424/75‑1 (MQ=255)

                                gTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGc                                      >  2:1770475/1‑75 (MQ=255)

                                gTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGc                                      >  2:654789/1‑75 (MQ=255)

                                   tCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaa                                   >  2:474516/1‑75 (MQ=255)

                                     tGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaaa                                 <  2:838834/75‑1 (MQ=255)

                                     tGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCCGCCaaaa                                 >  2:331109/1‑75 (MQ=255)

                                      gggTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaaaa                                >  2:1321464/1‑75 (MQ=255)

                                      gggTAGTGACTACGTAGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaaaa                                >  1:284616/1‑75 (MQ=255)

                                       ggTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCcaaaaa                               >  1:403339/1‑74 (MQ=255)

                                        gTAGTGTCTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAaga                              <  2:360297/75‑1 (MQ=255)

                                        gTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAaga                              <  1:217756/75‑1 (MQ=255)

                                         tAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAagag                             <  2:1375249/75‑1 (MQ=255)

                                         tAGTGACTACGTGGTGAAATCCAGCGTCGGTTAAATCCGCGGTTTGCGGACCAGGGGCTCGCCTGCCAACAcgag                             >  1:332125/1‑75 (MQ=255)

                                          aGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGt                            <  2:223596/75‑1 (MQ=255)

                                             gACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGcc                         <  1:1279653/75‑1 (MQ=255)

                                             gACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGcc                         >  1:689193/1‑75 (MQ=255)

                                             gACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGcc                         >  1:2140000/1‑75 (MQ=255)

                                              aCTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCg                        <  1:1994093/75‑1 (MQ=255)

                                                tACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGAc                      >  2:1343442/1‑75 (MQ=255)

                                                tACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGAc                      <  2:1070712/75‑1 (MQ=255)

                                                 aCGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGAct                     >  1:1218875/1‑75 (MQ=255)

                                                 aCGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGAct                     >  2:1883415/1‑75 (MQ=255)

                                                  cGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGActc                    <  2:525436/75‑1 (MQ=255)

                                                  cGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGActc                    <  2:1278267/75‑1 (MQ=255)

                                                   gtggggAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGActct                   <  1:393802/75‑1 (MQ=255)

                                                    tggtgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTa                  >  1:935381/1‑75 (MQ=255)

                                                       tgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCt               >  1:1930078/1‑75 (MQ=255)

                                                       tgAAAGCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCt               <  2:843916/75‑1 (MQ=255)

                                                        gAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTc              >  2:745819/1‑75 (MQ=255)

                                                        gAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTc              >  1:1011688/1‑75 (MQ=255)

                                                          aaTCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTcca            <  2:1035606/75‑1 (MQ=255)

                                                            tCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTccacc          <  1:881000/75‑1 (MQ=255)

                                                             ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTccaccc         >  1:509643/1‑74 (MQ=255)

                                                             ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTccacca         <  2:1789085/75‑1 (MQ=255)

                                                              cAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTccaccag        >  1:963817/1‑74 (MQ=255)

                                                               aGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCCCCAAg       >  2:1000808/1‑75 (MQ=255)

                                                                 cGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGAc     <  2:462199/75‑1 (MQ=255)

                                                                 cGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGAc     <  2:1867830/75‑1 (MQ=255)

                                                                 cGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGAc     <  1:1554048/75‑1 (MQ=255)

                                                                  gTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACg    >  2:2030560/1‑75 (MQ=255)

                                                                   tCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACgg   <  2:1270405/75‑1 (MQ=255)

                                                                    cGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGc  <  2:1420268/75‑1 (MQ=255)

GAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGC  >  REL606/1329444‑1329586

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