BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,422,099	C→T	19.3%	V79V (GTG→GTA)	wzxE ←	O‑antigen translocase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,422,099	0	C	T	19.3%	174.9 / 26.8	113	V79V (GTG→GTA)	wzxE	O‑antigen translocase
Reads supporting (aligned to +/- strand): ref base C (50/41); new base T (15/7); total (65/48)
Fisher's exact test for biased strand distribution p-value = 3.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.85e-01

Only 100 of 113 total aligned reads displayed.

GCCAGCACAAAAACCAGCGCCATCAGCGTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGCAACGTATTTGGTTACACCGTTAAA  >  minE/2422034‑2422161

gCCAGCACAAAAACCAGCGCCATCAGCGTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCAcg                                                             >  1:1360709/1‑69 (MQ=255)

   aGCACAAAAACCAGCGCCATTAGCGTAGAGAAGCCAAGTACCATCGCTGATGAAGTGCCGACTAcgc                                                            <  1:662715/67‑1 (MQ=255)

   aGCACAAAAACCAGCGCCATTAGCGTAGAGAAGCCAAGTACCATCGCTGATGAAGTGCCGACTAcgc                                                            <  1:2888025/67‑1 (MQ=255)

   aGCACAAAAACCAGCGCCATTAGCGTAGAGAAGCCAAGTACCATCGCTGATGAAGTGCCGACTAcgc                                                            <  1:2238634/67‑1 (MQ=255)

   aGCACAAAAACCAGCGCCATTAGCGTAGAGAAGCCAAGTACCATCGCTGATGAAGTGCCGACTAcgc                                                            <  1:1243116/67‑1 (MQ=255)

   aGCACAAAAACCAGCGCCATTAGCGTAGAGAAGCCAAGTACCATCGCTGATGAAGTGCCGACTAcgc                                                            <  1:3100722/67‑1 (MQ=255)

                      tCAGCGTAGAGAAACCAAGTACCATCGCTGATGAAGTGTCGACCACGCGGCGCAGCTGTTGCGGATTAt                                       <  1:617236/69‑1 (MQ=255)

                      tCAGCGTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTAt                                       <  1:2699713/69‑1 (MQ=255)

                      tCAGCGTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTAt                                       <  1:2831446/69‑1 (MQ=255)

                      tCAGCGTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTAt                                       <  1:3091915/69‑1 (MQ=255)

                      tCAGCGTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTAt                                       <  1:2792772/69‑1 (MQ=255)

                      tCAGCGTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTAt                                       <  1:1213460/69‑1 (MQ=255)

                      tCAGCGTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTAt                                       <  1:2431950/69‑1 (MQ=255)

                      tCAGCGTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTAt                                       <  1:35733/69‑1 (MQ=255)

                      tCAGCGTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTAt                                       <  1:2547074/69‑1 (MQ=255)

                        aGCGTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCa                                     >  1:2816499/1‑69 (MQ=255)

                        aGCGTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCa                                     >  1:655362/1‑69 (MQ=255)

                        aGCGTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCa                                     >  1:54189/1‑69 (MQ=255)

                        aGCGTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCa                                     >  1:956400/1‑69 (MQ=255)

                        aGCGTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCa                                     >  1:1807192/1‑69 (MQ=255)

                        aGCGTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCa                                     >  1:1775775/1‑69 (MQ=255)

                        aGCGTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCa                                     >  1:1085361/1‑69 (MQ=255)

                        aGCGTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCa                                     >  1:1059284/1‑69 (MQ=255)

                        aGCGTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCa                                     >  1:1210671/1‑69 (MQ=255)

                        aGCGTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCa                                     >  1:1504839/1‑69 (MQ=255)

                           gTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACTCGGCGCAGCTGTTGCGGATTATCATgg                                  <  1:1876899/69‑1 (MQ=255)

                           gTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATgg                                  <  1:2236135/69‑1 (MQ=255)

                           gTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATgg                                  <  1:603903/69‑1 (MQ=255)

                           gTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATgg                                  <  1:633088/69‑1 (MQ=255)

                           gTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATgg                                  <  1:2067056/69‑1 (MQ=255)

                           gTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATgg                                  <  1:2830972/69‑1 (MQ=255)

                           gTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATgg                                  <  1:3342516/69‑1 (MQ=255)

                           gTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATAATgg                                  <  1:920128/69‑1 (MQ=255)

                            tAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTa                                        >  1:1183697/1‑62 (MQ=255)

                            tAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGt                                 >  1:1561057/1‑69 (MQ=255)

                            tAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGt                                 >  1:37551/1‑69 (MQ=255)

                            tAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGt                                 >  1:1801210/1‑69 (MQ=255)

                            tAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGt                                 <  1:659009/69‑1 (MQ=255)

                            tAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGt                                 >  1:2529278/1‑69 (MQ=255)

                            tAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTCTTGCGGATTATCATGGt                                 >  1:534781/1‑69 (MQ=255)

                            tAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGAGGCGCAGCTGTTGCGGATTATCATGGt                                 <  1:2397762/69‑1 (MQ=255)

                              gagaAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTAc                               >  1:593003/1‑69 (MQ=255)

                              gagaAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTAc                               >  1:3373222/1‑69 (MQ=255)

                              gagaAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTAc                               >  1:3357753/1‑69 (MQ=255)

                              gagaAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTAc                               >  1:2963675/1‑69 (MQ=255)

                              gagaAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTAc                               >  1:2918332/1‑69 (MQ=255)

                              gagaAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTAc                               >  1:3344072/1‑69 (MQ=255)

                              gagaAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTAc                               >  1:2594611/1‑69 (MQ=255)

                              gagaAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTAc                               >  1:124775/1‑69 (MQ=255)

                              gagaAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTAc                               >  1:709834/1‑69 (MQ=255)

                              gagaAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTAc                               >  1:695167/1‑69 (MQ=255)

                              gagaAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTAc                               >  1:1870264/1‑69 (MQ=255)

                              gagaAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTAc                               >  1:2468445/1‑69 (MQ=255)

                              gagaAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTAc                               >  1:2514911/1‑69 (MQ=255)

                                   aCCAAGTACCATCGCTGATGAAGTGCCGACCAcgcg                                                           >  1:2891063/1‑36 (MQ=255)

                                   aCCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGc                          <  1:71069/69‑1 (MQ=255)

                                   aCCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGc                          <  1:3313598/69‑1 (MQ=255)

                                   aCCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGc                          <  1:3268988/69‑1 (MQ=255)

                                   aCCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGc                          <  1:3241647/69‑1 (MQ=255)

                                   aCCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGc                          <  1:3200417/69‑1 (MQ=255)

                                   aCCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGc                          <  1:3128076/69‑1 (MQ=255)

                                   aCCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGc                          <  1:3083478/69‑1 (MQ=255)

                                   aCCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGc                          <  1:2500540/69‑1 (MQ=255)

                                   aCCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGc                          <  1:173176/69‑1 (MQ=255)

                                   aCCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGc                          <  1:655474/69‑1 (MQ=255)

                                   aCCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGc                          <  1:2167055/69‑1 (MQ=255)

                                        gTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTAt                                       <  1:866030/51‑1 (MQ=255)

                                                gCTGATGAAGTGCCGACTACGCGGCGCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGtt             >  1:2762669/1‑69 (MQ=255)

                                                gCTGATGAAGTGCCGACTACGCGGCGCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGtt             >  1:414801/1‑69 (MQ=255)

                                                gCTGATGAAGTGCCGACTACGCGGCGCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGtt             >  1:334071/1‑69 (MQ=255)

                                                gCTGATGAAGTGCCGACTACGCGGCGCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGtt             >  1:659733/1‑69 (MQ=255)

                                                gCTGATGAAGTGCCGACTACGCGGCGCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGtt             >  1:3188032/1‑69 (MQ=255)

                                                gCTGATGAAGTGCCGACTACGCGGCGCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGtt             >  1:1007605/1‑69 (MQ=255)

                                                gCTGATGAAGTGCCGACTACGCGGCGCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGtt             >  1:1000292/1‑69 (MQ=255)

                                                gCTGATGAAGTGCCGACTACGCGGCGCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGtt             >  1:2750947/1‑69 (MQ=255)

                                                gCTGATGAAGTGCCGACTACGCGGCGCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGtt             >  1:2555261/1‑69 (MQ=255)

                                                gCTGATGAAGTGCCGACTACGCGGCGCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGtt             >  1:1721649/1‑69 (MQ=255)

                                                gCTGATGAAGTGCCGACTACGCGGCGCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGtt             >  1:1751319/1‑69 (MQ=255)

                                                gCTGATGAAGTGCCGACTACGCGGCGCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGtt             >  1:225686/1‑69 (MQ=255)

                                                gCTGATGAAGTGCCGACTACGCGGCGCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGtt             >  1:1993834/1‑69 (MQ=255)

                                                gCTGATGAAGTGCCGACTACGCGGCGCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGt              >  1:2916896/1‑68 (MQ=255)

                                                gCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGCAACGTATTTGGt              >  1:33423/1‑68 (MQ=255)

                                                gCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGCAACGTATTTGGt              >  1:1233826/1‑68 (MQ=255)

                                                gCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGCAACGTATTTGGt              >  1:1246493/1‑68 (MQ=255)

                                                gCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGCAACGTATTTGGt              >  1:707959/1‑68 (MQ=255)

                                                gCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGCAACGTATTTGGt              >  1:1966012/1‑68 (MQ=255)

                                                gCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGCAACGTATTTGGt              >  1:1997088/1‑68 (MQ=255)

                                                gCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGCAACGTATTTGGt              >  1:2306794/1‑68 (MQ=255)

                                                gCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGCAACGTATTTGGt              >  1:237337/1‑68 (MQ=255)

                                                gCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGCAACGTATTTGGt              >  1:2645982/1‑68 (MQ=255)

                                                gCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGCAACGTATTTGGt              >  1:2829720/1‑68 (MQ=255)

                                                      gaAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGCAACGTATTTGGTTacac         >  1:651833/1‑67 (MQ=255)

                                                           gCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGCAACGTATTTGGTTACACCGTTaaa  <  1:2029025/69‑1 (MQ=255)

                                                           gCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGCAACGTATTTGGTTACACCGTTaaa  <  1:2549463/69‑1 (MQ=255)

                                                           gCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGCAACGTATTTGGTTACACCGTTaaa  <  1:2782323/69‑1 (MQ=255)

                                                            ccGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGCAACGTATTTGGTTACACCg       >  1:1592881/1‑63 (MQ=255)

                                                            ccGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGCAACGTATTTGGTTACACCGTTaaa  <  1:1487972/68‑1 (MQ=255)

                                                            ccGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGCAACGTATTTGGTTACACCGTTaaa  <  1:817055/68‑1 (MQ=255)

                                                            ccGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGCAACGTATTTGGTTACACCGTTaaa  <  1:2717732/68‑1 (MQ=255)

                                                            ccGACCACGCGGCGCAGATGTTGCGGATTATCATGGt                                 <  1:2090355/37‑1 (MQ=255)

GCCAGCACAAAAACCAGCGCCATCAGCGTAGAGAAACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGCAACGTATTTGGTTACACCGTTAAA  >  minE/2422034‑2422161

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