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	99.4%	619.8 / NA	155	H33Y (CAC→TAC)	topA	DNA topoisomerase I
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (84/70); minor base G (0/1); total (84/71)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 179 total aligned reads displayed.

TGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCT  >  REL606/1329443‑1329587

 gAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTtaca                                                                      <  2:5174154/76‑1 (MQ=255)

   gTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATc                                                                    <  1:4804981/76‑1 (MQ=255)

   gTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATc                                                                    <  1:6393397/76‑1 (MQ=255)

    tCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATcc                                                                   >  2:7892000/1‑76 (MQ=255)

     ccccGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCcg                                                                  <  1:3904017/76‑1 (MQ=255)

      cccGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCcgc                                                                 <  1:498995/76‑1 (MQ=255)

       ccGGCAAAAGCCAAAACGATCAACAAGTATCTGAGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCcgcg                                                                <  1:1073121/76‑1 (MQ=255)

        cGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCCCGa                                                               >  2:2843937/1‑76 (MQ=255)

         ggCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAt                                                              <  1:1112339/76‑1 (MQ=255)

         ggCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAt                                                              >  1:5267368/1‑76 (MQ=255)

          gCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAtt                                                             <  2:7535981/76‑1 (MQ=255)

           cAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAttt                                                            <  1:1595981/76‑1 (MQ=255)

            aaaaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACCTCCGCGATTTg                                                           >  1:3397465/1‑76 (MQ=255)

            aaaaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGGTTTg                                                           >  2:4881902/1‑76 (MQ=255)

            aaaaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTg                                                           >  1:6721182/1‑76 (MQ=255)

            aaaaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTg                                                           >  1:5688983/1‑76 (MQ=255)

            aaaaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTg                                                           >  2:1913314/1‑76 (MQ=255)

            aaaaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTg                                                           >  2:4914637/1‑76 (MQ=255)

              aaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGcc                                                         >  2:6995830/1‑76 (MQ=255)

               aGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCg                                                        <  1:6052726/76‑1 (MQ=255)

                gCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGa                                                       <  2:2971230/76‑1 (MQ=255)

                gCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGa                                                       >  2:1403053/1‑76 (MQ=255)

                  cAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGAcc                                                     >  2:6103135/1‑76 (MQ=255)

                  cAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTAAATCCCCGATTTGGCGAcc                                                     >  2:2368199/1‑76 (MQ=255)

                   aaaaCGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCa                                                    <  2:105144/76‑1 (MQ=255)

                   aaaaCGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCa                                                    >  2:446530/1‑76 (MQ=255)

                   aaaaCGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCa                                                    >  2:1974258/1‑76 (MQ=255)

                    aaaCGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAg                                                   >  2:613964/1‑76 (MQ=255)

                        gATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGc                                               >  2:6790618/1‑76 (MQ=255)

                        gATCAACAAGTATCTGGGTAGTGACTACGTGGAGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGtagg                                               >  2:1710793/1‑73 (MQ=255)

                         aTCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCt                                              >  1:5486372/1‑76 (MQ=255)

                          tCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTc                                             >  2:1108055/1‑76 (MQ=255)

                           caacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCa                                            <  1:2736916/76‑1 (MQ=255)

                            aacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAg                                           <  2:1588893/76‑1 (MQ=255)

                             acaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGc                                          >  2:6564669/1‑76 (MQ=255)

                             acaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGc                                          <  1:6004486/76‑1 (MQ=255)

                              caaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCt                                         <  2:3166152/76‑1 (MQ=255)

                               aaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTg                                        <  2:5767082/76‑1 (MQ=255)

                               aaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTg                                        >  1:5226585/1‑76 (MQ=255)

                               aaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTg                                        >  1:2666170/1‑76 (MQ=255)

                                aGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGGGACTTGCCGACCAGTGGCTCAGCCGc                                       >  1:5392757/1‑76 (MQ=255)

                                aGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCCGc                                       >  2:3658178/1‑76 (MQ=255)

                                 gTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGcc                                      <  1:1271540/76‑1 (MQ=255)

                                 gTATCTGGGTAGTGACTACGTGGTGAAATACAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGcc                                      <  1:6306083/76‑1 (MQ=255)

                                  tATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCa                                     >  1:6759729/1‑76 (MQ=255)

                                   aTCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaa                                    >  1:1497447/1‑76 (MQ=255)

                                   aTCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaa                                    >  1:4804807/1‑76 (MQ=255)

                                   aTCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaa                                    <  2:6264834/76‑1 (MQ=255)

                                    tCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaa                                   <  1:2431059/76‑1 (MQ=255)

                                    tCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaa                                   <  1:1825976/76‑1 (MQ=255)

                                     cTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaaa                                  <  1:722537/76‑1 (MQ=255)

                                     cTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaaa                                  <  1:1967014/76‑1 (MQ=255)

                                     cTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaaa                                  >  2:6845047/1‑76 (MQ=255)

                                     cTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaaa                                  >  2:5216742/1‑76 (MQ=255)

                                         gTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAagag                              >  2:4157838/1‑76 (MQ=255)

                                         gTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAagag                              <  2:2183300/76‑1 (MQ=255)

                                          tAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGt                             <  2:5935538/76‑1 (MQ=255)

                                           aGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCCAAAAGAGTg                            >  1:4580756/1‑76 (MQ=255)

                                           aGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTg                            <  2:4740818/76‑1 (MQ=255)

                                           aGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTg                            >  1:7685043/1‑76 (MQ=255)

                                            gTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGc                           <  2:3170894/76‑1 (MQ=255)

                                             tGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGcc                          >  1:3159640/1‑76 (MQ=255)

                                              gACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCCAAAAGAGTGCCg                         >  1:7037936/1‑76 (MQ=255)

                                              gACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCg                         >  1:3292777/1‑76 (MQ=255)

                                               aCTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGa                        <  1:1808113/76‑1 (MQ=255)

                                                cTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGAc                       <  2:7291327/76‑1 (MQ=255)

                                                 tACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGAct                      <  1:2170309/76‑1 (MQ=255)

                                                  aCGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGActc                     >  2:6909582/1‑76 (MQ=255)

                                                   cGTGTTGAAACCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGActct                    <  2:806760/76‑1 (MQ=255)

                                                   cGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGActct                    <  2:901003/76‑1 (MQ=255)

                                                    gtggtgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGGGCCGACTCTa                   >  1:4754206/1‑76 (MQ=255)

                                                    gtggtaAAACCCAGCGTCGTTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTa                   <  2:6655985/76‑1 (MQ=255)

                                                     tggtgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTAc                  <  1:492623/76‑1 (MQ=255)

                                                     tggtgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTAc                  <  1:3012350/76‑1 (MQ=255)

                                                      ggtgAAGCCCAGCGTCGGTGACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTAcc                 <  2:3105865/76‑1 (MQ=255)

                                                        tgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTc               >  1:4206587/1‑76 (MQ=255)

                                                        tgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTc               <  2:2412546/76‑1 (MQ=255)

                                                         gAAATCCAGCGTCGTTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTcc              <  1:6032052/76‑1 (MQ=255)

                                                         gAAATCCAGCGTCGTTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTcc              <  1:4479742/76‑1 (MQ=255)

                                                         gAAACCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTcc              <  1:1398248/76‑1 (MQ=255)

                                                          aaaTCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTcca             >  1:6619960/1‑76 (MQ=255)

                                                          aaaTCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTcca             >  1:4850316/1‑76 (MQ=255)

                                                           aaTCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTccac            >  2:1893477/1‑76 (MQ=255)

                                                              ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCaa         >  2:3304953/1‑76 (MQ=255)

                                                              ccAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCaa         >  2:3965269/1‑76 (MQ=255)

                                                              ccAGCGTCGGTTACATCCGCGATTTGACGACCAGTGGCTCAGCTGCCCAAAAGAGGGCCGGCTCTACCTCCACCaa         >  2:5081549/1‑76 (MQ=255)

                                                               cAGCGTCGTTTACATCCGCTATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAg        <  2:2642284/76‑1 (MQ=255)

                                                               cAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAACAAGAGTGCCGACTCTACCTCCACCAAg        <  1:4132209/76‑1 (MQ=255)

                                                               cAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAg        <  1:3772993/76‑1 (MQ=255)

                                                               cAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAg        <  1:3935739/76‑1 (MQ=255)

                                                               cAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAg        <  2:3236885/76‑1 (MQ=255)

                                                               cAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAg        >  1:1377827/1‑76 (MQ=255)

                                                                aGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGa       >  1:7969270/1‑76 (MQ=255)

                                                                 gCGTCGGTTACATCCGCGATTTGCGGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGAc      <  1:1704781/76‑1 (MQ=255)

                                                                 gCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGAc      >  2:3762585/1‑76 (MQ=255)

                                                                 gCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGAc      <  2:7511168/76‑1 (MQ=255)

                                                                  cGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCGGACTCTACCTCCACCAAGACg     >  2:4781914/1‑76 (MQ=255)

                                                                   gTCGTTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACgg    <  1:6168738/76‑1 (MQ=255)

                                                                   gTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACgg    <  1:7413703/76‑1 (MQ=255)

                                                                     cGGTTACATCCGCGATTTTCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCt  >  1:1225776/1‑76 (MQ=255)

TGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCT  >  REL606/1329443‑1329587

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