BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	609,267	A→G	7.9%	S329S (TCT→TCC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	609,267	0	A	G	7.9%	285.2 / 11.1	101	S329S (TCT→TCC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (56/37); new base G (4/4); total (60/41)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 103 total aligned reads displayed.

GGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCA  >  Exported/609201‑609336

ggTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCagag                                                                    >  1:1969865/1‑70 (MQ=255)

 gtcaCCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGaa                                                                  <  1:4104932/68‑1 (MQ=255)

 gTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAgaga                                                                   <  1:2512797/70‑1 (MQ=255)

    aCCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAg                                                               >  1:845085/1‑71 (MQ=255)

    aCCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAg                                                               <  1:3508966/71‑1 (MQ=255)

      cAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGtt                                                             >  1:4221212/1‑71 (MQ=255)

      cAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGt                                                              >  1:3660695/1‑70 (MQ=255)

      cAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGt                                                              >  1:2647042/1‑70 (MQ=255)

         ggTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCAt                                                          <  1:1293650/71‑1 (MQ=255)

           tACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCg                                                        <  1:3889893/71‑1 (MQ=255)

           tACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCg                                                        <  1:2079290/71‑1 (MQ=255)

               aCCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTAtt                                                    >  1:1698968/1‑71 (MQ=255)

               aCCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTAt                                                     >  1:624818/1‑70 (MQ=255)

                ccGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc                                                   >  1:4124624/1‑71 (MQ=255)

                ccGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc                                                   >  1:572776/1‑71 (MQ=255)

                  gTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg                                                 <  1:1691527/71‑1 (MQ=255)

                   tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACa                                                                <  1:4986670/55‑1 (MQ=255)

                   tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACa                                                                <  1:3899602/55‑1 (MQ=255)

                   tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACa                                                                <  1:3307660/55‑1 (MQ=255)

                   tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg                                                 >  1:4205076/1‑70 (MQ=255)

                      aGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAg                                              >  1:4514107/1‑70 (MQ=255)

                          ccATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTGCGGAGTACg                                          >  1:3484677/1‑70 (MQ=255)

                           cATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAAc                                                                 <  1:4291754/46‑1 (MQ=255)

                           cATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAAc                                                                 <  1:2868013/46‑1 (MQ=255)

                            aTACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAAc                                                                 <  1:4092524/45‑1 (MQ=255)

                            aTACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAg                                       >  1:4629968/1‑71 (MQ=255)

                             tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAAc                                                                 >  1:2804630/1‑44 (MQ=255)

                             tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAg                                                               >  1:797551/1‑46 (MQ=255)

                             tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAgg                                      <  1:1712763/71‑1 (MQ=255)

                             tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAgg                                      <  1:3820930/71‑1 (MQ=255)

                               ccaccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGa                                     <  1:821195/70‑1 (MQ=255)

                               ccaccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGa                                     <  1:3361070/70‑1 (MQ=255)

                                 accGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg                                                 >  1:378396/1‑56 (MQ=255)

                                  ccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc                                                   >  1:5616054/1‑53 (MQ=255)

                                   cGATAGATTCGGTTGCCCAGATCGGGTCGCCAgaga                                                                   >  1:2719212/1‑36 (MQ=255)

                                   cGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGaa                                                                  >  1:4196245/1‑37 (MQ=255)

                                   cGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCg                                 <  1:1657633/70‑1 (MQ=255)

                                   cGATAGATTCGGTTGCCCAGATCGGGTAGCCAGAGaa                                                                  >  1:2396641/1‑37 (MQ=255)

                                    gATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATc                                                         >  1:503729/1‑45 (MQ=255)

                                     aTAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTc                                                            <  1:5490046/41‑1 (MQ=255)

                                     aTAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAc                              >  1:4952182/1‑71 (MQ=255)

                                      tAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc                                                   <  1:4114301/49‑1 (MQ=255)

                                      tAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGa                                <  1:2988351/68‑1 (MQ=255)

                                       aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGc                                  >  1:1460356/1‑65 (MQ=255)

                                       aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                             >  1:5282517/1‑70 (MQ=255)

                                       aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAc                              >  1:4756260/1‑69 (MQ=255)

                                       aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCa                            >  1:399609/1‑71 (MQ=255)

                                        gATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCa                            <  1:4702818/70‑1 (MQ=255)

                                        acttCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCa                            >  1:5649674/3‑70 (MQ=255)

                                            cGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACg                        >  1:4007674/1‑70 (MQ=255)

                                                 gCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGtt                   >  1:1803539/1‑70 (MQ=255)

                                                  cccAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTc                 <  1:3545860/71‑1 (MQ=255)

                                                  cccAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTc                 <  1:5186223/71‑1 (MQ=255)

                                                   ccAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCa                >  1:3512965/1‑71 (MQ=255)

                                                    cAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCAt               <  1:3130415/71‑1 (MQ=255)

                                                    cAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCAt               <  1:2266684/71‑1 (MQ=255)

                                                    cAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCg                                 <  1:3613448/53‑1 (MQ=255)

                                                    cAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCg                                 <  1:42153/53‑1 (MQ=255)

                                                    cAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCg                                 <  1:4897229/53‑1 (MQ=255)

                                                     aGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCAt               >  1:3455217/1‑70 (MQ=255)

                                                     aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt               >  1:4891136/1‑70 (MQ=255)

                                                     aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt               >  1:3794498/1‑70 (MQ=255)

                                                     aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt               >  1:123839/1‑70 (MQ=255)

                                                     aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt               >  1:2193704/1‑70 (MQ=255)

                                                     aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt               >  1:3589419/1‑70 (MQ=255)

                                                     aGATCGGGTCGCCAGAGAACAGTTCATCGTATTACGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt               >  1:4897799/1‑70 (MQ=255)

                                                      gATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGa             >  1:3407185/1‑71 (MQ=255)

                                                      gATCGGGTCGCCAGAGAACAGTTCATCGTATTTCGGAGTACGCAGGAAGCGa                                >  1:190192/1‑52 (MQ=255)

                                                      gATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGc                                         <  1:2008814/43‑1 (MQ=255)

                                                      gATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGaa                               >  1:5165173/1‑53 (MQ=255)

                                                      gATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGaa                               >  1:4145257/1‑53 (MQ=255)

                                                      gATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGaa                               >  1:2798119/1‑53 (MQ=255)

                                                       aTCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGa                                               >  1:440374/1‑36 (MQ=255)

                                                       aTCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGa                                               >  1:5447180/1‑36 (MQ=255)

                                                       aTCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAc            <  1:3607582/71‑1 (MQ=255)

                                                        ttgggTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGa             >  1:1024587/3‑69 (MQ=255)

                                                        tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCa                                        >  1:352507/1‑42 (MQ=255)

                                                        tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGa                                >  1:1816350/1‑50 (MQ=255)

                                                        tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATa                          >  1:23699/1‑56 (MQ=255)

                                                        tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATa                          >  1:3601350/1‑56 (MQ=255)

                                                        tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc           <  1:3888221/71‑1 (MQ=255)

                                                         cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                             >  1:5177806/1‑52 (MQ=255)

                                                         cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                             >  1:5637492/1‑52 (MQ=255)

                                                         cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                             >  1:2447193/1‑52 (MQ=255)

                                                         cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                             >  1:3362506/1‑52 (MQ=255)

                                                         cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                  >  1:722803/1‑63 (MQ=255)

                                                         cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCa                >  1:5626028/1‑65 (MQ=255)

                                                         cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc           <  1:2468755/70‑1 (MQ=255)

                                                         cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc           <  1:5527717/70‑1 (MQ=255)

                                                         cGGGTCGCCAGAGAACAGTTCATCGTATTACGGAGTACGCAGGAAGCg                                 >  1:5453611/1‑48 (MQ=255)

                                                           ggTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCa                            <  1:2814762/51‑1 (MQ=255)

                                                              cGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtg      >  1:4023666/1‑70 (MQ=255)

                                                               gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCg                                 >  1:3567643/1‑42 (MQ=255)

                                                               gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAc            >  1:4516847/1‑63 (MQ=255)

                                                               gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtggt    >  1:5303577/1‑71 (MQ=255)

                                                               gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtggt    >  1:3271539/1‑71 (MQ=255)

                                                                 cAGAGAACAGTTCATCGTATTCCTGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCa  <  1:5031927/71‑1 (MQ=255)

                                                                 cAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGt                    <  1:228596/53‑1 (MQ=255)

                                                                 cAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCa  <  1:278814/71‑1 (MQ=255)

                                                                  agagAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAg         >  1:1972603/1‑63 (MQ=255)

GGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCA  >  Exported/609201‑609336

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