BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,200,994	T→G	6.9%	V49V (GTT→GTG)	sdaA →	L‑serine deaminase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,200,994	0	T	G	6.9%	299.0 / 13.0	115	V49V (GTT→GTG)	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand): ref base T (42/65); new base G (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 = 1.00e+00

Only 100 of 126 total aligned reads displayed.

AAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGG  >  Exported/1200926‑1201062

aaaCAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTa                                                                    >  1:2306499/1‑71 (MQ=255)

aaaCAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTa                                                                    >  1:3959729/1‑71 (MQ=255)

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

  aCAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGttt                                                                     <  1:3069208/68‑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:3546404/71‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

               gatCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGtt                                                                      <  1:4535171/54‑1 (MQ=255)

               gatCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGtt                                                                      <  1:233003/54‑1 (MQ=255)

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

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

                atCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCt                                                     >  1:2479764/1‑70 (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:1062648/71‑1 (MQ=255)

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

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

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

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

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

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

                        gAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGttt                                                                     >  1:4812674/1‑46 (MQ=255)

                        gAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGttt                                                                     >  1:4484706/1‑46 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                                 ttACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTAt                                                                   <  1:569788/39‑1 (MQ=255)

                                  tACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCa                                                             <  1:954589/44‑1 (MQ=255)

                                   aCTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCa                                                             <  1:249137/43‑1 (MQ=255)

                                      ggATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCt                                                     <  1:773802/48‑1 (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: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:780093/1‑65 (MQ=255)

                                                aCTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGatat                          >  1:4445484/1‑65 (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:2844025/61‑1 (MQ=255)

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

                                                  tCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGat                            <  1:1827148/61‑1 (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)

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

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

                                                         gCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATggg           >  1:1900878/1‑71 (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:1442507/1‑64 (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:87060/38‑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: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)

AAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGG  >  Exported/1200926‑1201062

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