BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	190,956	C→T	49.1%	D85D (GAC→GAT)	yaeR →	predicted lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	190,956	0	C	T	49.1%	‑1.4 / 46.8	55	D85D (GAC→GAT)	yaeR	predicted lyase
Reads supporting (aligned to +/- strand): ref base C (2/26); new base T (12/15); total (14/41)
Fisher's exact test for biased strand distribution p-value = 1.88e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

CGCCGGAACGACCCAGCCGACCGGAAGCTTGCGGTCTGCGTCATCTGGCTTTTAGCGTTGATGACATCGATGCGGCAGTGGCGCACCTTGAAAGCCATAACGTGAAGTGTGAAACCATCCGTGTCGATCCA  >  minE/190892‑191022

cgccgGAGCGCCCCAGCCGACCGGAAGCTTGCGGCCTGCGTCATCTGGCTTTTAGCGTTGATGATATc                                                                 <  1:1195042/68‑1 (MQ=255)

cgccgGAGCGCCCCAGCCGACCGGAAGCTTGCGGCCTGCGTCATCTGGCTTTTAGCGTTGATGATATc                                                                 <  1:3226592/68‑1 (MQ=255)

cgccgGAGCGCCCCAGCCGACCGGAAGCTTGCGGCCTGCGTCATCTGGCTTTTAGCGTTGATGATATc                                                                 <  1:3280488/68‑1 (MQ=255)

cgccgGAGCGCCCCAGCCGACCGGAAGCTTGCGGCCTGCGTCATCTGGCTTTTAGCGTTGATGATATc                                                                 <  1:2840959/68‑1 (MQ=255)

cgccgGAGCGCCCCAGCCGACCGGAAGCTTGCGGCCTGCGTCATCTGGCTTTTAGCGTTGATGATATc                                                                 <  1:25280/68‑1 (MQ=255)

cgccgGAGCGCCCCAGCCGACCGGAAGCTTGCGGCCTGCGTCATCTGGCTTTTAGCGTTGATGATATc                                                                 <  1:65517/68‑1 (MQ=255)

cgccgGAGCGCCCCAGCCGACCGGAAGCTTGCGGCCTGCGTCATCTGGCTTTTAGCGTTGATGATATc                                                                 <  1:2189622/68‑1 (MQ=255)

cgccgGAGCGCCCCAGCCGACCGGAAGCTTGCGGCCTGCGTCATCTGGCTTTTAGCGTTGATGATATc                                                                 <  1:840893/68‑1 (MQ=255)

