BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,201,015	G→C	12.2%	T56T (ACG→ACC)	sdaA →	L‑serine deaminase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,201,015	0	G	C	12.2%	277.9 / 23.9	106	T56T (ACG→ACC)	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand): ref base G (38/54); new base C (9/4); total (47/59)
Fisher's exact test for biased strand distribution p-value = 7.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 108 total aligned reads displayed.

TGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTG  >  Exported/1200945‑1201081

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

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

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

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

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

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

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

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

     gAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTa                                                                <  1:3705772/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)

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

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

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

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

            gCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGccac                                                        <  1:4158358/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)

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

                     aTAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTa                                                                <  1:2541638/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: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:1827148/61‑1 (MQ=255)

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

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

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

                                 gcgTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCattat                                    <  1:205515/70‑1 (MQ=255)

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

                                         gttgACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCt                           >  1:557156/4‑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: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:5262900/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:347360/44‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

                                                              tCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAg                      <  1:1518438/55‑1 (MQ=255)

                                                              tCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCAc      >  1:1145933/1‑71 (MQ=255)

                                                                tcTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATg                                >  1:2541223/2‑43 (MQ=38)

                                                                tcTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATg                                >  1:3872182/2‑43 (MQ=38)

                                                                tcTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATg                                >  1:4313237/2‑43 (MQ=38)

                                                                  tGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTg  >  1:2724003/1‑71 (MQ=255)

                                                                  tGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAAc                 >  1:24917/1‑56 (MQ=255)

                                                                   gACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTg  >  1:2265245/1‑70 (MQ=255)

                                                                    aCGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCAACCGTg  <  1:2093796/69‑1 (MQ=255)

TGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTG  >  Exported/1200945‑1201081

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