BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,615,887	C→T	13.8%	G144G (GGC→GGT)	yhhS →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,615,887	0	C	T	13.8%	127.6 / 7.6	71	G144G (GGC→GGT)	yhhS	predicted transporter
Reads supporting (aligned to +/- strand): ref base C (18/43); new base T (6/4); total (24/47)
Fisher's exact test for biased strand distribution p-value = 7.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATCCTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTCACTTACGGGGCGATGGCGATG  >  minE/2615820‑2615950

tCATCCTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCt                                                                >  1:2095460/1‑69 (MQ=255)

tCATCCTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCt                                                                >  1:2760799/1‑69 (MQ=255)

tCATCCTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCt                                                                >  1:1156380/1‑69 (MQ=255)

  aTCCTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCg                                                              >  1:2750475/1‑69 (MQ=255)

  aTCCTTGGGATTGGACAAAGTTTTGCCGGAACGGGATCGACCCTGTGGGGTGTTGGCGTGGTTGGTTCg                                                              <  1:1065998/69‑1 (MQ=255)

  aTCCTTGGGATTGGACAAAGTTTTGCCGGAACGGGATCGACCCTGTGGGGTGTTGGCGTGGTTGGTTCg                                                              <  1:439600/69‑1 (MQ=255)

  aTCCTTGGGATTGGACAAAGTTTTGCCGGAACGGGATCGACCCTGTGGGGTGTTGGCGTGGTTGGTTCg                                                              <  1:1103167/69‑1 (MQ=255)

     cTTGGGATTGGGCAAAGTTTTGCCGGAACTGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTg                                                           <  1:780791/69‑1 (MQ=255)

     cTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTg                                                           <  1:602091/69‑1 (MQ=255)

     cTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTg                                                           <  1:73154/69‑1 (MQ=255)

     cTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTg                                                           <  1:3364955/69‑1 (MQ=255)

     cTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTg                                                           <  1:3094471/69‑1 (MQ=255)

     cTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTg                                                           <  1:1278747/69‑1 (MQ=255)

     cTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTg                                                           <  1:3019057/69‑1 (MQ=255)

     cTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTg                                                           <  1:3013725/69‑1 (MQ=255)

     cTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTg                                                           <  1:2867951/69‑1 (MQ=255)

     cTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTg                                                           <  1:2115302/69‑1 (MQ=255)

                   aaGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGt                                             >  1:82979/1‑69 (MQ=255)

                    aGTTTTGCCGGAACGGGATCGACCCTGTGGGGTGTTGGCGTGGTTGGTTCGCTACATATCGGGCGGGTg                                            >  1:1801836/1‑69 (MQ=255)

                    aGTTTTGCCGGAACGGGATCGACCCTGTGGGGTGTTGGCGTGGTTGGTTCGCTACATATCGGGCGGGTg                                            >  1:282370/1‑69 (MQ=255)

                    aGTTTTGCCGGAACGGGATCGACCCTGTGGGGTGTTGGCGTGGTTGGTTCGCTACATATCGGGCGGGTg                                            >  1:2633798/1‑69 (MQ=255)

                    aGTTTTGCCGGAACGGGATCGACCCTGTGGGGTGTTGGCGTGGTTGGTTCGCTACATATCGGGCGGGTg                                            >  1:257224/1‑69 (MQ=255)

                    aGTTTTGCCGGAACGGGATCGACCCTGTGGGGTGTTGGCGTGGTTGGTTCGCTACATATCGGGCGGGTg                                            >  1:2463124/1‑69 (MQ=255)

                    aGTTTTGCCGGAACGGGATCGACCCTGTGGGGTGTTGGCGTGGTTGGTTCGCTACATATCGGGCGGGTg                                            >  1:685470/1‑69 (MQ=255)

                              gAACGGGATCGACCCTGTGGGGTGTTGGCGTGGTTGGTTCGCTACATATCGGGCGGGTGATTTCGTGGa                                  <  1:2316723/69‑1 (MQ=255)

                                   ggATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGTAACGGc                             <  1:1180184/69‑1 (MQ=255)

                                   ggATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGc                             <  1:592903/69‑1 (MQ=255)

                                   ggATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGc                             <  1:548287/69‑1 (MQ=255)

                                   ggATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGc                             <  1:451382/69‑1 (MQ=255)

                                   ggATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGc                             <  1:1706325/69‑1 (MQ=255)

                                   ggATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGc                             <  1:2228810/69‑1 (MQ=255)

                                   ggATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGc                             <  1:3011885/69‑1 (MQ=255)

                                   ggATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGc                             <  1:1862796/69‑1 (MQ=255)

                                   ggATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGc                             <  1:2516071/69‑1 (MQ=255)

                                   ggATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGc                             <  1:112993/69‑1 (MQ=255)

                                      tCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTg                                                           <  1:683493/36‑1 (MQ=255)

                                        gACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGt                        <  1:386341/69‑1 (MQ=255)

                                        gACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGt                        <  1:1106214/69‑1 (MQ=255)

                                        gACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGt                        <  1:3234041/69‑1 (MQ=255)

                                        gACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGt                        <  1:3361016/69‑1 (MQ=255)

                                        gACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGt                        <  1:382207/69‑1 (MQ=255)

                                         aCCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGt                        <  1:1828730/68‑1 (MQ=255)

                                          cccTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTc                       <  1:1698719/68‑1 (MQ=255)

                                          cccTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTc                       <  1:1536501/68‑1 (MQ=255)

                                          cccTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTc                       <  1:238828/68‑1 (MQ=255)

                                          cccTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTc                       <  1:2771600/68‑1 (MQ=255)

                                            cTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTCACt                    <  1:2772372/69‑1 (MQ=255)

                                            cTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTCACt                    <  1:2028655/69‑1 (MQ=255)

                                            cTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTCACt                    <  1:127259/69‑1 (MQ=255)

                                            cTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTCACt                    <  1:2439204/69‑1 (MQ=255)

                                               tGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTc                       <  1:1222961/63‑1 (MQ=255)

                                               tGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTc                       <  1:1853432/63‑1 (MQ=255)

                                               tGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTc                       <  1:1873764/63‑1 (MQ=255)

                                               tGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTc                       <  1:438612/63‑1 (MQ=255)

                                               tGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTc                       <  1:2626004/63‑1 (MQ=255)

                                               tGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTc                       <  1:1430471/63‑1 (MQ=255)

                                               tGGGGCGTCGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTc                       <  1:2046949/63‑1 (MQ=255)

                                                            tGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTCACTTACGGggcgatggcga    >  1:1854688/1‑69 (MQ=255)

                                                            tGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTCACTTACGGggcgatggcga    >  1:1599505/1‑69 (MQ=255)

                                                            tGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTCACTTACGGggcgatggcga    >  1:1879307/1‑69 (MQ=255)

                                                            tGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTCACTTACGGggcgatggcga    >  1:2965952/1‑69 (MQ=255)

                                                            tGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTCACTTACGGggcgatggcga    >  1:2884667/1‑69 (MQ=255)

                                                            tGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTCACTTACGGggcgatggcga    >  1:1149445/1‑69 (MQ=255)

                                                             ggTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGt                        <  1:2582011/48‑1 (MQ=255)

                                                              gTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTCACTTACg                >  1:1622005/1‑55 (MQ=255)

                                                              gTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTCACTTACGGGgcgatggcgatg  >  1:1918806/1‑69 (MQ=255)

                                                              gTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTCACTTACGGGgcgatggcgatg  >  1:658647/1‑69 (MQ=255)

                                                              gTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTCACTTACGGGgcgatggcgatg  >  1:3110806/1‑69 (MQ=255)

                                                              gTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTCACTTACGGGgcgatggcgatg  >  1:3167697/1‑69 (MQ=255)

                                                              gTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTCACTTACGGGgcgatggcgatg  >  1:1618145/1‑69 (MQ=255)

                                                              gTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTCACTTACGGGgcgatggcgatg  >  1:1042250/1‑69 (MQ=255)

TCATCCTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTCACTTACGGGGCGATGGCGATG  >  minE/2615820‑2615950

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