cgccgGAGCGCCCCAGCCGACCGGAAGCTTGCGGCCTCCGTCATCTGGCTTTTAGCGTTGATGATATc                                                                 <  1:2626109/68‑1 (MQ=255)

 gccgGAACGACCCAGCCGACCGGAAGCTTGCGGTCTGCGTCATCTGGCTTTTAGCGTTGATGACATCGa                                                               <  1:2205422/69‑1 (MQ=255)

 gccgGAACGACCCAGCCGACCGGAAGCTTGCGGTCTGCGTCATCTGGCTTTTAGCGTTGATGACATCGa                                                               <  1:1019254/69‑1 (MQ=255)

 gccgGAACGACCCAGCCGACCGGAAGCTTGCGGTCTGCGTCATCTGGCTTTTAGCGTTGATGACATCGa                                                               <  1:389510/69‑1 (MQ=255)

 gccgGAACGACCCAGCCGACCGGAAGCTTGCGGTCTGCGTCATCTGGCTTTTAGCGTTGATGACATCGa                                                               <  1:691096/69‑1 (MQ=255)

 gccgGAACGACCCAGCCGACCGGAAGCTTGCGGTCTGCGTCATCTGGCTTTTAGCGTTGATGACATCGa                                                               <  1:2016495/69‑1 (MQ=255)

 gccgGAACGACCCAGCCGACCGGAAGCTTGCGGTCTGCGTCATCTGGCCTTTAGCGTTGATGACATCGa                                                               <  1:2997038/69‑1 (MQ=255)

       aCGACCCAGCCGACCGGAAGCTTGCGGTCTGCGTCATCTGGCTTTTAGCGTTGATGACATCGATGCGGc                                                         <  1:2098913/69‑1 (MQ=255)

       aCGACCCAGCCGACCGGAAGCTTGCGGTCTGCGTCATCTGGCTTTTAGCGTTGATGACATCGATGCGGc                                                         <  1:2129415/69‑1 (MQ=255)

       aCGACCCAGCCGACCGGAAGCTTGCGGTCTGCGTCATCTGGCTTTTAGCGTTGATGACATCGATGCGGc                                                         <  1:840925/69‑1 (MQ=255)

       aCGACCCAGCCGACCGGAAGCTTGCGGTCTGCGTCATCTGGCTTTTAGCGTTGATGACATCGATGCGGc                                                         <  1:2262509/69‑1 (MQ=255)

       aCGACCCAGCCGACCGGAAGCTTGCGGTCTGCGTCATCTGGCTTTTAGCGTTGATGACATCGATGCGGc                                                         <  1:946921/69‑1 (MQ=255)

       aCGACCCAGCCGACCGGAAGCTTGCGGTCTGCGTCATCTGGCTTTTAGCGTTGATGACATCGATGCGGc                                                         <  1:2569804/69‑1 (MQ=255)

       aCGACCCAGCCGACCGGAAGCTTGCGGTCTGCGTCATCTGGCTTTTAGCGTTGATGACATCGATGCGGc                                                         <  1:2624688/69‑1 (MQ=255)

       aCGACCCAGCCGACCGGAAGCTTGCGGTCTGCGTCATCTGGCTTTTAGCGTTGATGACATCGATGCGGc                                                         <  1:745444/69‑1 (MQ=255)

       aCGACCCAGCCGACCGGAAGCTTGCGGTCTGCGTCATCTGGCTTTTAGCGTTGATGACATCGATGCGGc                                                         <  1:2824731/69‑1 (MQ=255)

       aCGACCCAGCCGACCGGAAGCTTGCGGTCTGCGTCATCTGGCTTTTAGCGTTGATGACATCGATGCGGc                                                         <  1:1224589/69‑1 (MQ=255)

        cGACCCAGCCGACCGGAAGCTTGCGGTCTGCGTCATCTGGCTTTTAGCGTTGATGACATCGATGCGGc                                                         <  1:871671/68‑1 (MQ=255)

             cAGCCGACCGGAAGCTTGCGGCCTGCGTCATCTGGCTTTTAGCGTTGATGATTTCGATGCGGCAGTGgc                                                   <  1:2198116/69‑1 (MQ=255)

             cAGCCGACCGGAAGCTTGCGGCCTGCGTCATCTGGCTTTTAGCGTTGATGATATCGATGCGGCAGTGgc                                                   <  1:1996575/69‑1 (MQ=255)

             cAGCCGACCGGAAGCTTGCGGCCTGCGTCATCTGGCTTTTAGCGTTGATGATATCGATGCGGCAGTGgc                                                   <  1:372005/69‑1 (MQ=255)

             cAGCCGACCGGAAGCTTGCGGCCTGCGTCATCTGGCTTTTAGCGTTGATGATATCGATGCGGCAGTGgc                                                   <  1:2971666/69‑1 (MQ=255)

             cAGCCGACCGGAAGCTTGCGGCCTGCGTCATCTGGCTTTTAGCGTTGATGATATCGATGCGGCAGTGgc                                                   <  1:2898576/69‑1 (MQ=255)

              aGCCGACCGGAAGCTTGCGGTCTGCGTCATCTGGCTTTTAGCGTTGATGACATCGATGCGGCAGTGgcg                                                  <  1:587801/69‑1 (MQ=255)

              aGCCGACCGGAAGCTTGCGGTCTGCGTCATCTGGCTTTTAGCGTTGATGACATCGATGCGGCAGTGgcg                                                  <  1:955483/69‑1 (MQ=255)

               gCCGACCGGAAGCTTGCGGTCTGCGTCATCTGGCTTTTAGCGTTGATGACATCGATGCGGCAGTGgcgc                                                 <  1:3028594/69‑1 (MQ=255)

               gCCGACCGGAAGCTTGCGGTCTGCGTCATCTGGCTTTTAGCGTTGATGACATCGATGCGGCAGTGgcgc                                                 <  1:362786/69‑1 (MQ=255)

                  gACCGGAAGCTTGCGGCCTGCGTCATCTGGCTTTTAGCGTTGATGATATc                                                                 <  1:630122/50‑1 (MQ=255)

                        aaGCTTGCGGTCTGCGTCATCTGGCTTTTAGCGTTGATGACATCGATGCGGCAGTGGCGCACCTTGaaa                                        <  1:1254129/69‑1 (MQ=255)

                          gCTTGCGGTCTGCGTCATCTGGCTTTTAGCGTTGATGACATCGATGCGGCAGTGGCGCACCTTGAAAGc                                      <  1:565524/69‑1 (MQ=255)

                               cGGTCTGCGTCATCTGGCTTTTAGCGTTGATGACATCGATGCGGc                                                         <  1:206843/45‑1 (MQ=255)

                                                cTTTTAGCGTTGATGATATCGATGCGGGAGTGGCGCACCTTGAAAGCCATAACGTGAAGTGTGAAGCCa                >  1:1251550/1‑69 (MQ=255)

                                                cTTTTAGCGTTGATGATATCGATGCGGCAGTGGCGCACCTTGAAAGCCATaa                                 >  1:1975058/1‑52 (MQ=255)

                                                cTTTTAGCGTTGATGATATCGATGCGGCAGTGGCGCACCTTGAAAGCCATAACGTGAAGTGTGAAGCCa                >  1:3137874/1‑69 (MQ=255)

                                                cTTTTAGCGTTGATGATATCGATGCGGCAGTGGCGCACCTTGAAAGCCATAACGTGAAGTGTGAAGCCa                >  1:1524591/1‑69 (MQ=255)

                                                cTTTTAGCGTTGATGATATCGATGCGGCAGTGGCGCACCTTGAAAGCCATAACGTGAAGTGTGAAGCCa                >  1:1942527/1‑69 (MQ=255)

                                                cTTTTAGCGTTGATGATATCGATGCGGCAGTGGCGCACCTTGAAAGCCATAACGTGAAGTGTGAAGCCa                >  1:2045138/1‑69 (MQ=255)

                                                cTTTTAGCGTTGATGATATCGATGCGGCAGTGGCGCACCTTGAAAGCCATAACGTGAAGTGTGAAGCCa                >  1:2394013/1‑69 (MQ=255)

                                                cTTTTAGCGTTGATGATATCGATGCGGCAGTGGCGCACCTTGAAAGCCATAACGTGAAGTGTGAAGCCa                >  1:476454/1‑69 (MQ=255)

                                                cTTTTAGCGTTGATGATATCGATGCGGCAGTGGCGCACCTTGAAAGCCATAACGTGAAGTGTGAAGCCa                >  1:2955939/1‑69 (MQ=255)

                                                cTTTTAGCGTTGATGATATCGATGCGGCAGTGGCGCACCTTGAAAGCCATAACGTGAAGTGTGAAGCCa                >  1:2978161/1‑69 (MQ=255)

                                                cTTTTAGCGTTGATGATATCGATGCGGCAGTGGCGCACCTTGAAAGCCATAACGTGAAGTGTGAAGCCa                >  1:3405758/1‑69 (MQ=255)

                                                cTTTTAGCGTTGATGATATCGATGCGGCAGTGGCGCACCTTGAAAGCCATAACGTGAAGTGTGAAGCCa                >  1:3237859/1‑69 (MQ=255)

                                                           gatgaCATCGATGCGGCAGTGGCGCACCTTGAAAGc                                      <  1:836469/36‑1 (MQ=255)

                                                           gatgaCATCGATGCGGCAGTGGCGCACCTTGAAAGc                                      <  1:3219404/36‑1 (MQ=255)

                                                               aCATCGATGCGGCAGTGGCGCACCTTGAAAGCCATAACGTGAAGTGTGAAACCATCCGTGTCGATCCa  >  1:286218/1‑68 (MQ=255)

                                                               aCATCGATGCGGCAGTGGCGCACCTTGAAAGCCATAACGTGAAGTGTGAAACCATCCGTGTCGATCCa  >  1:492537/1‑68 (MQ=255)

CGCCGGAACGACCCAGCCGACCGGAAGCTTGCGGTCTGCGTCATCTGGCTTTTAGCGTTGATGACATCGATGCGGCAGTGGCGCACCTTGAAAGCCATAACGTGAAGTGTGAAACCATCCGTGTCGATCCA  >  minE/190892‑191022

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