BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	516,365	T→A	42.8%	S189S (TCT→TCA)	bioB →	biotin synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	516,365	0	T	A	42.8%	21.5 / 54.9	56	S189S (TCT→TCA)	bioB	biotin synthase
Reads supporting (aligned to +/- strand): ref base T (22/10); new base A (19/5); total (41/15)
Fisher's exact test for biased strand distribution p-value = 5.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCA  >  minE/516299‑516428

ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:487252/2‑69 (MQ=255)

ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:2171489/2‑69 (MQ=255)

ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:2394844/2‑69 (MQ=255)

ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:2480938/2‑69 (MQ=255)

ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:2141110/2‑69 (MQ=255)

ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:2611706/2‑69 (MQ=255)

ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:2108761/2‑69 (MQ=255)

ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:2690297/2‑69 (MQ=255)

ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:2740157/2‑69 (MQ=255)

ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:2911204/2‑69 (MQ=255)

ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:1871857/2‑69 (MQ=255)

ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:3267887/2‑69 (MQ=255)

ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:1580494/2‑69 (MQ=255)

ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:729876/2‑69 (MQ=255)

ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:1428991/2‑69 (MQ=255)

ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:83141/2‑69 (MQ=255)

ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:133067/2‑69 (MQ=255)

ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:1153245/2‑69 (MQ=255)

ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:835675/2‑69 (MQ=255)

       tATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCAtt                                                        >  1:497114/1‑69 (MQ=255)

       tATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCAtt                                                        >  1:1902532/1‑69 (MQ=255)

       tATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCAtt                                                        >  1:525333/1‑69 (MQ=255)

       tATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCAtt                                                        >  1:115452/1‑69 (MQ=255)

       tATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCAtt                                                        >  1:894929/1‑69 (MQ=255)

        aTCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCAtt                                                        >  1:1116876/1‑68 (MQ=255)

               aCGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCtt                                                >  1:943425/1‑69 (MQ=255)

               aCGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCtt                                                >  1:3424375/1‑69 (MQ=255)

               aCGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCtt                                                >  1:1974029/1‑69 (MQ=255)

                                 aaaaGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACTGt                                    >  1:952649/1‑63 (MQ=255)

                                 aaaaGTGCGAGAAGCTGGGATTAAAGTCTGTTCAGGTGGCATTGTGGGCTTAGGCGAAACGGTGAAAg                               <  1:3379258/68‑1 (MQ=255)

                                 aaaaGTGCGAGAAGCTGGGATTAAAGTCTGTTCAGGTGGCATTGTGGGCTTAGGCGAAACGGTGAAAg                               <  1:105884/68‑1 (MQ=255)

                                 aaaaGTGCGAGAAGCTGGGATTAAAGTCTGTTCAGGTGGCATTGTGGGCTTAGGCGAAACGGTGAAAg                               <  1:3346759/68‑1 (MQ=255)

                                 aaaaGTGCGAGAAGCTGGGATTAAAGTCTGTTCAGGTGGCATTGTGGGCTTAGGCGAAACGGTGAAAg                               <  1:1974904/68‑1 (MQ=255)

                                 aaaaGTGCGAGAAGCTGGGATTAAAGTCTGTTCAGGTGGCATTGTGGGCTTAGGCGAAACGGTGAAAg                               <  1:2160617/68‑1 (MQ=255)

                                              gCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGatta                 >  1:2959342/1‑69 (MQ=255)

                                              gCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGatta                 >  1:1613775/1‑69 (MQ=255)

                                              gCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGatta                 >  1:711984/1‑69 (MQ=255)

                                              gCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGatta                 >  1:7056/1‑69 (MQ=255)

                                                     tCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATtgctgc          >  1:11453/1‑69 (MQ=255)

                                                             tGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCa  >  1:2109822/1‑69 (MQ=255)

                                                             tGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCa  >  1:3181885/1‑69 (MQ=255)

                                                             tGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCa  >  1:1707278/1‑69 (MQ=255)

                                                             tGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCa  >  1:268721/1‑69 (MQ=255)

                                                             tGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCa  >  1:2285675/1‑69 (MQ=255)

                                                             tGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCa  >  1:2645870/1‑69 (MQ=255)

                                                             tGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGAGCCGGATTATTGCTGCAACTGGCa  >  1:1975528/1‑69 (MQ=255)

                                                                 cTGGCGGCATTGTGGGCTTAGGCTAAACGGTAAAAGATCGCGCCGGattat                <  1:429035/51‑1 (MQ=255)

                                                                 cTGGCGGCATTGTGGGCTTAGGCTAAACGGTAAAAGATCGCGCCGGattat                <  1:2723332/51‑1 (MQ=255)

                                                                 cTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGattat                <  1:607195/51‑1 (MQ=255)

                                                                 cTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGattat                <  1:1495902/51‑1 (MQ=255)

                                                                 cTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGattat                <  1:1527601/51‑1 (MQ=255)

                                                                 cTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGattat                <  1:1488136/51‑1 (MQ=255)

                                                                 cTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGattat                <  1:811334/51‑1 (MQ=255)

                                                                 cTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGattat                <  1:1382334/51‑1 (MQ=255)

                                                                 cTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGattat                <  1:3005915/51‑1 (MQ=255)

                                                                 cTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGattat                <  1:2153026/51‑1 (MQ=255)

ACGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCA  >  minE/516299‑516428

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