BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,615,866	A→G	19.5%	L137L (CTA→CTG)	yhhS →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,615,866	0	A	G	19.5%	88.0 / 13.3	61	L137L (CTA→CTG)	yhhS	predicted transporter
Reads supporting (aligned to +/- strand): ref base A (9/40); new base G (8/4); total (17/44)
Fisher's exact test for biased strand distribution p-value = 2.09e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

ATTACTTTGCCTGGGGCGCGTCATCCTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTCACT  >  minE/2615800‑2615932

attaCTTTGCCTGGGGCGCGTCATCCTTGGGTTTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATg                                                                  >  1:1392790/1‑69 (MQ=255)

attaCTTTGCCTGGGGCGCGTCATCCTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATg                                                                  >  1:2329291/1‑69 (MQ=255)

attaCTTTGCCTGGGGCGCGTCATCCTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATg                                                                  >  1:1266191/1‑69 (MQ=255)

            ggggCGCGTCATCCTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGt                                                      <  1:420052/69‑1 (MQ=255)

            ggggCGCGTCATCCTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGt                                                      <  1:2418081/69‑1 (MQ=255)

            ggggCGCGTCATCCTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGt                                                      <  1:1818937/69‑1 (MQ=255)

               gcgcGTCATCCTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTg                                                     <  1:2556552/67‑1 (MQ=255)

               gcgcGTCATCCTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTg                                                     <  1:1496594/67‑1 (MQ=255)

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

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

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

                    tCATCCTTGGGATTGGGCAAAGTTTAGCCGGAACGGGATCGACCCTATGGGGCGTTg                                                          >  1:949918/1‑57 (MQ=255)

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

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

                      aTCCTTGGGATTGGACAAAGTTTTGCCGGAACGGGATCGACCCTGTGGGGTGTTGGCGTGGTTGGTTCg                                            <  1:1065998/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:3364955/69‑1 (MQ=255)

                         cTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTg                                         <  1:3094471/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:73154/69‑1 (MQ=255)

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

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

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

                                        aGTTTTGCCGGAACGGGATCGACCCTGTGGGGTGTTGGCGTGGTTggtt                                              >  1:712405/1‑47 (MQ=255)

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

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

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

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

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

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

                                        aGTTTTGCCGGAACGGGATCGACCCTGTGGGGTGTTCGCGt                                                      >  1:1015621/1‑41 (MQ=37)

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

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

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

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

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

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

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

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

                                                       ggATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGc           <  1:2228810/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:3234041/69‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

ATTACTTTGCCTGGGGCGCGTCATCCTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTCACT  >  minE/2615800‑2615932

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