BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	150,743	T→G	55.1%	E27A (GAG→GCG)	CJ019_00156 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	150,743	0	T	G	55.1%	‑4.8 / 3.8	18	E27A (GAG→GCG)	CJ019_00156	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (3/4); new base G (8/2); total (11/7)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.11e-01

TCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGC  >  KT2440/150614‑150842

tCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCGCAAgg                                                                                                >  2:191875/1‑134 (MQ=255)

          tGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGGAGGTCGAGCCGGACGTCgg                                                                                      >  5:328788/1‑135 (MQ=255)

                 cACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGCCGCCGGCCGCGCCAGCGGCCGGTGGGGCAAggc                                                                               >  7:334051/1‑135 (MQ=255)

                  aCGCGGGGGTCCAGCTGTTGCTTGACGCGCGGGGGCCGCTTGCCCAGGCCCGCGGCGCGCCGGGCGTGGGAGGCTTTAGTTGGTTGGGGGCCCCGGCCGCCCGCCGGCAGCTCACGCGGCCGGGCGGCCAAggcg                                                                              <  4:52907/135‑1 (MQ=255)

                          gTCCAGCTGTTGGTTTCAGCGGGGGGGACGCTTGCCCTTGGCGTCGGCGAGCCGGGGTGGGGAGGTTTTAGTTGGTTGGGGGCCGAGGCCGTACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGcc                                                                      <  4:151691/135‑1 (MQ=255)

                                       ttGACGCGCGGGGGCCGCTGGCCCATGGCCTCGTCGAGCTGGGCTTGGGCGGCTTTAGTTTGTTGGGGGCCGCGGCCGGACGCCTGCAGCGCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAgg                                                         <  1:245124/135‑1 (MQ=255)

                                       ttGAAGCGCGGTGGCTGGTTGCCGTGGGCGTCGTCGAGCCGGGCTGGGGAGGGTTTAGTTTGTTGGGGGGCGCGGCCGCACGCCTGCAGCTCAAGGTCGTGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAgg                                                         <  4:99741/135‑1 (MQ=255)

                                       tggacgcgcGGGGGCAGGTTGCCCTTGGCGCCGTCGCGCCGGGCTTGGGAGGCTTTCGTTGTGTGGGGGGCGCGCCCGTCCGCCTGCAGCGCAAGCGGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAgg                                                         <  8:189423/130‑1 (MQ=255)

                                             cgcgGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCCGCGCCCGCGGCCGGGCGGCCCCCGCGGGCGTGGCGCGCAGCGGCCGGgcgccg                                                   >  5:335726/1‑133 (MQ=255)

                                               cgGTGGCAGCTTGCCCAGGGCCGCGTCGCGCCGGGCTTGGGCGGCTTTAGTTTGTTGGGGGGCGCGGCCGGACGCCGGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGcc                                                 <  8:92080/135‑1 (MQ=255)

                                                ggTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCGCAAGCGGGCGGGGGGCCCAGGCGGGCGGGCCGCGCCGCGGCAGGGCgccgccc                                                >  1:90115/1‑134 (MQ=255)

                                                                     tcgtcgAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCAGCAGCTCAAGCTGCAGGTCGGCCACGGCGGGCGTGCCGCGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACgcgc                           >  1:94046/1‑135 (MQ=255)

                                                                         cgAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGATGCAGGTCGGCCAAGGCGGGCGTGCCGCGCAGCGGAAGGGCGCAGCCAAGCAGCCCAGCCCGCCCGCGCTGca                       >  6:236400/1‑135 (MQ=255)

                                                                                     gggAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCGGCCGCGCAAGCGGGAGGGCGGCCCAGGCGGGCGGGCCGCGCCGCGTCCGGGCGCCGCCCCCCAGCCCCGCCCGCACGCGCGTCCCAGGGCGAGTAt           >  1:267248/1‑135 (MQ=255)

                                                                                       gAGGCTTTATTTGGTTGGGGGGCGCGGCCGCCCGCCTGAAGCAAACCCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGCCACGCACGCGCTTCACAGTGCGAGTATgg         <  6:163171/135‑1 (MQ=255)

                                                                                             ttAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCCGCGCAAGCTGCAGGTCGGCCCAGGCGGGCGTGCCCCGCACCGTCAGGGCGCCGCCCACCCGCCCCGCCCCCACCCGCGGCACCCTGCGAGTATTGGTTccg   >  8:176049/1‑133 (MQ=255)

                                                                                              tAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCGCAAGCTGCAGGTCGGCCCAGGCGGGCGTGCCGCCCCGCGGCAGGGCGCCGCCCAGCAGCCCGGCACCCCCGCGCGTGACAGTGCGCGTATGGCTTAAg   >  6:73945/1‑134 (MQ=255)

                                                                                               aGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCGCAAGCTGCCCGTCGGCCAACGCGGGCGTGCCGAGCACCGTCACGGCGCCGCCAAACCCCCCGGCACCCACCCGCTTCACAGTGCGAGTAATGCTTAAGc  >  4:135652/1‑134 (MQ=255)

TCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGC  >  KT2440/150614‑150842

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