BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,283,633	T→G	39.0%	V1804V (GTT→GTG)	yeeJ →	adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,283,633	0	T	G	39.0%	71.0 / 148.7	154	V1804V (GTT→GTG)	yeeJ	adhesin
Reads supporting (aligned to +/- strand): ref base T (44/50); new base G (29/31); total (73/81)
Fisher's exact test for biased strand distribution p-value = 8.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

Only 100 of 160 total aligned reads displayed.

TCAGGTAATGGTCGCACAAGACGTAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGCCAGAATACGGTCTC  >  minE/1283567‑1283701

tCAGGTAATGGTCGCACAAGACGTAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTTg                                                                     >  1:351678/1‑68 (MQ=255)

tCAGGTAATGGTCGCACAAGACGTAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTTg                                                                     >  1:3027910/1‑68 (MQ=255)

tCAGGTAATGGTCGCACAAGACGTAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTTg                                                                     >  1:838701/1‑68 (MQ=255)

tCAGGTAATGGTCGCACAAGACGTAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTTg                                                                     >  1:1411904/1‑68 (MQ=255)

tCAGGTAATGGTCGCACAAGACGTAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTTg                                                                     >  1:3128949/1‑68 (MQ=255)

tCAGGTAATGGTCGCACAAGACGTAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTTg                                                                     >  1:1330521/1‑68 (MQ=255)

            ctcacaAGACGTAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGt                                                        <  1:2628922/67‑1 (MQ=255)

            cGCACAAGACGTAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGt                                                        <  1:1878810/69‑1 (MQ=255)

            cGCACAAGACGTAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGt                                                        <  1:538715/69‑1 (MQ=255)

            cGCACAAGACGTAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGt                                                        <  1:172107/69‑1 (MQ=255)

            cGCACAAGACGTAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGt                                                        <  1:230795/69‑1 (MQ=255)

            cGCACAAGACGTAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGt                                                        <  1:2331235/69‑1 (MQ=255)

            cGCACAAGACGTAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGt                                                        <  1:1277279/69‑1 (MQ=255)

                       tAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTGACGCATCAACCCGCGACATTCAg                                               >  1:263726/1‑67 (MQ=255)

                       tAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTGACGCATCAACCCGCGACATTCAg                                               >  1:3143871/1‑67 (MQ=255)

                       tAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTGACGCATCAACCCGCGACATTCAg                                               >  1:1046275/1‑67 (MQ=255)

                       tAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTGACGCATCAACCCGCGACATTCAg                                               >  1:2455226/1‑67 (MQ=255)

                       tAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTGACGCATCAACCCGCGACATTCAg                                               >  1:1247586/1‑67 (MQ=255)

                       tAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTGACGCATCAACCCGCGACATTCAg                                               >  1:3374450/1‑67 (MQ=255)

                       tAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTGACGCATCAACCCGCGACATTCAg                                               >  1:2215375/1‑67 (MQ=255)

                       tAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTGACGCATCAACCCGCGACATTCAg                                               >  1:781646/1‑67 (MQ=255)

                       tAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTGACGCATCAACCCGCGACATTCAg                                               >  1:1156803/1‑67 (MQ=255)

                       tAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTGACGCATCAACCCGCGACATTCAg                                               >  1:966412/1‑67 (MQ=255)

                       tAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTGACGCATCAACCCGCGACATTCAg                                               >  1:2012720/1‑67 (MQ=255)

                              tAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGTGACATTCa                                                >  1:1544610/1‑59 (MQ=255)

                              tAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGcc                                      >  1:2844534/1‑69 (MQ=255)

                              tAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGcc                                      >  1:3082618/1‑69 (MQ=255)

                               aaaGCTCACGTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGcc                                      >  1:2420070/1‑68 (MQ=255)

                               aaaGCTCACGTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCa                                     >  1:457194/1‑69 (MQ=255)

                               aaaGCTCACGTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCa                                     >  1:2454369/1‑69 (MQ=255)

                               aaaGCTCACGTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCa                                     >  1:235274/1‑69 (MQ=255)

                               aaaGCTCACGTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCa                                     >  1:2077185/1‑69 (MQ=255)

                               aaaGCTCACGTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCa                                     >  1:1260548/1‑69 (MQ=255)

                               aaaGCTCACGTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCa                                     >  1:2761625/1‑69 (MQ=255)

                                        gTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCa                             <  1:3222699/68‑1 (MQ=255)

                                        gTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCa                             <  1:2497275/68‑1 (MQ=255)

                                        gTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCa                             <  1:1440732/68‑1 (MQ=255)

                                        gTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCa                             <  1:2664338/68‑1 (MQ=255)

                                        gTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCa                             <  1:1355326/68‑1 (MQ=255)

                                        gTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCa                             <  1:1632077/68‑1 (MQ=255)

                                        gTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCa                             <  1:2030673/68‑1 (MQ=255)

                                        gTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCa                             <  1:2888990/68‑1 (MQ=255)

                                        gTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCa                             <  1:2921562/68‑1 (MQ=255)

                                        gTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCa                             <  1:2987925/68‑1 (MQ=255)

                                        gTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCa                             <  1:1850185/68‑1 (MQ=255)

                                              gACCAGTTTGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGatc                      <  1:741137/69‑1 (MQ=255)

                                                  aGTTTGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGAtcatca                  <  1:87668/69‑1 (MQ=255)

                                                  aGTTTGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGAtcatca                  <  1:3353212/69‑1 (MQ=255)

                                                  aGTTTGGCAACCCGGTGACGCATCAACCCGCGACATTCAGTGCAGCGCCATCCTCGCAAATGAtcatca                  <  1:2745746/69‑1 (MQ=255)

                                                  aGTTTGGCAACCCGGTGACGCATCAACCCGCGACATTCAGTGCAGCGCCATCCTCGCAAATGAtcatca                  <  1:2797582/69‑1 (MQ=255)

                                                  aGTTTGGCAACCCGGTGACGCATCAACCCGCGACATTCAGTGCAGCGCCATCCTCGCAAATGAtcatca                  <  1:265572/69‑1 (MQ=255)

                                                  aGTTTGGCAACCCGGTGACGCATCAACCCGCGACATTCAGTGCAGCGCCATCCTCGCAAATGAtcatca                  <  1:5016/69‑1 (MQ=255)

                                                  aGTTTGGCAACCCGGTGACGCATCAACCCGCGACATTCAGTGCAGCGCCATCCTCGCAAATGAtcatca                  <  1:197313/69‑1 (MQ=255)

                                                  aGTTTGGCAACCCGGTGACGCATCAACCCGCGACATTCAGTGCAGCGCCATCCTCGCAAATGAtcatca                  <  1:1993109/69‑1 (MQ=255)

                                                  aGTTTGGCAACCCGGTGACGCATCAACCCGCGACATTCAGTGCAGCGCCATCCTCGCAAATGAtcatca                  <  1:651460/69‑1 (MQ=255)

                                                  aGTTTGGCAACCCGGTGACGCATCAACCCGCGACATTCAGTGCAGCGCCATCCTCGCAAATGAtcatca                  <  1:611915/69‑1 (MQ=255)

                                                  aGTTTGGCAACCCGGTGACGCATCAACCCGCGACATTCAGTGCAGCGCCATCCTCGCAAATGAtcatca                  <  1:1784725/69‑1 (MQ=255)

                                                  aGTTTGGCAACCCGGTGACGCATCAACCCGCGACATTCAGTGCAGCGCCATCCTCGCAAATGAtcatca                  <  1:2234343/69‑1 (MQ=255)

                                                     ttGGCACCCCGGTGACGCATCAACCCGCGACATTCAGTGCAGCGCCATCCTCGCAAATGatcat                    >  1:2785270/1‑64 (MQ=255)

                                                     ttGGCAACCCGGTGACGCATCAACCCGCGACATTCAGTGCAGCGCCATCCTCGCAAATGAtcatca                  >  1:677780/1‑66 (MQ=255)

                                                     ttGGCAACCCGGTGACGCATCAACCCGCGACATTCAGTGCAGCGCCATCCTCGCAAATGAtcatca                  >  1:440205/1‑66 (MQ=255)

                                                     ttGGCAACCCGGTGACGCATCAACCCGCGACATTCAGTGCAGCGCCATCCTCGCAAATGAtcatca                  >  1:1728508/1‑66 (MQ=255)

                                                     ttGGCAACCCGGTGACGCATCAACCCGCGACATTCAGTGCAGCGCCATCCTCGCAAATGAtcatca                  >  1:1284713/1‑66 (MQ=255)

                                                     ttGGCAACCCGGTGACGCATCAACCCGCGACATTCAGTGCAGCGCCATCCTCGCAAATGAtcatca                  >  1:3392176/1‑66 (MQ=255)

                                                     ttGGCAACCCGGTGACGCATCAACCCGCGACATTCAGTGCAGCGCCATCCTCGCAAATGAtcatca                  >  1:2774322/1‑66 (MQ=255)

                                                     ttGGCAACCCGGTGACGCATCAACCCGCGACATTCAGTGCAGCGCCATCCTCGCAAATGAtcatca                  >  1:590036/1‑66 (MQ=255)

                                                     ttGGCAACCCGGTGACGCATCAACCCGCGACATTCAGTGCAGCGCCATCCTCGCAAATGAtcatca                  >  1:2524789/1‑66 (MQ=255)

                                                     ttGGCAACCCGGTGACGCATCAACCCGCGACATTCAGTGCAGCGCCATCCTCGCAAATGAtcatca                  >  1:2651329/1‑66 (MQ=255)

                                                     ttGGCAACCCGGTGACGCATCAACCCGCGACATTCAGTGCAGCGCCATCCTCGCAAATGAtcatca                  >  1:705968/1‑66 (MQ=255)

                                                      tGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATg                         >  1:2551830/1‑58 (MQ=255)

                                                      tGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATg                         >  1:1795270/1‑58 (MQ=255)

                                                      tGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGcc               >  1:1232607/1‑68 (MQ=255)

                                                      tGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGcc               >  1:570494/1‑68 (MQ=255)

                                                      tGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGcc               >  1:1927183/1‑68 (MQ=255)

                                                        gCAACCCGGTTGCGCATCAACCCGTTACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGCCAg             <  1:310945/68‑1 (MQ=255)

                                                        gCAACCCGGTTGCGCATCAACCCGTTACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGCCAg             <  1:3087680/68‑1 (MQ=255)

                                                        gCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCa                             <  1:1230116/52‑1 (MQ=255)

                                                        gCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGCCAg             <  1:893808/68‑1 (MQ=255)

                                                        gCAACCCGGTGACGCATCAACCCGCGACATTCAGTGCAGCGCCATCCTCGCAAATGATCATCAGCCAg             <  1:1767154/68‑1 (MQ=255)

                                                        gCAACCCGGTGACGCATCAACCCGCGACATTCAGTGCAGCGCCATCCTCGCAAATGATCATCAGCCAg             <  1:1101434/68‑1 (MQ=255)

                                                        gCAACCCGGTGACGCATCAACCCGCGACATTCAGTGCAGCGCCATCCTCGCAAATGATCATCAGCCAg             <  1:3087929/68‑1 (MQ=255)

                                                        gCAACCCGGTGACGCATCAACCCGCGACATTCAGTGCAGCGCCATCCTCGCAAATGATCATCAGCCAg             <  1:2273663/68‑1 (MQ=255)

                                                          aaCCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGCCAg             <  1:2475685/66‑1 (MQ=255)

                                                              cGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGCCAGAATACgg      >  1:1505943/1‑69 (MQ=255)

                                                              cGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGCCAGAATACgg      >  1:1166542/1‑69 (MQ=255)

                                                              cGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGCCAGAATACgg      >  1:841837/1‑69 (MQ=255)

                                                              cGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGCCAGAATACgg      >  1:2717074/1‑69 (MQ=255)

                                                              cGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGCCAGAATACgg      >  1:2563652/1‑69 (MQ=255)

                                                              cGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGCCAGAATACgg      >  1:14359/1‑69 (MQ=255)

                                                              cGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGCCAGAATACgg      >  1:2113080/1‑69 (MQ=255)

                                                              cGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGCCAGAATACgg      >  1:1598744/1‑69 (MQ=255)

                                                              cGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGCCAGAATACgg      >  1:1631440/1‑69 (MQ=255)

                                                              cGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGCCAGAATACgg      >  1:3321552/1‑69 (MQ=255)

                                                              cGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGCCAGAATACgg      >  1:1956498/1‑69 (MQ=255)

                                                              cGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGCCAGAATACgg      >  1:390221/1‑69 (MQ=255)

                                                                  tGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGCCAGAATACGGtctc  <  1:603184/69‑1 (MQ=255)

                                                                  tGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGCCAGAATACGGtctc  <  1:13987/69‑1 (MQ=255)

                                                                  tGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGCCAGAATACGGtctc  <  1:1362620/69‑1 (MQ=255)

                                                                  tGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGCCAGAATACGGtctc  <  1:215684/69‑1 (MQ=255)

                                                                  tGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGCCAGAATACGGtctc  <  1:2819224/69‑1 (MQ=255)

TCAGGTAATGGTCGCACAAGACGTAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGCCAGAATACGGTCTC  >  minE/1283567‑1283701

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