BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,201,074	C→A	16.0%	A76D (GCC→GAC) ‡	sdaA →	L‑serine deaminase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,201,074	0	C	A	16.0%	292.0 / 29.2	106	A76D (GCC→GAC) ‡	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand): ref base C (45/44); new base A (6/11); total (51/55)
Fisher's exact test for biased strand distribution p-value = 2.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

Only 100 of 112 total aligned reads displayed.

TGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTCTGCTGCTGGCA  >  Exported/1201005‑1201144

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

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

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

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

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

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

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

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

            gTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATc                                                           <  1:3545152/71‑1 (MQ=255)

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

               aaGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGa                                                         >  1:797924/1‑70 (MQ=255)

               aaGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGa                                                         >  1:4609480/1‑70 (MQ=255)

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

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

                aGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGatt                                                       <  1:1721752/71‑3 (MQ=255)

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

                 ggCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc                                                      <  1:3683581/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)

                   ccaccaCACCGATATCGCCATTTTTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTa                                                    >  1:5507829/1‑71 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                            tgATATCGCCATTATTATGGGCCTTGCGGGTAACCTGCCGGATACCGTGGATATCGat                                                        <  1:2910886/57‑2 (MQ=255)

                              atatCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGAc                                                        <  1:5570109/56‑1 (MQ=255)

                              atatCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGAc                                                        <  1:2376697/56‑1 (MQ=255)

                              atatCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGttt                                          >  1:4485513/1‑70 (MQ=255)

                               tatCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTa                                        <  1:14276/71‑1 (MQ=255)

                               tatCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTa                                        <  1:2652223/71‑1 (MQ=255)

                               tatCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTa                                        <  1:5695205/71‑1 (MQ=255)

                                atCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCAc                                                                     >  1:5268274/1‑41 (MQ=255)

                                atCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAt                                       >  1:2619947/1‑71 (MQ=255)

                                 tCGCCATTATTATGGGTCTTGCAGGTAAGGAACCTGCCACCGTGGatat                                                            >  1:8664/1‑49 (MQ=255)

                                 tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCg                                                                   >  1:181865/1‑42 (MQ=255)

                                 tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCg                                                                   >  1:799848/1‑42 (MQ=255)

                                 tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGatat                                                            >  1:3658793/1‑49 (MQ=255)

                                 tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAtt                                      >  1:4179570/1‑71 (MQ=255)

                                 tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAt                                       <  1:349867/70‑1 (MQ=255)

                                 tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAt                                       <  1:5274320/70‑1 (MQ=255)

                                 tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAt                                       <  1:5637959/70‑1 (MQ=255)

                                 tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAt                                       <  1:1733972/70‑1 (MQ=255)

                                 tCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTAt                                       >  1:2231524/1‑70 (MQ=255)

                                  cGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTa                                                    >  1:42484/1‑56 (MQ=255)

                                    ccATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGt                                                                  <  1:334685/40‑1 (MQ=255)

                                    ccATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGt                                                                  <  1:1404212/40‑1 (MQ=255)

                                    ccATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTcgc                                   <  1:438116/71‑1 (MQ=255)

                                     cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                                      <  1:4998195/51‑1 (MQ=255)

                                     cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                                      <  1:3015882/51‑1 (MQ=255)

                                     cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                                      <  1:166450/51‑1 (MQ=255)

                                     cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                                      <  1:3332616/51‑1 (MQ=255)

                                     cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                                      <  1:712724/51‑1 (MQ=255)

                                     cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTcgcg                                  <  1:4439293/71‑1 (MQ=255)

                                     cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTcgc                                   >  1:458109/1‑70 (MQ=255)

                                      attattATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                                      <  1:4008774/50‑1 (MQ=255)

                                      attattATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGa                                 >  1:3070030/1‑71 (MQ=255)

                                        tattatGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACg                               >  1:5206861/1‑71 (MQ=255)

                                        tattatGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACg                               >  1:5009563/1‑71 (MQ=255)

                                        tattatGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTAt                                       <  1:918092/63‑1 (MQ=255)

                                         attatGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTc                                     >  1:2963194/1‑64 (MQ=255)

                                          ttatGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTa                             <  1:4607497/71‑1 (MQ=255)

                                            atGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTag                            <  1:1330065/70‑1 (MQ=255)

                                             tGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTagaa                          <  1:2686352/71‑1 (MQ=255)

                                               ggTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTagaaga                        <  1:5659484/71‑1 (MQ=255)

                                                 tCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAgc                      >  1:4706057/1‑71 (MQ=255)

                                                 tCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAgc                      <  1:1277218/71‑1 (MQ=255)

                                                 tCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAgc                      <  1:5406685/71‑1 (MQ=255)

                                                     gCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAtt                                      >  1:4062546/1‑51 (MQ=255)

                                                     gCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTc                                     >  1:3193066/1‑52 (MQ=255)

                                                     gCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTc                                     >  1:1314832/1‑52 (MQ=255)

                                                     gCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGcgcg                   <  1:3664178/70‑1 (MQ=255)

                                                     gCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGcgcg                   <  1:3789326/70‑1 (MQ=255)

                                                        ggTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACg               <  1:4591619/71‑1 (MQ=255)

                                                         gTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGt              <  1:83741/71‑1 (MQ=255)

                                                          tAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGcgcg                   <  1:2144097/65‑1 (MQ=255)

                                                           aaCGAACCTGCCACCGTGGATATCGACAGTATTCCCGGtt                                           >  1:2542514/1‑40 (MQ=255)

                                                           aaCGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAgcgc                    >  1:5675277/1‑63 (MQ=255)

                                                           aaCGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAgcgc                    >  1:2158501/1‑63 (MQ=255)

                                                            aCGAACCTGCCACCGTGGATATCGACAGTATTCCCGGttt                                          <  1:608601/40‑1 (MQ=255)

                                                            aCGAACCTGCCACCGTGGATATCGACAGTATTCCCGGttt                                          <  1:5580462/40‑1 (MQ=255)

                                                            aCGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAAc                >  1:1699940/1‑66 (MQ=255)

                                                            aCGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAAc                >  1:2893734/1‑66 (MQ=255)

                                                            aCGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTctg           <  1:1717414/71‑1 (MQ=255)

                                                            aCGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTctg           >  1:5413966/1‑71 (MQ=255)

                                                             cGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTagaa                          >  1:543348/1‑55 (MQ=255)

                                                               aaCCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTctgctg        >  1:4854142/1‑71 (MQ=255)

                                                               aaCCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTctgatg        >  1:1298047/1‑71 (MQ=255)

                                                                 ccTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTctgctgct      <  1:2169312/71‑1 (MQ=255)

                                                                  cTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAgc                      <  1:2971960/54‑1 (MQ=255)

                                                                     ccACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTCTGCTGCTGGca  <  1:2405849/71‑1 (MQ=255)

                                                                     ccACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTCTGCTGCTGGca  <  1:648721/71‑1 (MQ=255)

TGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTCTGCTGCTGGCA  >  Exported/1201005‑1201144

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 21 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`