BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,201,009	G→T	7.4%	S54S (TCG→TCT)	sdaA →	L‑serine deaminase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,201,009	0	G	T	7.4%	347.4 / 7.2	107	S54S (TCG→TCT)	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand): ref base G (40/59); new base T (5/3); total (45/62)
Fisher's exact test for biased strand distribution p-value = 2.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

ATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCAC  >  Exported/1200939‑1201077

atgatCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCg                                                                      <  1:1894508/71‑1 (MQ=255)

atgatCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCg                                                                      >  1:2254368/1‑71 (MQ=255)

atgatCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCg                                                                      <  1:1719046/71‑1 (MQ=255)

   atCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCt                                                                    <  1:1738992/70‑1 (MQ=255)

   atCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCt                                                                    >  1:2479764/1‑70 (MQ=255)

      tGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCGCTGTCGCTGACg                                                                >  1:5159125/1‑71 (MQ=255)

      tGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACg                                                                >  1:139326/1‑71 (MQ=255)

        gTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGAc                                                                 <  1:4536826/68‑1 (MQ=255)

        gTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGAc                                                                 <  1:5312343/68‑1 (MQ=255)

         tCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGt                                                             <  1:1062648/71‑1 (MQ=255)

         tCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGt                                                             <  1:138073/71‑1 (MQ=255)

          cGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGt                                                             >  1:1280185/1‑70 (MQ=255)

          cGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGt                                                             >  1:7122/1‑70 (MQ=255)

           gAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTaa                                                           <  1:3658153/71‑1 (MQ=255)

           gAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTa                                                            <  1:3107155/70‑1 (MQ=255)

           gAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTa                                                            <  1:1808172/70‑1 (MQ=255)

           gAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTa                                                            <  1:3705772/70‑1 (MQ=255)

            aaaaaGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTaaa                                                          <  1:3949212/71‑1 (MQ=255)

              aaaGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAgg                                                        >  1:3083624/1‑71 (MQ=255)

              aaaGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAgg                                                        >  1:4885595/1‑71 (MQ=255)

              aaaGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAgg                                                        >  1:1916355/1‑71 (MQ=255)

               aaGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGc                                                       <  1:3214937/71‑1 (MQ=255)

                  gCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGccac                                                    <  1:4158358/71‑1 (MQ=255)

                  gCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGccac                                                    <  1:1785162/71‑1 (MQ=255)

                  gCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGccac                                                    <  1:4456703/71‑1 (MQ=255)

                  gCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGccac                                                    <  1:4463620/71‑1 (MQ=255)

                         ggATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCt                                                                    <  1:773802/48‑1 (MQ=255)

                         ggATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACAcc                                              >  1:4572350/1‑70 (MQ=255)

                          gATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACTGGTAAAGGCCACCACAcc                                              <  1:5236837/69‑1 (MQ=255)

                          gATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACAcc                                              <  1:3455469/69‑1 (MQ=255)

                          gATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGa                                            >  1:5047711/1‑71 (MQ=255)

                           aTAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTa                                                            <  1:3828588/54‑1 (MQ=255)

                           aTAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTa                                                            <  1:5176416/54‑1 (MQ=255)

                           aTAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTa                                                            <  1:2541638/54‑1 (MQ=255)

                              acGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATc                                        >  1:2310198/2‑71 (MQ=255)

                                  tACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCAGACGGGTAAAGGCCACCACACCGATATCGCCa                                    >  1:460753/1‑71 (MQ=255)

                                   aCTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGatat                                         >  1:780093/1‑65 (MQ=255)

                                   aCTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGatat                                         >  1:4445484/1‑65 (MQ=255)

                                   aCTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCat                                   >  1:4016157/1‑71 (MQ=255)

                                   aCTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCat                                   >  1:3678227/1‑71 (MQ=255)

                                   aCTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCat                                   >  1:3298186/1‑71 (MQ=255)

                                    cTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCat                                   <  1:2910008/70‑1 (MQ=255)

                                     tCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGat                                           <  1:5104891/61‑1 (MQ=255)

                                     tCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGat                                           <  1:1827148/61‑1 (MQ=255)

                                     tCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGat                                           <  1:2844025/61‑1 (MQ=255)

                                     tCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGat                                           <  1:5310405/61‑1 (MQ=255)

                                       gcgTTGCCGTGGACGTTTATGGTTCACTTTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCattat                                >  1:2162360/1‑70 (MQ=255)

                                       gcgTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCg                                       >  1:5673414/1‑63 (MQ=255)

                                       gcgTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCattatt                               <  1:4326109/71‑1 (MQ=255)

                                       gcgTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCattat                                <  1:2320579/70‑1 (MQ=255)

                                       gcgTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCattat                                <  1:205515/70‑1 (MQ=255)

                                       gcgTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCattat                                <  1:2738794/70‑1 (MQ=255)

                                         gTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCAttatta                              <  1:3583519/70‑1 (MQ=255)

                                         gTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATtattat                             <  1:3305787/71‑1 (MQ=255)

                                           tGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATgg                           >  1:1163418/1‑71 (MQ=255)

                                            gCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATggg                          >  1:1900878/1‑71 (MQ=255)

                                            gCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATggg                          >  1:1241847/1‑71 (MQ=255)

                                            gCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATgg                           <  1:4498710/70‑1 (MQ=255)

                                             ccGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCacac                                               >  1:2770013/1‑49 (MQ=255)

                                               gttgACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCt                       >  1:557156/4‑71 (MQ=255)

                                                tGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACCCCTATATCGCCATTATTATGGGTCtt                      >  1:4107399/1‑71 (MQ=255)

                                                tGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCtt                      <  1:4504614/71‑1 (MQ=255)

                                                 ggACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATg                            >  1:435331/1‑64 (MQ=255)

                                                 ggACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATg                            >  1:1442507/1‑64 (MQ=255)

                                                  gACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTg                     <  1:5210681/70‑1 (MQ=255)

                                                  gACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTg                     >  1:5487791/1‑70 (MQ=255)

                                                  gACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGc                    <  1:3776526/71‑1 (MQ=255)

                                                  gACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACcaca                                                >  1:2326386/1‑43 (MQ=12)

                                                  gACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCAttatta                              >  1:5683594/1‑61 (MQ=255)

                                                  gACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCtgg                     <  1:2125531/70‑3 (MQ=255)

                                                   aCGTTTATGGTTCACTGTCGCTGACGGGTAAAGGccac                                                    <  1:5252144/38‑1 (MQ=255)

                                                   aCGTTTATGGTTCACTGTCGCTGACGGGTAAAGGccac                                                    <  1:5262900/38‑1 (MQ=255)

                                                   aCGTTTATGGTTCACTGTCGCTGACGGGTAAAGGccac                                                    <  1:87060/38‑1 (MQ=255)

                                                   aCGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGc                    <  1:2226196/70‑1 (MQ=255)

                                                    cGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAg                  <  1:4697721/71‑1 (MQ=255)

                                                    cGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAg                  <  1:1742101/71‑1 (MQ=255)

                                                     gTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGa                                            <  1:347360/44‑1 (MQ=255)

                                                     gTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGa                                            <  1:3319133/44‑1 (MQ=255)

                                                      tttATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAgg                 >  1:3736191/1‑70 (MQ=255)

                                                      tttATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAgg                 >  1:3126734/1‑70 (MQ=255)

                                                      tttATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTACTATg                            <  1:4664236/59‑1 (MQ=255)

                                                        tATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGa                                            >  1:374146/1‑41 (MQ=255)

                                                        tATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGt                         >  1:5135927/1‑60 (MQ=255)

                                                        tATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTa               <  1:3053035/70‑1 (MQ=255)

                                                        tATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTa               <  1:2150604/70‑1 (MQ=255)

                                                             ttCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAc         <  1:481076/71‑1 (MQ=255)

                                                               cACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCt       >  1:1876963/1‑71 (MQ=255)

                                                                aCTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAcc        <  1:1358070/69‑1 (MQ=255)

                                                                aCTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCt       <  1:106621/70‑1 (MQ=255)

                                                                aCTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTg      >  1:4793423/1‑71 (MQ=255)

                                                                  tGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGc     >  1:5539730/1‑70 (MQ=255)

                                                                  tGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGc     >  1:2059994/1‑70 (MQ=255)

                                                                  tGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGc     >  1:5175328/1‑70 (MQ=255)

                                                                  tGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGc     >  1:852152/1‑70 (MQ=255)

                                                                   gTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGc                    <  1:2941518/54‑1 (MQ=255)

                                                                   gTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGc                    <  1:2980099/54‑1 (MQ=255)

                                                                   gTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGc                    <  1:3867789/54‑1 (MQ=255)

                                                                   gTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGa           <  1:2673851/63‑1 (MQ=255)

                                                                   gTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAc         >  1:2793122/1‑65 (MQ=255)

                                                                    tCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCAc  >  1:1145933/1‑71 (MQ=255)

ATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCAC  >  Exported/1200939‑1201077

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`