BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	454,254	A→G	75.8%	S277S (AGT→AGC)	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	454,254	0	A	G	75.8%	98.8 / 24.8	78	S277S (AGT→AGC)	ybgH	predicted transporter
Reads supporting (aligned to +/- strand): ref base A (10/9); new base G (56/3); total (66/12)
Fisher's exact test for biased strand distribution p-value = 7.29e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GATATAAAGGCTAATCGAGCTACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTT  >  minE/454200‑454315

gATATAAAGGCTAATCGAGCTACCCTCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                                 >  1:3045005/1‑69 (MQ=255)

gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                                 >  1:1135104/1‑69 (MQ=255)

gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                                 >  1:913099/1‑69 (MQ=255)

gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                                 >  1:747532/1‑69 (MQ=255)

gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                                 >  1:707075/1‑69 (MQ=255)

gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                                 >  1:670355/1‑69 (MQ=255)

gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                                 >  1:645792/1‑69 (MQ=255)

gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                                 >  1:522320/1‑69 (MQ=255)

gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                                 >  1:488157/1‑69 (MQ=255)

gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                                 >  1:3027596/1‑69 (MQ=255)

gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                                 >  1:281568/1‑69 (MQ=255)

gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                                 >  1:2783203/1‑69 (MQ=255)

gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                                 >  1:2656722/1‑69 (MQ=255)

gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                                 >  1:2622523/1‑69 (MQ=255)

gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                                 >  1:1692809/1‑69 (MQ=255)

gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                                 >  1:129256/1‑69 (MQ=255)

gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                                 >  1:1534273/1‑69 (MQ=255)

gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                                 >  1:1648013/1‑69 (MQ=255)

gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                                 >  1:11688/1‑69 (MQ=255)

gATATAAAGGCTAATCGAGATACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                                 >  1:3081202/1‑69 (MQ=255)

      aaGGCTAATCGAGCTACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTcac                                           >  1:720240/1‑69 (MQ=255)

      aaGGCTAATCGAGCTACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTcac                                           >  1:265833/1‑69 (MQ=255)

      aaGGCTAATCGAGCTACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTcac                                           >  1:2105492/1‑69 (MQ=255)

               cGAGCTACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTg                                                   >  1:825206/1‑52 (MQ=255)

               cGAGCTACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTg                                                   >  1:2933455/1‑52 (MQ=255)

               cGAGCTACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTg                                                   >  1:2335684/1‑52 (MQ=255)

               cGAGCTACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTg                                                   >  1:32057/1‑52 (MQ=255)

               cGAGCTACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTg                                                   >  1:1965174/1‑52 (MQ=255)

               cGAGCTACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTCACAATCAGc                                    <  1:182124/67‑1 (MQ=255)

               cGAGCTACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTCACAATCAGc                                    <  1:1222627/67‑1 (MQ=255)

               cGAGCTACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTCACAATCAGc                                    <  1:1414855/67‑1 (MQ=255)

                  gCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCa                                <  1:286840/68‑1 (MQ=255)

                  gCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCa                                <  1:3367426/68‑1 (MQ=255)

                  gCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCa                                <  1:3187811/68‑1 (MQ=255)

                    tACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCt                             <  1:811201/69‑1 (MQ=255)

                    tACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCt                             <  1:1993122/69‑1 (MQ=255)

                    tACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCt                             <  1:1251521/69‑1 (MQ=255)

                    tACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCt                             <  1:1581935/69‑1 (MQ=255)

                              ttgtGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTg                   >  1:1747834/1‑69 (MQ=255)

                              ttgtGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTg                   >  1:58892/1‑69 (MQ=255)

                                          ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGtttt        >  1:2305692/1‑68 (MQ=255)

                                          ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGtttt        >  1:1790977/1‑68 (MQ=255)

                                          ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTc       >  1:395850/1‑69 (MQ=255)

                                          ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTc       >  1:1169824/1‑69 (MQ=255)

                                          ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTc       >  1:311570/1‑69 (MQ=255)

                                          ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTc       >  1:2927921/1‑69 (MQ=255)

                                          ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTc       >  1:772244/1‑69 (MQ=255)

                                          ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTc       >  1:1662344/1‑69 (MQ=255)

                                          ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTc       >  1:2382367/1‑69 (MQ=255)

                                          ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTc       >  1:2266002/1‑69 (MQ=255)

                                          ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTc       >  1:2110715/1‑69 (MQ=255)

                                          ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTc       >  1:911619/1‑69 (MQ=255)

                                          ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTc       >  1:11234/1‑69 (MQ=255)

                                          ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTc       >  1:969270/1‑69 (MQ=255)

                                          ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTc       >  1:1940833/1‑69 (MQ=255)

                                          ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTc       >  1:1872910/1‑69 (MQ=255)

                                          ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTAAGCTGTTTTTGGtttt        >  1:633953/1‑68 (MQ=255)

                                            aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGtt          >  1:321731/1‑64 (MQ=255)

                                            aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg     >  1:89835/1‑69 (MQ=255)

                                            aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg     >  1:166451/1‑69 (MQ=255)

                                            aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg     >  1:971492/1‑69 (MQ=255)

                                            aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg     >  1:2652927/1‑69 (MQ=255)

                                            aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg     >  1:1724399/1‑69 (MQ=255)

                                            aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg     >  1:1947926/1‑69 (MQ=255)

                                            aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg     >  1:2022828/1‑69 (MQ=255)

                                            aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg     >  1:2050359/1‑69 (MQ=255)

                                            aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg     >  1:2116057/1‑69 (MQ=255)

                                            aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg     >  1:2211205/1‑69 (MQ=255)

                                            aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg     >  1:2632299/1‑69 (MQ=255)

                                            aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg     >  1:354016/1‑69 (MQ=255)

                                            aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg     >  1:2989567/1‑69 (MQ=255)

                                            aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg     >  1:3020108/1‑69 (MQ=255)

                                            aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg     >  1:1511849/1‑69 (MQ=255)

                                            aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg     >  1:3140766/1‑69 (MQ=255)

                                            aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg     >  1:3383827/1‑69 (MQ=255)

                                            aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCACCt                             >  1:236556/1‑45 (MQ=255)

                                               acaacaTACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCtt  <  1:277675/69‑1 (MQ=255)

                                               acaacaTACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCtt  <  1:1179325/69‑1 (MQ=255)

GATATAAAGGCTAATCGAGCTACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTT  >  minE/454200‑454315

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