BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,388,901	G→A	34.9%	L181L (CTC→CTT)	pldB ←	lysophospholipase L(2)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,388,901	0	G	A	34.9%	88.0 / 66.8	117	L181L (CTC→CTT)	pldB	lysophospholipase L(2)
Reads supporting (aligned to +/- strand): ref base G (26/50); new base A (33/8); total (59/58)
Fisher's exact test for biased strand distribution p-value = 2.06e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 120 total aligned reads displayed.

CGCGCCAGCGCCCGGTGCCTATTGCATAGCCATCACGGAAACGTGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAGCGCAAT  >  minE/2388837‑2388969

         gcCCGGTGCCTATTGCATAGCCATCACGGAAACGTGGATGCGCTTCGGCCCAGTTGAGGATCTCCCGTg                                                         <  1:3153257/69‑1 (MQ=255)

              gTGCCTATTGCATAGCCATCACGGAAACGTGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATa                                                    <  1:950979/69‑1 (MQ=255)

              gTGCCTATTGCATAGCCATCACGGAAACGTGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATa                                                    <  1:923888/69‑1 (MQ=255)

              gTGCCTATTGCATAGCCATCACGGAAACGTGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATa                                                    <  1:2386073/69‑1 (MQ=255)

              gTGCCTATTGCATAGCCATCACGGAAACGTGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATa                                                    <  1:1998623/69‑1 (MQ=255)

              gTGCCTATTGCATAGCCATCACGGAAACGTGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATa                                                    <  1:1785217/69‑1 (MQ=255)

              gTGCCTATTGCATAGCCATCACGGAAACGTGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATa                                                    <  1:356765/69‑1 (MQ=255)

                 ccTATTGCATAGCCATCACGGAAACGTGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATAAAt                                                 <  1:2066324/69‑1 (MQ=255)

                    atcgCATAGCCATCACGGAAACGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGAc                                              <  1:2343399/66‑1 (MQ=255)

                    atcgCATAGCCATCACGGAAACGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGAc                                              <  1:2523937/66‑1 (MQ=255)

                    atcgCATAGCCATCACGGAAACGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGAc                                              <  1:1605799/66‑1 (MQ=255)

                    atcgCATAGCCATCACGGAAACGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGAc                                              <  1:1574049/66‑1 (MQ=255)

                    atcgCATAGCCATCACGGAAACGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGAc                                              <  1:1543456/66‑1 (MQ=255)

                    atcgCATAGCCATCACGGAAACGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGAc                                              <  1:872299/66‑1 (MQ=255)

                    atcgCATAGCCATCACGGAAACGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGAc                                              <  1:1318744/66‑1 (MQ=255)

                          tAGCCATCACGGAAACGTGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATAAATGACGGCa                                          >  1:2950841/1‑67 (MQ=255)

                          tAGCCATCACGGAAACGTGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATAAATGACGGCAt                                         >  1:740487/1‑68 (MQ=255)

                          tAGCCATCACGGAAACGTGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATAAATGACGGCAt                                         >  1:2805838/1‑68 (MQ=255)

                            gCCATCACGGAAACGTGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATAAATGACGGCATACg                                      >  1:2833810/1‑69 (MQ=255)

                            gCCATCACGGAAACGTGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATAAATGACGGCATACg                                      >  1:1261013/1‑69 (MQ=255)

                                  aCGGAAACGTGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATc                                  >  1:1356118/1‑67 (MQ=255)

                                  aCGGAAACGTGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCAc                                >  1:2985304/1‑69 (MQ=255)

                                  aCGGAAACGTGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCAc                                >  1:2302876/1‑69 (MQ=255)

                                  aCGGAAACGTGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCAc                                >  1:2714622/1‑69 (MQ=255)

                                  aCGGAAACGTGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCAc                                >  1:2848349/1‑69 (MQ=255)

                                  aCGGAAACGTGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCAc                                >  1:3366097/1‑69 (MQ=255)

                                  aCGGAAACGTGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCAc                                >  1:611406/1‑69 (MQ=255)

                                  aCGGAAACGTGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGACATAAATGACGGCATACGAATAAc                                >  1:2717597/1‑69 (MQ=255)

                                        aCGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAAt                                                 >  1:3227908/1‑46 (MQ=255)

                                        aCGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATg                                                >  1:2773519/1‑47 (MQ=255)

                                        aCGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACAATcc                           >  1:687560/1‑68 (MQ=255)

                                        aCGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACAATcc                           >  1:990772/1‑68 (MQ=255)

                                        aCGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACAATcc                           >  1:3436420/1‑68 (MQ=255)

                                        aCGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACAATcc                           >  1:2488900/1‑68 (MQ=255)

                                        aCGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACAATcc                           >  1:856470/1‑68 (MQ=255)

                                        aCGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACAATcc                           >  1:3336825/1‑68 (MQ=255)

                                        aCGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACAATcc                           >  1:465051/1‑68 (MQ=255)

                                        aCGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACAATcc                           >  1:3032605/1‑68 (MQ=255)

                                        aCGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACAATcc                           >  1:561641/1‑68 (MQ=255)

                                        aCGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACAATcc                           >  1:3253924/1‑68 (MQ=255)

                                        aCGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACAATcc                           >  1:423222/1‑68 (MQ=255)

                                        aCGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACAATcc                           >  1:3122756/1‑68 (MQ=255)

                                        aCGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACAATcc                           >  1:3015647/1‑68 (MQ=255)

                                        aCGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACAATcc                           >  1:1701464/1‑68 (MQ=255)

                                        aCGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACAATcc                           >  1:1423945/1‑68 (MQ=255)

                                        aCGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACAATcc                           >  1:1377941/1‑68 (MQ=255)

                                        aCGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACAATcc                           >  1:1840834/1‑68 (MQ=255)

                                        aCGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACAATcc                           >  1:1352996/1‑68 (MQ=255)

                                        aCGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACAATcc                           >  1:1816051/1‑68 (MQ=255)

                                        aCGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACAATcc                           >  1:1763143/1‑68 (MQ=255)

                                        aCGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACAATcc                           >  1:1710655/1‑68 (MQ=255)

                                        aCGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACAATcc                           >  1:2122756/1‑68 (MQ=255)

                                        aCGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACAATcc                           >  1:1646669/1‑68 (MQ=255)

                                        aCGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACAATcc                           >  1:2211294/1‑68 (MQ=255)

                                        aCGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACAATcc                           >  1:1666636/1‑68 (MQ=255)

                                        aCGTGGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACAATcc                           >  1:2293462/1‑68 (MQ=255)

                                        aCGTGGATGCGCTTCGGACAAGTTAAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACAATc                            >  1:3253965/1‑67 (MQ=255)

                                        aCGTCGATGCGCTTCGGCCCAGTTAAGGATCTGCCGTGCCATAAAt                                                 >  1:963505/1‑46 (MQ=255)

                                           tGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATAAATGACGGCATa                                        >  1:1174714/1‑52 (MQ=255)

                                           tGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCaaa                       >  1:2988434/1‑69 (MQ=255)

                                           tGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCaaa                       >  1:1319249/1‑69 (MQ=255)

                                           tGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCaaa                       >  1:1285452/1‑69 (MQ=255)

                                           tGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCaaa                       <  1:1274247/69‑1 (MQ=255)

                                           tGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCaaa                       <  1:1066052/69‑1 (MQ=255)

                                           tGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCaaa                       >  1:2385664/1‑69 (MQ=255)

                                           tGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCaa                        >  1:500135/1‑68 (MQ=255)

                                           tGGATGCACTTCGGCCCAGTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCaaa                       >  1:1296819/1‑69 (MQ=255)

                                                 cgcTTCGGCCCAGTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTgg                 >  1:2303249/1‑69 (MQ=255)

                                                        gCCCAGTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGt           <  1:759720/68‑1 (MQ=255)

                                                        gCCCAGTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGt           <  1:1984336/68‑1 (MQ=255)

                                                             gTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAgcgc     <  1:1877154/69‑1 (MQ=255)

                                                             gTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAgcgc     <  1:1036174/69‑1 (MQ=255)

                                                             gTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAgcgc     <  1:952352/69‑1 (MQ=255)

                                                             gTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAgcgc     <  1:2518972/69‑1 (MQ=255)

                                                             gTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAgcgc     <  1:1105659/69‑1 (MQ=255)

                                                             gTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAgcgc     <  1:1107894/69‑1 (MQ=255)

                                                             gTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAgcgc     <  1:1111887/69‑1 (MQ=255)

                                                             gTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAgcgc     <  1:833639/69‑1 (MQ=255)

                                                             gTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAgcgc     <  1:2673139/69‑1 (MQ=255)

                                                             gTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAgcgc     <  1:690736/69‑1 (MQ=255)

                                                             gTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAgcgc     <  1:1210254/69‑1 (MQ=255)

                                                             gTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAgcgc     <  1:2166448/69‑1 (MQ=255)

                                                             gTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAgcgc     <  1:2790661/69‑1 (MQ=255)

                                                             gTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAgcgc     <  1:3129494/69‑1 (MQ=255)

                                                             gTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAgcgc     <  1:442005/69‑1 (MQ=255)

                                                             gTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAgcgc     <  1:435753/69‑1 (MQ=255)

                                                             gTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAgcgc     <  1:1847144/69‑1 (MQ=255)

                                                             gTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAgcgc     <  1:4224/69‑1 (MQ=255)

                                                             gTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAgcgc     <  1:368075/69‑1 (MQ=255)

                                                             gTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAgcgc     <  1:2840499/69‑1 (MQ=255)

                                                             gTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAgcgc     <  1:3380194/69‑1 (MQ=255)

                                                             gTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAgcgc     <  1:1380669/69‑1 (MQ=255)

                                                             gTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAgcgc     <  1:2886874/69‑1 (MQ=255)

                                                             gTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAgcgc     <  1:1816685/69‑1 (MQ=255)

                                                             gTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAgcgc     <  1:30662/69‑1 (MQ=255)

                                                             gTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAgcgc     <  1:3067440/69‑1 (MQ=255)

                                                             gTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAgcgc     <  1:3181373/69‑1 (MQ=255)

                                                             gTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAgcgc     <  1:1699246/69‑1 (MQ=255)

                                                                gAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAGCGCaa   >  1:286637/1‑68 (MQ=255)

                                                                gAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAGCGCAAt  >  1:248618/1‑69 (MQ=255)

CGCGCCAGCGCCCGGTGCCTATTGCATAGCCATCACGGAAACGTGGATGCGCTTCGGCCCAGTTGAGGATCTGCCGTGCCATAAATGACGGCATACGAATCACGATCCCAAACATTGGCGCAGTTAGCGCAAT  >  minE/2388837‑2388969

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