BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	422,926	A→C	46.1%	P234P (CCT→CCG)	nagD ←	UMP phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	422,926	0	A	C	46.1%	0.6 / 154.2	128	P234P (CCT→CCG)	nagD	UMP phosphatase
Reads supporting (aligned to +/- strand): ref base A (39/30); new base C (36/23); total (75/53)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

Only 100 of 128 total aligned reads displayed.

ACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAG  >  minE/422862‑422985

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

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

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

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

     gTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTg                                                    >  1:1292949/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)

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

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

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

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

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

       tGCTATTCAGATAACGTCGATTTCAGCGACCGATGGGTAAATCCAACTTGGGCGGAACGGCATACTGTc                                                  >  1:2893230/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)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                          aTTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGcgacga                                <  1:2205991/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)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                       ggggCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCa   <  1:1398767/68‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                           cGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAaca                            <  1:3132230/39‑1 (MQ=255)

ACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAG  >  minE/422862‑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: `‑`