BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,201,033	C→T	19.6%	T62T (ACC→ACT)	sdaA →	L‑serine deaminase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,201,033	0	C	T	19.6%	253.4 / 51.4	122	T62T (ACC→ACT)	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand): ref base C (51/47); new base T (11/13); total (62/60)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.17e-01

Only 100 of 127 total aligned reads displayed.

GGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGT  >  Exported/1200964‑1201102

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)

         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)

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

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

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

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

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

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

              gcgTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCattatt                                                        <  1:4326109/71‑1 (MQ=255)

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

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

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

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

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

                         gACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGc                                             <  1:3776526/71‑1 (MQ=255)

                         gACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCAttatta                                                       >  1:5683594/1‑61 (MQ=255)

                         gACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCtgg                                              <  1:2125531/70‑3 (MQ=255)

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

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

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

                           cGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAg                                           <  1:1742101/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)

                               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)

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

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

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

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

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

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

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

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

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

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

                                               tGACCGGTAAAGGCCACCACACTGATATCGCCattatt                                                        <  1:4687686/38‑1 (MQ=38)

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

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

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

                                                    ggTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGatat                  >  1:5640380/1‑71 (MQ=255)

                                                     gTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGatat                  >  1:16107/1‑70 (MQ=255)

                                                      tAAAGGCCACCACACCGATATCGCCATTATTATgg                                                    >  1:4947558/1‑35 (MQ=255)

                                                      tAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTaa                                       >  1:1116551/1‑48 (MQ=255)

                                                      tAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTaa                                       >  1:1601581/1‑48 (MQ=255)

                                                      tAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTaa                                       >  1:2002513/1‑48 (MQ=255)

                                                      tAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACg                                     <  1:568434/50‑1 (MQ=255)

                                                      tAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACg                                     <  1:3137000/50‑1 (MQ=255)

                                                      tAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACg                                     <  1:4425788/50‑1 (MQ=255)

                                                      tAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACg                                     <  1:2645510/50‑1 (MQ=255)

                                                       aaaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGata                   >  1:2191146/1‑67 (MQ=255)

                                                        aaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGa               <  1:3802516/70‑1 (MQ=255)

                                                          ggCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc            <  1:3683581/71‑4 (MQ=255)

                                                          ggCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc            <  1:822517/71‑4 (MQ=255)

                                                          ggCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc            <  1:2724485/71‑4 (MQ=255)

                                                          ggCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCg                         >  1:3985719/1‑58 (MQ=255)

                                                           gCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc            >  1:3360804/1‑67 (MQ=255)

                                                           gCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGa                     >  1:1748441/1‑61 (MQ=255)

                                                            ccaccaCACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGatt             <  1:5359276/68‑3 (MQ=255)

                                                              accacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGa                                    <  1:2920708/43‑1 (MQ=255)

                                                              accacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGa                                    <  1:520764/43‑1 (MQ=255)

                                                               ccacACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATc                 <  1:1903575/61‑1 (MQ=255)

                                                               ccacACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc            <  1:5044630/66‑4 (MQ=255)

                                                               ccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAcc                                 <  1:4402897/45‑1 (MQ=255)

                                                               ccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAcc                                 <  1:5312026/45‑1 (MQ=255)

                                                               ccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTg                               <  1:1413433/47‑1 (MQ=255)

                                                                cacaCTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATc                 <  1:4232881/60‑1 (MQ=255)

                                                                cacaCCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAc                                  >  1:4101675/1‑43 (MQ=255)

                                                                cacaCCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAc                                  >  1:1883555/1‑43 (MQ=255)

                                                                 acacCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGa                     >  1:1579627/1‑55 (MQ=255)

                                                                 acacCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTccc     >  1:2749234/1‑71 (MQ=255)

                                                                  cacCGATATCGCCATTATTATGGTTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCg    >  1:688421/1‑71 (MQ=255)

                                                                  cacCGATATCGCCATTATTATGGGTCTTGCAGGTaa                                       >  1:523184/1‑36 (MQ=255)

                                                                  cacCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCg    >  1:4176217/1‑71 (MQ=255)

                                                                  cacCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCg    >  1:1683737/1‑71 (MQ=255)

                                                                   acCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCgg   >  1:5515268/1‑71 (MQ=255)

                                                                   acCGATATCGCCATTATTATGGGACTTCCACGTaa                                       >  1:4293903/1‑35 (MQ=21)

                                                                    ccGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGt  >  1:1125878/1‑71 (MQ=255)

                                                                    ccGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGt  >  1:3155375/1‑71 (MQ=255)

GGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGT  >  Exported/1200964‑1201102

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