BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	240,000	G→A	39.0%	N317N (AAC→AAT)	araJ ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	240,000	0	G	A	39.0%	17.9 / 24.4	41	N317N (AAC→AAT)	araJ	predicted transporter
Reads supporting (aligned to +/- strand): ref base G (5/20); new base A (14/2); total (19/22)
Fisher's exact test for biased strand distribution p-value = 2.76e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAGGGCAAATAATCCCGC  >  minE/239932‑240050

ccGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCg                                                    <  1:2443484/69‑1 (MQ=255)

ccGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCg                                                    <  1:1295024/69‑1 (MQ=255)

ccGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCg                                                    <  1:2158927/69‑1 (MQ=255)

ccGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCg                                                    <  1:2153798/69‑1 (MQ=255)

ccGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCg                                                    <  1:3048536/69‑1 (MQ=255)

ccGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCg                                                    <  1:3294598/69‑1 (MQ=255)

ccGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCg                                                    <  1:365147/69‑1 (MQ=255)

ccGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCg                                                    <  1:329774/69‑1 (MQ=255)

 cGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCg                                                    <  1:2693588/68‑1 (MQ=255)

 cGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCg                                                    <  1:1783868/68‑1 (MQ=255)

     ggCGCTACCGAGTTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTg                                               <  1:1854504/69‑1 (MQ=255)

     ggCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTg                                               <  1:2294620/69‑1 (MQ=255)

     ggCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTg                                               <  1:100852/69‑1 (MQ=255)

     ggCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTg                                               <  1:2734435/69‑1 (MQ=255)

     ggCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTg                                               <  1:2490298/69‑1 (MQ=255)

     ggCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTg                                               <  1:2470395/69‑1 (MQ=255)

     ggCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTg                                               <  1:1675882/69‑1 (MQ=255)

     ggCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTg                                               <  1:288826/69‑1 (MQ=255)

     ggCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTg                                               <  1:1392321/69‑1 (MQ=255)

     ggCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTg                                               <  1:1867127/69‑1 (MQ=255)

     ggCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTg                                               <  1:1741006/69‑1 (MQ=255)

      atccTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTg                                               <  1:462791/65‑1 (MQ=255)

      atccTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTg                                               <  1:1347726/65‑1 (MQ=255)

      atccTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTg                                               <  1:2440933/65‑1 (MQ=255)

      atccTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTg                                               <  1:2476625/65‑1 (MQ=255)

      atccTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTg                                               <  1:458776/65‑1 (MQ=255)

      atccTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTg                                               <  1:2520404/65‑1 (MQ=255)

      atccTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTg                                               <  1:2852222/65‑1 (MQ=255)

      atccTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTg                                               <  1:1612992/65‑1 (MQ=255)

                    aaaCGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTg                                               <  1:628037/54‑1 (MQ=255)

                    aaaCGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCGTTGGCg                                                    <  1:1013738/49‑1 (MQ=255)

                       cGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCATTTTGCAGTAACAATATTTGCAg                             <  1:2188125/69‑1 (MQ=255)

                       cGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCATTTTGCAGTAACAATATTTGCAg                             <  1:1659183/69‑1 (MQ=255)

                                       gCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAGGGCa             >  1:2836595/1‑69 (MQ=255)

                                       gCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAGGGCa             >  1:3350643/1‑69 (MQ=255)

                                          gCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAGGGCAaat          >  1:1824074/1‑69 (MQ=255)

                                           cACCTAATAACTCTCCGCCTTTGGCATTTTGCAGTAACaa                                      >  1:2596149/1‑40 (MQ=38)

                                           cACCTAATAACTCTCCGCCTTTGGCATTTTGCAGTAACAATATTTGCGGTGGTGCTGAAAGGGCAaat          >  1:2559719/1‑68 (MQ=255)

                                           cACCTAATAACTCTCCGCCTTTGGCATTTTGCAGTAACAATATTTGCAGTGGTGCTGAAAGGGGAaata         >  1:3210773/1‑69 (MQ=255)

                                           cACCTAATAACTCTCCGCCTTTGGCATTTTGCAGTAACAATATTTGCAGTGGTGCTGAAAGGGCAaata         >  1:1725982/1‑69 (MQ=255)

                                           cACCTAATAACTCTCCGCCTTTGGCATTTTGCAGTAACAATATTTGCAGTGGTGCTGAAAGGGCAaata         >  1:834407/1‑69 (MQ=255)

                                           cACCTAATAACTCTCCGCCTTTGGCATTTTGCAGTAACAATATTTGCAGTGGTGCTGAAAGGGCAaata         >  1:101816/1‑69 (MQ=255)

                                           cACCTAATAACTCTCCGCCTTTGGCATTTTGCAGTAACAATATTTGCAGTGGTGCTGAAAGGGCAaata         >  1:1386159/1‑69 (MQ=255)

                                           cACCTAATAACTCTCCGCCTTTGGCATTTTGCAGTAACAATATTTGCAGTGGTGCTGAAAGGGCAaata         >  1:1415425/1‑69 (MQ=255)

                                           cACCTAATAACTCTCCGCCTTTGGCATTTTGCAGTAACAATATTTGCAGTGGTGCTGAAAGGGCAaata         >  1:399514/1‑69 (MQ=255)

                                           cACCTAATAACTCTCCGCCTTTGGCATTTTGCAGTAACAATATTTGCAGTGGTGCTGAAAGGGCAaata         >  1:3306532/1‑69 (MQ=255)

                                           cACCTAATAACTCTCCGCCTTTGGCATTTTGCAGTAACAATATTTGCAGTGGTGCTGAAAGGGCAaata         >  1:2147928/1‑69 (MQ=255)

                                           cACCTAATAACTCTCCGCCTTTGGCATTTTGCAGTAACAATATTTGCAGTGGTGCTGAAAGGGCAaata         >  1:3044471/1‑69 (MQ=255)

                                           cACCTAATAACTCTCCGCCTTTGGCATTTTGCAGTAACAATATTTGCAGTGGTGCTGAAAGGGCAaat          >  1:2696679/1‑68 (MQ=255)

                                           cACCTAATAACTCTCCGCCTTTGGCATTTTGCAGTAACAATATTTGCAGTGGTGCTGAAAGGGCAaat          >  1:895045/1‑68 (MQ=255)

                                                   aaCTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAGGGCAAATAATCCcgc  >  1:827470/1‑68 (MQ=255)

                                                   aaCTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAGGGCAAATAATCCcgc  >  1:2865873/1‑68 (MQ=255)

CCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAGGGCAAATAATCCCGC  >  minE/239932‑240050

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