BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,892	T→C	9.4%	K454K (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,892	0	T	C	9.4%	264.7 / 16.6	96	K454K (AAA→AAG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base T (35/52); new base C (4/5); total (39/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

Only 100 of 108 total aligned reads displayed.

CATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACT  >  minE/608823‑608962

cATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGtt                                                                       >  1:1882162/1‑71 (MQ=255)

cATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGtt                                                                       <  1:1769378/71‑1 (MQ=255)

 aTCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGttt                                                                      <  1:2184695/71‑1 (MQ=255)

 aTCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGttt                                                                      <  1:3865584/71‑1 (MQ=255)

 aTCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGttt                                                                      <  1:1845000/71‑1 (MQ=255)

  tCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGtttt                                                                     <  1:1003753/71‑1 (MQ=255)

   cACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGttttt                                                                    >  1:1600743/1‑71 (MQ=255)

    aCTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGttttt                                                                    >  1:704688/1‑70 (MQ=255)

    aCTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTc                                                                   >  1:1941100/1‑71 (MQ=255)

    aCTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTc                                                                   >  1:2894769/1‑71 (MQ=255)

    aCTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTc                                                                   >  1:3898571/1‑71 (MQ=255)

     cTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCg                                                                  <  1:1599451/71‑1 (MQ=255)

     cTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCg                                                                  <  1:1325117/71‑1 (MQ=255)

        catcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGt                                                                 <  1:2748844/69‑1 (MQ=255)

        catcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTc                                                                <  1:654914/70‑1 (MQ=255)

        catcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTc                                                                >  1:2463515/1‑70 (MQ=255)

        catcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTc                                                                <  1:3666135/70‑1 (MQ=255)

        catcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTc                                                                <  1:600650/70‑1 (MQ=255)

        catcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTc                                                                >  1:1806489/1‑70 (MQ=255)

            atAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGttttt                                                                    >  1:2886941/1‑62 (MQ=255)

            atAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAAcgc                                                           >  1:1091706/1‑71 (MQ=255)

                 tCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACgccgcc                                                       <  1:1676283/70‑1 (MQ=255)

                  cAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTc                                                                   >  1:3849779/1‑57 (MQ=255)

                  cAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGt                                                     >  1:3295514/1‑71 (MQ=255)

                  cAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGt                                                     >  1:3908735/1‑71 (MQ=255)

                   aGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACgccgcc                                                       <  1:2494603/68‑1 (MQ=255)

                   aGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGtt                                                    <  1:368580/71‑1 (MQ=255)

                   aGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGtt                                                    <  1:1245158/71‑1 (MQ=255)

                    ggACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTg                                                   <  1:3280851/71‑1 (MQ=255)

                    ggACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTg                                                   <  1:355522/71‑1 (MQ=255)

                     gACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACgccgcc                                                       >  1:2014325/1‑66 (MQ=255)

                      aCATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAt                                                 >  1:3088825/1‑71 (MQ=255)

                       cATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAtt                                                >  1:3686774/1‑71 (MQ=255)

                         tCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCg                                                                  <  1:2366769/51‑1 (MQ=255)

                         tCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTg                                                   <  1:1961296/66‑1 (MQ=255)

                           gCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTc                                                                <  1:498474/51‑1 (MQ=255)

                           gCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTc                                                                <  1:3384560/51‑1 (MQ=255)

                           gCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAAc                                                             >  1:2949624/1‑54 (MQ=255)

                           gCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCg                                             <  1:2592926/70‑1 (MQ=255)

                           gCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCg                                             >  1:3643924/1‑70 (MQ=255)

                            ccTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTa                                           <  1:3931434/71‑1 (MQ=255)

                             cTTTGATCGGTTCAGACTTCGGACCATCCTGCATTTTCAGTTTTTCGTCAACgccgcc                                                       <  1:2740702/58‑1 (MQ=255)

                             cTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACgccgcc                                                       <  1:2102384/58‑1 (MQ=255)

                             cTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACgccgcc                                                       <  1:3753135/58‑1 (MQ=255)

                             cTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTAc                                          <  1:1548144/71‑1 (MQ=255)

                             cTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTAc                                          <  1:255345/71‑1 (MQ=255)

                              tttGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCaa                                                              <  1:1231372/50‑1 (MQ=255)

                                tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                        >  1:344849/1‑70 (MQ=255)

                                tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                        >  1:271456/1‑70 (MQ=255)

                                tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                        >  1:447005/1‑70 (MQ=255)

                                tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                        >  1:3777530/1‑70 (MQ=255)

                                  aTCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCAt                                     <  1:3374726/71‑1 (MQ=255)

                                   tCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATg                                    >  1:3559816/1‑71 (MQ=255)

                                    cGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATgg                                   >  1:3071503/1‑71 (MQ=255)

                                     ggTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGtt                                                    >  1:1746406/1‑53 (MQ=255)

                                     ggTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGt                                  <  1:1785360/71‑1 (MQ=255)

                                      gTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAAcgccgc                                                        <  1:4080758/48‑1 (MQ=255)

                                       ttCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCg                                                                  <  1:3979470/37‑1 (MQ=255)

                                       ttCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCg                                                                  <  1:1304130/37‑1 (MQ=255)

                                        tCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGccgccg                                                      >  1:3150346/1‑48 (MQ=255)

                                         cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTc                              >  1:962584/1‑71 (MQ=255)

                                         cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTc                              >  1:508069/1‑71 (MQ=255)

                                          aGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCg                                             <  1:2526462/55‑1 (MQ=255)

                                             cTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATg                                    >  1:1254041/1‑61 (MQ=255)

                                              ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAAc                                                             <  1:2367990/35‑1 (MQ=255)

                                              ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTg                                                   <  1:2825671/45‑1 (MQ=255)

                                               tCGGCCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAgg                        >  1:1510627/1‑71 (MQ=255)

                                               tCGGACCAACCTGCATTTTCAGTTTTTCGTCAACTCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAgg                        >  1:2552676/1‑71 (MQ=255)

                                                   aCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                    <  1:2792016/71‑1 (MQ=255)

                                                   aCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                    <  1:1038274/71‑1 (MQ=255)

                                                   aCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                    <  1:1971130/71‑1 (MQ=255)

                                                   aCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                    <  1:2491041/71‑1 (MQ=255)

                                                    ccgaCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAgg                        >  1:3826218/4‑66 (MQ=255)

                                                     cgaCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGtta                     >  1:1234575/3‑67 (MQ=255)

                                                      aaCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                 <  1:3040962/71‑1 (MQ=255)

                                                      aaCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                 <  1:1016534/71‑1 (MQ=255)

                                                      aaCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                 <  1:2910320/71‑1 (MQ=255)

                                                      aaCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                 <  1:2681248/71‑1 (MQ=255)

                                                       aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcg                >  1:1329070/1‑71 (MQ=255)

                                                        ccTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgc               >  1:2043102/1‑71 (MQ=255)

                                                          tGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                        <  1:1021816/44‑1 (MQ=255)

                                                          tGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                        <  1:2522047/44‑1 (MQ=255)

                                                          tGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                        <  1:1945260/44‑1 (MQ=255)

                                                          tGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgcg              >  1:24167/1‑70 (MQ=255)

                                                            cATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACgcgc             <  1:1895898/69‑1 (MQ=255)

                                                                ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg         >  1:3654931/1‑69 (MQ=255)

                                                                ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg         >  1:2381497/1‑69 (MQ=255)

                                                                ttCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAcca         >  1:1494808/1‑68 (MQ=255)

                                                                  cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag      <  1:27440/70‑1 (MQ=255)

                                                                  cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag      <  1:102482/70‑1 (MQ=255)

                                                                  cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAccaa        <  1:1916967/68‑3 (MQ=255)

                                                                  cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAccaa        <  1:3865116/68‑3 (MQ=255)

                                                                  cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAcca         <  1:3618961/67‑2 (MQ=255)

                                                                  cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag      <  1:1703505/70‑1 (MQ=255)

                                                                  cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag      <  1:1401060/70‑1 (MQ=255)

                                                                   aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa    >  1:595649/1‑71 (MQ=255)

                                                                   aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaaga     <  1:2728383/70‑1 (MQ=255)

                                                                   aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaaga     <  1:3160660/70‑1 (MQ=255)

                                                                    gTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAAc   <  1:2532648/71‑1 (MQ=255)

                                                                     tttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGAACGCGCACCGAAGAACt  >  1:2008479/1‑71 (MQ=255)

CATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACT  >  minE/608823‑608962

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