BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	115,752	C→T	39.9%	A83A (GCG→GCA)	yacH ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	115,752	0	C	T	39.9%	28.3 / 37.0	60	A83A (GCG→GCA)	yacH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (9/27); new base T (4/20); total (13/47)
Fisher's exact test for biased strand distribution p-value = 5.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACGG  >  minE/115689‑115813

tGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGa                                                           >  1:1185824/1‑68 (MQ=255)

tGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGa                                                           >  1:880100/1‑68 (MQ=255)

tGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGa                                                           >  1:1913621/1‑68 (MQ=255)

tGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGa                                                           >  1:2380867/1‑68 (MQ=255)

tGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGa                                                           >  1:126973/1‑68 (MQ=255)

tGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGa                                                           >  1:2296664/1‑68 (MQ=255)

tGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGa                                                           >  1:2019721/1‑68 (MQ=255)

                   aTCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGACCCATTGATCTATTTGTGc                                        <  1:2292634/68‑1 (MQ=255)

                   aTCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGACCCATTGATCTATTTGTGc                                        <  1:1927186/68‑1 (MQ=255)

                   aTCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGACCCATTGATCTATTTGTGc                                        <  1:4942/68‑1 (MQ=255)

                             gCGAAAGTAGGGCGTCGGGGTATATCTCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                             <  1:2668343/69‑1 (MQ=255)

                             gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCATTGCTGAATg                             <  1:3387176/69‑1 (MQ=255)

                             gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                             <  1:2187023/69‑1 (MQ=255)

                             gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                             <  1:2282317/69‑1 (MQ=255)

                             gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                             <  1:1974296/69‑1 (MQ=255)

                             gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                             <  1:1945015/69‑1 (MQ=255)

                             gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                             <  1:2169835/69‑1 (MQ=255)

                             gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                             <  1:2178574/69‑1 (MQ=255)

                             gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                             <  1:2819282/69‑1 (MQ=255)

                             gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                             <  1:1840732/69‑1 (MQ=255)

                             gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                             <  1:1839019/69‑1 (MQ=255)

                             gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                             <  1:1766218/69‑1 (MQ=255)

                             gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                             <  1:288219/69‑1 (MQ=255)

                             gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                             <  1:1539974/69‑1 (MQ=255)

                             gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                             <  1:3051442/69‑1 (MQ=255)

                             gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                             <  1:3081135/69‑1 (MQ=255)

                             gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                             <  1:3182131/69‑1 (MQ=255)

                             gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                             <  1:654405/69‑1 (MQ=255)

                             gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                             <  1:7026/69‑1 (MQ=255)

                               gAAAGTAGGGCGTCGGGGTATAGCGCGAGAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                             <  1:973161/67‑1 (MQ=255)

                                    tAGGGCGTCGGGATACAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGAttt                      >  1:225097/1‑69 (MQ=255)

                                     aGGGCGTCGGGATACAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTg                     >  1:2525818/1‑69 (MQ=255)

                                             ggggTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATGCGGATTTGACGAcggc             >  1:2349546/1‑69 (MQ=255)

                                             ggggTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATGCGGATTTGACGAcggc             >  1:1752683/1‑69 (MQ=255)

                                             ggggTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATGCGGATTTGACGAcggc             >  1:2646863/1‑69 (MQ=255)

                                             ggggTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGAGCTGAATGCGGATTTGACGAcggc             >  1:1924490/1‑69 (MQ=255)

                                             gggATACAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                         <  1:1161205/41‑1 (MQ=255)

                                             gggATACAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                         <  1:911778/41‑1 (MQ=255)

                                             gggATACAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                         <  1:2646368/41‑1 (MQ=255)

                                             gggATACAGCGCGA‑GGGCGCGACCCATTGATCTATTtgtg                                         <  1:3398690/40‑1 (MQ=39)

                                               gATACAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                         <  1:2672119/39‑1 (MQ=255)

                                                  aCAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg        <  1:2446344/69‑1 (MQ=255)

                                                  aCAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg        <  1:1914786/69‑1 (MQ=255)

                                                  aCAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg        <  1:1878231/69‑1 (MQ=255)

                                                  aCAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg        <  1:2362157/69‑1 (MQ=255)

                                                  aCAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg        <  1:144692/69‑1 (MQ=255)

                                                  aCAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg        <  1:14412/69‑1 (MQ=255)

                                                  aCAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg        <  1:1437560/69‑1 (MQ=255)

                                                  aCAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg        <  1:551754/69‑1 (MQ=255)

                                                  aCAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg        <  1:1349907/69‑1 (MQ=255)

                                                   cAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                         <  1:2086310/35‑1 (MQ=255)

                                                   cAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                         <  1:1193323/35‑1 (MQ=255)

                                                   cAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                         <  1:1208716/35‑1 (MQ=255)

                                                   cAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                         <  1:1219932/35‑1 (MQ=255)

                                                   cAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                         <  1:393848/35‑1 (MQ=255)

                                                   cAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                         <  1:3407067/35‑1 (MQ=255)

                                                   cAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                         <  1:1902935/35‑1 (MQ=255)

                                                   cAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg        <  1:1400908/68‑1 (MQ=255)

                                                        cgACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACgg  <  1:2037405/69‑1 (MQ=255)

                                                        cgACGGGCGAGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg        <  1:2496546/63‑1 (MQ=255)

TGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACGG  >  minE/115689‑115813

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