BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	422,899	T→C	51.9%	S243S (TCA→TCG)	nagD ←	UMP phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	422,899	0	T	C	51.9%	7.9 / 100.9	115	S243S (TCA→TCG)	nagD	UMP phosphatase
Reads supporting (aligned to +/- strand): ref base T (26/29); new base C (39/20); total (65/50)
Fisher's exact test for biased strand distribution p-value = 5.82e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.78e-01

Only 100 of 115 total aligned reads displayed.

ATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACA  >  minE/422839‑422959

atcaGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTa                                                      <  1:919764/69‑1 (MQ=255)

atcaGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTa                                                      <  1:1999319/69‑1 (MQ=255)

atcaGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTa                                                      <  1:581537/69‑1 (MQ=255)

atcaGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTa                                                      <  1:2080260/69‑1 (MQ=255)

atcaGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTa                                                      <  1:3378532/69‑1 (MQ=255)

atcaGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTa                                                      <  1:2742735/69‑1 (MQ=255)

      gCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCa                                                <  1:2631788/69‑1 (MQ=255)

            cACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTggg                                          >  1:2849684/1‑69 (MQ=255)

                       aCGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCa                               <  1:353470/69‑1 (MQ=255)

                       aCGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCa                               <  1:1951311/69‑1 (MQ=255)

                            gTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTg                          >  1:3102068/1‑69 (MQ=255)

                            gTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTg                          >  1:1292949/1‑69 (MQ=255)

                            gTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTg                          >  1:1447138/1‑69 (MQ=255)

                            gTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTg                          >  1:235521/1‑69 (MQ=255)

                            gTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTg                          >  1:638337/1‑69 (MQ=255)

                              tGCTATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATcc                                                 >  1:2566841/1‑44 (MQ=255)

                              tGCTATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCa                                                >  1:203814/1‑45 (MQ=255)

                              tGCTATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCa                               >  1:3199427/1‑62 (MQ=255)

                              tGCTATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGt                         >  1:94083/1‑68 (MQ=255)

                              tGCTATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGt                         >  1:2477595/1‑68 (MQ=255)

                              tGCTATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTc                        >  1:2247079/1‑69 (MQ=255)

                              tGCTATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTc                        >  1:2893230/1‑69 (MQ=255)

                              tGCTATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTc                        >  1:3240141/1‑69 (MQ=255)

                              tGCTATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTc                        >  1:2580649/1‑69 (MQ=255)

                              tGCTATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTc                        >  1:2274349/1‑69 (MQ=255)

                              tGCTATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTc                        >  1:453198/1‑69 (MQ=255)

                              tGCTATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTc                        >  1:22435/1‑69 (MQ=255)

                              tGCTATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTc                        >  1:2215228/1‑69 (MQ=255)

                              tGCTATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTc                        >  1:2183672/1‑69 (MQ=255)

                              tGCTATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTc                        >  1:2160900/1‑69 (MQ=255)

                              tGCTATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTc                        >  1:708074/1‑69 (MQ=255)

                              tGCTATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTc                        >  1:1682245/1‑69 (MQ=255)

                              tGCTATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTc                        >  1:1939339/1‑69 (MQ=255)

                              tGCTATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTc                        >  1:919863/1‑69 (MQ=255)

                                 taatCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGat                     >  1:933782/4‑69 (MQ=255)

                                 tATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACttg                                           >  1:2348423/1‑45 (MQ=255)

                                 tATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGat                     >  1:2316592/1‑69 (MQ=255)

                                 tATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGat                     >  1:1020460/1‑69 (MQ=255)

                                 tATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGat                     >  1:2461212/1‑69 (MQ=255)

                                 tATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGat                     >  1:409871/1‑69 (MQ=255)

                                 tATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGat                     >  1:166873/1‑69 (MQ=255)

                                 tATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGat                     >  1:1895533/1‑69 (MQ=255)

                                 tATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGat                     >  1:752830/1‑69 (MQ=255)

                                 tATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGat                     >  1:479102/1‑69 (MQ=255)

                                 tATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGat                     >  1:1046303/1‑69 (MQ=255)

                                 tATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGat                     >  1:1029898/1‑69 (MQ=255)

                                 tATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGat                     >  1:910958/1‑69 (MQ=255)

                                 tATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGa                      >  1:1792632/1‑68 (MQ=255)

                                 tATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGa                      >  1:1179572/1‑68 (MQ=255)

                                 tATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGa                      >  1:2342581/1‑68 (MQ=255)

                                 tATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGa                      >  1:3233329/1‑68 (MQ=255)

                                 tATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTAGGCGGAACGGCATACTGTCGat                     >  1:2964281/1‑69 (MQ=255)

                                     cAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATAtcg                 >  1:1693121/1‑69 (MQ=255)

                                     cAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATAtcg                 >  1:940484/1‑69 (MQ=255)

                                     cAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATAtcg                 >  1:1097380/1‑69 (MQ=255)

                                     cAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATAtcg                 >  1:2437889/1‑69 (MQ=255)

                                     cAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATAtcg                 >  1:559150/1‑69 (MQ=255)

                                     cAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATAtcg                 >  1:2311015/1‑69 (MQ=255)

                                     cAGATAACGTCGATTTCAGCGACTGACGGGAAAATACAGCTgg                                           >  1:2637350/1‑43 (MQ=37)

                                                gATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGATATCGTCGAGtgacga      <  1:343366/69‑1 (MQ=255)

                                                gATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGATATCGTCGAGtgacga      <  1:1026522/69‑1 (MQ=255)

                                                gATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGATATCGTCGAGtgacga      <  1:420301/69‑1 (MQ=255)

                                                gATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGATATCGTCGAGtgacga      <  1:2933608/69‑1 (MQ=255)

                                                gATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGATATCGTCGAGtgacga      <  1:1244821/69‑1 (MQ=255)

                                                gATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGATATCGTCGAGtgacga      <  1:1247779/69‑1 (MQ=255)

                                                gATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGATATCGTCGAGtgacga      <  1:2271787/69‑1 (MQ=255)

                                                gATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGATATCGTCGAGtgacga      <  1:2142286/69‑1 (MQ=255)

                                                gATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGATATCGTCGAGtgacga      <  1:2035679/69‑1 (MQ=255)

                                                gATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGATATCGTCGAGtgacga      <  1:2682260/69‑1 (MQ=255)

                                                gATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGATATCGTCGAGtgacga      <  1:1725857/69‑1 (MQ=255)

                                                gATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGATATCGTCGAGtgacga      <  1:1885375/69‑1 (MQ=255)

                                                gATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGATATCGTCGAGtgacga      <  1:182216/69‑1 (MQ=255)

                                                gATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGATATCGTCGAGtgacga      <  1:1242390/69‑1 (MQ=255)

                                                gATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGATATCGTCGAGtgacga      <  1:179569/69‑1 (MQ=255)

                                                gATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGATATCGTCGAGtcacga      <  1:792820/69‑1 (MQ=255)

                                                gATTTCAGCGACAGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGATATCGTCGAGtgacga      <  1:2013887/69‑1 (MQ=255)

                                                 aTTTCAGCGACTGACGTGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGcgacga      <  1:409220/68‑1 (MQ=255)

                                                 aTTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGcgacga      <  1:495190/68‑1 (MQ=255)

                                                 aTTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGcgacga      <  1:831857/68‑1 (MQ=255)

                                                 aTTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGcgacga      <  1:2173910/68‑1 (MQ=255)

                                                 aTTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGcgacga      <  1:2205991/68‑1 (MQ=255)

                                                 aTTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGcgacga      <  1:2892907/68‑1 (MQ=255)

                                                 aTTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGcgacga      <  1:1222609/68‑1 (MQ=255)

                                                 aTTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGcgacga      <  1:3161703/68‑1 (MQ=255)

                                                 aTTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGcgacga      <  1:1618867/68‑1 (MQ=255)

                                                 aTTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGcgacga      <  1:2642336/68‑1 (MQ=255)

                                                 aTTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGcgacga      <  1:2308887/68‑1 (MQ=255)

                                                 aTTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGATATCGTCGAGtgacga      <  1:1255585/68‑1 (MQ=255)

                                                   ttCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGcga         >  1:107039/1‑63 (MQ=255)

                                                   ttCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGaaa    >  1:215031/1‑68 (MQ=255)

                                                   ttCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAac   >  1:783118/1‑69 (MQ=255)

                                                   ttCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAac   >  1:3152988/1‑69 (MQ=255)

                                                   ttCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAac   >  1:1336560/1‑69 (MQ=255)

                                                   ttCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAac   >  1:743327/1‑69 (MQ=255)

                                                   ttCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAac   >  1:584370/1‑69 (MQ=255)

                                                   ttCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAac   >  1:2198169/1‑69 (MQ=255)

                                                    tCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAaca  <  1:834601/69‑1 (MQ=255)

                                                    tCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAaca  <  1:3300609/69‑1 (MQ=255)

                                                    tCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAaca  <  1:2202372/69‑1 (MQ=255)

                                                    tCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAaca  <  1:2909090/69‑1 (MQ=255)

ATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACA  >  minE/422839‑422959

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