BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,506,604	C→T	62.1%	R245R (CGG→CGA)	yfbS ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,506,604	0	C	T	62.1%	56.4 / 71.3	108	R245R (CGG→CGA)	yfbS	predicted transporter
Reads supporting (aligned to +/- strand): ref base C (24/17); new base T (45/22); total (69/39)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

Only 100 of 108 total aligned reads displayed.

AGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCGATA  >  minE/1506546‑1506669

aGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCaa                                                         <  1:2450911/69‑1 (MQ=255)

aGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCaa                                                         <  1:3325288/69‑1 (MQ=255)

aGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCaa                                                         <  1:1944691/69‑1 (MQ=255)

aGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCaa                                                         <  1:2846184/69‑1 (MQ=255)

aGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCaa                                                         <  1:1562906/69‑1 (MQ=255)

aGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCaa                                                         <  1:1353016/69‑1 (MQ=255)

 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACTTAATTTGAGATCGTCCAGTCGTTGACCAAt                                                        >  1:2404824/1‑69 (MQ=255)

 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTg                                                              >  1:94951/1‑63 (MQ=255)

 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGCCCAAt                                                        >  1:1492457/1‑69 (MQ=255)

 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCaa                                                         >  1:1265655/1‑68 (MQ=255)

 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                                        >  1:22446/1‑69 (MQ=255)

 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                                        >  1:3224354/1‑69 (MQ=255)

 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                                        >  1:3059620/1‑69 (MQ=255)

 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                                        >  1:201513/1‑69 (MQ=255)

 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                                        >  1:1996005/1‑69 (MQ=255)

 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                                        >  1:246128/1‑69 (MQ=255)

 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                                        >  1:345296/1‑69 (MQ=255)

 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                                        >  1:411442/1‑69 (MQ=255)

 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                                        >  1:1902827/1‑69 (MQ=255)

 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                                        >  1:891384/1‑69 (MQ=255)

 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                                        >  1:3018104/1‑69 (MQ=255)

 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                                        >  1:3355321/1‑69 (MQ=255)

 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                                        >  1:624061/1‑69 (MQ=255)

 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                                        >  1:1664365/1‑69 (MQ=255)

 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                                        >  1:2589399/1‑69 (MQ=255)

 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                                        >  1:2403023/1‑69 (MQ=255)

 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                                        >  1:818378/1‑69 (MQ=255)

 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                                        >  1:687173/1‑69 (MQ=255)

 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                                        >  1:757971/1‑69 (MQ=255)

 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                                        >  1:770109/1‑69 (MQ=255)

 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                                        >  1:814195/1‑69 (MQ=255)

 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                                        >  1:895347/1‑69 (MQ=255)

 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                                        >  1:2546606/1‑69 (MQ=255)

 gCGTTCAACACCGATGACGTTAGCACCATAACGCTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                                        >  1:2421059/1‑69 (MQ=255)

 gCGTTCAACACAGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                                        >  1:1746644/1‑69 (MQ=255)

          acCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGa                                               <  1:134659/69‑1 (MQ=255)

             gATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATcc                                            <  1:1590121/69‑1 (MQ=255)

              aTGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATcc                                            <  1:626901/68‑1 (MQ=255)

                gACGTTTGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCgg                                          >  1:441444/1‑68 (MQ=255)

                gACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCaa                                                         <  1:1243018/53‑1 (MQ=255)

                gACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCggg                                         >  1:2054223/1‑69 (MQ=255)

                gACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCggg                                         >  1:597353/1‑69 (MQ=255)

                     tAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGaa                                    <  1:2928920/69‑1 (MQ=255)

                     tAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGaa                                    <  1:3046915/69‑1 (MQ=255)

                     tAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGaa                                    <  1:2924991/69‑1 (MQ=255)

                     tAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGaa                                    <  1:2425349/69‑1 (MQ=255)

                     tAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGaa                                    <  1:236819/69‑1 (MQ=255)

                     tAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGaa                                    <  1:2084769/69‑1 (MQ=255)

                     tAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGaa                                    <  1:2052085/69‑1 (MQ=255)

                     tAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGaa                                    <  1:3380156/69‑1 (MQ=255)

                     tAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGaa                                    <  1:526215/69‑1 (MQ=255)

                     tAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGaa                                    <  1:1138439/69‑1 (MQ=255)

                          ccATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCg                                                >  1:86768/1‑52 (MQ=255)

                          ccATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATa                                  >  1:3336402/1‑66 (MQ=255)

                          ccATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATa                                  >  1:3269526/1‑66 (MQ=255)

                          ccATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATa                                  >  1:3372213/1‑66 (MQ=255)

                          ccATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATa                                  >  1:314851/1‑66 (MQ=255)

                          ccATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATa                                  >  1:2678373/1‑66 (MQ=255)

                          ccATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATa                                  >  1:1022097/1‑66 (MQ=255)

                          ccATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATa                                  >  1:277893/1‑66 (MQ=255)

                          ccATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATa                                  >  1:1228404/1‑66 (MQ=255)

                          ccATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATa                                  >  1:1377213/1‑66 (MQ=255)

                          ccATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATa                                  >  1:1482271/1‑66 (MQ=255)

                          ccATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATa                                  >  1:1000889/1‑66 (MQ=255)

                          ccATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATa                                  >  1:2946549/1‑66 (MQ=255)

                           cATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCa                              <  1:428290/69‑1 (MQ=255)

                           cATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCa                              <  1:10652/69‑1 (MQ=255)

                           cATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCa                              <  1:458030/69‑1 (MQ=255)

                            aTAACGCTCACGTAATTTGAGATCATCCAGCCTTTGACCAATCATTGGCGATCCGGGGCGAATAGCCa                              <  1:741905/68‑1 (MQ=255)

                               aCGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGAc                           >  1:1483738/1‑68 (MQ=255)

                               aCGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACg                          >  1:2418544/1‑69 (MQ=255)

                               aCGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACg                          >  1:2832073/1‑69 (MQ=255)

                               aCGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACg                          >  1:1218260/1‑69 (MQ=255)

                               aCGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACg                          >  1:801023/1‑69 (MQ=255)

                               aCGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACg                          >  1:1701211/1‑69 (MQ=255)

                               aCGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACg                          >  1:450307/1‑69 (MQ=255)

                                       gTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGTGGCGAATAGCCAGACGTCGCGCa                   <  1:3386470/68‑1 (MQ=255)

                                       gTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCa                   <  1:2125490/68‑1 (MQ=255)

                                       gTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCa                   <  1:103334/68‑1 (MQ=255)

                                       gTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCa                   <  1:1615824/68‑1 (MQ=255)

                                       gTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCa                   <  1:586569/68‑1 (MQ=255)

                                       gTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCa                   <  1:1774913/68‑1 (MQ=255)

                                       gTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCa                   <  1:1954198/68‑1 (MQ=255)

                                       gTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCa                   <  1:394952/68‑1 (MQ=255)

                                       gTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCa                   <  1:379499/68‑1 (MQ=255)

                                       gTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCa                   <  1:3396368/68‑1 (MQ=255)

                                             tGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACg                 >  1:1875655/1‑64 (MQ=255)

                                             tGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACg                 >  1:2072737/1‑64 (MQ=255)

                                             tGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACg                 >  1:2718964/1‑64 (MQ=255)

                                             tGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCg             >  1:2593758/1‑68 (MQ=255)

                                                    atcCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCg     <  1:2877431/69‑1 (MQ=255)

                                                    atcCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCg     <  1:823289/69‑1 (MQ=255)

                                                    atcCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCg     <  1:3204901/69‑1 (MQ=255)

                                                       cAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCGata  >  1:1976367/1‑69 (MQ=255)

                                                       cAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCGata  >  1:1984726/1‑69 (MQ=255)

                                                       cAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCGata  >  1:1671586/1‑69 (MQ=255)

                                                       cAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCGata  >  1:2648942/1‑69 (MQ=255)

                                                       cAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCGat   >  1:3189581/1‑68 (MQ=255)

                                                       cAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCGa    >  1:1361532/1‑67 (MQ=255)

                                                       cAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCAGTCAGTCGa    >  1:33868/1‑67 (MQ=255)

AGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCGATA  >  minE/1506546‑1506669

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