BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	609,000	A→G	26.0%	A418A (GCT→GCC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	609,000	0	A	G	26.0%	151.1 / 57.0	96	A418A (GCT→GCC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (31/40); new base G (14/11); total (45/51)
Fisher's exact test for biased strand distribution p-value = 3.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 141 total aligned reads displayed.

TTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAG  >  Exported/608930‑609069

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

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

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

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

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

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

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

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

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

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

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

        aGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGc                                                               <  1:4923163/71‑1 (MQ=255)

         ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                                                     >  1:4642752/1‑64 (MQ=255)

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

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

                   cgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgt                                                    <  1:1980355/71‑1 (MQ=255)

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

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

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

                     cACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc                                                               <  1:3759162/58‑1 (MQ=255)

                     cACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc                                                               <  1:4355125/58‑1 (MQ=255)

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

                     cACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAAtc                                                                 <  1:4774095/56‑2 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                          aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                               >  1:4654861/1‑69 (MQ=255)

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

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

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

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

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

                             aaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                                                     <  1:5574477/44‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                     ttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg                                        <  1:659474/65‑1 (MQ=255)

                                     ttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg                                        <  1:5176578/65‑1 (MQ=255)

                                      tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAt                                 >  1:268613/1‑71 (MQ=255)

                                       gTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGa                                      >  1:5244188/1‑65 (MQ=255)

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

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

                                          tACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAg                              <  1:1626106/70‑1 (MQ=255)

                                           aCCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAgg                             <  1:3253912/70‑1 (MQ=255)

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

                                              aaCGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCAt                         <  1:2096741/71‑1 (MQ=255)

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

                                                cGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg                       >  1:822877/1‑71 (MQ=255)

                                                 gATCATCGGTCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGt                      >  1:3346541/1‑71 (MQ=255)

                                                  accattGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCAtt                     <  1:2958484/69‑1 (MQ=255)

                                                   tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                           >  1:1007848/1‑64 (MQ=255)

                                                   tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                           >  1:1557215/1‑64 (MQ=255)

                                                   tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                           >  1:4006540/1‑64 (MQ=255)

                                                   tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtc                    <  1:3016969/71‑1 (MQ=255)

                                                    catcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                   >  1:5542797/1‑71 (MQ=255)

                                                    catcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                   >  1:1482523/1‑71 (MQ=255)

                                                     atcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAg                              <  1:2042020/59‑1 (MQ=255)

                                                     atcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                   <  1:439156/70‑1 (MQ=255)

                                                      ttgggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAg                              <  1:5126323/56‑1 (MQ=255)

                                                      ttgggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCa                 >  1:5684195/3‑71 (MQ=255)

                                                      ttgggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCa                 >  1:3290511/3‑71 (MQ=255)

                                                      ttgggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCa                 >  1:4556021/3‑71 (MQ=255)

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

                                                        gggCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATa               >  1:230515/1‑71 (MQ=255)

                                                        gggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTca                                         >  1:1596959/1‑44 (MQ=255)

                                                            ttACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTg            <  1:2953341/70‑1 (MQ=255)

                                                             tACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTg            >  1:683461/1‑69 (MQ=255)

                                                              aCGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGt                                           <  1:3564746/37‑1 (MQ=255)

                                                              aCGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGt                                           >  1:3345073/1‑37 (MQ=255)

                                                              aCGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTg            >  1:4535093/1‑68 (MQ=255)

                                                                gcagcaGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCaagg       <  1:3895082/71‑4 (MQ=255)

                                                                 cagcaAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAgaga      <  1:2118396/71‑1 (MQ=255)

                                                                   gcaggcaaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAggaag    <  1:2366258/71‑1 (MQ=255)

                                                                   gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                 <  1:553177/58‑1 (MQ=255)

                                                                   gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                 <  1:1710426/58‑1 (MQ=255)

                                                                   gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag    <  1:2007039/71‑1 (MQ=255)

                                                                    caagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga        <  1:312338/66‑1 (MQ=255)

                                                                     aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGa                                      >  1:1423125/1‑35 (MQ=255)

                                                                     aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGa                                      >  1:3979835/1‑35 (MQ=255)

                                                                     aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGa                                      >  1:4328136/1‑35 (MQ=255)

                                                                     aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa     >  1:3710429/1‑68 (MQ=255)

                                                                     aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa     >  1:1216229/1‑68 (MQ=255)

                                                                     aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAagag  <  1:316184/71‑1 (MQ=255)

                                                                     aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAagag  >  1:3078264/1‑71 (MQ=255)

TTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAG  >  Exported/608930‑609069

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