BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,928	G→C	25.6%	T442T (ACC→ACG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,928	0	G	C	25.6%	214.3 / 77.7	125	T442T (ACC→ACG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base G (53/40); new base C (18/14); total (71/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.28e-01

Only 100 of 130 total aligned reads displayed.

GGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAG  >  Exported/608860‑608997

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)

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

       aCTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGcc                                                              <  1:1821597/71‑1 (MQ=255)

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

         ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATg                                                                       <  1:4109264/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)

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

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

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

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

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

                cgaCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcg                                                     >  1:5607572/3‑70 (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:2942887/71‑1 (MQ=255)

                  aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcg                                                   <  1:1813025/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)

                      gCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCa                                               >  1:763935/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)

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

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

                             cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgcgtg                                                 >  1:5664726/1‑60 (MQ=255)

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

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

                             cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                                         <  1:3515728/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)

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

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

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

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

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

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

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

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

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

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

                                    tCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                       >  1:3078734/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)

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

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

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

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

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

                                           ccccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCa                                  <  1:2393552/61‑1 (MQ=255)

                                           ccccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCa                                  <  1:1155144/61‑1 (MQ=255)

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

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

                                             ccgccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTAc                         >  1:5557805/1‑70 (MQ=255)

                                              cgccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                                         <  1:4665827/53‑1 (MQ=255)

                                              cgccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                                         <  1:1526179/53‑1 (MQ=255)

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

                                              cgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTgct                            >  1:3426410/1‑64 (MQ=255)

                                              cgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTgct                            >  1:1081659/1‑64 (MQ=255)

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

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

                                                ccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGtt                                                         <  1:3790739/35‑1 (MQ=16)

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

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

                                                    taaTTGCTTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCa                                  <  1:438023/52‑1 (MQ=255)

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

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

                                                     aaTTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAt                                 >  1:4177436/2‑54 (MQ=255)

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

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

                                                       ttGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                                         <  1:3908913/44‑1 (MQ=12)

                                                         gCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAcca                                            >  1:4149446/1‑38 (MQ=9)

                                                         gCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaagaa                                       >  1:4326718/1‑44 (MQ=12)

                                                         gCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg            >  1:5340523/1‑71 (MQ=255)

                                                          catACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                       >  1:2913451/3‑43 (MQ=255)

                                                          cGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                       >  1:2067710/1‑43 (MQ=255)

                                                              cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAAc                                      >  1:5458941/1‑40 (MQ=255)

                                                              cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAAc                                      >  1:136434/1‑40 (MQ=255)

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

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

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

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

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

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

                                                                 cATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg     >  1:2583007/1‑70 (MQ=255)

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

                                                                   tGGTTTTCGGCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg     <  1:4912499/68‑1 (MQ=255)

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

GGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAG  >  Exported/608860‑608997

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