BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,894,340	C→T	25.0%	W32* (TGG→TGA)	ygcQ ←	predicted flavoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,894,340	0	C	T	25.0%	102.2 / 33.6	56	W32* (TGG→TGA)	ygcQ	predicted flavoprotein
Reads supporting (aligned to +/- strand): ref base C (37/5); new base T (12/2); total (49/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.58e-01

TCGCTGCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTC  >  W3110S.gb/2894274‑2894403

tCGCTGCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGt                                                             >  1:465318/1‑71 (MQ=255)

 cGCTGCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTg                                                            >  1:1166818/1‑71 (MQ=255)

 cGCTGCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTg                                                            >  1:624173/1‑71 (MQ=255)

 cGCTGCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTg                                                            >  1:2525152/1‑71 (MQ=255)

       cACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAg                                                      >  1:1773444/1‑71 (MQ=255)

          tGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAggg                                                    >  1:3739720/1‑70 (MQ=255)

          tGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGAGCCAggg                                                    >  1:2708343/1‑70 (MQ=255)

                   aatgCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCt                                          >  1:1983063/4‑71 (MQ=255)

                   aaCGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCt                                          >  1:2183398/1‑71 (MQ=255)

                           cAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCt                                  >  1:13235/1‑71 (MQ=255)

                             gCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTa                                                  >  1:3892497/1‑53 (MQ=255)

                             gCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTa                                                  >  1:3000023/1‑53 (MQ=255)

                              cACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTgc                                <  1:314410/70‑1 (MQ=255)

                               aCCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTgcgc                              >  1:3810721/1‑71 (MQ=255)

                               aCCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTgcgc                              >  1:3767750/1‑71 (MQ=255)

                               aCCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTgcgc                              >  1:2189511/1‑71 (MQ=255)

                               aCCTGTTCCGCCACTAAAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTgcgc                              >  1:705108/1‑71 (MQ=255)

                                ccTGTTCCGCCACTACAGGCTGAGGTTCAATCTGTCa                                                               >  1:452814/1‑37 (MQ=255)

                                ccTGTTCCGCCACTACAGGCTGAGGTTCAATCTGTCa                                                               >  1:3580764/1‑37 (MQ=255)

                                ccTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGcgcc                             >  1:2283987/1‑71 (MQ=255)

                                ccTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGcgcc                             >  1:3449883/1‑71 (MQ=255)

                                ccTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGcgcc                             >  1:3106508/1‑71 (MQ=255)

                                ccTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGcgcc                             >  1:3094737/1‑71 (MQ=255)

                                      ccGCCACTACAGGCTGAGGTTCAATCTGTCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCgccagc                       >  1:2689525/1‑71 (MQ=255)

                                      ccGCCACTACAGGCTGAGGTTCAATCTGTCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCgccagc                       >  1:2840444/1‑71 (MQ=255)

                                      ccGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCgccagc                       >  1:962481/1‑71 (MQ=255)

                                      ccGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCgccagc                       >  1:1291552/1‑71 (MQ=255)

                                            cTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCtt                 >  1:720553/1‑71 (MQ=255)

                                            cTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCtt                 <  1:653223/71‑1 (MQ=255)

                                              aCAGGCTGAGGTTCAATCTGTCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGc               >  1:2449020/1‑71 (MQ=255)

                                              aCAGGCTGAGGTTCAATCTGTCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGc               >  1:3197256/1‑71 (MQ=255)

                                              aCAGGCTGAGGTTCAATCTGTCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGc               >  1:2928617/1‑71 (MQ=255)

                                              aCAGGCTGAGGTTCAATCTGTCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGc               >  1:1424865/1‑71 (MQ=255)

                                              aCAGGCTGAGGTTCAATCTGTCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGc               >  1:2147828/1‑71 (MQ=255)

                                              aCAGGCTGAGGTTCAATCTGTCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGc               >  1:438777/1‑71 (MQ=255)

                                               cAGGCTGAGGTTCAATCTGTCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGc               <  1:736341/70‑1 (MQ=255)

                                               cAGGCTGAGGTTCAATCTGTCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGc               <  1:2533761/70‑1 (MQ=255)

                                               cAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGaaa                                      >  1:4024851/1‑47 (MQ=255)

                                               cAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGc               >  1:4007543/1‑70 (MQ=255)

                                               cAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGc               >  1:3600611/1‑70 (MQ=255)

                                               cAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGc               <  1:849219/70‑1 (MQ=255)

                                               cAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTAAAACCGCTGaaa                                      >  1:911261/1‑47 (MQ=255)

                                               cAGGCTGAGGTTCAATCTGCCAGTGAGCCAGGGTACAACCGCTGaaa                                      >  1:2579877/1‑47 (MQ=255)

                                                 ggCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGa                                        <  1:917689/43‑1 (MQ=255)

                                                      aGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCgg       >  1:3509956/1‑71 (MQ=255)

                                                        gTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGcgccgc                           <  1:2522804/49‑1 (MQ=255)

                                                         ttCAATCTGTCAGTGCGCCAGGGTACAACAGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCAt    >  1:1693393/1‑71 (MQ=255)

                                                         ttCAATATGTCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCAt    >  1:524777/1‑71 (MQ=255)

                                                           cAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATtc  >  1:3520189/1‑71 (MQ=255)

                                                           cAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATtc  >  1:2814042/1‑71 (MQ=255)

                                                           cAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCAGCCAGCAAGCTTGCTATTGCGGCAtt   >  1:3235206/1‑70 (MQ=255)

                                                                tGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCa     >  1:249895/1‑63 (MQ=255)

                                                                tGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCa     >  1:495710/1‑63 (MQ=255)

                                                                tGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCa     >  1:2220440/1‑63 (MQ=255)

                                                                tGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCa     >  1:3038178/1‑63 (MQ=255)

                                                                tGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCa     >  1:1157877/1‑63 (MQ=255)

TCGCTGCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTC  >  W3110S.gb/2894274‑2894403

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