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	96.2%	77.6 / ‑5.5	26	H33Y (CAC→TAC)	topA	DNA topoisomerase I
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (10/15); minor base G (0/1); total (10/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

TGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCT  >  REL606/1329443‑1329587

tGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTtac                                                                       >  2:707442/1‑76 (MQ=255)

     ccccGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCcg                                                                  >  2:863379/1‑76 (MQ=255)

        cGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGa                                                               >  2:57267/1‑76 (MQ=255)

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

             aaaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTGCATCCGCGATTTGc                                                          >  2:676159/1‑76 (MQ=255)

             aaaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGc                                                          >  1:51711/1‑76 (MQ=255)

                 ccAAAACGATCAACACGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTTCATCCTCGATTTGCCggc                                                      <  2:987037/76‑3 (MQ=255)

                          tCAACAAGTATCCGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGGTTTGCCGACCAGTGGCTc                                             <  1:557413/76‑1 (MQ=255)

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

                               aaGTATCGTGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTg                                        <  2:545850/76‑1 (MQ=255)

                                  tATCTNGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCa                                     <  2:869281/76‑1 (MQ=255)

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

                                           aGTGANTACGTGGTGAAATNCAGCGTCGGNTACATCCGCGATATGCCGACCAGTGGCTCGGCNGCCANAAAGAGTg                            <  2:167695/76‑1 (MQ=255)

                                           aGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGATTg                            >  1:518418/1‑76 (MQ=255)

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

                                              gACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCg                         <  2:1022923/76‑1 (MQ=255)

                                                cTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGAc                       >  1:460400/1‑76 (MQ=255)

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

                                                   ngtggtgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGActct                    <  1:920019/75‑1 (MQ=255)

                                                    gtggtgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTa                   <  1:490752/76‑1 (MQ=255)

                                                      ggtgAAAACCAGCGTCGGTTACATCCGGGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTAcc                 <  2:556448/76‑1 (MQ=255)

                                                       gtgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCt                <  2:690535/76‑1 (MQ=255)

                                                        tgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGGTGCCAAAAAGAGTGCCGACTCTACCTc               <  1:99650/76‑1 (MQ=255)

                                                        tgAAATCCAGCGTCGGTTACATCCGCGATCTGCTGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTc               <  1:975944/76‑1 (MQ=255)

                                                               cAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCAAAAAAGAGTGCCGACTCTACCTCCACCAAg        >  2:10021/1‑76 (MQ=255)

                                                                  cGTCGGTGACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACg     <  2:1016473/76‑1 (MQ=255)

                                                                   gTCGGTTACATCCGCGATTTGCCGACCAGTGGCGCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACgg    <  1:863379/76‑1 (MQ=255)

                                                                     cGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCAACCAAGACGGCt  >  1:978604/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` `< 32 ≤` `ATCG/ATCG` `< 38 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`