BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	647,932	T→C	37.6%	R204R (CGT→CGC)	pepN →	aminopeptidase N

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	647,932	0	T	C	37.6%	34.6 / 43.2	69	R204R (CGT→CGC)	pepN	aminopeptidase N
Reads supporting (aligned to +/- strand): ref base T (16/27); new base C (13/13); total (29/40)
Fisher's exact test for biased strand distribution p-value = 3.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.16e-01

TGTTTGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCTT  >  minE/647877‑647980

tGTTTGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGa                                       >  1:1225786/1‑67 (MQ=255)

tGTTTGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGa                                       >  1:1310765/1‑67 (MQ=255)

tGTTTGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGa                                       >  1:289107/1‑67 (MQ=255)

    tGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAg                                  >  1:2184568/1‑68 (MQ=255)

    tGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGc                                 <  1:2900456/69‑1 (MQ=255)

    tGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGc                                 <  1:3322482/69‑1 (MQ=255)

    tGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGc                                 <  1:2151529/69‑1 (MQ=255)

    tGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGc                                 <  1:3264087/69‑1 (MQ=255)

    tGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGc                                 <  1:2392034/69‑1 (MQ=255)

    tGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGc                                 <  1:1737607/69‑1 (MQ=255)

    tGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGc                                 >  1:606134/1‑69 (MQ=255)

    tGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGc                                 <  1:1555614/69‑1 (MQ=255)

    tGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGc                                 <  1:922613/69‑1 (MQ=255)

    tGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGc                                 <  1:1145068/69‑1 (MQ=255)

    tGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGc                                 <  1:1149406/69‑1 (MQ=255)

    tGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGc                                 <  1:2850298/69‑1 (MQ=255)

    tGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGc                                 <  1:1522304/69‑1 (MQ=255)

        ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                              <  1:597660/67‑1 (MQ=255)

        ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                              <  1:2864050/67‑1 (MQ=255)

        ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                              <  1:3386729/67‑1 (MQ=255)

        ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                              <  1:2906822/67‑1 (MQ=255)

        ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                              <  1:1861711/67‑1 (MQ=255)

        ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                              <  1:1761963/67‑1 (MQ=255)

        ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                              <  1:1297621/67‑1 (MQ=255)

        ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                              <  1:701460/67‑1 (MQ=255)

        ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                              <  1:742985/67‑1 (MQ=255)

        ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                              <  1:1537979/67‑1 (MQ=255)

        ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTAACCGCTCTGGCCGTGAAGTGGCGCt                              <  1:275477/67‑1 (MQ=255)

        gtggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                              <  1:342843/67‑1 (MQ=255)

          ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTgg                            >  1:2622762/1‑68 (MQ=255)

          ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                           >  1:258749/1‑69 (MQ=255)

          ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                           >  1:2484557/1‑69 (MQ=255)

          ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                           >  1:451210/1‑69 (MQ=255)

          ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                           >  1:2030282/1‑69 (MQ=255)

          ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                           >  1:2010221/1‑69 (MQ=255)

          ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                           >  1:343062/1‑69 (MQ=255)

          ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                           >  1:1650299/1‑69 (MQ=255)

          ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                           >  1:1467933/1‑69 (MQ=255)

          ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                           >  1:1248359/1‑69 (MQ=255)

          ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                           >  1:1102925/1‑69 (MQ=255)

                     attttGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                              <  1:756278/53‑1 (MQ=255)

                                  gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc    >  1:3330291/1‑68 (MQ=255)

                                  gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc    >  1:2915786/1‑68 (MQ=255)

                                  gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc    >  1:448900/1‑68 (MQ=255)

                                  gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc    >  1:2323823/1‑68 (MQ=255)

                                  gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc    >  1:707418/1‑68 (MQ=255)

                                  gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc    >  1:243221/1‑68 (MQ=255)

                                  gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc    >  1:1549525/1‑68 (MQ=255)

                                  gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc    >  1:1944714/1‑68 (MQ=255)

                                  gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc    >  1:1668907/1‑68 (MQ=255)

                                  gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc    >  1:1869964/1‑68 (MQ=255)

                                  gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc    >  1:1712456/1‑68 (MQ=255)

                                  gcgcGATACCTTCACTACCCGCTCTGGCCGTGAACTGGCGCTGGAGTTGTATGTCGATCGTGGCaa      >  1:2596225/1‑66 (MQ=255)

                                  gcgcGATACCTTCACTACCCGCTATGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc    >  1:2069169/1‑68 (MQ=255)

                                         aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCTCGGCAACCtt  <  1:236491/63‑1 (MQ=255)

                                         aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt  <  1:986495/63‑1 (MQ=255)

                                         aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt  <  1:2703997/63‑1 (MQ=255)

                                         aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt  <  1:1132963/63‑1 (MQ=255)

                                         aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt  <  1:1383609/63‑1 (MQ=255)

                                         aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt  <  1:1470369/63‑1 (MQ=255)

                                         aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt  <  1:612270/63‑1 (MQ=255)

                                         aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt  <  1:1867879/63‑1 (MQ=255)

                                         aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt  <  1:1898983/63‑1 (MQ=255)

                                         aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt  <  1:391927/63‑1 (MQ=255)

                                         aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt  <  1:3387716/63‑1 (MQ=255)

                                         aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt  <  1:2214016/63‑1 (MQ=255)

                                         aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt  <  1:2219895/63‑1 (MQ=255)

                                         aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt  <  1:1083147/63‑1 (MQ=255)

                                         aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGGCGATCGCGGCAACCtt  <  1:2221256/63‑1 (MQ=255)

TGTTTGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCTT  >  minE/647877‑647980

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