BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	512,666	C→T	47.0%	V357V (GTG→GTA)	ybhC ←	predicted pectinesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	512,666	0	C	T	47.0%	6.6 / 56.3	68	V357V (GTG→GTA)	ybhC	predicted pectinesterase
Reads supporting (aligned to +/- strand): ref base C (23/13); new base T (18/14); total (41/27)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AGATCACCGCATCGGCCCACGGTTTAGCCGTGTTAAAACCTTCGTTGATGGCGCTATCACGGATCACCACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGCCCAGTTGCGCCACACCATC  >  minE/512599‑512725

agaTCACCGCATCGGCCCACGGTTTAGCCGTGTTAAAACCTTCGTTGATGGCGCTATCACGGATcacc                                                             <  1:3267468/68‑1 (MQ=255)

            cGGCCCACGGTTTAGCCGTGTTAAAACCTTCGTTGATGGCGCTATCACGGATCACCACCTGACCGTTgg                                                <  1:1969827/69‑1 (MQ=255)

                          gCCGTGTTAAAACCTTCGTTGATGGCGCTATCACGGATCACCACCTGACCGTTGGTATTgg                                          >  1:3182492/1‑61 (MQ=255)

                          gCCGTGTTAAAACCTTCGTTGATGGCGCTATCACGGATCACCACCTGACCGTTGGTATTGGCATCAACa                                  >  1:96018/1‑69 (MQ=255)

                          gCCGTGTTAAAACCTTCGTTGATGGCGCTATCACGGATCACCACCTGACCGTTGGTATTGGCATCAACa                                  >  1:1393359/1‑69 (MQ=255)

                          gCCGTGTTAAAACCTTCGTTGATGGCGCTATCACGGATCACCACCTGACCGTTGGTATTGGCATCAACa                                  >  1:2509818/1‑69 (MQ=255)

                          gCCGTGTTAAAACCTTCGTTGATGGCGCTATCACGGATCACCACCTGACCGTTGGTATTGGCATCAACa                                  >  1:1414426/1‑69 (MQ=255)

                                ttAAAACCTTCGTTGATGGCGCTATCACGGATCACCACCTGACCGTTgg                                                <  1:2841542/49‑1 (MQ=255)

                                ttAAAACCTTCGTTGATGGCGCTATCACGGATCACCACCTGACCGTTgg                                                <  1:2524747/49‑1 (MQ=255)

                                ttAAAACCTTCGTTGATGGCGCTATCACGGATCACCACCTGACCGTTgg                                                <  1:245025/49‑1 (MQ=255)

                                ttAAAACCTTCGTTGATGGCGCTATCACGGATCACCACCTGACCGTTgg                                                <  1:1489969/49‑1 (MQ=255)

                                ttAAAACCTTCGTTGATGGCGCTATCACGGATCACCACCTGACCGTTgg                                                <  1:1645030/49‑1 (MQ=255)

                                ttAAAACCTTCGTTGATGGCGCTATCACGGATCACCACCTGACCGTTgg                                                <  1:1583606/49‑1 (MQ=255)

                                ttAAAACCTTCGTTGATGGCGCTATCACGGATCACCACCTGACCGTTgg                                                <  1:1586017/49‑1 (MQ=255)

                                  aaaaCCTTCGTTGATGGCGCTATGACGGATCACTACCTGACCGTTGGTATTGGCATCAACATCCagcga                          >  1:3142215/1‑69 (MQ=255)

                                  aaaaCCTTCGTTGATGGCGCTATCACGGATCAGTACCTGACCGTTGGTATTGGCATCAACATCCagcga                          >  1:404508/1‑69 (MQ=255)

                                  aaaaCCTTCGTTGATGGCGCTATCACGGATCACTACCTGACCGTTGGTATTGGCa                                        >  1:507519/1‑55 (MQ=255)

                                  aaaaCCTTCGTTGATGGCGCTATCACGGATCACTACCTGACCGTTGGTATTGGCATCaa                                    >  1:233733/1‑59 (MQ=255)

                                  aaaaCCTTCGTTGATGGCGCTATCACGGATCACTACCTGACCGTTGGTATTGGCATCAACATCCagcga                          >  1:98478/1‑69 (MQ=255)

                                  aaaaCCTTCGTTGATGGCGCTATCACGGATCACTACCTGACCGTTGGTATTGGCATCAACATCCagcga                          >  1:3423157/1‑69 (MQ=255)

                                  aaaaCCTTCGTTGATGGCGCTATCACGGATCACTACCTGACCGTTGGTATTGGCATCAACATCCagcga                          >  1:2512100/1‑69 (MQ=255)

                                  aaaaCCTTCGTTGATGGCGCTATCACGGATCACTACCTGACCGTTGGTATTGGCATCAACATCCagcga                          >  1:1401109/1‑69 (MQ=255)

                                  aaaaCCTTCGTTGATGGCGCTATCACGGATCACTACCTGACCGTTGGTATTGGCATCAACATCCagcga                          >  1:367888/1‑69 (MQ=255)

                                  aaaaCCTTCGTTGATGGCGCTATCACGGATCACTACCTGACCGTTGGTATTGGCATCAACATCCagcga                          >  1:249602/1‑69 (MQ=255)

                                  aaaaCCTTCGTTGATGGCGCTATCACGGATCACTACCTGACCGTTGGTATTGGCATCAACATCCagcga                          >  1:2933277/1‑69 (MQ=255)

                                  aaaaCCTTCGTTGATGGCGCTATCACGGATCACTACCTGACCGTTGGTATTGGCATCAACATCCagcga                          >  1:2962057/1‑69 (MQ=255)

                                  aaaaCCTTCGTTGATGGCGCTATCACGGATCACTACCTGACCGTTGGTATTGGCATCAACATCCagcga                          >  1:1807201/1‑69 (MQ=255)

                                  aaaaCCTTCGTTGATGGCGCTATCACGGATCACTACCTGACCGTTGGTATTGGCATCAACATCCagcga                          >  1:628636/1‑69 (MQ=255)

                                  aaaaCCTTCGTTGATGGCGCTATCACGGATCACTACCTGACCGTTGGTATTGGCATCAACATCCagcga                          >  1:3078428/1‑69 (MQ=255)

                                  aaaaCCTTCGTTGATGGCGCTATCACGGATCACTACCTGACCGTTGGTATTGGCATCAACATCCagcga                          >  1:1585818/1‑69 (MQ=255)

                                  aaaaCCTTCGTTGATGGCGCTATCACGGATCACTACCTGACCGTTGGTATTGGCATCAACATCCagcg                           >  1:67750/1‑68 (MQ=255)

                                  aaaaCCTTCGTTGATGGCGCTATCACGGATCACTACCAGACCGTTGGTATTGGCATCAACATCCagcga                          >  1:974012/1‑69 (MQ=255)

                                        ttCGTTGATGGCGCTATCACGGATCACTACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGcc                    <  1:2958628/69‑1 (MQ=255)

                                        ttCGTTGATGGCGCTATCACGGATCACTACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGcc                    <  1:1115924/69‑1 (MQ=255)

                                        ttCGTTGATGGCGCTATCACGGATCACTACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGcc                    <  1:1016250/69‑1 (MQ=255)

                                        ttCGTTGATGGCGCTATCACGGATCACTACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGcc                    <  1:2871089/69‑1 (MQ=255)

                                        ttCGTTGATGGCGCTATCACGGATCACTACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGcc                    <  1:1632037/69‑1 (MQ=255)

                                        ttCGTTGATGGCGCTATCACGGATCACTACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGcc                    <  1:3272159/69‑1 (MQ=255)

                                        ttCGTTGATGGCGCTATCACGGATCACTACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGcc                    <  1:354914/69‑1 (MQ=255)

                                        ttCGTTGATGGCGCTATCACGGATCACTACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGcc                    <  1:2487616/69‑1 (MQ=255)

                                        ttCGTTGATGGCGCTATCACGGATCACTACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGcc                    <  1:1403709/69‑1 (MQ=255)

                                        ttCGTTGATGGCGCTATCACGGATCACTACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGcc                    <  1:2263767/69‑1 (MQ=255)

                                        ttCGTTGATGGCGCTATCACGGATCACTACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGcc                    <  1:1504829/69‑1 (MQ=255)

                                        ttCGTTGATGGCGCTATCACCGATCACTACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGcc                    <  1:1764993/69‑1 (MQ=255)

                                         tCGTTGATGGCGCTATCACGGATCACTACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGcc                    <  1:929118/68‑1 (MQ=255)

                                         tCGTTGATGGCGCTATCACGGATCACCACCTGACCGTTGGTATTGGCATCa                                     >  1:3129444/1‑51 (MQ=255)

                                         tCGTTGATGGCGCTATCACGGATCACCACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGccc                   >  1:1084554/1‑69 (MQ=255)

                                         tCGTTGATGGCGCTATCACGGATCACCACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGccc                   >  1:1154530/1‑69 (MQ=255)

                                         tCGTTGATGGCGCTATCACGGATCACCACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGc                     >  1:720843/1‑67 (MQ=255)

                                            ttGATGGCGCTATCACGGATCACTACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGcc                    <  1:3135126/65‑1 (MQ=255)

                                                 ggCGCTATCACGGATCACCACCTGACCGGTGGTATTGGCATCAACATCCAGCGAGCGGCCCAGTTgcgc           <  1:2379445/69‑1 (MQ=255)

                                                       aTCACGGATCACCACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGCCCAGTTGCGCCACAcc     <  1:193246/69‑1 (MQ=255)

                                                          aCGGATCACCACCTGACCGTTGGTATTGGCATCAACATCCagcgag                         >  1:1664532/1‑46 (MQ=255)

                                                          aCGGATCACCACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGCCCAGt                <  1:1789804/55‑1 (MQ=255)

                                                          aCGGATCACCACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGCCCAGt                <  1:3022594/55‑1 (MQ=255)

                                                          aCGGATCACCACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGCCCAGt                >  1:1305472/1‑55 (MQ=255)

                                                          aCGGATCACCACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGCCCAGt                >  1:1501655/1‑55 (MQ=255)

                                                          aCGGATCACCACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGCCCAGt                >  1:752569/1‑55 (MQ=255)

                                                          aCGGATCACCACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGCCCAGt                <  1:724414/55‑1 (MQ=255)

                                                          aCGGATCACCACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGCCCAGt                >  1:2532312/1‑55 (MQ=255)

                                                          aCGGATCACCACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGCCCAGTTGCGCCACACCATc  >  1:2932800/1‑69 (MQ=255)

                                                          aCGGATCACCACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGCCCAGTTGCGCCACACCATc  >  1:2825947/1‑69 (MQ=255)

                                                          aCGGATCACCACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGCCCAGTTGCGCCACACCATc  >  1:2509711/1‑69 (MQ=255)

                                                          aCGGATCACCACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGCCCAGTTGCGCCACACCATc  >  1:3249862/1‑69 (MQ=255)

                                                          aCGGATCACCACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGCCCAGTTGCGCCACACCATc  >  1:2571982/1‑69 (MQ=255)

                                                          aCGGATCACCACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGCCCAGTTGCGCCACACCATc  >  1:3418353/1‑69 (MQ=255)

                                                          aCGGATCACCACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGCCCAGTTGCGCCACACCATc  >  1:3348499/1‑69 (MQ=255)

                                                          aCGGATCACCACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGCCCAGTTGCGCCACACCATc  >  1:1285952/1‑69 (MQ=255)

                                                          aCGGATCACCACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGCCCAGTTGCGCCACACCATc  >  1:3289823/1‑69 (MQ=255)

AGATCACCGCATCGGCCCACGGTTTAGCCGTGTTAAAACCTTCGTTGATGGCGCTATCACGGATCACCACCTGACCGTTGGTATTGGCATCAACATCCAGCGAGCGGCCCAGTTGCGCCACACCATC  >  minE/512599‑512725

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