BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	423,289	G→A	55.1%	F113F (TTC→TTT)	nagD ←	UMP phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	423,289	0	G	A	55.1%	22.9 / 104.9	126	F113F (TTC→TTT)	nagD	UMP phosphatase
Reads supporting (aligned to +/- strand): ref base G (18/38); new base A (31/38); total (50/76)
Fisher's exact test for biased strand distribution p-value = 1.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

Only 100 of 126 total aligned reads displayed.

CGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATGAATCAGTGCGCCTTC  >  minE/423228‑423354

cgaAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAg                                                            <  1:3215182/69‑1 (MQ=255)

cgaAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAg                                                            <  1:1494237/69‑1 (MQ=255)

cgaAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAg                                                            <  1:1749224/69‑1 (MQ=255)

cgaAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAg                                                            <  1:406439/69‑1 (MQ=255)

cgaAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAg                                                            <  1:814027/69‑1 (MQ=255)

       ggCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTcac                                                     <  1:382251/69‑1 (MQ=255)

       ggCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTcac                                                     <  1:2344567/69‑1 (MQ=255)

       ggCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTcac                                                     <  1:1560165/69‑1 (MQ=255)

       ggCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTcac                                                     <  1:1517442/69‑1 (MQ=255)

       ggCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTcac                                                     <  1:347774/69‑1 (MQ=255)

       ggCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTcac                                                     <  1:2375629/69‑1 (MQ=255)

       ggCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTcac                                                     <  1:2168091/69‑1 (MQ=255)

       ggCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTcac                                                     <  1:1638748/69‑1 (MQ=255)

       ggCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTcac                                                     <  1:2691409/69‑1 (MQ=255)

       ggCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTcac                                                     <  1:1833212/69‑1 (MQ=255)

        gctgctTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTcac                                                     <  1:1406270/68‑1 (MQ=255)

        gctgctTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTcac                                                     <  1:2944378/68‑1 (MQ=255)

                       atGTCCCAGTTTTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACAt                                                   <  1:3120795/55‑1 (MQ=39)

                       atGTCCCAGTTGTAGGAACTCGTTTCGCCAACAATCACAAAATCAGGGTTCACAt                                                   <  1:1980264/55‑1 (MQ=39)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGaaa                                     >  1:1868548/1‑69 (MQ=255)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGaaa                                     >  1:1767121/1‑69 (MQ=255)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACATAATCAGGGTTCACAt                                                   <  1:137026/55‑1 (MQ=255)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACAt                                                   <  1:3423538/55‑1 (MQ=255)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACAt                                                   <  1:2406193/55‑1 (MQ=255)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACAt                                                   <  1:357472/55‑1 (MQ=255)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACAt                                                   <  1:3363122/55‑1 (MQ=255)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACAt                                                   <  1:3318782/55‑1 (MQ=255)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACAt                                                   <  1:2453458/55‑1 (MQ=255)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACAt                                                   <  1:2482380/55‑1 (MQ=255)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACAt                                                   <  1:1929326/55‑1 (MQ=255)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACAt                                                   <  1:2491652/55‑1 (MQ=255)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACAt                                                   <  1:2763570/55‑1 (MQ=255)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACAt                                                   <  1:2769702/55‑1 (MQ=255)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACAt                                                   <  1:175754/55‑1 (MQ=255)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACAt                                                   <  1:1218406/55‑1 (MQ=255)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACAt                                                   <  1:877186/55‑1 (MQ=255)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACAt                                                   <  1:1153612/55‑1 (MQ=255)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACAt                                                   <  1:1250021/55‑1 (MQ=255)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACAt                                                   <  1:1867163/55‑1 (MQ=255)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACAt                                                   <  1:862988/55‑1 (MQ=255)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACAt                                                   <  1:1411891/55‑1 (MQ=255)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACAt                                                   <  1:1507905/55‑1 (MQ=255)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACAt                                                   <  1:794173/55‑1 (MQ=255)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACAt                                                   <  1:648286/55‑1 (MQ=255)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACAt                                                   <  1:1648442/55‑1 (MQ=255)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACAt                                                   <  1:1656409/55‑1 (MQ=255)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACAt                                                   <  1:580784/55‑1 (MQ=255)

                       atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACAt                                                   <  1:481888/55‑1 (MQ=255)

                         gTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAAcc                                   >  1:212491/1‑69 (MQ=255)

                          tCCCAGTTGTAGGACCGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCg                                  >  1:3187824/1‑69 (MQ=255)

                          tCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACATc                                                  >  1:2660256/1‑53 (MQ=255)

                          tCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACATc                                                  >  1:2623093/1‑53 (MQ=255)

                          tCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACATc                                                  >  1:1027706/1‑53 (MQ=255)

                          tCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACATc                                                  >  1:2968247/1‑53 (MQ=255)

                          tCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACATc                                                  >  1:2589538/1‑53 (MQ=255)

                          tCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACATc                                                  >  1:3268402/1‑53 (MQ=255)

                          tCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACATc                                                  >  1:3305760/1‑53 (MQ=255)

                          tCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACATc                                                  >  1:1001752/1‑53 (MQ=255)

                          tCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACATc                                                  >  1:2198354/1‑53 (MQ=255)

                          tCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACATc                                                  >  1:195714/1‑53 (MQ=255)

                             cAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTc                                                       >  1:1595542/1‑45 (MQ=255)

                             cAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACATCgg                                                >  1:125181/1‑52 (MQ=255)

                             cAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACATCGGt                                               >  1:95062/1‑53 (MQ=255)

                             cAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACATCGGt                                               >  1:924110/1‑53 (MQ=255)

                             cAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACATCGGt                                               >  1:804189/1‑53 (MQ=255)

                             cAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACATCGGt                                               >  1:981222/1‑53 (MQ=255)

                             cAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACATCGGt                                               >  1:1278419/1‑53 (MQ=255)

                             cAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACATCGGt                                               >  1:1740335/1‑53 (MQ=255)

                             cAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACATCGGt                                               >  1:1574987/1‑53 (MQ=255)

                             cAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACATCGGt                                               >  1:2357326/1‑53 (MQ=255)

                             cAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACATCGGt                                               >  1:2418908/1‑53 (MQ=255)

                             cAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACATCGGt                                               >  1:1573429/1‑53 (MQ=255)

                             cAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACATCGGt                                               >  1:2571081/1‑53 (MQ=255)

                             cAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACATCGGt                                               >  1:1221023/1‑53 (MQ=255)

                             cAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACATCGGt                                               >  1:2724142/1‑53 (MQ=255)

                             cAGTTGTAGGAACGCGTTTCGCCAACAATCACAAAATCAGGGTTCACATCGGt                                               >  1:2823663/1‑53 (MQ=255)

                                  gTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTaa                                             <  1:1033862/50‑1 (MQ=255)

                                  gTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTaa                                             <  1:725887/50‑1 (MQ=255)

                                  gTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTaa                                             <  1:196618/50‑1 (MQ=255)

                                  gTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAAc                                    >  1:3322642/1‑59 (MQ=255)

                                  gTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAAc                                    >  1:1683188/1‑59 (MQ=255)

                                  gTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTa                          <  1:2615726/69‑1 (MQ=255)

                                  gTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTa                          <  1:770411/69‑1 (MQ=255)

                                  gTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTa                          <  1:1005438/69‑1 (MQ=255)

                                       aaCGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAg                                          >  1:2687047/1‑48 (MQ=255)

                                       aaCGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTg                            >  1:1267258/1‑62 (MQ=255)

                                       aaCGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTg                            >  1:2849156/1‑62 (MQ=255)

                                       aaCGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTg                            >  1:2240478/1‑62 (MQ=255)

                                       aaCGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTg                            >  1:2480083/1‑62 (MQ=255)

                                       aaCGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTg                            >  1:3123959/1‑62 (MQ=255)

                                                       aaTCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATGAATCAGTGCGcc     <  1:20808/69‑1 (MQ=255)

                                                       aaTCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATGAATCAGTGCGcc     <  1:2406895/69‑1 (MQ=255)

                                                       aaTCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATGAATCAGTGCGcc     <  1:1338629/69‑1 (MQ=255)

                                                       aaTCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATGAATCAGTGCGcc     <  1:3290261/69‑1 (MQ=255)

                                                       aaTCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATGAATCAGTGCGcc     <  1:3260914/69‑1 (MQ=255)

                                                       aaTCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATGAATCAGTGCGcc     <  1:2666238/69‑1 (MQ=255)

                                                       aaTCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATGAATCAGTGCGcc     <  1:1115539/69‑1 (MQ=255)

                                                       aaTCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATGAATCAGTGCGcc     <  1:2714674/69‑1 (MQ=255)

                                                          cACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATGAATCAGTGCGCCTTc  >  1:515512/1‑69 (MQ=255)

                                                          cACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATGAATCAGTGCGCCTTc  >  1:1994416/1‑69 (MQ=255)

CGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATGAATCAGTGCGCCTTC  >  minE/423228‑423354

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