BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,201,003	A→G	8.5%	S52S (TCA→TCG)	sdaA →	L‑serine deaminase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,201,003	0	A	G	8.5%	297.0 / 13.1	106	S52S (TCA→TCG)	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand): ref base A (42/55); new base G (6/3); total (48/58)
Fisher's exact test for biased strand distribution p-value = 2.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 113 total aligned reads displayed.

TCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTG  >  Exported/1200936‑1201073

tcgATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTg                                                                     <  1:561487/71‑1 (MQ=255)

tcgATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTg                                                                     <  1:3501568/71‑1 (MQ=255)

 cgATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGt                                                                    >  1:3661490/1‑71 (MQ=255)

  gatgatCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTc                                                                   >  1:5333594/1‑71 (MQ=255)

  gatgatCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTc                                                                   >  1:5316686/1‑71 (MQ=255)

  gatgatCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTc                                                                   >  1:5263452/1‑71 (MQ=255)

  gatgatCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTc                                                                   >  1:712552/1‑71 (MQ=255)

   atgatCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTc                                                                   >  1:3320730/1‑70 (MQ=255)

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

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

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

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

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

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

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

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

             cGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTgg                                                                    >  1:1304183/1‑58 (MQ=255)

             cGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGt                                                                    >  1:845967/1‑59 (MQ=255)

             cGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGt                                                                    >  1:774716/1‑59 (MQ=255)

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

              gAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGa                                                              >  1:4514924/1‑64 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                         aCTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCa                                                                        <  1:249137/43‑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)

                              aTAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTc                                                                   <  1:5162358/43‑1 (MQ=255)

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

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

                               agacGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCacac                                           >  1:5109075/4‑66 (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:3298186/1‑71 (MQ=255)

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

                                      aCTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCat                               >  1:3678227/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:5310405/61‑1 (MQ=255)

                                        tCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGat                                       <  1:2844025/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)

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

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

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

                                               gCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATggg                      >  1:1900878/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:5487791/1‑70 (MQ=255)

                                                     gACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTg                 <  1:5210681/70‑1 (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:87060/38‑1 (MQ=255)

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

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

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

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

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

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

                                                        gTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGa                                        <  1:347360/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)

                                                           tATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTaa          >  1:5435198/1‑71 (MQ=255)

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

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

                                                               gTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGa       <  1:912467/70‑1 (MQ=255)

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

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

                                                                   aCTGTCGCTGACGGGTAAAGTCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAcc    <  1:3951882/69‑1 (MQ=255)

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

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

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

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

TCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTG  >  Exported/1200936‑1201073

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