BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	437,386	G→A	48.0%	R176R (AGG→AGA)	pgm →	phosphoglucomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	437,386	0	G	A	48.0%	4.0 / 66.0	77	R176R (AGG→AGA)	pgm	phosphoglucomutase
Reads supporting (aligned to +/- strand): ref base G (9/31); new base A (20/17); total (29/48)
Fisher's exact test for biased strand distribution p-value = 5.25e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAATACAATCCGCCAAATGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGTATCTCCCTCGACGAAGCGATGGCA  >  minE/437323‑437452

cAAATACAATCCGCCAAATGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGGGCCaa                                                               >  1:864496/1‑69 (MQ=255)

             ccAAATGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGct                                                  <  1:1038849/69‑2 (MQ=255)

             ccAAATGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGct                                                  <  1:2278366/69‑2 (MQ=255)

             ccAAATGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGct                                                  <  1:2414921/69‑2 (MQ=255)

             ccAAATGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGct                                                  <  1:1913071/69‑2 (MQ=255)

             ccAAATGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGct                                                  <  1:2513125/69‑2 (MQ=255)

             ccAAATGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGct                                                  <  1:2839400/69‑2 (MQ=255)

             ccAAATGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGct                                                  <  1:305345/69‑2 (MQ=255)

             ccAAATGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGct                                                  <  1:1603096/69‑2 (MQ=255)

             ccAAATGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGct                                                  <  1:585713/69‑2 (MQ=255)

             ccAAATGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGct                                                  <  1:861805/69‑2 (MQ=255)

             ccAAATGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGct                                                  <  1:1552589/69‑2 (MQ=255)

             ccAAATGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGct                                                  <  1:92887/69‑2 (MQ=255)

             ccAAATGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGca                                                  <  1:2742295/69‑2 (MQ=255)

             ccAAATGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAAAGCACTGCTGGct                                                  <  1:2746957/69‑2 (MQ=255)

               aaaTGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCATTGCTGGCTGa                                                >  1:2211822/1‑69 (MQ=255)

               aaaTGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGCTGa                                                >  1:2715999/1‑69 (MQ=255)

               aaaTGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGCTGa                                                >  1:2255243/1‑69 (MQ=255)

               aaaTGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGCTGa                                                >  1:1069562/1‑69 (MQ=255)

               aaaTGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGCTGa                                                >  1:2296462/1‑69 (MQ=255)

               aaaTGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGCTGa                                                >  1:932394/1‑69 (MQ=255)

               aaaTGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGCTGa                                                >  1:916935/1‑69 (MQ=255)

               aaaTGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGCTGa                                                >  1:254685/1‑69 (MQ=255)

               aaaTGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGCTGa                                                >  1:870623/1‑69 (MQ=255)

               aaaTGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGCTGa                                                >  1:2847680/1‑69 (MQ=255)

               aaaTGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGCTGa                                                >  1:841200/1‑69 (MQ=255)

               aaaTGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGCTGa                                                >  1:311329/1‑69 (MQ=255)

               aaaTGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGCTGa                                                >  1:59234/1‑69 (MQ=255)

               aaaTGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGCTGa                                                >  1:1499805/1‑69 (MQ=255)

               aaaTGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGCTGa                                                >  1:2195457/1‑69 (MQ=255)

               aaaTGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGCTGa                                                >  1:1525661/1‑69 (MQ=255)

               aaaTGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGCTGa                                                >  1:2181765/1‑69 (MQ=255)

               aaaTGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGCTGa                                                >  1:1743034/1‑69 (MQ=255)

               aaaTGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGCTGa                                                >  1:1090930/1‑69 (MQ=255)

               aaaTGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGCTGa                                                >  1:1706783/1‑69 (MQ=255)

                 aTGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGGGCCAACGCACTGCTGGCCGATg                                              <  1:236595/69‑1 (MQ=255)

                 aTGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGGGCCAACGCACTGCTGGCCGATg                                              <  1:1886511/69‑1 (MQ=255)

                 aTGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGGGCCAACGCACTGCTGGCCGATg                                              <  1:3078218/69‑1 (MQ=255)

                 aTGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGGGCCAACGCACTGCTGGCCGATg                                              <  1:2162990/69‑1 (MQ=255)

                  tggtggCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGGGCCAACGCACTGCTGGCCGATg                                              <  1:1546294/68‑1 (MQ=255)

                  tggtggCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGGGCCAACGCACTGCTGGCCGATg                                              <  1:2201773/68‑1 (MQ=255)

                    gtggCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGct                                                  <  1:2719592/62‑2 (MQ=255)

                       gCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGct                                                  <  1:823523/59‑2 (MQ=255)

                                  aCCAACGTCACTAAAGTGGTGGAAGACAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAg                             <  1:1449313/69‑1 (MQ=255)

                                  aCCAACGTCACTAAAGTGGTGGAAGACAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAg                             <  1:823041/69‑1 (MQ=255)

                                  aCCAACGTCACTAAAGTGGTGGAAGACAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAg                             <  1:2670131/69‑1 (MQ=255)

                                     aaCGTCACTAAAGTGGTGGAAGACAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGt                          >  1:1262256/1‑69 (MQ=255)

                                     aaCGTCACTAAAGTGGTGGAAGACAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGt                          >  1:1258203/1‑69 (MQ=255)

                                     aaCGTCACTAAAGTGGTGGAAGACAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGt                          >  1:1882268/1‑69 (MQ=255)

                                     aaCGTCACTAAAGTGGTGGAAGACAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGt                          >  1:1646172/1‑69 (MQ=255)

                                     aaCGTCACTAAAGTGGTGGAAGACAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGt                          >  1:2933990/1‑69 (MQ=255)

                                     aaCGTCACTAAAGTGGTGGAAGACAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGt                          >  1:1624745/1‑69 (MQ=255)

                                     aaCGTCACTAAAGTGGTGGAAGACAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGt                          >  1:3384669/1‑69 (MQ=255)

                                        gTCACTAAAGTGGTGGAAGACAGAGCCAACGCACTGCTGGct                                                  <  1:2297104/42‑2 (MQ=255)

                                              aaaGTGGTGGAAGACAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGTATCTCCCTc                 <  1:2230812/69‑1 (MQ=255)

                                                      ggAAGACAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGTATCTCCCTCGACGAAGc         <  1:2209400/69‑1 (MQ=255)

                                                      ggAAGACAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGTATCTCCCTCGACGAAGc         <  1:1178970/69‑1 (MQ=255)

                                                      ggAAGACAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGTATCTCCCTCGACGAAGc         <  1:429485/69‑1 (MQ=255)

                                                      ggAAGACAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGTATCTCCCTCGACGAAGc         <  1:2112591/69‑1 (MQ=255)

                                                      ggAAGACAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGTATCTCCCTCGACGAAGc         <  1:3148445/69‑1 (MQ=255)

                                                      ggAAGACAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGTATCTCCCTCGACGAAGc         <  1:2101958/69‑1 (MQ=255)

                                                      ggAAGACAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGTATCTCCCTCGACGAAGc         <  1:1743515/69‑1 (MQ=255)

                                                      ggAAGACAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGTATCTCCCTCGACGAAGc         <  1:2609382/69‑1 (MQ=255)

                                                       gAAGACAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGTATCTCCCTCGACGAAGc         <  1:2816461/68‑1 (MQ=255)

                                                         aGACAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGTATCTCCCTCGACGAAGCGAt      <  1:321065/69‑1 (MQ=255)

                                                         aGACAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGTATCTCCCTCGACGAAGCGAt      <  1:1560430/69‑1 (MQ=255)

                                                         aGACAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGTATCTCCCTCGACGAAGCGAt      <  1:66455/69‑1 (MQ=255)

                                                           aCAGGGCCAACGCACTGCTTGCCGATGGCCTGAAAGGCGTGAAGCGTATCTCCCTCGACGAAGCGATgg    <  1:1584895/69‑1 (MQ=255)

                                                           aCAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGTATCTCCCTCGACGAAGc         <  1:825647/64‑1 (MQ=255)

                                                           aCAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGTATCTCCCTCGACGAAGCGATgg    <  1:94482/69‑1 (MQ=255)

                                                           aCAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGTATCTCCCTCGACGAAGCGATgg    <  1:2559355/69‑1 (MQ=255)

                                                           aCAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGTATCTCCCTCGACGAAGCGATgg    <  1:1871602/69‑1 (MQ=255)

                                                           aCAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGTATCTCCCTCGACGAAGCGATgg    <  1:2737270/69‑1 (MQ=255)

                                                           aCAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGTATCTCCCTCGACGAAGCGATgg    <  1:3348371/69‑1 (MQ=255)

                                                           aCAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGTATCTCCCTCGACGAAGCGATgg    <  1:3327740/69‑1 (MQ=255)

                                                           aCAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGCATCTCCCTCGACGAAGCGATgg    <  1:3123059/69‑1 (MQ=255)

                                                             aGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGTATCTCCCTCGACGAAGCGATGGCa  >  1:2703327/1‑69 (MQ=255)

CAAATACAATCCGCCAAATGGTGGCCCGGCTGATACCAACGTCACTAAAGTGGTGGAAGACAGGGCCAACGCACTGCTGGCCGATGGCCTGAAAGGCGTGAAGCGTATCTCCCTCGACGAAGCGATGGCA  >  minE/437323‑437452

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