BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	158,626	G→A	13.0%	S167S (AGC→AGT)	pcnB ←	poly(A) polymerase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	158,626	0	G	A	13.0%	127.0 / 10.0	46	S167S (AGC→AGT)	pcnB	poly(A) polymerase I
Reads supporting (aligned to +/- strand): ref base G (7/33); new base A (4/2); total (11/35)
Fisher's exact test for biased strand distribution p-value = 2.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

TCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAA  >  W3110S.gb/158561‑158694

tCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGtt                                                                 <  1:1599618/71‑1 (MQ=255)

tCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGtt                                                                 <  1:24175/71‑1 (MQ=255)

tCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGtt                                                                 <  1:1889085/71‑1 (MQ=255)

tCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGtt                                                                 <  1:367740/71‑1 (MQ=255)

tCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGtt                                                                 <  1:128242/71‑1 (MQ=255)

tCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGtt                                                                 <  1:910509/71‑1 (MQ=255)

 ccTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGtt                                                                 <  1:874065/70‑1 (MQ=255)

  cTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGa                                                               <  1:828538/71‑1 (MQ=255)

  cTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGa                                                               <  1:882995/71‑1 (MQ=255)

         tCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGa                                                         >  1:1274893/1‑70 (MQ=255)

           cTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGAAAt                                                      >  1:345186/1‑71 (MQ=255)

           cTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGAAAt                                                      >  1:664579/1‑71 (MQ=255)

           cTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGAAAt                                                      >  1:1144120/1‑71 (MQ=255)

           cTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGAAAt                                                      >  1:1716981/1‑71 (MQ=255)

                          gTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGaaa                                                       >  1:823074/1‑55 (MQ=255)

                          gTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGaaa                                                       >  1:1073185/1‑55 (MQ=255)

                          gTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGaaa                                                       >  1:1723104/1‑55 (MQ=255)

                          gTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGaaa                                                       >  1:1380477/1‑55 (MQ=255)

                              tCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGtcttc                                   <  1:198299/71‑1 (MQ=255)

                              tCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGtcttc                                   <  1:225593/71‑1 (MQ=255)

                              tCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGtcttc                                   <  1:1189163/71‑1 (MQ=255)

                                 cggacggTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCttcttc                                <  1:1554865/71‑1 (MQ=255)

                                 cggacggTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCttcttc                                <  1:770228/71‑1 (MQ=255)

                                  ggacggTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCttcttc                                <  1:1156820/70‑1 (MQ=255)

                                   tacggTAAAATCCGCTACGCTGTAATAGAGGCTGTTGTTAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGa                              <  1:1595824/70‑1 (MQ=255)

                                   gacggTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGa                              <  1:646854/71‑1 (MQ=255)

                                   gacggTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGa                              <  1:1098075/71‑1 (MQ=255)

                                   gacggTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGa                              <  1:282701/71‑1 (MQ=255)

                                   gacggTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGa                              <  1:771362/71‑1 (MQ=255)

                                   gacggTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGa                              <  1:1806528/71‑1 (MQ=255)

                                   gacggTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGa                              <  1:1801619/71‑1 (MQ=255)

                                   gacggTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGa                              <  1:1727537/71‑1 (MQ=255)

                                   gacggTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGa                              <  1:1615167/71‑1 (MQ=255)

                                   gacggTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGa                              <  1:888818/71‑1 (MQ=255)

                                   gacggTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGa                              <  1:1450760/71‑1 (MQ=255)

                                   gacggTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGa                              <  1:1127319/71‑1 (MQ=255)

                                    acggTAAAATCCGCTACGCTGTAATACAGGCTGtt                                                                 <  1:749703/35‑1 (MQ=255)

                                               cGCTACGCTGTAATACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGaaaa                  <  1:1780096/71‑1 (MQ=255)

                                               cGCTACGCTGTAATACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGaaaa                  <  1:740770/71‑1 (MQ=255)

                                                              cAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCgcagc    >  1:185248/1‑70 (MQ=255)

                                                              cAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCgcagc    >  1:1606905/1‑70 (MQ=255)

                                                               aGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCaa  <  1:532210/71‑1 (MQ=255)

                                                               aGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCaa  <  1:1881520/71‑1 (MQ=255)

                                                               aGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCaa  <  1:1468379/71‑1 (MQ=255)

                                                               aGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCaa  <  1:1168972/71‑1 (MQ=255)

                                                               aGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCaa  <  1:983941/71‑1 (MQ=255)

TCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAA  >  W3110S.gb/158561‑158694

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