BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AP009048	4,372,381	G→C	7.0%	I209M (ATC→ATG) ‡	aspA ←	aspartate ammonia‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AP009048	4,372,381	0	G	C	7.0%	338.9 / 10.9	114	I209M (ATC→ATG) ‡	aspA	aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand): ref base G (60/45); new base C (4/4); total (65/49)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.83e-01

Only 100 of 120 total aligned reads displayed.

TGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGA  >  AP009048/4372312‑4372444

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

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

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

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

 gCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGa                                                                <  1:7152575/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)

    ccAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCt                                                            <  1:1535195/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)

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

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

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

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

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

             aCTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGcgcg                                                   >  1:3594286/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)

              cTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGcgcg                                                   >  1:4283487/1‑70 (MQ=255)

                 ccAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAgg                                                                 >  1:231384/1‑53 (MQ=255)

                 ccAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAgg                                                                 >  1:4708766/1‑53 (MQ=255)

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

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

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

                   agcagcagTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAAt                                               <  1:5792350/69‑1 (MQ=255)

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

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

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

                      agcagTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGa                                                                >  1:5479446/1‑49 (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:2181014/1‑71 (MQ=255)

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

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

                           ttCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAgg                                                                 <  1:7931557/43‑1 (MQ=255)

                           ttCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAgg                                                                 <  1:4151386/43‑1 (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:1974921/46‑1 (MQ=25)

                                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:5362153/1‑55 (MQ=255)

                                     cGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCt                                           >  1:2568497/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:7819872/1‑71 (MQ=255)

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

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

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

                                       ttGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTa                         <  1:3990543/71‑1 (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)

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

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

                                                tttCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCt                >  1:2806819/1‑71 (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:5521118/60‑1 (MQ=255)

                                                         ttcttTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTc                  <  1:4328521/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:2869484/1‑70 (MQ=255)

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

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

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

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

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

                                                              tCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGa  >  1:7681023/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)

                                                                    ggATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTgg        >  1:5957336/1‑59 (MQ=255)

TGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGA  >  AP009048/4372312‑4372444

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