BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	422,914	G→A	52.2%	S238S (AGC→AGT)	nagD ←	UMP phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	422,914	0	G	A	52.2%	2.8 / 104.2	115	S238S (AGC→AGT)	nagD	UMP phosphatase
Reads supporting (aligned to +/- strand): ref base G (34/21); new base A (37/23); total (71/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

Only 100 of 115 total aligned reads displayed.

CACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGT  >  minE/422851‑422979

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:1447138/1‑69 (MQ=255)

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

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

                gTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTg                                              >  1:235521/1‑69 (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:1939339/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:2580649/1‑69 (MQ=255)

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

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

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

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

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

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

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

                  tGCTATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTc                                            >  1:708074/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:752830/1‑69 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                    gATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTCGATATCGTCGAGtgacga                          <  1:2271787/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:1618867/68‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                        tCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAaca                      <  1:1790179/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:678412/1‑60 (MQ=255)

                                                       tAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGaa                >  1:2901715/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:286012/1‑69 (MQ=255)

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

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

CACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGT  >  minE/422851‑422979

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