BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	303,157	C→A	58.6%	L158L (CTG→CTT)	tesB ←	acyl‑CoA thioesterase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	303,157	0	C	A	58.6%	16.9 / 55.7	58	L158L (CTG→CTT)	tesB	acyl‑CoA thioesterase II
Reads supporting (aligned to +/- strand): ref base C (9/15); new base A (8/26); total (17/41)
Fisher's exact test for biased strand distribution p-value = 3.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGCCAGCGATTGGGCG  >  minE/303099‑303223

gATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACg                                                            <  1:621875/67‑1 (MQ=255)

      gTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCa                                                      >  1:18862/1‑67 (MQ=255)

      gTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAg                                                     >  1:440437/1‑68 (MQ=255)

      gTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAg                                                     >  1:1992219/1‑68 (MQ=255)

      gTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAg                                                     >  1:2162492/1‑68 (MQ=255)

      gTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAg                                                     >  1:2197203/1‑68 (MQ=255)

      gTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAg                                                     >  1:2713285/1‑68 (MQ=255)

      gTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAg                                                     >  1:2939933/1‑68 (MQ=255)

         cTGCGACGTGCCCTTTCAGTGGGTTATTAAACTCCACCGGACGGACTTCAAGCGGACGATCGCtgatga                                                 <  1:2464234/69‑1 (MQ=255)

         cTGCGACGTGCCCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCAAGCGGACGATCGCtgatga                                                 <  1:2899351/69‑1 (MQ=255)

         cTGCGACGTGCCCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCAAGCGGACGATCGCtgatga                                                 <  1:2698076/69‑1 (MQ=255)

         cTGCGACGTGCCCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCAAGCGGACGATCGCtgatga                                                 <  1:1024596/69‑1 (MQ=255)

         cTGCGACGTGCCCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCAAGCGGACGATCGCtgatga                                                 <  1:580510/69‑1 (MQ=255)

         cTGCGACGTGCCCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCAAGCGGACGATCGCtgatga                                                 <  1:2227090/69‑1 (MQ=255)

         cTGCGACGTGCCCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCAAGCGGACGATCGCtgatga                                                 <  1:600358/69‑1 (MQ=255)

         cTGCGACGTGCCCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCAAGCGGACGATCGCtgatga                                                 <  1:3022640/69‑1 (MQ=255)

         cTGCGACGTGCCCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCAAGCGGACGATCGCtgatga                                                 <  1:3058459/69‑1 (MQ=255)

         cTGCGACGTGCCCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCAAGCGGACGATCGCtgatga                                                 <  1:3135025/69‑1 (MQ=255)

         cTGCGACGTGCCCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCAAGCGGACGATCGCtgatga                                                 <  1:2010637/69‑1 (MQ=255)

         cTGCGACGTGCCCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCAAGCGGACGATCGCtgatga                                                 <  1:2009699/69‑1 (MQ=255)

         cTGCGACGTGCCCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCAAGCGGACGATCGCtgatga                                                 <  1:3395823/69‑1 (MQ=255)

         cTGCGACGTGCCCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCAAGCGGACGATCGCtgatga                                                 <  1:1945307/69‑1 (MQ=255)

         cTGCGACGTGCCCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCAAGCGGACGATCGCtgatga                                                 <  1:190695/69‑1 (MQ=255)

         cTGCGACGTGCCCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCAAGCGGACGATCGCtgatga                                                 <  1:1786706/69‑1 (MQ=255)

         cTGCGACGTGCCCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCAAGCGGACGATCGCtgatga                                                 <  1:1326619/69‑1 (MQ=255)

         cTGCGACGTGCCCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCAAGCGGACGATCGCtgatga                                                 <  1:1962215/69‑1 (MQ=255)

         cTGCGACGTGCCCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCAAGCGGACGATCGCtgatga                                                 <  1:69640/69‑1 (MQ=255)

         cTGCGACGTGCCCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCAAGCGGACGATCGCtgatga                                                 <  1:966729/69‑1 (MQ=255)

         cTGCGACGTGCCCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCAAGCGGACGATCGCtgatga                                                 <  1:1597548/69‑1 (MQ=255)

         cTGCGACGTGCCCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCAAGCGGACGATCGCtgatga                                                 <  1:1295102/69‑1 (MQ=255)

         cTGCGACGTGCCCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCAAGCGGACGATCGCtgatga                                                 <  1:1418079/69‑1 (MQ=255)

         cTGCGACGTGCCCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCAAGCGGACGATCGCtgatga                                                 <  1:99604/69‑1 (MQ=255)

         cTGCGACGTGCCCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCAAGCGGACGATCGCtgatga                                                 <  1:1482255/69‑1 (MQ=255)

                       ttCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGc                                   <  1:1516111/69‑1 (MQ=255)

                                  aTGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTggcggc                          >  1:2347122/1‑67 (MQ=255)

                                  aTGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTggcggc                          >  1:2771926/1‑67 (MQ=255)

                                   tGAAACTCCACCGGACGGACTTCAAGCGGACGATCGCtgatga                                                 <  1:2140388/43‑1 (MQ=38)

                                       aCTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGc                                   <  1:1298239/53‑1 (MQ=255)

                                             ccGGACGGACTTCAAGCGGACGATCGCTGATGAATTTATCTTTCAGCACTGGCGgcag                        >  1:2220326/1‑58 (MQ=39)

                                             ccGGACGGACTTCAAGCGGACGATCGCTGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGcc              >  1:3047373/1‑68 (MQ=255)

                                             ccGGACGGACTTCAAGCGGACGATCGCTGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGCCa             >  1:1117397/1‑69 (MQ=255)

                                             ccGGACGGACTTCAAGCGGACGATCGCTGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGCCa             >  1:1089763/1‑69 (MQ=255)

                                             ccGGACGGACTTCAAGCGGACGATCGCTGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGCCa             >  1:325049/1‑69 (MQ=255)

                                             ccGGACGGACTTCAAGCGGACGATCGCTGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGCCa             >  1:3391330/1‑69 (MQ=255)

                                             ccGGACGGACTTCAAGCGGACGATCGCTGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGCCa             >  1:2788906/1‑69 (MQ=255)

                                             ccGGACGGACTTCAAGCGGACGATCGCTGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGCCa             >  1:2573516/1‑69 (MQ=255)

                                                gactgactTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGCCAGCg          <  1:1693241/69‑1 (MQ=255)

                                                gacggacTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGCCAGCg          <  1:1041422/69‑1 (MQ=255)

                                                gacggacTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGCCAGCg          <  1:719462/69‑1 (MQ=255)

                                                gacggacTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGCCAGCg          <  1:3235781/69‑1 (MQ=255)

                                                gacggacTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGCCAGCg          <  1:2139773/69‑1 (MQ=255)

                                                gacggacTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGCCAGCg          <  1:1680467/69‑1 (MQ=255)

                                                gacggacTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGCCAGCg          <  1:1509822/69‑1 (MQ=255)

                                                gacggacTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGCCAGCg          <  1:1333404/69‑1 (MQ=255)

                                                gacggacTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGCCAGCg          <  1:1270593/69‑1 (MQ=255)

                                                        tCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGCCAGCGATTGGGCg  <  1:380533/69‑1 (MQ=255)

                                                        tCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGCCAGCGATTGGGCg  <  1:2390739/69‑1 (MQ=255)

                                                        tCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGCCAGCGATTGGGCg  <  1:2223750/69‑1 (MQ=255)

                                                         ccAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGCCAGCGATTGGGCg  <  1:2416919/68‑1 (MQ=255)

GATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGCCAGCGATTGGGCG  >  minE/303099‑303223

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