BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	422,917	G→T	51.7%	P237P (CCC→CCA)	nagD ←	UMP phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	422,917	0	G	T	51.7%	0.3 / 143.4	114	P237P (CCC→CCA)	nagD	UMP phosphatase
Reads supporting (aligned to +/- strand): ref base G (32/23); new base T (36/23); total (68/46)
Fisher's exact test for biased strand distribution p-value = 8.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 129 total aligned reads displayed.

CACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAG  >  minE/422851‑422985

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

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

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

                gTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAg                                                            >  1:1368319/1‑61 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                     ggTAAATCCAACTTGGGCGGAACGGCATACTGTCGATATCGTCGAGtgacga                                <  1:3268478/52‑1 (MQ=37)

                                                      gTAAATCCAACTTGGGCGGAACGGCATACTGTCGATATCGTCGAGtgacga                                <  1:2276037/51‑1 (MQ=37)

                                                       tAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAg                        >  1:1694501/1‑58 (MQ=255)

                                                       tAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGaa                      >  1:2901715/1‑60 (MQ=255)

                                                       tAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGaa                      >  1:678412/1‑60 (MQ=255)

                                                            ccAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGaaa                     >  1:1491602/1‑56 (MQ=255)

                                                            ccAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGt        >  1:1498721/1‑69 (MQ=255)

                                                            ccAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGt        >  1:1541869/1‑69 (MQ=255)

                                                            ccAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGt        >  1:2808265/1‑69 (MQ=255)

                                                             cAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAacac                           <  1:316204/49‑1 (MQ=255)

                                                                 tGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCa   <  1:3116854/69‑1 (MQ=255)

                                                                 tGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCa   <  1:2155418/69‑1 (MQ=255)

                                                                  ggggCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAg  <  1:2489885/69‑1 (MQ=255)

                                                                  ggggCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAg  >  1:2627625/1‑69 (MQ=255)

                                                                  ggggCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAg  >  1:527968/1‑69 (MQ=255)

                                                                  ggggCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAg  <  1:1655132/69‑1 (MQ=255)

                                                                  ggggCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAg  >  1:2803609/1‑69 (MQ=255)

                                                                  ggggCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAg  >  1:370946/1‑69 (MQ=255)

                                                                  ggggCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAg  >  1:899794/1‑69 (MQ=255)

                                                                  ggggCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAg  <  1:3394878/69‑1 (MQ=255)

                                                                  ggggCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAg  >  1:180756/1‑69 (MQ=255)

                                                                  ggggCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAg  >  1:1012896/1‑69 (MQ=255)

CACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAG  >  minE/422851‑422985

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