BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,906,931	A→G	46.4%	V464V (GTT→GTC)	pyrG ←	CTP synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,906,931	0	A	G	46.4%	7.7 / 99.8	69	V464V (GTT→GTC)	pyrG	CTP synthetase
Reads supporting (aligned to +/- strand): ref base A (11/26); new base G (10/22); total (21/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

CAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGGCTATCGTCAACCAACTGGC  >  W3110S.gb/2906865‑2906986

cAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTcaac                                                       <  1:1215431/69‑1 (MQ=255)

cAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTcaac                                                       <  1:2487204/69‑1 (MQ=255)

cAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTcaac                                                       <  1:2647252/69‑1 (MQ=255)

    ccTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACTATGACGCTCGACAAtt                                                   >  1:550081/1‑69 (MQ=255)

       gCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTcaac                                                       <  1:1222967/62‑1 (MQ=255)

       gCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTcaac                                                       <  1:1478780/62‑1 (MQ=255)

       gCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAAAGGTGACGATGACGCTCAACAATTGTc                                                <  1:3346040/69‑1 (MQ=255)

        cATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCGACAATTGTCg                                               >  1:658018/1‑69 (MQ=255)

        cATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCGACAATTGTCg                                               >  1:568959/1‑69 (MQ=255)

        cATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCGACAATTGTCg                                               >  1:3728980/1‑69 (MQ=255)

            ttCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCGACAATTGTCGgcgc                                           >  1:2405572/1‑69 (MQ=255)

            ttCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAAtt                                                   >  1:543367/1‑61 (MQ=255)

              cAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCGACAATTGTCGGCGCAt                                         >  1:3090298/1‑69 (MQ=255)

              cAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCGACAATTGTCGGCGCAt                                         >  1:2680137/1‑69 (MQ=255)

              cAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCAt                                         >  1:17006/1‑69 (MQ=255)

              cAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCAt                                         >  1:1572124/1‑69 (MQ=255)

                 tCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGCCGATGACGCTCAACAATTGTCGGCGCATTg                                       >  1:1794629/1‑68 (MQ=255)

                 tCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCGACAATTGTCGGCGCATTg                                       <  1:2669027/68‑1 (MQ=255)

                 tCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCGACAATTGTCGGCGCATTg                                       <  1:1770018/68‑1 (MQ=255)

                 tCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCGACAATTGTCGGCGCATTg                                       <  1:2664846/68‑1 (MQ=255)

                 tCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCGACAATTGTCGGCGCATTg                                       <  1:1974768/68‑1 (MQ=255)

                 tCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCGACAATTGTCGGCGCATTg                                       <  1:202245/68‑1 (MQ=255)

                 tCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCGACAATTGTCGGCGCATTg                                       <  1:2536395/68‑1 (MQ=255)

                 tCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTg                                       >  1:1798907/1‑68 (MQ=255)

                 tCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTg                                       >  1:3021133/1‑68 (MQ=255)

                       tCAAAAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCt                                >  1:2310828/1‑69 (MQ=255)

                           cAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCg                            <  1:2621507/69‑1 (MQ=255)

                           cAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCg                            <  1:201223/69‑1 (MQ=255)

                           cAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCg                            <  1:3683997/69‑1 (MQ=255)

                              cATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCa                                          >  1:104559/1‑52 (MQ=255)

                              cATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTg                                       <  1:3339306/55‑1 (MQ=255)

                              cATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTg                                       <  1:3195041/55‑1 (MQ=255)

                              cATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTg                                       <  1:2750130/55‑1 (MQ=255)

                              cATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTg                                       <  1:1158803/55‑1 (MQ=255)

                              cATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCGAAc                         <  1:2313936/69‑1 (MQ=255)

                               aTGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCGAAcc                        <  1:1884411/69‑1 (MQ=255)

                                tgttgtTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCg                            <  1:1477852/64‑1 (MQ=255)

                                  ttgttgACTTCGTAACGGTGACGATGACGCTCGACAATTGTCGGCGCATTg                                       <  1:1752037/51‑1 (MQ=255)

                                  ttgttgACTTCGTAACGGTGACGATGACGCTCGACAATTGTCGGCGCATTg                                       <  1:54256/51‑1 (MQ=255)

                                  ttgttgACTTCGTAACGGTGACGATGACGCTCGACAATTGTCGGCGCATTg                                       <  1:1578265/51‑1 (MQ=255)

                                  ttgttgACTTCGTAACGGTGACGATGACGCTCGACAATTGTCGGCGCATTg                                       <  1:3656284/51‑1 (MQ=255)

                                  ttgttgACTTCGTAACGGTGACGATGACGCTCGACAATTGTCGGCGCATTg                                       <  1:1260157/51‑1 (MQ=255)

                                  ttgttgACTTCGTAACGGTGACGATGACGCTCGACAATTGTCGGCGCATTg                                       <  1:2343887/51‑1 (MQ=255)

                                   tgttgACTTCGTAACGGTGACGATGACGCTCGACAATTGTCGGCGCATTg                                       <  1:1064808/50‑1 (MQ=255)

                                   tgttgACTTCGTAACGGTGACGATGACGCTCGACAATTGTCGGCGCATTg                                       <  1:3339532/50‑1 (MQ=255)

                                   tgttgACTTCGTAACGGTGACGATGACGCTCGACAATTGTCGGCGCATTg                                       <  1:2252189/50‑1 (MQ=255)

                                    gttgACTTCGTAACGGTGACGATGACGCTCGACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGGCt                   >  1:2066177/1‑69 (MQ=255)

                                      tgACTTCGTAACGGTGACGATGACGCTCAACAAtt                                                   <  1:2040076/35‑1 (MQ=255)

                                      tgACTTCGTAACGGTGACGATGACGCTCAACAAtt                                                   <  1:3474982/35‑1 (MQ=255)

                                      tgACTTCGTAACGGTGACGATGACGCTCAACAAtt                                                   <  1:3287464/35‑1 (MQ=255)

                                      tgACTTCGTAACGGTGACGATGACGCTCAACAAtt                                                   <  1:3067957/35‑1 (MQ=255)

                                      tgACTTCGTAACGGTGACGATGACGCTCAACAAtt                                                   >  1:2893175/1‑35 (MQ=255)

                                      tgACTTCGTAACGGTGACGATGACGCTCAACAAtt                                                   <  1:37589/35‑1 (MQ=255)

                                      tgACTTCGTAACGGTGACGATGACGCTCAACAAtt                                                   <  1:677311/35‑1 (MQ=255)

                                      tgACTTCGTAACGGTGACGATGACGCTCAACAAtt                                                   <  1:1938230/35‑1 (MQ=255)

                                          ttCGTAACGGTGACGATGACGCTCAACAATTGTCTGCGCATTGTACAGCTGGCGAACCAGGCTATCGTc             <  1:79684/69‑1 (MQ=255)

                                           tCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGGCTATCGTca            >  1:663170/1‑69 (MQ=255)

                                           tCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCGAACCACGCTATCGTc             >  1:1039172/1‑68 (MQ=255)

                                               aaCGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGGCTATCGTcaacc         <  1:1883982/68‑1 (MQ=255)

                                               aaCGGTGACGATGACGCTAGACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGGCTATCGTcaacca        <  1:2076269/69‑1 (MQ=255)

                                                aCGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGGCTATCGTcaacca        <  1:2677085/68‑1 (MQ=255)

                                                    tGACGATGACGCTCGACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGGCTATCGTCAACCAACTgg   <  1:3225986/69‑1 (MQ=255)

                                                    tGACGATGACGCTCGACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGGCTATCGTCAACCAACTgg   <  1:2867192/69‑1 (MQ=255)

                                                    tGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGGCTATCGTCAACCAACTgg   <  1:39067/69‑1 (MQ=255)

                                                     gaggATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGGCTATCGTCAACCAACTGGc  <  1:1508185/66‑1 (MQ=255)

                                                     gACGATGACGCTCGACATTTGTCGGCGCATTGTACAGCTGGCGAACCAGGCTATCGTCAACGAACTgg   <  1:479103/68‑1 (MQ=255)

                                                     gACGATGACGCTCGACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGGCTATCGTCAACCAACTgg   <  1:1395600/68‑1 (MQ=255)

                                                     gACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGGCTATCGTCAACCAACTGGc  <  1:423836/69‑1 (MQ=255)

                                                     gACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGGCTATCGTCAACCAACTGGc  <  1:1061089/69‑1 (MQ=255)

                                                        gATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGGCTATCGTCAACCAACTgg   <  1:2126421/65‑1 (MQ=255)

                                                             cGCTCGACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGGCTATCGTCAACCAACTgg   >  1:1610611/1‑60 (MQ=255)

                                                             cGCTCGACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGGCTATCGTCAACCAACTgg   >  1:567332/1‑60 (MQ=255)

                                                             cGCTCGACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGGCTATCGTCAACCAACTgg   <  1:3342341/60‑1 (MQ=255)

                                                             cGCTCGACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGGCTATCGTCAACCAACTgg   <  1:3591259/60‑1 (MQ=255)

CAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGGCTATCGTCAACCAACTGGC  >  W3110S.gb/2906865‑2906986

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