BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	62,818	T→G	55.1%	N214T (AAC→ACC) ‡	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	62,818	0	T	G	55.1%	18.0 / 46.5	49	N214T (AAC→ACC) ‡	polB	DNA polymerase II
Reads supporting (aligned to +/- strand): ref base T (18/4); new base G (19/8); total (37/12)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACG  >  minE/62752‑62868

tGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGtt                                                   >  1:2533482/1‑68 (MQ=255)

tGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGtt                                                   >  1:1105027/1‑68 (MQ=255)

tGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGtt                                                   >  1:1410492/1‑68 (MQ=255)

tGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGtt                                                   >  1:1735555/1‑68 (MQ=255)

  tttttGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTgg                                                 >  1:2904841/1‑68 (MQ=255)

  tttttGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTgg                                                 >  1:1281197/1‑68 (MQ=255)

  tttttGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTgg                                                 >  1:163195/1‑68 (MQ=255)

          gcaTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAg                                         >  1:550040/1‑68 (MQ=255)

           caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                       >  1:2657457/1‑69 (MQ=255)

           caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                       >  1:2765880/1‑69 (MQ=255)

           caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                       >  1:3132254/1‑69 (MQ=255)

           caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                       >  1:2291475/1‑69 (MQ=255)

           caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                       >  1:385853/1‑69 (MQ=255)

           caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                       >  1:2173530/1‑69 (MQ=255)

           caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                       >  1:2135923/1‑69 (MQ=255)

           caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                       >  1:2133030/1‑69 (MQ=255)

           caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                       >  1:2092325/1‑69 (MQ=255)

           caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                       >  1:53269/1‑69 (MQ=255)

           caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                       >  1:714925/1‑69 (MQ=255)

           caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                       >  1:1917570/1‑69 (MQ=255)

           caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                       >  1:732424/1‑69 (MQ=255)

           caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                       >  1:1156732/1‑69 (MQ=255)

           caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                       >  1:1129254/1‑69 (MQ=255)

           caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                       >  1:955319/1‑69 (MQ=255)

           caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                       >  1:981789/1‑69 (MQ=255)

           caTTCGCAGATCGAACTGCACCACGTTCCAACAGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                       >  1:82231/1‑69 (MQ=255)

           caTTCGCAGATAGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                       >  1:2148703/1‑69 (MQ=255)

                     tCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGtttt                             >  1:1044803/1‑69 (MQ=255)

                     tCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGtttt                             >  1:1055302/1‑69 (MQ=255)

                      cGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGTCAAACCAGGCGTTGAGtttt                             >  1:2711010/1‑68 (MQ=255)

                      cGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGtttt                             >  1:1028482/1‑68 (MQ=255)

                      cGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGtttt                             >  1:624519/1‑68 (MQ=255)

                      cGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGtttt                             >  1:2062702/1‑68 (MQ=255)

                      cGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGtttt                             >  1:2269854/1‑68 (MQ=255)

                      cGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGtttt                             >  1:2344173/1‑68 (MQ=255)

                      cGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGtttt                             >  1:2497394/1‑68 (MQ=255)

                      cGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGtttt                             >  1:2774838/1‑68 (MQ=255)

                      cGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGttt                              >  1:2578415/1‑67 (MQ=255)

                           tGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGc                       <  1:877215/69‑1 (MQ=255)

                           tGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGc                       <  1:737671/69‑1 (MQ=255)

                                 acGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTgc                 >  1:353229/1‑69 (MQ=255)

                                 acGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTgc                 >  1:2013876/1‑69 (MQ=255)

                                             atgatCACATCAGGATCATGAGTCGCAAACCAGGCGTTAAGTTTTTCCAGTAGCAGCGGGCGGCTGGCg     <  1:2671661/69‑1 (MQ=255)

                                             atgatCACATCAGGATCATGAGTCGCAAACCAGGCGTTAAGTTTTTCCAGTAGCAGCGGGCGGCTGGCg     <  1:3021665/69‑1 (MQ=255)

                                             atgatCACATCAGGATCATGAGTCGCAAACCAGGCGTTAAGTTTTTCCAGTAGCAGCGGGCGGCTGGCg     <  1:901390/69‑1 (MQ=255)

                                             atgatCACATCAGGATCATGAGTCGCAAACCAGGCGTTAAGTTTTTCCAGTAGCAGCGGGCGGCTGGCg     <  1:1490062/69‑1 (MQ=255)

                                             atgatCACATCAGGATCATGAGTCGCAAACCAGGCGTTAAGTTTTTCCAGTAGCAGCGGGCGGCTGGCg     <  1:1516117/69‑1 (MQ=255)

                                             atgatCACATCAGGATCATGAGTCGCAAACCAGGCGTTAAGTTTTTCCAGTAGCAGCGGGCGGCTGGCg     <  1:2060673/69‑1 (MQ=255)

                                             atgatCACATCAGGATCATGAGTCGCAAACCAGGCGTTAAGTTTTTCCAGTAGCAGCGGGCGGCTGGCg     <  1:28236/69‑1 (MQ=255)

                                             atgatCACATCAGGATCATGAGTCGCAAACCAGGCGTTAAGTTTTTCCAGTAGCAGCGGGCGGCTGGCg     <  1:3284486/69‑1 (MQ=255)

                                                atCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACg  <  1:3272295/69‑1 (MQ=255)

                                                atCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACg  <  1:3270849/69‑1 (MQ=255)

                                                  cacaTCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTgcggg              >  1:324677/1‑55 (MQ=255)

TGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACG  >  minE/62752‑62868

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