BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,991	0	T	G	24.6%	203.1 / 98.4	142	V421V (GTA→GTC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base T (48/59); new base G (20/15); total (68/74)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 146 total aligned reads displayed.

CAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTG  >  Exported/608922‑609059

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

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

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

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

   cATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                                                                     <  1:4089256/68‑1 (MQ=255)

   cATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                                                                     <  1:5202882/68‑1 (MQ=255)

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

     tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACAATCATCGGGCTTACg                                                                   <  1:508464/68‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                  gTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGt                                                   <  1:2925223/71‑1 (MQ=255)

                      gCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                                                                     <  1:3726161/49‑1 (MQ=255)

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

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

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

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

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

                            gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgt                                          >  1:140475/1‑70 (MQ=255)

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

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

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

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

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

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

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

                               ccGAAGAACTGCATTTGTTTACCAACGATCATCGGCCTTACgcag                                                                <  1:771054/45‑1 (MQ=255)

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

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

                                cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgca                                                                 <  1:3427250/43‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                  aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                     >  1:4654861/1‑69 (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)

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

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

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

                                     aaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTc                                <  1:1569837/71‑1 (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:762743/53‑1 (MQ=255)

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

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

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

                                           aTTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACg                          <  1:1878427/71‑1 (MQ=255)

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

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

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

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

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

                                                   aCCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAgg                   <  1:3253912/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)

                                                           tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                 >  1:1007848/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)

                                                            catcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct         >  1:4750025/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:3290511/3‑71 (MQ=255)

                                                              ttgggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCa       >  1:5684195/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:642556/1‑44 (MQ=255)

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

CAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTG  >  Exported/608922‑609059

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