BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	125,387	A→C	26.5%	N193K (AAT→AAG)	gcd ←	glucose dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	125,387	0	A	C	26.5%	149.9 / 65.7	117	N193K (AAT→AAG)	gcd	glucose dehydrogenase
Reads supporting (aligned to +/- strand): ref base A (44/42); new base C (16/15); total (60/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 117 total aligned reads displayed.

TGGTGATTTCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTAC  >  minE/125320‑125448

tGGTGATTTCACCCGGATCGTTCTGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGAt                                                              <  1:2773703/69‑1 (MQ=255)

tGGTGATTTCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGAt                                                              <  1:2074351/69‑1 (MQ=255)

 ggTGATTTCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGAtt                                                             >  1:2559424/1‑69 (MQ=255)

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

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

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

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

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

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

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

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

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

        tCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGAc                                                       >  1:2454404/1‑68 (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:286724/68‑1 (MQ=255)

          aCCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGt                                                     <  1:3303751/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)

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

          aCCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGt                                                     <  1:1117799/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:716020/1‑51 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                           aCGGAACCCCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                   >  1:882678/1‑69 (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:2576589/1‑69 (MQ=255)

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

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

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

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

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

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

                                           aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                   >  1:1045577/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:2983539/1‑68 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                            cttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg   <  1:3214108/68‑2 (MQ=255)

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

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

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

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

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

TGGTGATTTCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTAC  >  minE/125320‑125448

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