BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,095,825	A→G	37.2%	L372L (TTG→CTG)	garK ←	glycerate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,095,825	0	A	G	37.2%	52.4 / 60.1	94	L372L (TTG→CTG)	garK	glycerate kinase I
Reads supporting (aligned to +/- strand): ref base A (20/39); new base G (22/13); total (42/52)
Fisher's exact test for biased strand distribution p-value = 9.70e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 108 total aligned reads displayed.

ACGTGAAGCACGGCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCAACATCATCGGTCAGGCTACCCGC  >  minE/2095757‑2095891

aCGTGAAGCACGGCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCaa                                                                    >  1:1775864/1‑69 (MQ=255)

aCGTGAAGCACGGCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCaa                                                                    >  1:3006178/1‑69 (MQ=255)

aCGTGAAGCACGGCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCaa                                                                    >  1:629285/1‑69 (MQ=255)

 cGTGAAGCACGGCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCaa                                                                    <  1:1120450/68‑1 (MQ=255)

 cGTGAAGCACGGCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCaa                                                                    <  1:26477/68‑1 (MQ=255)

  gTGAAGCACGGCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCaa                                                                    <  1:936626/67‑1 (MQ=255)

  gTGAAGCACGGCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCaa                                                                    <  1:3087738/67‑1 (MQ=255)

   tGAAGCACGGCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCaa                                                                    <  1:3186851/66‑1 (MQ=255)

   tGAAGCACGGCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCaa                                                                    <  1:1425372/66‑1 (MQ=255)

   tGAAGCACGGCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCaa                                                                    <  1:864824/66‑1 (MQ=255)

       gCACGGCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTg                                                             <  1:2354849/69‑1 (MQ=255)

       gCACGGCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTg                                                             <  1:1762779/69‑1 (MQ=255)

       gCACGGCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTg                                                             <  1:3036649/69‑1 (MQ=255)

             cAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGAc                                                        >  1:3088602/1‑68 (MQ=255)

             cAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGAc                                                        >  1:917541/1‑68 (MQ=255)

                 tatTGTCATAAGCACCG‑GGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGACCGCATc                                                  <  1:3086184/69‑1 (MQ=255)

                  atTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGACCGCATc                                                  <  1:2701400/69‑1 (MQ=255)

                  atTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGACCGCATc                                                  <  1:1047626/69‑1 (MQ=255)

                  atTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGACCGCATc                                                  <  1:3407294/69‑1 (MQ=255)

                  atTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGACCGCATc                                                  <  1:352172/69‑1 (MQ=255)

                  atTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGACCGCATc                                                  <  1:1082203/69‑1 (MQ=255)

                  atTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGACCGCATc                                                  <  1:3058201/69‑1 (MQ=255)

                  atTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGACCGCATc                                                  <  1:160056/69‑1 (MQ=255)

                  atTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGACCGCATc                                                  <  1:510203/69‑1 (MQ=255)

                  atTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGACCGCATc                                                  >  1:1478354/1‑69 (MQ=255)

                  atTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGACCGCATc                                                  <  1:862674/69‑1 (MQ=255)

                  atTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGACCGCATc                                                  <  1:2896247/69‑1 (MQ=255)

                  atTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGACCGCATc                                                  <  1:1351055/69‑1 (MQ=255)

                  atTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGACCGCATc                                                  <  1:12948/69‑1 (MQ=255)

                  atTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGACCGCATc                                                  <  1:2679765/69‑1 (MQ=255)

                   ttGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCaa                                                                    <  1:2219361/50‑1 (MQ=255)

                   ttGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCaa                                                                    <  1:1745044/50‑1 (MQ=255)

                   ttGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCaa                                                                    <  1:1619022/50‑1 (MQ=255)

                   ttGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGACCGCATc                                                  <  1:3369437/68‑1 (MQ=255)

                   ttGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGACCGCATc                                                  <  1:3375969/68‑1 (MQ=255)

                   ttGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGACCGCATc                                                  <  1:1005621/68‑1 (MQ=255)

                   ttGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGACCGCATc                                                  <  1:3145625/68‑1 (MQ=255)

                   ttGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGACCGCATc                                                  <  1:1502434/68‑1 (MQ=255)

                   ttGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGACCGCATc                                                  <  1:391571/68‑1 (MQ=255)

                    tGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTg                                                             <  1:247989/56‑1 (MQ=255)

                     gTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGACCGCATCAAt                                               >  1:764753/1‑69 (MQ=255)

                             ccccGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGACCGCATCAATGCCATGCt                                       <  1:164195/69‑1 (MQ=255)

                             ccccGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGACCGCATCAATGCCATGCt                                       <  1:2172863/69‑1 (MQ=255)

                                    gAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAg                                                          <  1:3377812/43‑1 (MQ=255)

                                    gAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGACCGCATCAATGCCATGCTGATGTAc                                >  1:1168077/1‑69 (MQ=255)

                                         cTTCGTCTAACGTACCAATGCTGGTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAAc                             >  1:1962422/1‑67 (MQ=255)

                                         cTTCGTCTAACGTACCAATGCTGGTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACg                            >  1:528357/1‑68 (MQ=255)

                                         cTTCGTCTAACGTACCAATGCTGGTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACg                            >  1:2389056/1‑68 (MQ=255)

                                         cTTCGTCTAACGTACCAATGCTGGTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACg                            >  1:2954473/1‑68 (MQ=255)

                                         cTTCGTCTAACGTACCAATGCTGGTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACg                            >  1:594916/1‑68 (MQ=255)

                                         cTTCGTCTAACGTACCAATGCTGGTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACg                            >  1:704359/1‑68 (MQ=255)

                                         cTTCGTCTAACGTACCAATGCTGGTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACg                            >  1:720013/1‑68 (MQ=255)

                                         cTTCGTCTAACGTACCAATGCTGGTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACg                            >  1:721623/1‑68 (MQ=255)

                                         cTTCGTCTAACGTACCAATGCTGGTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACg                            >  1:2668252/1‑68 (MQ=255)

                                         cTTCGTCTAACGTACCAATGCTGGTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACg                            >  1:3317439/1‑68 (MQ=255)

                                         cTTCGTCTAACGTACCAATGCTGGTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACg                            >  1:1409991/1‑68 (MQ=255)

                                         cTTCGTCTAACGTACCAATGCTGGTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACg                            >  1:2595639/1‑68 (MQ=255)

                                         cTTCGTCTAACGTACCAATGCTGGTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACg                            >  1:2477465/1‑68 (MQ=255)

                                         cTTCGTCTAACGTACCAATGCTGGTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACg                            >  1:1754568/1‑68 (MQ=255)

                                         cTTCGTCTAACGTACCAATGCTGGTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACg                            >  1:1180985/1‑68 (MQ=255)

                                         cTTCGTCTAACGTACCAATGCTGGTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACg                            >  1:1497668/1‑68 (MQ=255)

                                         cTTCGTCTAACGTACCAATGCTGGTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACg                            >  1:1277488/1‑68 (MQ=255)

                                         cTTCGTCTAACGTACCAATGCTGGTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACg                            >  1:1265491/1‑68 (MQ=255)

                                         cTTCGTCTAACGTACCAATGCTGGTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACg                            >  1:1848413/1‑68 (MQ=255)

                                         cTTCGTCTAACGTACCAATGCTGGTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACg                            >  1:2071420/1‑68 (MQ=255)

                                         cTTCGTCTAACGTACCAATGCTGGTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACg                            >  1:1958453/1‑68 (MQ=255)

                                                    gTACCTATGCTGGTCAATACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCAACATCATCg                >  1:1576896/1‑69 (MQ=255)

                                                             cTGGTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGT‑CAACGCCAACATCATTGGTCAGGCTAc      <  1:1885838/69‑1 (MQ=255)

                                                             cTGGTCAATACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCAACATCATCGGTCAGGCTa       >  1:1578544/1‑69 (MQ=255)

                                                             atGGTCAATACGCTGAAGACCGCATCAATGCCATGCTg                                      >  1:2048673/2‑38 (MQ=255)

                                                               ggTCAGTACGCTGAAGACCGCATCAATGCCATTCTGATGTACAACGCCAACATCATTGGTCAGGCTAc      <  1:2471975/68‑1 (MQ=255)

                                                               ggTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCAACATCATTGGTCAGGCTAc      <  1:1972250/68‑1 (MQ=255)

                                                               ggTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCAACATCATTGGTCAGGCTAc      <  1:1911767/68‑1 (MQ=255)

                                                               ggTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCAACATCATTGGTCAGGCTAc      <  1:3172782/68‑1 (MQ=255)

                                                               ggTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCAACATCATTGGTCAGGCTAc      <  1:797888/68‑1 (MQ=255)

                                                               ggTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCAACATCATTGGTCAGGCTAc      <  1:1260430/68‑1 (MQ=255)

                                                               ggTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCAACATCATTGGTCAGGCTAc      <  1:2204099/68‑1 (MQ=255)

                                                               ggTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCAACATCATTGGTCAGGCTAc      <  1:2286790/68‑1 (MQ=255)

                                                               ggTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCAACATCATTGGTCAGGCTAc      <  1:2450826/68‑1 (MQ=255)

                                                               ggTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCAACATCATTGGTCAGGCTAc      <  1:45967/68‑1 (MQ=255)

                                                               ggTCAGTACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCAACATCATTGGTCAGGCTAc      <  1:399881/68‑1 (MQ=255)

                                                                gTCAATACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCaa                        <  1:711748/49‑1 (MQ=255)

                                                                gTCAATACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCaa                        <  1:1720617/49‑1 (MQ=255)

                                                                gTCAATACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCaa                        <  1:3419301/49‑1 (MQ=255)

                                                                  cAATACGCTGAAGCCCGCATCAATGCCCTGCTGATGTACAACGCCAACATCATCGGTCAGGCTACCCGc  >  1:2614234/1‑69 (MQ=255)

                                                                  cAATACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCaa                        <  1:1683351/47‑1 (MQ=255)

                                                                  cAATACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCaa                        <  1:2180764/47‑1 (MQ=255)

                                                                  cAATACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCaa                        >  1:64344/1‑47 (MQ=255)

                                                                  cAATACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCaa                        <  1:1399190/47‑1 (MQ=255)

                                                                  cAATACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCAACATCATCGGTCAg           <  1:1284970/60‑1 (MQ=255)

                                                                  cAATACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCAACATCATCGGTCAg           >  1:1231288/1‑60 (MQ=255)

                                                                  cAATACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCAACATCATCGGTCAg           >  1:1029209/1‑60 (MQ=255)

                                                                  cAATACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCAACATCATCGGTCAg           >  1:1786326/1‑60 (MQ=255)

                                                                  cAATACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCAACATCATCGGTCAg           >  1:1790232/1‑60 (MQ=255)

                                                                  cAATACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCAACATCATCGGTCAg           >  1:1795433/1‑60 (MQ=255)

                                                                  cAATACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCAACATCATCGGTCAg           >  1:3215777/1‑60 (MQ=255)

                                                                  cAATACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCAACATCATCGGTCAg           <  1:2670330/60‑1 (MQ=255)

                                                                  cAATACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCAACATCATCGGTCAGGCTACCCGc  >  1:2178755/1‑69 (MQ=255)

                                                                  cAATACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCAACATCATCGGTCAGGCTACCCGc  >  1:2618352/1‑69 (MQ=255)

                                                                  cAATACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCAACATCATCGGTAAg           >  1:2358443/1‑60 (MQ=255)

ACGTGAAGCACGGCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAATACGCTGAAGACCGCATCAATGCCATGCTGATGTACAACGCCAACATCATCGGTCAGGCTACCCGC  >  minE/2095757‑2095891

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