BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,201,000	T→C	7.2%	G51G (GGT→GGC)	sdaA →	L‑serine deaminase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,201,000	0	T	C	7.2%	295.6 / 10.3	110	G51G (GGT→GGC)	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand): ref base T (42/60); new base C (5/3); total (47/63)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

Only 100 of 113 total aligned reads displayed.

GTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAC  >  Exported/1200931‑1201070

gTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGtt                                                                       <  1:3099586/71‑1 (MQ=255)

 ttCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTc                                                                      <  1:2124328/71‑1 (MQ=255)

 ttCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTc                                                                      <  1:1355089/71‑1 (MQ=255)

 ttCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTc                                                                      <  1:5388239/71‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

           atCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCt                                                             <  1:2025536/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:5312343/68‑1 (MQ=255)

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

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

                 tCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGt                                                      <  1:1062648/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:4885595/1‑71 (MQ=255)

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

                       aaGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGc                                                <  1:3214937/71‑1 (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)

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

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

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

                                    agacGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCacac                                        >  1:5109075/4‑66 (MQ=255)

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

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

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

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

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

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

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

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

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

                                             tCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGat                                    <  1:5104891/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:205515/70‑1 (MQ=255)

                                               gcgTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCattat                         <  1:2320579/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:1241847/1‑71 (MQ=255)

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

                                                    gCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATgg                    <  1:4498710/70‑1 (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)

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

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

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

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

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

                                                           aCGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGc             <  1:2226196/70‑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:3126734/1‑70 (MQ=255)

                                                              tttATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAgg          >  1:3736191/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:2150604/70‑1 (MQ=255)

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

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

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

GTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAC  >  Exported/1200931‑1201070

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