BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,991	T→G	15.0%	V421V (GTA→GTC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,991	0	T	G	15.0%	248.0 / 51.5	127	V421V (GTA→GTC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base T (40/68); new base G (10/9); total (50/77)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

Only 100 of 127 total aligned reads displayed.

CAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTG  >  minE/608922‑609059

     tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcag                                                                >  1:2136562/1‑71 (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:3293049/1‑70 (MQ=255)

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

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

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

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

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

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

            cGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcag                                                                <  1:351852/64‑1 (MQ=255)

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

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

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

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

                aGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGc                                                     >  1:1523287/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)

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

                       cACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCa                                              >  1:1887477/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:3810445/71‑1 (MQ=255)

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

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

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

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

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

                                caaagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCgctgtt                                      >  1:2861963/3‑70 (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:3363081/71‑1 (MQ=255)

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

                                 gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                     >  1:2659473/1‑70 (MQ=255)

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

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

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

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

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

                                   agaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa                                   >  1:3119929/1‑70 (MQ=255)

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

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

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

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

                                    gaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGt                                 <  1:3451581/71‑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:1416918/63‑1 (MQ=255)

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

                                       cTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg                              >  1:1524361/1‑71 (MQ=255)

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

                                          cATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCgg                             <  1:785074/69‑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)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                       aCGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGa                            >  1:2558266/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)

                                                         gATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt                                                <  1:3776659/35‑1 (MQ=255)

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

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

                                                             atcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct         >  1:39150/1‑70 (MQ=255)

                                                             atcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct         >  1:2943661/1‑70 (MQ=255)

                                                              tcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt                                         <  1:2213414/37‑1 (MQ=255)

                                                              tcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg                              <  1:1884846/48‑1 (MQ=255)

                                                              tcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGGCATCGTTCTCa       <  1:1346632/71‑1 (MQ=255)

                                                               cGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa       >  1:3156580/1‑70 (MQ=255)

                                                                  gCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACt   <  1:3340462/71‑1 (MQ=255)

                                                                   cTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATa     <  1:884056/68‑1 (MQ=255)

                                                                   cTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTg  >  1:1478820/1‑71 (MQ=255)

                                                                   cTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTAATACTg  <  1:2620261/71‑1 (MQ=255)

CAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTG  >  minE/608922‑609059

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