BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,970	0	T	C	23.3%	292.6 / 87.4	159	K428K (AAA→AAG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base T (47/75); new base C (21/16); total (68/91)
Fisher's exact test for biased strand distribution p-value = 5.89e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 166 total aligned reads displayed.

AACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAT  >  Exported/608901‑609038

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

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

 aCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                                                     >  1:3459613/1‑70 (MQ=255)

   cccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAc                                                                  >  1:3756992/2‑71 (MQ=255)

     cgccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCa                                                                <  1:2065457/71‑1 (MQ=255)

     cgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAcc                                                                 >  1:5186119/1‑70 (MQ=255)

       ccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTa                                                                   >  1:5516307/1‑66 (MQ=255)

       ccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCaat                                                              <  1:2375314/71‑2 (MQ=255)

        cgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACg                                                             <  1:2771617/71‑1 (MQ=255)

        cgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACg                                                             <  1:775996/71‑1 (MQ=255)

          ttGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGat                                                           >  1:4718237/1‑71 (MQ=255)

            gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                                          <  1:2818379/70‑1 (MQ=255)

            gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                                          <  1:4810051/70‑1 (MQ=255)

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

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

               tGCGTACAGCATGGTTTTCCCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTa                                                                   >  1:3366110/1‑58 (MQ=255)

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

                     cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                                                >  1:2240311/1‑71 (MQ=255)

                     cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                                                >  1:5572364/1‑71 (MQ=255)

                      aGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTAc                                               <  1:196137/71‑1 (MQ=255)

                         aTGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                                                >  1:3132667/1‑67 (MQ=255)

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

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

                            gTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcat                                                        <  1:2207997/56‑1 (MQ=255)

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

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

                             ttttCGCCAGGTTTGCACGCGCACCGATGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                                                <  1:4923260/63‑1 (MQ=255)

                             ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                                         <  1:399597/54‑1 (MQ=255)

                             ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                                         <  1:576990/54‑1 (MQ=255)

                             ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCtt                                                 <  1:1943784/62‑1 (MQ=255)

                             ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaa                                        <  1:5109958/71‑1 (MQ=255)

                              tttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                                          >  1:5134384/1‑52 (MQ=255)

                              tttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaag                                       >  1:5591309/1‑71 (MQ=255)

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

                               ttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagcagg                                       >  1:5184458/1‑68 (MQ=255)

                               ttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagcagg                                       >  1:32789/1‑68 (MQ=255)

                               ttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagcagg                                       >  1:5207214/1‑68 (MQ=255)

                               ttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagcagg                                       >  1:5332640/1‑68 (MQ=255)

                                tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt                                                  <  1:5587630/58‑1 (MQ=255)

                                tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                      <  1:2032018/70‑1 (MQ=255)

                                tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                      <  1:3043363/70‑1 (MQ=255)

                                tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                      <  1:5274704/70‑1 (MQ=255)

                                tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                      <  1:5547664/70‑1 (MQ=255)

                                tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                      <  1:202592/70‑1 (MQ=255)

                                tCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCagcaggc                                      <  1:4334380/70‑1 (MQ=255)

                                 cGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagca                                     <  1:2937922/70‑1 (MQ=255)

                                  gCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                    <  1:4224293/70‑1 (MQ=255)

                                  gCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                    <  1:4653418/70‑1 (MQ=255)

                                    cAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCaatc                                                             <  1:1905524/43‑3 (MQ=19)

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

                                      ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                                          <  1:4240566/44‑1 (MQ=255)

                                      ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                                          >  1:1341037/1‑44 (MQ=255)

                                       gTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt                                                  <  1:4772410/51‑1 (MQ=255)

                                        tttGCACGCGCACCGAAGAACTGCATTTGTTTACCAAc                                                              <  1:1089070/38‑1 (MQ=255)

                                        tttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcat                                                        <  1:1834487/44‑1 (MQ=255)

                                           gCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                                          <  1:2935523/39‑1 (MQ=255)

                                           gCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                                     <  1:2511788/44‑1 (MQ=255)

                                           gCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                                     <  1:4919366/44‑1 (MQ=255)

                                            cACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                      <  1:2138586/58‑1 (MQ=255)

                                             aCGCGCACCGAAGAACTGCATTTGTTTACCAACGatcat                                                        >  1:4914727/1‑39 (MQ=255)

                                             aCGCGCACCGAAGAACTGCATTTGTTTACCAACGatcat                                                        >  1:1657053/1‑39 (MQ=255)

                                             aCGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                                         <  1:252889/38‑1 (MQ=255)

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

                                                 gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt                    >  1:3664825/1‑71 (MQ=255)

                                                 gcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAgcaggcaa                                    <  1:4796127/55‑1 (MQ=255)

                                                  cACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttg                   <  1:2154030/71‑1 (MQ=255)

                                                  cACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCgctg                   >  1:3185884/1‑71 (MQ=255)

                                                   aCCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagc                                          <  1:5172893/47‑1 (MQ=16)

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

                                                     caaagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGctgttg                <  1:5360228/69‑1 (MQ=255)

                                                     caaagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGctgttg                <  1:3782391/69‑1 (MQ=255)

                                                     cGAAGATCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                <  1:2340665/71‑1 (MQ=255)

                                                     cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                               <  1:192719/56‑1 (MQ=255)

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

                                                      gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACgcagca                                         >  1:4493100/1‑45 (MQ=255)

                                                      gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACgcagca                                         >  1:2612485/1‑45 (MQ=255)

                                                      aaagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGctgttg                >  1:4544652/2‑70 (MQ=255)

                                                      aaagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGctgttg                >  1:633249/2‑70 (MQ=255)

                                                       aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATa                                  >  1:1469284/1‑51 (MQ=255)

                                                       aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATa                                  >  1:4730602/1‑51 (MQ=255)

                                                       aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATa                                  >  1:4163056/1‑51 (MQ=255)

                                                       aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc                                <  1:3781532/53‑1 (MQ=255)

                                                       aagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACgcagc                                          >  1:1492323/1‑43 (MQ=16)

                                                       aagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAAtc                                  >  1:4440321/1‑50 (MQ=255)

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

                                                        agaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCagcaggc                                      <  1:3828394/46‑1 (MQ=16)

                                                        agaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCagcaggc                                      <  1:3829407/46‑1 (MQ=16)

                                                         gaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGt            >  1:24494/1‑71 (MQ=255)

                                                          aaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtcgc                     <  1:3580464/61‑2 (MQ=255)

                                                          aaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTc           <  1:1569837/71‑1 (MQ=255)

                                                            cTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                    >  1:5648645/1‑44 (MQ=255)

                                                            cTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                <  1:2567882/64‑1 (MQ=255)

                                                            cTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa              >  1:5081927/1‑66 (MQ=255)

                                                             tGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTc                          <  1:762743/53‑1 (MQ=255)

                                                             tGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTc                          <  1:2918683/53‑1 (MQ=255)

                                                              gCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAAt                                   >  1:1195933/1‑43 (MQ=255)

                                                               cATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                    >  1:4831019/1‑41 (MQ=255)

                                                                  ttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTc                          <  1:3881398/48‑1 (MQ=255)

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

                                                                     tttACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAt  >  1:4548540/1‑69 (MQ=255)

AACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAT  >  Exported/608901‑609038

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`