BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	125,390	A→G	27.6%	H192H (CAT→CAC)	gcd ←	glucose dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	125,390	0	A	G	27.6%	116.8 / 47.3	112	H192H (CAT→CAC)	gcd	glucose dehydrogenase
Reads supporting (aligned to +/- strand): ref base A (42/39); new base G (16/15); total (58/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

Only 100 of 112 total aligned reads displayed.

TCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTAC  >  minE/125328‑125448

tCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGAc                                                       >  1:2622750/1‑68 (MQ=255)

tCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGAc                                                       >  1:2454404/1‑68 (MQ=255)

tCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGAc                                                       >  1:2421276/1‑68 (MQ=255)

tCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGAc                                                       >  1:3072961/1‑68 (MQ=255)

tCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGAc                                                       >  1:2312628/1‑68 (MQ=255)

tCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGAc                                                       >  1:1650022/1‑68 (MQ=255)

tCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGAc                                                       >  1:1194412/1‑68 (MQ=255)

tCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGAc                                                       >  1:838609/1‑68 (MQ=255)

  aCCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGt                                                     <  1:1117799/68‑1 (MQ=255)

  aCCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGt                                                     <  1:286724/68‑1 (MQ=255)

  aCCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGt                                                     <  1:2292375/68‑1 (MQ=255)

  aCCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGt                                                     <  1:215269/68‑1 (MQ=255)

  aCCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGt                                                     <  1:3238851/68‑1 (MQ=255)

  aCCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGt                                                     <  1:3303751/68‑1 (MQ=255)

  aCCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGt                                                     <  1:115778/68‑1 (MQ=255)

  aCCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGt                                                     <  1:3381895/68‑1 (MQ=255)

  aCCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGt                                                     <  1:148951/68‑1 (MQ=255)

  aCCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGt                                                     <  1:1486294/68‑1 (MQ=255)

  aCCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGt                                                     <  1:592047/68‑1 (MQ=255)

   cccGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGtt                                                    <  1:2603887/68‑1 (MQ=255)

    ccGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGt                                                     <  1:3308856/66‑1 (MQ=255)

                gctgctTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAAtt                                       <  1:1368783/68‑1 (MQ=255)

                       cacaTCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATc                                                 >  1:3099726/1‑51 (MQ=255)

                       cacaTCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATc                                                 >  1:1109638/1‑51 (MQ=255)

                       cacaTCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATc                                                 >  1:1763777/1‑51 (MQ=255)

                       cacaTCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATc                                                 >  1:716020/1‑51 (MQ=255)

                       cacaTCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATc                                                 >  1:2023254/1‑51 (MQ=255)

                       cacaTCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATc                                                 >  1:2452610/1‑51 (MQ=255)

                       cacaTCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATc                                                 >  1:3053474/1‑51 (MQ=255)

                       cacaTCGCCAGTACGGAACACCCAGGCTTCTTT‑‑GATTATGGACGTTAt                                                  >  1:3406218/1‑48 (MQ=37)

                                 gTACGGAACACCCAGGCTTCTTTCAGATTATTGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAg                     <  1:2765523/69‑1 (MQ=255)

                                 gTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAg                     <  1:2740856/69‑1 (MQ=255)

                                 gTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAg                     <  1:1970540/69‑1 (MQ=255)

                                 gTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAg                     <  1:3212579/69‑1 (MQ=255)

                                 gTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAg                     <  1:237929/69‑1 (MQ=255)

                                 gTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAg                     <  1:536265/69‑1 (MQ=255)

                                 gTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAg                     <  1:1615944/69‑1 (MQ=255)

                                 gTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAg                     <  1:2487930/69‑1 (MQ=255)

                                 gTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAg                     <  1:358125/69‑1 (MQ=255)

                                 gTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAg                     <  1:1448472/69‑1 (MQ=255)

                                  tACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATc                                                 >  1:1784138/1‑40 (MQ=255)

                                  tACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATc                                                 >  1:109427/1‑40 (MQ=255)

                                  tACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATc                                                 >  1:3215073/1‑40 (MQ=255)

                                  tACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATc                                                 >  1:2773897/1‑40 (MQ=255)

                                  tACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATc                                                 >  1:55523/1‑40 (MQ=255)

                                   aCGGAACCCCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                   >  1:882678/1‑69 (MQ=255)

                                   aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTa                                          >  1:128164/1‑46 (MQ=255)

                                   aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCGGTGGTGAAAAGCg                   >  1:2630930/1‑69 (MQ=255)

                                   aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAgt                              >  1:2902718/1‑57 (MQ=255)

                                   aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGc                    >  1:859148/1‑68 (MQ=255)

                                   aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                   >  1:28759/1‑69 (MQ=255)

                                   aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                   >  1:3065593/1‑69 (MQ=255)

                                   aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                   >  1:2826107/1‑69 (MQ=255)

                                   aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                   >  1:3380948/1‑69 (MQ=255)

                                   aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                   >  1:460851/1‑69 (MQ=255)

                                   aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                   >  1:1056445/1‑69 (MQ=255)

                                   aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                   >  1:1045577/1‑69 (MQ=255)

                                   aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                   >  1:2377523/1‑69 (MQ=255)

                                   aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                   >  1:2494485/1‑69 (MQ=255)

                                   aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                   >  1:1119219/1‑69 (MQ=255)

                                   aCGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCa                                >  1:1259497/1‑56 (MQ=255)

                                   aCGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCg                   >  1:792122/1‑69 (MQ=255)

                                   aCGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCg                   >  1:2426253/1‑69 (MQ=255)

                                    cGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCg                   >  1:1387996/1‑68 (MQ=255)

                                    cGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCg                   >  1:1499356/1‑68 (MQ=255)

                                    cGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCg                   >  1:1904489/1‑68 (MQ=255)

                                    cGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCg                   >  1:2983539/1‑68 (MQ=255)

                                    cGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCg                   >  1:2305292/1‑68 (MQ=255)

                                       aaCACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGa               <  1:810253/69‑1 (MQ=255)

                                       aaCACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGa               <  1:66677/69‑1 (MQ=255)

                                       aaCACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGa               <  1:1971037/69‑1 (MQ=255)

                                       aaCACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGa               <  1:2680212/69‑1 (MQ=255)

                                       aaCACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGa               <  1:2271238/69‑1 (MQ=255)

                                       aaCACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGa               <  1:2625045/69‑1 (MQ=255)

                                       aaCACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGa               <  1:2358052/69‑1 (MQ=255)

                                       aaCACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGa               <  1:3146023/69‑1 (MQ=255)

                                                 cttcttTCAGATTATGGACGTTATCGGCGTTAATTTGt                                    >  1:2565070/1‑38 (MQ=255)

                                                 cttcttTCAGATTATGGACGTTATCGGCGTTAATTTGTTTc                                 >  1:1076942/1‑41 (MQ=255)

                                                 cttcttTCAGATTATGGACGTTATCGGCGTTAATTTGTTTc                                 >  1:1220988/1‑41 (MQ=255)

                                                 cttcttTCAGATTATGGACGTTATCGGCGTTAATTTGTTTc                                 >  1:2947691/1‑41 (MQ=255)

                                                 cttcttTCAGATTATGGACGTTATCGGCGTTAATTTGTTTc                                 >  1:68826/1‑41 (MQ=255)

                                                 cttcttTCAGATTATGGACGTTATCGGCGTTAATTTGTTTc                                 >  1:2825725/1‑41 (MQ=255)

                                                 cttcttTCAGATTATGGACGTTATCGGCGTTAATTTGTTTc                                 >  1:1477887/1‑41 (MQ=255)

                                                 cttcttTCAGATTATGGACGTTATCGGCGTTAATTTGTTTc                                 >  1:1251358/1‑41 (MQ=255)

                                                   tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg   <  1:2499633/69‑2 (MQ=255)

                                                   tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg   <  1:2592203/69‑2 (MQ=255)

                                                   tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg   <  1:2588433/69‑2 (MQ=255)

                                                   tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg   <  1:689254/69‑2 (MQ=255)

                                                   tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg   <  1:547953/69‑2 (MQ=255)

                                                   tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg   <  1:1402113/69‑2 (MQ=255)

                                                   tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg   <  1:2434310/69‑2 (MQ=255)

                                                   tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg   <  1:1521373/69‑2 (MQ=255)

                                                   tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg   <  1:1942959/69‑2 (MQ=255)

                                                   tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg   <  1:325195/69‑2 (MQ=255)

                                                   tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg   <  1:299409/69‑2 (MQ=255)

                                                   tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg   <  1:3221688/69‑2 (MQ=255)

                                                      ttCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg   <  1:1353458/66‑2 (MQ=255)

                                                      ttCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTAc  <  1:3369575/67‑1 (MQ=255)

                                                      ttCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTAc  <  1:2620457/67‑1 (MQ=255)

                                                      ttCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTAc  <  1:2330705/67‑1 (MQ=255)

TCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTAC  >  minE/125328‑125448

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