BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,904	G→C	13.9%	G450G (GGC→GGG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,904	0	G	C	13.9%	272.6 / 30.3	108	G450G (GGC→GGG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base G (51/42); new base C (7/8); total (58/50)
Fisher's exact test for biased strand distribution p-value = 5.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 110 total aligned reads displayed.

ATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTT  >  Exported/608835‑608971

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

       aGGATATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGtt                                                             >  1:4215052/1‑71 (MQ=255)

       aGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGtt                                                             >  1:808107/1‑71 (MQ=255)

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

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

        ggACATAGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGtt                                                             >  1:2870897/1‑70 (MQ=255)

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

            aTCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGa                                                           <  1:2356034/68‑1 (MQ=255)

             tCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGc                                                       <  1:2574888/71‑1 (MQ=255)

             tCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGc                                                       <  1:4092283/71‑1 (MQ=255)

              cGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCg                                                      >  1:2707379/1‑71 (MQ=255)

              cGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCg                                                      >  1:3047628/1‑71 (MQ=255)

                ccTTTGATCGGTTCAGACTTCTGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTa                                                    <  1:2881538/71‑1 (MQ=255)

                ccTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAt                                                          >  1:4518497/1‑65 (MQ=255)

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

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

                  tttGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTa                                                    >  1:1517222/1‑69 (MQ=255)

                     gATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGa                                                           >  1:1843127/1‑59 (MQ=255)

                     gATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGa                                                           >  1:2231155/1‑59 (MQ=255)

                         ggTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAAAGCCGCCGTTGATTGCGTACAGCATGGt                                           <  1:3284944/71‑1 (MQ=255)

                          gTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGtt                                          <  1:65963/71‑1 (MQ=255)

                           ttCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGttt                                         <  1:4326348/71‑1 (MQ=255)

                           ttCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGttt                                         <  1:3358694/71‑1 (MQ=255)

                           ttCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGttt                                         <  1:5431602/71‑1 (MQ=255)

                             cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGtttt                                        >  1:5268579/1‑70 (MQ=255)

                             cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGtttt                                        <  1:2553028/70‑1 (MQ=255)

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

                               gACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGc                                     >  1:1264034/1‑71 (MQ=255)

                               gACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGc                                     >  1:1343391/1‑71 (MQ=255)

                                  ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAAcgc                                                                    >  1:545681/1‑37 (MQ=255)

                                  ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCa                                               <  1:1190732/58‑1 (MQ=255)

                                  ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCa                                               <  1:3575453/58‑1 (MQ=255)

                                  ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCa                                               <  1:4165177/58‑1 (MQ=255)

                                  ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCa                                               <  1:2527513/58‑1 (MQ=255)

                                  ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATg                                             <  1:4109264/60‑1 (MQ=255)

                                  ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATg                                             <  1:4931906/60‑1 (MQ=255)

                                  ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATg                                             <  1:1342268/60‑1 (MQ=255)

                                  ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGtt                                          >  1:2264402/1‑63 (MQ=255)

                                    cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATg                                             >  1:3448639/1‑58 (MQ=255)

                                    cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATg                                             >  1:4370299/1‑58 (MQ=255)

                                    cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGc                                     >  1:2451882/1‑66 (MQ=255)

                                    cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAgg                                 <  1:1832260/70‑1 (MQ=255)

                                    cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAg                                  <  1:808652/69‑1 (MQ=255)

                                    cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGt                                >  1:3961435/1‑71 (MQ=255)

                                    cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGt                                >  1:4920491/1‑71 (MQ=255)

                                         cgaCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcg                           >  1:5607572/3‑70 (MQ=255)

                                         cAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCa                           >  1:2086234/1‑71 (MQ=255)

                                          gaCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgc                          >  1:4024147/2‑69 (MQ=255)

                                           aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGcc                                    >  1:2503781/1‑60 (MQ=255)

                                           aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcg                         <  1:4610772/71‑1 (MQ=255)

                                           aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcg                         <  1:1813025/71‑1 (MQ=255)

                                           aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcg                         <  1:2942887/71‑1 (MQ=255)

                                           aCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgcg                         <  1:62444/71‑1 (MQ=255)

                                             cTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgcg                       <  1:3899536/71‑1 (MQ=255)

                                             cTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgcg                       >  1:2418760/1‑71 (MQ=255)

                                             cTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgcgtg                       <  1:604767/71‑3 (MQ=255)

                                               gCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTa                                                    >  1:2643916/1‑40 (MQ=255)

                                               gCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTa                                                    >  1:2273886/1‑40 (MQ=255)

                                               gCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCa                     >  1:763935/1‑71 (MQ=255)

                                               gCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCa                     >  1:2833316/1‑71 (MQ=255)

                                                cATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCAc                    >  1:664257/1‑71 (MQ=255)

                                                  ttttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCAcc                   >  1:3300909/1‑70 (MQ=255)

                                                  ttttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCAcc                   >  1:3640368/1‑70 (MQ=255)

                                                    ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                 <  1:3972192/38‑1 (MQ=255)

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

                                                     tCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                 <  1:5693198/37‑1 (MQ=255)

                                                     tCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                  >  1:4919861/1‑68 (MQ=255)

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

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

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

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

                                                       aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa             <  1:2485728/71‑1 (MQ=255)

                                                         tttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACt           <  1:3820837/71‑1 (MQ=255)

                                                          ttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGtttt                                        >  1:1781154/1‑41 (MQ=255)

                                                          ttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGtttt                                        >  1:899933/1‑41 (MQ=255)

                                                          ttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg          >  1:4708645/1‑71 (MQ=255)

                                                          ttctCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTg          >  1:5586913/4‑71 (MQ=255)

                                                           tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAg                                  >  1:2433046/1‑46 (MQ=255)

                                                           tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAg                                  >  1:2557281/1‑46 (MQ=255)

                                                           tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAg                                  >  1:42513/1‑46 (MQ=255)

                                                           tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGtt                               >  1:3815029/1‑49 (MQ=255)

                                                           tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg          >  1:2813492/1‑70 (MQ=255)

                                                           tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg          >  1:3031028/1‑70 (MQ=255)

                                                           tttCGTCAACGCCACCGTTGATTGCGTACAGCATGGTTTTCGCCAGGtt                               >  1:1673985/1‑49 (MQ=255)

                                                            ctCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGc         <  1:5301396/69‑1 (MQ=255)

                                                            ctCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGc         <  1:4795152/69‑1 (MQ=255)

                                                             tCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAt       >  1:2057723/1‑71 (MQ=255)

                                                             tCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAt       >  1:1163239/1‑71 (MQ=255)

                                                             tCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa             >  1:3078734/1‑65 (MQ=255)

                                                             tCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa             >  1:4202879/1‑65 (MQ=255)

                                                             tCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa             >  1:4648764/1‑65 (MQ=255)

                                                             tCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCa        <  1:2517823/70‑1 (MQ=255)

                                                              cGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAtt      >  1:462022/1‑71 (MQ=255)

                                                              cGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCAt       <  1:3273245/70‑1 (MQ=255)

                                                                tCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTg    <  1:5268998/71‑1 (MQ=255)

                                                                tCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTg    <  1:1694310/71‑1 (MQ=255)

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

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

                                                                   aCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACTCGCACCGAAGAACTGCATTTGtt  <  1:2118390/70‑1 (MQ=255)

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

ATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTT  >  Exported/608835‑608971

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