BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,981	T→C	15.8%	I425V (ATC→GTC) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,981	0	T	C	15.8%	280.8 / 42.2	120	I425V (ATC→GTC) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base T (39/62); new base C (10/9); total (49/71)
Fisher's exact test for biased strand distribution p-value = 3.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 130 total aligned reads displayed.

GATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGT  >  minE/608913‑609049

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

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

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

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

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

       tACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt                                                             <  1:1903957/71‑1 (MQ=255)

         cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                                                           <  1:3093794/71‑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)

                gTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcagca                                                    <  1:256287/71‑1 (MQ=255)

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

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

                 ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcagca                                                    <  1:453039/70‑1 (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)

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

                       ccAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATa                                             <  1:2004439/71‑1 (MQ=255)

                       ccAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAAtc                                             <  1:2743812/71‑2 (MQ=255)

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

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

                        cAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                                                    >  1:437652/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)

                            tttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                                          <  1:3466562/69‑1 (MQ=255)

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

                              tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAg                                       <  1:1345190/70‑1 (MQ=255)

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

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

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

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

                                  cgtgcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTc                                  <  1:3689070/68‑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)

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

                                         caaagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCgctgtt                            >  1:2861963/3‑70 (MQ=255)

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

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

                                         cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGAAAGCAATAGCGTAGTCATCGttgttg                           <  1:3363081/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)

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

                                           aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa                         <  1:2336403/71‑1 (MQ=255)

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

                                            agaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa                         >  1:3119929/1‑70 (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:2700263/71‑1 (MQ=255)

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

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

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

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

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

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

                                               aCTGCATTTGTTTAGCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                             <  1:1304360/63‑1 (MQ=255)

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

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

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

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

                                               aCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                             <  1:7253/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)

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

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

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

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

                                                       tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCa              <  1:3882253/70‑1 (MQ=255)

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

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

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

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

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

                                                           tACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa                         >  1:1184445/1‑55 (MQ=255)

                                                            aCCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGt        >  1:1691041/1‑71 (MQ=255)

                                                             ccAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGt        <  1:2234179/70‑1 (MQ=255)

                                                             ccAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc       <  1:2214462/71‑1 (MQ=255)

                                                              cAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg                    >  1:1200391/1‑57 (MQ=255)

                                                                aCGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGa                  >  1:2496742/1‑57 (MQ=255)

                                                                aCGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCAt     >  1:1839732/1‑70 (MQ=255)

                                                                 cGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAgg         <  1:596952/65‑1 (MQ=255)

                                                                 cGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg   <  1:1121057/71‑1 (MQ=255)

                                                                 cGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg   <  1:3905242/71‑1 (MQ=255)

                                                                  gATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGt  <  1:3827410/71‑1 (MQ=255)

                                                                    tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCa                                    >  1:469559/1‑35 (MQ=255)

GATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGT  >  minE/608913‑609049

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