BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,513,416	T→C	50.0%	F203S (TTT→TCT)	gltS →	glutamate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,513,416	0	T	C	50.0%	‑2.4 / 75.2	50	F203S (TTT→TCT)	gltS	glutamate transporter
Reads supporting (aligned to +/- strand): ref base T (15/10); new base C (18/7); total (33/17)
Fisher's exact test for biased strand distribution p-value = 5.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCGCTATCTGGTGAAACACTCCACCACGCCGAACGGTATTCCGGATGACCAGGAAGTCCCGACGGCGTTTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCTGATTGAAACTATCGCGCTGATTGCT  >  minE/2513346‑2513483

gcgcgcTATCTGGTGAAACACTCCACCACGCCGAACGGTATTCCGGATGACCAGGAAGTCCCGACGGCGtt                                                                     <  1:2954781/71‑1 (MQ=255)

        tCTGGTGAAACACTCCACCACGCCGAACGGTATTCCGGATGACCAGGAAGTCCCGACGGCGTTTGaaaa                                                               <  1:2941673/69‑1 (MQ=255)

        tCTGGTGAAACACTCCACCACGCCGAACGGTATTCCGGATGACCAGGAAGTCCCGACGGCGTTTGAAAAg                                                              <  1:2248915/70‑1 (MQ=255)

          tGGTGAAACACTCCACCACGCCGAACGGTATTCCGGATGACCAGGAAGTCCCGACGGCGTTTGAAAAGCCg                                                           >  1:3585084/1‑71 (MQ=255)

          tGGTGAAACACTCCACCACGCCGAACGGTATTCCGGATGACCAGGAAGTCCCGACGGCGTCTGAAAAGCCg                                                           >  1:1401167/1‑71 (MQ=255)

          tGGTGAAACACTCCACCACGCCGAACGGTATTCCGGATGACCAGGAAGTCCCGACGGCGTCTGAAAAGCCg                                                           >  1:2503630/1‑71 (MQ=255)

          tGGTGAAACACTCCACCACGCCGAACGGTATTCCGGATGACCAGGAAGTCCCGACGGCGTCTGAAAAGCCg                                                           >  1:3091952/1‑71 (MQ=255)

            gTGAAACACTCCACCACGCCGAACGGTATTCCGGATGACCAGGAAGTCCCGACGGCGTTTGAAAAGCCg                                                           <  1:2636900/69‑1 (MQ=255)

            gTGAAACACTCCACCACGCCGAACGGTATTCCGGATGACCAGGAAGTCCCGACGGCGTTTGAAAAGCCg                                                           <  1:2019014/69‑1 (MQ=255)

                           acGCCGAACGGTATTCCGGATGACCAGGAAGTCCCGACGGCGTTTGAAAAGCCGGATGTGGGACGCATGAt                                          <  1:2208741/71‑1 (MQ=255)

                            cGCCGAACGGTATTCCGGATGACCAGGAAGTCCCGACGGCGTTTGAAAAGCCGGATGTGGGACGCATGATc                                         <  1:841966/71‑1 (MQ=255)

                            cGCCGAACGGTATTCCGGATGACCAGGAAGTCCCGACGGCGTCTGAAAAGCCGGATGTGGGACGCATGATc                                         <  1:2473447/71‑1 (MQ=255)

                               cGAACGGTATTCCGGATGACCAGGAAGTCCCGACGGCGTTTGAAAAGCCGGATGTGGGACGCATGATc                                         <  1:3721744/68‑1 (MQ=255)

                                 aaCGGTATTCCGGATGACCAGGAAGTCCCGACGGCGTTTGAAAAGCCGGATGTGGGACGCATGATCACCTc                                    <  1:3749560/71‑1 (MQ=255)

                                   cGGTATTCCGGATGACCAGGAAGTCCCGACGGCGTCTGAAAAGCCGGATGTGGGACGCATGATCACCTCg                                   >  1:622554/1‑70 (MQ=255)

                                          ccGGATGACCAGGAAGTCCCGACGGCGTCTGAAAAGCCGGATGTGGGACGCATGATCAc                                       >  1:3378883/1‑59 (MQ=255)

                                          ccGGATGACCAGGAAGTCCCGACGGCGTCTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCt                           >  1:2320918/1‑71 (MQ=255)

                                           cGGATGACCAGGAAGTCCCGACGGCGTTTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCt                           <  1:3616250/70‑1 (MQ=255)

                                           cGGATGACCAGGAAGTCCCGACGGCGTCTGAAAAGCCGGATGTGGGTCGCATGATCACCTCGTTGGTGCt                           <  1:129363/70‑1 (MQ=255)

                                           cGGATGACCAGGAAGTCCCGACGGCGTCTGAAAAGCCGGATGTGGGACGCATg                                            >  1:868454/1‑53 (MQ=255)

                                           cGGATGACCAGGAAGTCCCGACGGCGTCTGAAAAGCCGGATGTGGGACGCATGa                                           >  1:1310330/1‑54 (MQ=255)

                                           cGGATGACCAGGAAGTCCCGACGGCGTCTGAAAAGCCGGATGTGGGACGCATGa                                           >  1:312025/1‑54 (MQ=255)

                                           cGGATGACCAGGAAGTCCCGACGGCGTCTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCt                           <  1:649528/70‑1 (MQ=255)

                                               tGACCAGGAAGTCCCGACGGCGTCTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCTGAtt                       >  1:927954/1‑70 (MQ=255)

                                               tGACCAGGAAGTCCCGACGGCGTCTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCTGAtt                       >  1:2377918/1‑70 (MQ=255)

                                                   cAGGAAGTCCCGACGGCGTTTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCTGATTGAAAc                  >  1:1694270/1‑71 (MQ=255)

                                                   cAGGAAGTCCCGACGGCGTTTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCTGATTGAAAc                  >  1:1762047/1‑71 (MQ=255)

                                                   cAGGAAGTCCCGACGGCGTTTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCTGATTGAAAc                  >  1:4064101/1‑71 (MQ=255)

                                                   cAGGAAGTCCCGACGGCGTTTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCTGATTGAAAc                  >  1:3372296/1‑71 (MQ=255)

                                                   cAGGAAGTCCCGAAGGCGTTTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCTGATTGAAAc                  >  1:2713238/1‑71 (MQ=255)

                                                       aaGTCCCGACGGCGTTTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCTGATTGAAACTAt               >  1:2667155/1‑70 (MQ=255)

                                                       aaGTCCCGACGGCGTTTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCTGATTGAAACTATc              >  1:3169539/1‑71 (MQ=255)

                                                       aaGTCCCGACGGCGTTTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCTGATTGAAACTATc              >  1:1943216/1‑71 (MQ=255)

                                                        aGTCCCGACGGCGTTTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCTGATTGAAACTATcg             >  1:1714031/1‑71 (MQ=255)

                                                        aGTCCCGACGGCGTTTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCTGATTGAAACTATcg             >  1:4088491/1‑71 (MQ=255)

                                                        aGTCCCGACGGCGTTTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCTGATTGAAACTATcg             >  1:3094358/1‑71 (MQ=255)

                                                        aGTCCCGACGGCGTCTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCTGATTGAAAc                  >  1:2876496/1‑66 (MQ=255)

                                                        aGTCCCGACGGCGTCTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCTGATTGAAAc                  >  1:3699962/1‑66 (MQ=255)

                                                        aGTCCCGACGGCGTCTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCTGATTGAAAc                  >  1:3196154/1‑66 (MQ=255)

                                                        aGTCCCGACGGCGTCTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTAGTGCTGATTGAAAc                  >  1:4077456/1‑66 (MQ=255)

                                                         gTCCCGACGGCGTTTGAAAAGCCGGATGTTGGACGCATGATCACCTCGTTg                                >  1:2200895/1‑51 (MQ=255)

                                                         gTCCCGACGGCGTTTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTg                                >  1:3479308/1‑51 (MQ=255)

                                                         gTCCCGACGGCGTTTGAAAAGCCGGATGTGGGACGCATGAT‑ACCTCGTTg                                >  1:69410/1‑50 (MQ=255)

                                                           cccGACGGCGTCTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCt                           <  1:916379/54‑1 (MQ=255)

                                                           cccGACGGCGTCTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCTGATTGAAACTATCgcgc          >  1:4011770/1‑71 (MQ=255)

                                                           cccGACGGCGTCTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCTGATTGAAACTATCgcgc          >  1:2077042/1‑71 (MQ=255)

                                                           cccGACGGCGTCTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCTGATTGAAACTATCgcgc          >  1:2328548/1‑71 (MQ=255)

                                                            ccGACGGCGTCTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCTGATTGAAACTa                <  1:2455504/64‑1 (MQ=255)

                                                            ccGACGGCGTCTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCTGATTGAAACTa                <  1:3095146/64‑1 (MQ=255)

                                                               aCGGCGTCTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCTGATTGAAACTATCGCGCTGAt      <  1:1994106/71‑1 (MQ=255)

                                                                   cGTTTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCTGATTGAAACTATCGCGCTGATTGCt  <  1:3345981/71‑1 (MQ=255)

GCGCGCTATCTGGTGAAACACTCCACCACGCCGAACGGTATTCCGGATGACCAGGAAGTCCCGACGGCGTTTGAAAAGCCGGATGTGGGACGCATGATCACCTCGTTGGTGCTGATTGAAACTATCGCGCTGATTGCT  >  minE/2513346‑2513483

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