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	98.2%	197.6 / ‑6.1	56	H33Y (CAC→TAC)	topA	DNA topoisomerase I
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (24/31); total (24/32)
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.97e-01

GAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAA  >  REL606/1329444‑1329590

gAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTtaca                                                                         >  2:830472/1‑76 (MQ=255)

  gTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATc                                                                       >  2:1963624/1‑76 (MQ=255)

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

  gTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATc                                                                       >  2:1325643/1‑76 (MQ=255)

  gTCCCCGGCAAAAGCCAAAAAGATCAACAAGTATCTGGGTAGTGACTACGTGGCGAAAGCCAGCGTCGGTTACATc                                                                       >  1:94810/1‑76 (MQ=255)

    ccacGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCcg                                                                     >  2:1866098/4‑76 (MQ=255)

     cccGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCcgc                                                                    >  2:1570162/1‑76 (MQ=255)

          cAAAAGCCAAAACGAACAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGAttt                                                               >  2:552740/1‑76 (MQ=255)

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

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

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

               gCCAAAACGATCNACNAGTATCTGGGTAGTGACTNCGTGGTTAAATCCAGCGTCGGTTACATCCGCGATTTGCCGa                                                          <  1:1658681/76‑1 (MQ=255)

                 cAAAACGATCAACANGTATCTGGGTAGTGACTACGTGGTGAAGTCCAGCGTCGGTTACATTCGCGATTTGCCGAcc                                                        <  1:1742663/76‑1 (MQ=255)

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

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

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

                    aaCGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGt                                                     >  2:1006827/1‑76 (MQ=255)

                     aCGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTg                                                    <  2:138929/76‑1 (MQ=255)

                        aTCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGCCCAGTGGCt                                                 >  1:1249740/1‑76 (MQ=255)

                         tCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTc                                                <  2:654583/76‑1 (MQ=255)

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

                         tCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTc                                                <  2:149070/76‑1 (MQ=255)

                         tCAACAAGTATCAGAGTAGAGACTACGTGCTGAAATCCAGGGTCGGTTACCTCCGCGATTTGCCGACGAGTGGCTc                                                <  2:40853/76‑1 (MQ=255)

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

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

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

                              aaGTATCTGGGTAGTGACTACGTGGTAAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTg                                           <  1:1097992/76‑1 (MQ=255)

                               aGTATCTGGGTAGTGCCTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGc                                          <  1:1690715/76‑1 (MQ=255)

                               aGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGc                                          >  2:1537292/1‑76 (MQ=255)

                               aGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGc                                          >  2:836845/1‑76 (MQ=255)

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

                                   tCTGGGTATTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaa                                      <  2:537538/76‑1 (MQ=255)

                                     tGGGTAGTGACTGCGGGGTGAAATCCAGCATCGCTTACATCCGCGATCCGCCGACCAGTGGCTCAGCGGCCaaaaa                                    >  1:1973366/1‑76 (MQ=255)

                                     tGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaaaa                                    <  2:1145308/76‑1 (MQ=255)

                                     tGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaaaa                                    >  2:1020233/1‑76 (MQ=255)

                                     tGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaaaa                                    >  2:1511924/1‑76 (MQ=255)

                                       ggTAGTGACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAaga                                  <  1:88786/76‑1 (MQ=255)

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

                                          aGTGACTACGTGGTGAAATCCAGCGTCGGTTACATGTGCGATTTGCTGACCAGTGGTTCAGTTCCCAAAAATAGTg                               >  2:43586/1‑76 (MQ=255)

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

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

                                             gACTACGTGGTGAAATCCAGCGTCGGTTACATCCGCGATTTACCGACCCGTAGCTCAGCTGCCCAAANGAGTGCCg                            >  2:346684/1‑76 (MQ=255)

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

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

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

                                               cTACGTGGTGAAATCCAGCGNCGGTTACATGCGCGANTTGCCGACCAGTGACACAGCTCCCAAAAAGAGTGCCGAc                          <  2:2033624/76‑1 (MQ=255)

                                                tACGTGGGGATATCCAGCATCGGTTACATCCGCGATTTGNNGACCAGTNGCTCAGCTGCCAAAAAGAGTGCCGAct                         <  1:2004966/76‑1 (MQ=255)

                                                   gtggtgAAATTCAGCGTGGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAATGAGTGCCGACTCTa                      >  1:169311/1‑76 (MQ=255)

                                                   gtggtgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTa                      >  1:826549/1‑76 (MQ=255)

                                                    tggtgAAATCCAGCNTCGGTTACATCCGCGATTTGCAGACCAGTGGCTCAGCTGCCAAACAGAGTGCAGACTCTAc                     >  1:503/1‑76 (MQ=255)

                                                     ggtgAAATCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTAcc                    <  2:488522/76‑1 (MQ=255)

                                                           aTCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTTCCTccacc              <  1:1901608/76‑1 (MQ=255)

                                                            tCCAGCGTCGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTccacca             <  2:1071058/76‑1 (MQ=255)

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

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

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

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

                                                                    cGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCt     >  1:529545/1‑76 (MQ=255)

                                                                    cGGTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCt     <  1:1943623/76‑1 (MQ=255)

                                                                     ggTTACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTa    >  2:1212267/1‑76 (MQ=255)

                                                                       tCACATCCGCGATTTGCCGACCAGTGGCGCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTaaa  <  1:776426/76‑1 (MQ=255)

GAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAA  >  REL606/1329444‑1329590

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: `‑`