BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	4,372,384	A→C	7.6%	I209S (ATC→AGC) ‡	aspA ←	aspartate ammonia‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,372,384	0	A	C	7.6%	270.5 / 17.9	106	I209S (ATC→AGC) ‡	aspA	aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand): ref base A (57/41); new base C (4/4); total (61/45)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.07e-01

Only 100 of 116 total aligned reads displayed.

TGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGA  >  W3110S.gb/4372314‑4372446

tGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGa                                                                >  1:3321856/1‑71 (MQ=255)

tGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGa                                                                >  1:3525264/1‑71 (MQ=255)

 gCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGa                                                                <  1:3357126/70‑1 (MQ=255)

 gCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGa                                                                <  1:5691808/70‑1 (MQ=255)

 gCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGa                                                                <  1:5927738/70‑1 (MQ=255)

 gCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGa                                                                <  1:7743350/70‑1 (MQ=255)

 gCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGAt                                                               >  1:12686/1‑71 (MQ=255)

 gCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGAt                                                               >  1:2080844/1‑71 (MQ=255)

     cAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTg                                                           <  1:7032690/71‑1 (MQ=255)

     cAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTg                                                           <  1:7574707/71‑1 (MQ=255)

      aaGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTg                                                           >  1:3063841/1‑70 (MQ=255)

       aGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGaa                                                         <  1:2334834/71‑1 (MQ=255)

        ggTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGaaa                                                        <  1:7812721/71‑1 (MQ=255)

        ggTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGaaa                                                        <  1:6128985/71‑1 (MQ=255)

         gTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCATCCTGCTGAAAg                                                       <  1:2269192/71‑1 (MQ=39)

         gTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAg                                                       <  1:815366/71‑1 (MQ=255)

           tAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAgcg                                                     >  1:7037643/1‑71 (MQ=255)

            aaCTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAgcgc                                                    >  1:6452254/1‑71 (MQ=255)

            aaCTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAgcgc                                                    >  1:4930888/1‑71 (MQ=255)

             aCTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGcgcg                                                   >  1:4400198/1‑71 (MQ=255)

             aCTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACATCTTCTTTCAGCAGGATGCTGaaa                                                        <  1:6707198/66‑1 (MQ=255)

              cTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAgc                                                      <  1:1898698/67‑1 (MQ=255)

              cTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCgg                                                  <  1:1399948/71‑1 (MQ=255)

                 ccAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCG‑GGAAt                                               <  1:3248142/70‑1 (MQ=255)

                  cagcagcagTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAAt                                               <  1:816556/70‑1 (MQ=255)

                  cagcagcagTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAAt                                               <  1:6926868/70‑1 (MQ=255)

                    gcagcagTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGa                                                          <  1:5802025/57‑1 (MQ=255)

                    gcagcagTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGa                                                          <  1:2190559/57‑1 (MQ=255)

                     cagcagTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCt                                           <  1:1574993/71‑1 (MQ=255)

                      agcagTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGa                                                                >  1:5479446/1‑49 (MQ=255)

                      agcagTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGa                                                                >  1:2516694/1‑49 (MQ=255)

                      agcagTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGa                                                                >  1:4310047/1‑49 (MQ=255)

                      agcagTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGaaa                                                        >  1:5525471/1‑57 (MQ=255)

                      agcagTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGaaa                                                        >  1:7667994/1‑57 (MQ=255)

                      agcagTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGaaa                                                        >  1:6879509/1‑57 (MQ=255)

                      agcagTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTc                                             >  1:6407386/1‑68 (MQ=255)

                      agcagTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCt                                           >  1:509515/1‑70 (MQ=255)

                         agTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGTATGCTGAAAGCGCGGAATTCCTGAcc                                       >  1:7835927/1‑71 (MQ=255)

                         agTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGAcc                                       >  1:3606126/1‑71 (MQ=255)

                         agTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGAcc                                       >  1:2181014/1‑71 (MQ=255)

                         agTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGAcc                                       >  1:6617374/1‑71 (MQ=255)

                           ttCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCATCCTGCTGaa                                                         <  1:3947840/51‑1 (MQ=37)

                           ttCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGAcc                                       >  1:4814568/1‑69 (MQ=255)

                            tCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGAcc                                       >  1:3664653/1‑68 (MQ=255)

                             cAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCATCCTGCTGAAAg                                                       >  1:399184/1‑51 (MQ=25)

                             cAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCATCCTGCTGAAAg                                                       >  1:2826664/1‑51 (MQ=25)

                             cAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAgg                                   >  1:5874773/1‑71 (MQ=255)

                             cAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAgg                                   >  1:3723860/1‑71 (MQ=255)

                             cAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAgg                                   >  1:1316387/1‑71 (MQ=255)

                                cGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCATCCTGCTGaa                                                         <  1:4509031/46‑1 (MQ=25)

                                cGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGa                                                          <  1:4431619/45‑1 (MQ=255)

                                cGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTc                                <  1:2417919/71‑1 (MQ=255)

                                 ggTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAgg                                   <  1:816872/67‑1 (MQ=255)

                                 ggTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCa                               >  1:6596330/1‑71 (MQ=255)

                                  gTACGTTGGATGTTTTTCACTTCTTCTTTCAGCATCCTGCTGAAAg                                                       >  1:5159642/1‑46 (MQ=25)

                                  gTACGTTGGATGTTTTTCACTTCTTCTTTCAGCATCCTGCTGAAAg                                                       >  1:1674472/1‑46 (MQ=25)

                                   tACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGATCGAGGGt                                 >  1:4388545/1‑67 (MQ=255)

                                    aCGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAgg                                   >  1:5898437/1‑64 (MQ=255)

                                     cGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCt                                           >  1:5846187/1‑55 (MQ=255)

                                     cGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCt                                           >  1:2568497/1‑55 (MQ=255)

                                     cGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCt                                           >  1:5362153/1‑55 (MQ=255)

                                     cGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCgg                           <  1:413099/71‑1 (MQ=255)

                                     cGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCgg                           >  1:3081470/1‑71 (MQ=255)

                                      gTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCgg                           >  1:6416359/1‑70 (MQ=255)

                                       ttGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTa                         >  1:1119571/1‑71 (MQ=255)

                                       ttGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTa                         >  1:7212826/1‑71 (MQ=255)

                                        tGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAGTTCCTGACCGAGGGTCATCgg                           >  1:7204313/1‑68 (MQ=255)

                                        tGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTg                                          >  1:1115830/1‑53 (MQ=255)

                                        tGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCgg                           >  1:5544864/1‑68 (MQ=255)

                                         ggATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCt                                           >  1:1341073/1‑51 (MQ=255)

                                         ggATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACt                       >  1:4855067/1‑71 (MQ=255)

                                            tGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCg                    <  1:5887028/71‑1 (MQ=255)

                                             gTTTTTCACTTCTTCTTTCAGCTGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGt                   >  1:4124634/1‑71 (MQ=255)

                                             gTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGt                   >  1:7296642/1‑71 (MQ=255)

                                              tttttCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGGACTGCGt                   <  1:7253049/70‑1 (MQ=255)

                                               ttttCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGt                          <  1:6586656/62‑1 (MQ=255)

                                               ttttCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTcc                 >  1:2933053/1‑71 (MQ=255)

                                                tttCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCa                               <  1:400095/56‑1 (MQ=255)

                                                tttCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCgg                           <  1:3723362/60‑1 (MQ=255)

                                                tttCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCt                >  1:2806819/1‑71 (MQ=255)

                                                  tCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGa                                     <  1:1733692/48‑1 (MQ=255)

                                                  tCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTg               >  1:3691422/1‑70 (MQ=255)

                                                    aCTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAg            <  1:800261/71‑1 (MQ=255)

                                                      ttcttcttTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCt          >  1:1436256/1‑71 (MQ=255)

                                                         ttcttTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTc                  <  1:4328521/60‑1 (MQ=255)

                                                         ttcttTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTc                  <  1:5521118/60‑1 (MQ=255)

                                                         ttcttTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTggg       >  1:7189350/1‑71 (MQ=255)

                                                         ttcttTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTggg       >  1:1700291/1‑71 (MQ=255)

                                                           cttTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTgg        <  1:7561862/68‑1 (MQ=255)

                                                             ttCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTAc    >  1:1796332/1‑70 (MQ=255)

                                                             ttCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTAc    >  1:2869484/1‑70 (MQ=255)

                                                              tCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTAc    <  1:4000533/69‑1 (MQ=255)

                                                              tCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTAc    <  1:7896649/69‑1 (MQ=255)

                                                              tCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACg   <  1:7307676/70‑1 (MQ=255)

                                                              tCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGa  >  1:5146422/1‑71 (MQ=255)

                                                              tCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGa  >  1:7681023/1‑71 (MQ=255)

                                                              tCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGa  >  1:5716055/1‑71 (MQ=255)

                                                               cagcagGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTAc    <  1:4248660/68‑1 (MQ=255)

                                                                 gcagGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCgg                           <  1:6229561/43‑1 (MQ=255)

                                                                    ggATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGc                     <  1:3144998/46‑1 (MQ=255)

TGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGA  >  W3110S.gb/4372314‑4372446

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