BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,201,006	G→C	7.1%	L53L (CTG→CTC)	sdaA →	L‑serine deaminase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,201,006	0	G	C	7.1%	334.3 / 10.1	112	L53L (CTG→CTC)	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand): ref base G (46/57); new base C (5/3); total (52/60)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

TCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGC  >  Exported/1200936‑1201074

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

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

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

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

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

   atgatCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCg                                                                   <  1:1719046/71‑1 (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: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)

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

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

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

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

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

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

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

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

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

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

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

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

                     gCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGccac                                                 <  1:1785162/71‑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)

                              aTAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTa                                                         <  1:3828588/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:3298186/1‑71 (MQ=255)

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

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

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

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

                                          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)

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

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

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

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

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

                                                    ggACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATg                         >  1:435331/1‑64 (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)

                                                      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)

                                                         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: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)

                                                                   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)

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

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

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

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

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

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

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

TCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGC  >  Exported/1200936‑1201074

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