BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	511,757	G→A	43.1%	G600G (GGG→GGA)	ybhJ →	predicted hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	511,757	0	G	A	43.1%	10.0 / 28.2	44	G600G (GGG→GGA)	ybhJ	predicted hydratase
Reads supporting (aligned to +/- strand): ref base G (6/18); new base A (16/2); total (22/22)
Fisher's exact test for biased strand distribution p-value = 4.48e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

CGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCA  >  minE/511692‑511816

cgcgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCgggg                                                             <  1:1100151/66‑1 (MQ=255)

cgcgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCgggg                                                             <  1:962960/66‑1 (MQ=255)

cgcgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCgggg                                                             <  1:411278/66‑1 (MQ=255)

cgcgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCgggg                                                             <  1:3237569/66‑1 (MQ=255)

cgcgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCgggg                                                             <  1:2948825/66‑1 (MQ=255)

cgcgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCgggg                                                             <  1:2935992/66‑1 (MQ=255)

cgcgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCgggg                                                             <  1:2447715/66‑1 (MQ=255)

cgcgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCgggg                                                             <  1:2330311/66‑1 (MQ=255)

cgcgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCgggg                                                             <  1:225575/66‑1 (MQ=255)

cgcgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCgggg                                                             <  1:1759093/66‑1 (MQ=255)

cgcgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCgggg                                                             <  1:1674695/66‑1 (MQ=255)

cgcgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCgggg                                                             <  1:1095621/66‑1 (MQ=255)

cgcgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCgggg                                                             <  1:1411773/66‑1 (MQ=255)

 gcgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCgggg                                                             <  1:1886061/65‑1 (MQ=255)

  cgATCCAGGTTATGTTGGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGAAATGt                                                        <  1:141210/69‑1 (MQ=255)

   gATCCAGGTTATGTTGGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGAAATGt                                                        <  1:2568481/68‑1 (MQ=255)

                  aGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCgggg                                                             <  1:792294/48‑1 (MQ=255)

                  aGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCgggg                                                             <  1:631278/48‑1 (MQ=255)

                  aGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCgggg                                                             <  1:3219616/48‑1 (MQ=255)

                  aGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCgggg                                                             <  1:1764168/48‑1 (MQ=255)

                  aGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCgggg                                                             <  1:243281/48‑1 (MQ=255)

                            aaGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGTAATGTCAGCGAGCTGACAGAGgt                                      <  1:1447633/61‑1 (MQ=255)

                            aaGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGgt                                      <  1:1238025/61‑1 (MQ=255)

                            aaGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGgt                                      <  1:1588610/61‑1 (MQ=255)

                            aaGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGgt                                      <  1:538138/61‑1 (MQ=255)

                            aaGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGgt                                      <  1:1898343/61‑1 (MQ=255)

                            aaGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGgt                                      <  1:1302467/61‑1 (MQ=255)

                            aaGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGgt                                      <  1:2528560/61‑1 (MQ=255)

                            aaGCGACTGCTGAGCTGAAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGgt                                      <  1:3220981/61‑1 (MQ=255)

                               cGACTGCTGAGCTGGAAAATCAGCGTCTGGCgggg                                                             <  1:1766053/35‑1 (MQ=255)

                                gACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGAAATGTCAGCGAGCTGACAGAGGTGTTTgcgcgc                            >  1:2889171/1‑67 (MQ=255)

                                gACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGAAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCa                           >  1:3328395/1‑68 (MQ=255)

                                gACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGAAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCa                           >  1:219150/1‑68 (MQ=255)

                                gACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGAAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCa                           >  1:455890/1‑68 (MQ=255)

                                gACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGAAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCa                           >  1:450397/1‑68 (MQ=255)

                                gACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGAAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCa                           >  1:394770/1‑68 (MQ=255)

                                gACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGAAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCa                           >  1:3374280/1‑68 (MQ=255)

                                gACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGAAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCa                           >  1:250663/1‑68 (MQ=255)

                                gACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGAAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCa                           >  1:3309143/1‑68 (MQ=255)

                                gACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGAAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCa                           >  1:3272171/1‑68 (MQ=255)

                                gACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGAAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCa                           >  1:3229031/1‑68 (MQ=255)

                                gACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGAAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCa                           >  1:1246488/1‑68 (MQ=255)

                                gACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGAAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCa                           >  1:1377644/1‑68 (MQ=255)

                                gACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGAAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCa                           >  1:265731/1‑68 (MQ=255)

                                gACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGAAAAGTCAGCGAGCTGACAGAGGAGTTTGCGCGCa                           >  1:798582/1‑68 (MQ=255)

                                gACTGCTGAGCTGGAAAATCAGCGTCTGG‑GGGAAATGTCAGCGAGCTGACAGAGGTGTTTgcgcgc                            >  1:115555/1‑66 (MQ=255)

                                   tgctgAGCTGGAAAATCAGCGTCTGGCGGGTAATGTCAGCGAGCTGACAGAGgt                                      <  1:449398/54‑1 (MQ=255)

                                                  tCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTc         >  1:953530/1‑68 (MQ=255)

                                                  tCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTc         >  1:2650039/1‑68 (MQ=255)

                                                  tCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTc         >  1:859072/1‑68 (MQ=255)

                                                  tCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTc         >  1:1319003/1‑68 (MQ=255)

                                                  tCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTc         >  1:1613033/1‑68 (MQ=255)

                                                  tCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTc         >  1:1315150/1‑68 (MQ=255)

                                                        tCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCTGGAGCa  <  1:1755852/69‑1 (MQ=255)

                                                        tCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCa  <  1:2844054/69‑1 (MQ=255)

                                                        tCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCa  <  1:2514650/69‑1 (MQ=255)

                                                        tCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCa  <  1:728378/69‑1 (MQ=255)

                                                        tCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCa  <  1:1803993/69‑1 (MQ=255)

                                                        tCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCa  <  1:3427543/69‑1 (MQ=255)

                                                         cTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCa  <  1:1373512/68‑1 (MQ=255)

                                                         cTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCa  <  1:334754/68‑1 (MQ=255)

                                                          tGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAg    <  1:2559661/65‑1 (MQ=255)

                                                          tGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAg    <  1:2844299/65‑1 (MQ=255)

                                                          tGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAg    <  1:961038/65‑1 (MQ=255)

CGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCA  >  minE/511692‑511816

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