BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,388,189	A→G	46.4%	G75G (GGT→GGC)	yigL ←	predicted hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,388,189	0	A	G	46.4%	1.8 / 45.9	56	G75G (GGT→GGC)	yigL	predicted hydrolase
Reads supporting (aligned to +/- strand): ref base A (11/19); new base G (22/4); total (33/23)
Fisher's exact test for biased strand distribution p-value = 3.60e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

TTGTCGTTGACTACGCCAAACAGATCGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACTTAATCTCCAGATTATCGCG  >  minE/2388121‑2388257

ttGTCGTTGACTACGCCAAACAGATCGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTa                                                                      <  1:1713551/69‑1 (MQ=255)

ttGTCGTTGACTACGCCAAACAGATCGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTa                                                                      <  1:3412243/69‑1 (MQ=255)

ttGTCGTTGACTACGCCAAACAGATCGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTa                                                                      <  1:14611/69‑1 (MQ=255)

ttGTCGTTGACTACGCCAAACAGATCGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTa                                                                      <  1:2332429/69‑1 (MQ=255)

ttGTCGTTGACTACGCCAAACAGATCGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTa                                                                      <  1:2247117/69‑1 (MQ=255)

ttGTCGTTGACTACGCCAAACAGATCGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTa                                                                      <  1:666444/69‑1 (MQ=255)

ttGTCGTTGACTACGCCAAACAGATCGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTa                                                                      <  1:3066533/69‑1 (MQ=255)

                         cGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCAc                                             <  1:1735913/69‑1 (MQ=255)

                         cGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCAc                                             <  1:1571577/69‑1 (MQ=255)

                         cGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCAc                                             <  1:1455187/69‑1 (MQ=255)

                         cGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCAc                                             <  1:2372789/69‑1 (MQ=255)

                         cGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCAc                                             <  1:971682/69‑1 (MQ=255)

                          gCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCAcc                                            <  1:1436329/69‑1 (MQ=255)

                          gCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCAcc                                            <  1:3008474/69‑1 (MQ=255)

                          gCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCAc                                             <  1:1896685/68‑1 (MQ=255)

                                     cgcgATCCAGGTTATGATCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTa                                  >  1:2578235/1‑68 (MQ=255)

                                     cgcgATCCAGGTTATGAGCAAATATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTa                                  >  1:1190348/1‑68 (MQ=255)

                                     cgcgATCCAGGTTATGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTa                                  >  1:1880586/1‑68 (MQ=255)

                                     cgcgATCCAGGTTATGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTa                                  >  1:1134455/1‑68 (MQ=255)

                                     cgcgATCAAGGTTATGAGCAAAAATCAGGTTGCCATCAAGATCGTGAACGCGCGCGCCATTGGAGGTa                                  >  1:276227/1‑68 (MQ=255)

                                      gcgATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAAt                                <  1:1874949/69‑1 (MQ=255)

                                      gcgATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAAt                                >  1:2995231/1‑69 (MQ=255)

                                      gcgATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAAt                                >  1:41172/1‑69 (MQ=255)

                                      gcgATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAAt                                <  1:1945759/69‑1 (MQ=255)

                                      gcgATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAAt                                >  1:1017582/1‑69 (MQ=255)

                                           ccAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCAt                             >  1:3142476/1‑67 (MQ=255)

                                              ggTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTaa                         <  1:2906638/68‑1 (MQ=255)

                                              ggTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTaa                         <  1:2985614/68‑1 (MQ=255)

                                              ggTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTaa                         <  1:2640553/68‑1 (MQ=255)

                                              ggTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTaa                         <  1:865891/68‑1 (MQ=255)

                                              ggTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTaa                         <  1:2484541/68‑1 (MQ=255)

                                              ggTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTaa                         <  1:3064321/68‑1 (MQ=255)

                                              ggTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTaa                         <  1:1355763/68‑1 (MQ=255)

                                              ggTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTaa                         <  1:1722473/68‑1 (MQ=255)

                                                ttATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGa                       >  1:1484006/1‑68 (MQ=255)

                                                ttATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGa                       >  1:2604691/1‑68 (MQ=255)

                                                ttATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGa                       >  1:1732665/1‑68 (MQ=255)

                                                 tATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCa                              >  1:234902/1‑60 (MQ=255)

                                                 tATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGAc                      >  1:1152082/1‑68 (MQ=255)

                                                 tATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGAc                      >  1:1853521/1‑68 (MQ=255)

                                                  aTGAGCAAGAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                              >  1:1341656/1‑59 (MQ=255)

                                                  aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAAt                                >  1:1419630/1‑57 (MQ=255)

                                                  aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATc                               >  1:1795305/1‑58 (MQ=255)

                                                  aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                              >  1:3143561/1‑59 (MQ=255)

                                                  aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                              >  1:2253128/1‑59 (MQ=255)

                                                  aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                              >  1:429727/1‑59 (MQ=255)

                                                  aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                              >  1:2968155/1‑59 (MQ=255)

                                                  aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                              >  1:521874/1‑59 (MQ=255)

                                                  aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                              >  1:1115152/1‑59 (MQ=255)

                                                  aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                              >  1:255217/1‑59 (MQ=255)

                                                  aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                              >  1:137106/1‑59 (MQ=255)

                                                  aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                              >  1:2418287/1‑59 (MQ=255)

                                                  aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                              >  1:1738687/1‑59 (MQ=255)

                                                  aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                              >  1:1761502/1‑59 (MQ=255)

                                                  aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                              >  1:2750752/1‑59 (MQ=255)

                                                  aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                              >  1:1901317/1‑59 (MQ=255)

                                                  aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                              >  1:1948842/1‑59 (MQ=255)

                                                  aTGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACtt                    <  1:1985291/69‑1 (MQ=255)

                                                     aGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCATGTAAGACTTAAt                 <  1:3297449/69‑1 (MQ=255)

                                                     aGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCATGTAAGACTTAAt                 <  1:1119772/69‑1 (MQ=255)

                                                     aGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCATGTAAGACTTAAt                 <  1:3309369/69‑1 (MQ=255)

                                                     aGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCATGTAAGACTTAAt                 <  1:3063921/69‑1 (MQ=255)

                                                              cAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACTTAATCTCCAGAtt        >  1:2107349/1‑69 (MQ=255)

                                                                    aCCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACTTAATCTCCAGATTATcgcg  >  1:2056260/1‑69 (MQ=255)

                                                                    aCCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACTTAATCTCCAGATTATcgcg  >  1:2492610/1‑69 (MQ=255)

                                                                    aCCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACTTAATCTCCAGATTATcgcg  >  1:2745946/1‑69 (MQ=255)

TTGTCGTTGACTACGCCAAACAGATCGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACTTAATCTCCAGATTATCGCG  >  minE/2388121‑2388257

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