BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,505,865	C→A	8.6%	G107G (GGC→GGA)	yicJ →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,505,865	0	C	A	8.6%	269.6 / 9.5	117	G107G (GGC→GGA)	yicJ	predicted transporter
Reads supporting (aligned to +/- strand): ref base C (54/52); new base A (3/7); total (57/60)
Fisher's exact test for biased strand distribution p-value = 3.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 117 total aligned reads displayed.

GCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCAATATC  >  minE/2505803‑2505928

gcACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATg                                                           <  1:2061248/69‑1 (MQ=255)

gcACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATg                                                           <  1:1855641/69‑1 (MQ=255)

gcACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATg                                                           <  1:1585036/69‑1 (MQ=255)

   cTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGCCGGATCTCAGTATGAACGGCAAAATGATc                                                        <  1:608431/69‑1 (MQ=255)

   cTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATc                                                        <  1:2446071/69‑1 (MQ=255)

   cTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATc                                                        <  1:2521193/69‑1 (MQ=255)

   cTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATc                                                        <  1:2525123/69‑1 (MQ=255)

   cTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATc                                                        <  1:6494/69‑1 (MQ=255)

   cTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATc                                                        <  1:1758167/69‑1 (MQ=255)

   cTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATc                                                        >  1:3354679/1‑69 (MQ=255)

    tGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCt                                                       <  1:2792443/69‑1 (MQ=255)

     gCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTa                                                      >  1:1000073/1‑69 (MQ=255)

     gCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTa                                                      >  1:1864074/1‑69 (MQ=255)

     gCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTa                                                      >  1:1573311/1‑69 (MQ=255)

     gCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTa                                                      >  1:3398244/1‑69 (MQ=255)

     gCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTa                                                      >  1:1046318/1‑69 (MQ=255)

     gCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTa                                                      >  1:1428270/1‑69 (MQ=255)

     gCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTa                                                      >  1:451319/1‑69 (MQ=255)

           cGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATgcagc                                                >  1:2379145/1‑69 (MQ=255)

           cGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATgcagc                                                >  1:1983094/1‑69 (MQ=255)

           cGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATgcagc                                                >  1:703103/1‑69 (MQ=255)

           cGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATgcagc                                                >  1:2819015/1‑69 (MQ=255)

           cGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATgcagc                                                >  1:435313/1‑69 (MQ=255)

           cGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATgcagc                                                >  1:1685434/1‑69 (MQ=255)

           cGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATgcagc                                                >  1:1689628/1‑69 (MQ=255)

           cGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATgcagc                                                >  1:1690934/1‑69 (MQ=255)

           cGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATgcagc                                                >  1:2946534/1‑69 (MQ=255)

              gATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATg                                                           <  1:150125/55‑1 (MQ=255)

              gATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGAAAAATGATTTATGCAGCAAt                                             >  1:1484798/1‑69 (MQ=255)

              gATCGTCTGTGTACTGGCCTATAGCACGCCAGATCCCAGTATGAACGGAAAAATGATTTATGCAGCAAt                                             >  1:2447248/1‑69 (MQ=255)

                    cTGTTTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAAttactta                                       <  1:1137668/69‑1 (MQ=255)

                    cTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAAttactta                                       <  1:700442/69‑1 (MQ=255)

                    cTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAAttactta                                       <  1:2135895/69‑1 (MQ=255)

                    cTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAAttactta                                       <  1:715524/69‑1 (MQ=255)

                    cTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAAttactta                                       <  1:1534132/69‑1 (MQ=255)

                    cTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAAttactta                                       <  1:2274384/69‑1 (MQ=255)

                    cTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAAttactta                                       <  1:1410994/69‑1 (MQ=255)

                    cTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAAttactta                                       <  1:2279137/69‑1 (MQ=255)

                    cTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAAttactta                                       <  1:3049881/69‑1 (MQ=255)

                     tgtgTACTGGCCTATATCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAAttactta                                       <  1:2913430/68‑1 (MQ=255)

                       tgtACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATg                                                           <  1:2895067/46‑1 (MQ=255)

                       tgtACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATg                                                           <  1:3373228/46‑1 (MQ=255)

                       tgtACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATg                                                           <  1:2233300/46‑1 (MQ=255)

                       tgtACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATg                                                           <  1:1279330/46‑1 (MQ=255)

                       tgtACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATg                                                           <  1:1495757/46‑1 (MQ=255)

                       tgtACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATg                                                           <  1:2094809/46‑1 (MQ=255)

                       tgtACTGGCCTATAGCACGCCAGATCTCAGTATAAACGGCAAAATg                                                           <  1:1452284/46‑1 (MQ=255)

                           cTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTa                                >  1:2777242/1‑69 (MQ=255)

                           cTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTa                                >  1:685335/1‑69 (MQ=255)

                                   taGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtacctt                        >  1:1801813/1‑69 (MQ=255)

                                     gCACGCCCGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtacctta                       <  1:3250288/68‑1 (MQ=255)

                                     gCACGCCAGATCTCAGTATGAACGGTAAAATGATCTATGCAGCAATTACTTACACCCTActtacctta                       <  1:3342271/68‑1 (MQ=255)

                                     gCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtacctta                       <  1:2192390/68‑1 (MQ=255)

                                      cACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtaccttac                      >  1:737720/1‑68 (MQ=255)

                                      cACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtaccttac                      >  1:1853974/1‑68 (MQ=255)

                                      cACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtaccttac                      >  1:1276140/1‑68 (MQ=255)

                                        cGCCAGATCTCAGTATGAACGGAAAAATGATTTATGCAGCAATTACTTACACCCTACTTACCTTACttt                   <  1:1818093/69‑1 (MQ=255)

                                        cGCCAGATCTCAGTATGAACGGAAAAATGATTTATGCAGCAATTACTTACACCCTACTTACCTTACttt                   <  1:2579151/69‑1 (MQ=255)

                                        cGCCAGATCTCAGCATGAACGGAAAAATGATTTATGCAGCAATTACTTACACCCTACTTACCTTACttt                   <  1:3163103/69‑1 (MQ=255)

                                         gCCAGATCTCAGTATGAACGGAAAAATGATTTATGCAGCAATTACTTACACCCTACTTACCTTACttt                   <  1:1223696/68‑1 (MQ=255)

                                         gCCAGATCTCAGTATGAACGGAAAAATGATTTATGCAGCAATTACTTACACCCTACTTACCTTACttt                   <  1:1817400/68‑1 (MQ=255)

                                             gATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATAcc              <  1:1500972/69‑1 (MQ=255)

                                             gATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATAcc              <  1:1264697/69‑1 (MQ=255)

                                             gATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATAcc              <  1:1886783/69‑1 (MQ=255)

                                                ctcAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtacct                         >  1:3346835/1‑55 (MQ=255)

                                                ctcAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtacct                         >  1:1748320/1‑55 (MQ=255)

                                                ctcAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtacct                         >  1:646086/1‑55 (MQ=255)

                                                ctcAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtacct                         >  1:1813401/1‑55 (MQ=255)

                                                ctcAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtacct                         >  1:655843/1‑55 (MQ=255)

                                                ctcAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtacct                         >  1:123091/1‑55 (MQ=255)

                                                ctcAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtacct                         >  1:1144602/1‑55 (MQ=255)

                                                ctcAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtacct                         >  1:1991660/1‑55 (MQ=255)

                                                ctcAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtacct                         >  1:1045373/1‑55 (MQ=255)

                                                ctcAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtacct                         >  1:2086836/1‑55 (MQ=255)

                                                ctcAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtacct                         >  1:245586/1‑55 (MQ=255)

                                                     tATGAACTGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCaa      <  1:2533934/69‑1 (MQ=255)

                                                     tATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCgtcgtc        <  1:2043669/67‑1 (MQ=255)

                                                     tATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCgtcgtc        <  1:1127346/67‑1 (MQ=255)

                                                     tATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCgtcgtc        <  1:2803944/67‑1 (MQ=255)

                                                     tATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCaa      <  1:3397028/69‑1 (MQ=255)

                                                     tATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCaa      >  1:2215162/1‑69 (MQ=255)

                                                     tATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCaa      <  1:1249005/69‑1 (MQ=255)

                                                     tATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCaa      <  1:2518170/69‑1 (MQ=255)

                                                     tATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCaa      <  1:574500/69‑1 (MQ=255)

                                                     tATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCaa      <  1:3403318/69‑1 (MQ=255)

                                                     tATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCaa      >  1:3320386/1‑69 (MQ=255)

                                                     tATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCaa      <  1:736224/69‑1 (MQ=255)

                                                     tATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCaa      >  1:3172243/1‑69 (MQ=255)

                                                     tATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCaa      <  1:1766498/69‑1 (MQ=255)

                                                      aTGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCaa      <  1:2041920/68‑1 (MQ=255)

                                                      aTGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCaa      <  1:1719040/68‑1 (MQ=255)

                                                       tGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCAata    >  1:1705361/1‑69 (MQ=255)

                                                       tGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCAata    >  1:2420141/1‑69 (MQ=255)

                                                       tGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCAata    >  1:2738735/1‑69 (MQ=255)

                                                       tGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTAACCTTACTTTATACCGTCGTCAata    >  1:2768219/1‑69 (MQ=255)

                                                         aaCGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTGAaaata  >  1:1392156/1‑68 (MQ=255)

                                                         aaCGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCAatat   >  1:3234015/1‑68 (MQ=255)

                                                         aaCGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCAATATc  >  1:1553512/1‑69 (MQ=255)

                                                         aaCGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCAATATc  >  1:2638863/1‑69 (MQ=255)

                                                         aaCGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCAATATc  >  1:2444491/1‑69 (MQ=255)

GCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCAATATC  >  minE/2505803‑2505928

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