BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,283,647	T→C	43.2%	V1809A (GTG→GCG)	yeeJ →	adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,283,647	0	T	C	43.2%	37.9 / 113.0	141	V1809A (GTG→GCG)	yeeJ	adhesin
Reads supporting (aligned to +/- strand): ref base T (35/45); new base C (29/32); total (64/77)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 154 total aligned reads displayed.

CGCACAAGACGTAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGCCAGAATACGGTCTCTACTAATACGCAG  >  minE/1283579‑1283714

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

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

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

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

cGCACAAGACGTAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGt                                                                     <  1:1954896/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)

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

cGCACAAGACGTAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGCTGCGCATCAACCCGt                                                                     <  1:743804/69‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                       gTTTGGCAACCCGGTGACGCATCAACCCGCGACATTCAGTGCAGCGCCATCCTCGCAAATGAtcatca                               <  1:1991760/68‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                 ccGCGACATTCAGTGCAGCGCCATCCTCGCAAATGAtcatca                               <  1:3145058/42‑1 (MQ=38)

                                                                   gTGACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGCCAGAATACGGTCTCTACTAATACGCAg  <  1:503152/69‑1 (MQ=255)

CGCACAAGACGTAGCAGTTAAAGCTCACGTTAACGACCAGTTTGGCAACCCGGTTGCGCATCAACCCGTGACATTCAGTGCAGAGCCATCCTCGCAAATGATCATCAGCCAGAATACGGTCTCTACTAATACGCAG  >  minE/1283579‑1283714

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