BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,970	T→C	17.2%	K428K (AAA→AAG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,970	0	T	C	17.2%	263.8 / 45.7	116	K428K (AAA→AAG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base T (42/54); new base C (11/9); total (53/63)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 129 total aligned reads displayed.

CAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATC  >  minE/608900‑609039

 aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                                                      <  1:2912248/71‑1 (MQ=255)

 aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                                                      >  1:3542483/1‑71 (MQ=255)

 aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                                                      >  1:2708089/1‑71 (MQ=255)

 aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGt                                                                       <  1:2048483/70‑1 (MQ=255)

 aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGt                                                                       <  1:2327653/70‑1 (MQ=255)

  aCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGttt                                                                     <  1:2576364/71‑1 (MQ=255)

       gccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCaa                                                                >  1:2542308/1‑71 (MQ=255)

        ccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAAc                                                               <  1:3946145/71‑1 (MQ=255)

          gTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGa                                                             >  1:2423550/1‑71 (MQ=255)

          gTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGa                                                             >  1:3502839/1‑71 (MQ=255)

             gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGat                                                            >  1:1525207/1‑69 (MQ=255)

              aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcat                                                         <  1:3354044/71‑1 (MQ=255)

              aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                                          <  1:3074143/70‑1 (MQ=255)

              aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                                          <  1:3315875/70‑1 (MQ=255)

              aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                                          <  1:2226859/70‑1 (MQ=255)

                 gCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCg                                                       <  1:3931339/70‑1 (MQ=255)

                   gTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGc                                                    >  1:1335388/1‑71 (MQ=255)

                      cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                                                 <  1:3093794/71‑1 (MQ=255)

                        gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                                                      <  1:2632525/48‑1 (MQ=255)

                        gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCggg                                                     <  1:1882474/65‑1 (MQ=255)

                        gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                                      >  1:3197307/1‑64 (MQ=255)

                           tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcag                                            >  1:2136562/1‑71 (MQ=255)

                           gcgTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcag                                            <  1:574363/69‑1 (MQ=255)

                            cgTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcag                                            >  1:3598461/2‑70 (MQ=255)

                            cgTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcag                                            <  1:1296101/69‑1 (MQ=255)

                             gTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTAcgccgc                                           <  1:1695153/70‑1 (MQ=255)

                             gTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcagc                                           <  1:2876279/70‑1 (MQ=255)

                             gTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcagc                                           >  1:3293049/1‑70 (MQ=255)

                              ttttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagcag                                         >  1:3433081/1‑70 (MQ=255)

                               tttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagcagg                                        >  1:3199830/1‑69 (MQ=255)

                                ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAAGGATCATCGGGCTTACGCagcaag                                        >  1:1778673/1‑70 (MQ=255)

                                ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                                      >  1:3470676/1‑56 (MQ=255)

                                ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaag                                        >  1:527550/1‑70 (MQ=255)

                                ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaag                                        >  1:1305649/1‑70 (MQ=255)

                                 tCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCagcaggc                                       >  1:3820084/1‑70 (MQ=255)

                                     cAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                                          >  1:437652/1‑47 (MQ=255)

                                     cAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                                          >  1:2690320/1‑47 (MQ=255)

                                     cAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAAtcg                                  <  1:2406376/71‑3 (MQ=255)

                                      aGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc                                 >  1:746943/1‑71 (MQ=255)

                                      aGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGc                                 >  1:1523287/1‑71 (MQ=255)

                                       ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCATCAAGCAATAGCg                                >  1:3409611/1‑71 (MQ=255)

                                       ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                                >  1:1617089/1‑71 (MQ=255)

                                       ggTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATcgcg                                >  1:584689/1‑71 (MQ=255)

                                        gTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAAc                                                               <  1:2842939/39‑1 (MQ=255)

                                          ttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAg                             >  1:3045171/1‑71 (MQ=255)

                                           tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt                            >  1:1450037/1‑71 (MQ=255)

                                             cACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCa                          >  1:1887477/1‑71 (MQ=255)

                                             cACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCa                          >  1:1403558/1‑71 (MQ=255)

                                              aCGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                     <  1:430575/59‑1 (MQ=255)

                                              aCGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                     <  1:435288/59‑1 (MQ=255)

                                               cgtgcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTc                        <  1:3689070/68‑1 (MQ=255)

                                               cgcgCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATc                        <  1:2905324/71‑1 (MQ=255)

                                               cgcgCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATc                        <  1:3810445/71‑1 (MQ=255)

                                                gtgcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTaa                             >  1:3461263/3‑64 (MQ=255)

                                                gcgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCg                       <  1:2264441/71‑1 (MQ=255)

                                                 cgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt                                                   >  1:238756/1‑42 (MQ=255)

                                                  gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt                     <  1:2275509/71‑1 (MQ=255)

                                                   cACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCgctg                    <  1:444757/71‑1 (MQ=255)

                                                    aCCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCgctgt                   <  1:1484372/71‑1 (MQ=255)

                                                      caaagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATcg                          >  1:1257362/3‑61 (MQ=255)

                                                      cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcag                                            >  1:3632553/1‑44 (MQ=255)

                                                      cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCTTCgt                      <  1:203859/66‑1 (MQ=255)

                                                      cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                 <  1:404777/71‑1 (MQ=255)

                                                      cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGAAAGCAATAGCGTAGTCATCGttgttg                 <  1:76445/71‑1 (MQ=255)

                                                       gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGa                <  1:2835979/71‑1 (MQ=255)

                                                        aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACg                                               >  1:2513884/1‑39 (MQ=255)

                                                        aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACg                                               >  1:3085477/1‑39 (MQ=255)

                                                        aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagca                                      <  1:2171720/48‑1 (MQ=255)

                                                        aagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGctgttg                 >  1:2354153/1‑69 (MQ=255)

                                                          gaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                   <  1:909458/65‑1 (MQ=255)

                                                          gaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                   <  1:1567467/65‑1 (MQ=255)

                                                          gaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGt             >  1:2804005/1‑71 (MQ=255)

                                                          gaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGt             <  1:3471068/71‑1 (MQ=255)

                                                          gaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGt             <  1:2700263/71‑1 (MQ=255)

                                                          gaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGt             <  1:3833286/71‑1 (MQ=255)

                                                          gaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGt             <  1:3451581/71‑1 (MQ=255)

                                                           aaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                                <  1:264776/51‑1 (MQ=255)

                                                           aaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                                <  1:2019029/51‑1 (MQ=255)

                                                           aaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                                <  1:2014626/51‑1 (MQ=255)

                                                           aaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTc            <  1:1842783/71‑1 (MQ=255)

                                                            aCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                   <  1:7253/63‑1 (MQ=255)

                                                            aCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                   <  1:1412339/63‑1 (MQ=255)

                                                            aCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                   <  1:2532501/63‑1 (MQ=255)

                                                            aCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                   <  1:1416918/63‑1 (MQ=255)

                                                                cATTTGTTTACCAACGATCATCGGGCTTACgcagc                                           <  1:1246759/35‑1 (MQ=255)

                                                                cATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCgg         <  1:908338/69‑1 (MQ=255)

                                                                cATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGAc       >  1:1919654/1‑71 (MQ=255)

                                                                cATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGCTGTTGAAGTCCGGAc       <  1:796670/71‑1 (MQ=255)

                                                                  tttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                                <  1:1631342/44‑1 (MQ=255)

                                                                  tttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                 <  1:1793417/59‑1 (MQ=255)

                                                                  tttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACg      >  1:1357/1‑70 (MQ=255)

                                                                   ttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGa        >  1:2575094/1‑67 (MQ=255)

                                                                   ttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGa        >  1:3136309/1‑67 (MQ=255)

                                                                   ttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCa    >  1:1980019/1‑71 (MQ=255)

                                                                    tGTTTACCAACGATCATCGGGCTTACGCAGCGAGCAATAGCGTa                              >  1:592778/1‑44 (MQ=255)

                                                                     gTTTACCAACGATCATCGGGCTTACGCagcaagca                                      <  1:185198/35‑1 (MQ=255)

                                                                     gTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAt   >  1:1297535/1‑70 (MQ=255)

                                                                      tttACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttg                    <  1:1340539/52‑1 (MQ=255)

                                                                      tttACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATc  <  1:1992040/70‑1 (MQ=255)

                                                                      tttACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATc  <  1:2813371/70‑1 (MQ=255)

CAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATC  >  minE/608900‑609039

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