BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,872,697	T→C	35.1%	F440F (TTT→TTC)	ygcE →	predicted kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,872,697	0	T	C	35.1%	86.1 / 76.9	128	F440F (TTT→TTC)	ygcE	predicted kinase
Reads supporting (aligned to +/- strand): ref base T (43/40); new base C (32/13); total (75/53)
Fisher's exact test for biased strand distribution p-value = 3.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.02e-01

Only 100 of 133 total aligned reads displayed.

GGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAGACCT  >  minE/1872632‑1872759

ggTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTtata                                                             >  1:908499/1‑69 (MQ=255)

ggTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTtata                                                             >  1:2231693/1‑69 (MQ=255)

ggTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTtata                                                             >  1:1947710/1‑69 (MQ=255)

ggTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTtata                                                             >  1:2515628/1‑69 (MQ=255)

  tCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATAcc                                                           <  1:1783625/69‑1 (MQ=255)

    tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCg                                                          >  1:2233681/1‑68 (MQ=255)

    tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCg                                                          >  1:230725/1‑68 (MQ=255)

    tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa                                                         >  1:1663494/1‑69 (MQ=255)

    tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa                                                         >  1:2553891/1‑69 (MQ=255)

    tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa                                                         >  1:3238365/1‑69 (MQ=255)

    tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa                                                         >  1:3326747/1‑69 (MQ=255)

    tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa                                                         >  1:62038/1‑69 (MQ=255)

    tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa                                                         >  1:3371099/1‑69 (MQ=255)

    tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa                                                         >  1:1335231/1‑69 (MQ=255)

    tGGGGCTATACCCGCATTACACAACGGCTGTCGATAATATGGTTCGCGTGAAAGATATCTTCATACCGa                                                         >  1:2431865/1‑69 (MQ=255)

    tGGGGCGATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa                                                         >  1:3020/1‑69 (MQ=255)

      gggCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGa                                                         >  1:2714913/1‑67 (MQ=255)

              cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc                                                 <  1:3327540/67‑1 (MQ=255)

              cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc                                                 <  1:388596/67‑1 (MQ=255)

              cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc                                                 <  1:2139977/67‑1 (MQ=255)

              cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc                                                 <  1:2761953/67‑1 (MQ=255)

              cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc                                                 <  1:2844417/67‑1 (MQ=255)

              cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc                                                 <  1:726139/67‑1 (MQ=255)

              cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc                                                 <  1:1512264/67‑1 (MQ=255)

              cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc                                                 <  1:1425078/67‑1 (MQ=255)

              cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa                                               >  1:2826826/1‑69 (MQ=255)

              cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa                                               >  1:221201/1‑69 (MQ=255)

              cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa                                               >  1:1675949/1‑69 (MQ=255)

              cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa                                               >  1:1524093/1‑69 (MQ=255)

              cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa                                               >  1:268767/1‑69 (MQ=255)

              cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa                                               >  1:1360916/1‑69 (MQ=255)

              cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa                                               >  1:3381351/1‑69 (MQ=255)

              cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa                                               >  1:1481573/1‑69 (MQ=255)

              cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa                                               >  1:205960/1‑69 (MQ=255)

              cccGGATTACGCAACGGCTGTCGATAACATGGATCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa                                               >  1:2132856/1‑69 (MQ=255)

              ccaggATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa                                               >  1:2277287/4‑69 (MQ=255)

              cacGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc                                                 <  1:931381/65‑1 (MQ=255)

                cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg                                             <  1:200735/69‑1 (MQ=255)

                cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg                                             <  1:349137/69‑1 (MQ=255)

                cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg                                             <  1:1833466/69‑1 (MQ=255)

                cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg                                             <  1:1616353/69‑1 (MQ=255)

                           aCGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTac                                  <  1:2418151/69‑1 (MQ=255)

                           aCGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTac                                  <  1:267360/69‑1 (MQ=255)

                            cGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTac                                  <  1:312903/68‑1 (MQ=255)

                             ggCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTacga                                <  1:2289315/69‑1 (MQ=255)

                             ggCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTacga                                <  1:1059939/69‑1 (MQ=255)

                                 gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCt                                    >  1:1707379/1‑61 (MQ=255)

                                 gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGActcg                            >  1:2638911/1‑69 (MQ=255)

                                 gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg                            >  1:429752/1‑69 (MQ=255)

                                 gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg                            >  1:429925/1‑69 (MQ=255)

                                 gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg                            >  1:469108/1‑69 (MQ=255)

                                 gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg                            >  1:1894122/1‑69 (MQ=255)

                                 gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg                            >  1:243214/1‑69 (MQ=255)

                                 gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg                            >  1:2509417/1‑69 (MQ=255)

                                 gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg                            >  1:857445/1‑69 (MQ=255)

                                 gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg                            >  1:2990893/1‑69 (MQ=255)

                                 gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg                            >  1:1703275/1‑69 (MQ=255)

                                 gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg                            >  1:3016647/1‑69 (MQ=255)

                                 gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg                            >  1:1248946/1‑69 (MQ=255)

                                 gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg                            >  1:632220/1‑69 (MQ=255)

                                 gTCAATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg                            >  1:122519/1‑69 (MQ=255)

                                       aaCATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAAt                      <  1:1151296/69‑1 (MQ=255)

                                       aaCATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAAt                      <  1:763810/69‑1 (MQ=255)

                                       aaCATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAAt                      <  1:1358426/69‑1 (MQ=255)

                                       aaCATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAAt                      <  1:1943788/69‑1 (MQ=255)

                                       aaCATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAAt                      <  1:3074233/69‑1 (MQ=255)

                                             gTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGAcgc                             <  1:1497708/56‑1 (MQ=255)

                                             gTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGAcgc                             <  1:2132518/56‑1 (MQ=255)

                                             gTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGAcgc                             <  1:1654844/56‑1 (MQ=255)

                                             gTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGAcgc                             <  1:1796865/56‑1 (MQ=255)

                                             gTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGAcgc                             <  1:1076624/56‑1 (MQ=255)

                                             gTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGAt                          >  1:637478/1‑59 (MQ=255)

                                             gTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGAt                          >  1:113316/1‑59 (MQ=255)

                                             gTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGAt                          >  1:2033154/1‑59 (MQ=255)

                                             gTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGAt                          >  1:3251220/1‑59 (MQ=255)

                                             gTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGAt                          >  1:2979071/1‑59 (MQ=255)

                                             gTTCGCGTGAAAGATATCTTTATACAGATTGAGAGCAATGCCAAACGCTACGACGCGAt                          >  1:2655439/1‑59 (MQ=255)

                                                        aGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAg      <  1:853531/68‑1 (MQ=255)

                                                        aGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAg      <  1:1315552/68‑1 (MQ=255)

                                                        aGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAg      <  1:2643744/68‑1 (MQ=255)

                                                        aGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAg      <  1:2655556/68‑1 (MQ=255)

                                                        aGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAg      <  1:1779973/68‑1 (MQ=255)

                                                        aGATATCTTCATACCGATTGAGAGCAATCCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAg      <  1:1361286/68‑1 (MQ=255)

                                                           tatCTTTATCCCGATAGAGAGCAATGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAGACCt  >  1:2998009/1‑69 (MQ=255)

                                                           tatCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAGACCt  >  1:512466/1‑69 (MQ=255)

                                                           tatCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAGACCt  >  1:1969878/1‑69 (MQ=255)

                                                           tatCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAGACCt  >  1:1869415/1‑69 (MQ=255)

                                                           tatCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAGACCt  >  1:176356/1‑69 (MQ=255)

                                                           tatCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAGACCt  >  1:581070/1‑69 (MQ=255)

                                                           tatCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAGACCt  >  1:1469878/1‑69 (MQ=255)

                                                           tatCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAGACCt  >  1:1626504/1‑69 (MQ=255)

                                                           tatCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAGACCt  >  1:464023/1‑69 (MQ=255)

                                                           tatCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAGACCt  >  1:2852609/1‑69 (MQ=255)

                                                           tatCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAGACCt  >  1:2861998/1‑69 (MQ=255)

                                                           tatCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAGACCt  >  1:3372197/1‑69 (MQ=255)

                                                           tatCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAGACCt  >  1:1339196/1‑69 (MQ=255)

                                                           tatCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAGACCt  >  1:1605689/1‑69 (MQ=255)

                                                            atCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAATAAAgg                 >  1:864019/1‑53 (MQ=255)

                                                            atCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAATAAAgg                 >  1:1047769/1‑53 (MQ=255)

                                                            atCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGAGATGAATAAAgg                 >  1:969113/1‑53 (MQ=255)

GGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAGACCT  >  minE/1872632‑1872759

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`