BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,200,997	T→C	6.9%	Y50Y (TAT→TAC)	sdaA →	L‑serine deaminase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,200,997	0	T	C	6.9%	329.7 / 10.5	115	Y50Y (TAT→TAC)	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand): ref base T (42/65); new base C (5/3); total (47/68)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

Only 100 of 121 total aligned reads displayed.

AACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAA  >  Exported/1200927‑1201065

aaCAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCGGTGGACGTTTAt                                                                      <  1:2383997/71‑1 (MQ=255)

aaCAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTAt                                                                      <  1:2589193/71‑1 (MQ=255)

 aCAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATg                                                                     <  1:887002/71‑1 (MQ=255)

 aCAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATg                                                                     <  1:3819702/71‑1 (MQ=255)

 aCAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATg                                                                     <  1:1663697/71‑1 (MQ=255)

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

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

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

     ttCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTc                                                                 <  1:2124328/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:5316686/1‑71 (MQ=255)

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

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

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

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

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

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

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

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

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

                    gTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGAc                                                     <  1:5312343/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:774716/1‑59 (MQ=255)

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

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

                       gAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTAt                                                                      >  1:225609/1‑48 (MQ=255)

                       gAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTAt                                                                      >  1:2492291/1‑48 (MQ=255)

                       gAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTAt                                                                      >  1:850558/1‑48 (MQ=255)

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

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

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

                       gAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTa                                                <  1:3107155/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:1916355/1‑71 (MQ=255)

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

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

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

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

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

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

                                 tACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCa                                                                <  1:954589/44‑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)

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

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

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

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

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

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

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

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

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

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

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

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

                                                   gcgTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCattat                    <  1:205515/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)

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

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

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

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

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

AACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAA  >  Exported/1200927‑1201065

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