BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,979	G→C	14.8%	I425M (ATC→ATG) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,979	0	G	C	14.8%	316.1 / 40.5	128	I425M (ATC→ATG) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base G (46/63); new base C (10/9); total (56/72)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.64e-01

Only 100 of 131 total aligned reads displayed.

GTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGT  >  minE/608910‑609049

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)

        cGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCggg                                                               >  1:1808979/1‑71 (MQ=255)

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

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

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

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

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

                 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)

                    ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcagca                                                    <  1:453039/70‑1 (MQ=255)

                    ttttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagcag                                                   >  1:3433081/1‑70 (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:2448292/1‑47 (MQ=255)

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

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

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

                                 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)

                                    aCGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                               <  1:430575/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)

                                        gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt                               <  1:2275509/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)

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

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

                                              aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACg                                                         >  1:3085477/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)

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

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

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

                                                gaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGt                       >  1:2804005/1‑71 (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:264776/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:2532501/63‑1 (MQ=255)

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

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

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

                                                         ttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGa                  >  1:3136309/1‑67 (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:1340539/52‑1 (MQ=255)

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

                                                            tttACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATc            <  1:2813371/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: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)

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

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

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

GTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGT  >  minE/608910‑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: `‑`