BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,201,069	A→G	14.3%	E74E (GAA→GAG) ‡	sdaA →	L‑serine deaminase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,201,069	0	A	G	14.3%	276.9 / 34.2	119	E74E (GAA→GAG) ‡	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand): ref base A (52/50); new base G (6/11); total (58/61)
Fisher's exact test for biased strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01

Only 100 of 121 total aligned reads displayed.

TTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTCTGCTG  >  Exported/1201000‑1201138

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)

   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)

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

           tGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTg                                                           >  1:2724003/1‑71 (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)

                 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)

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

                      ggCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc                                                <  1:3683581/71‑4 (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)

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

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

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

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

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

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

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

                              cacCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCg                                        >  1:1683737/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:2376697/56‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                      tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAt                                 <  1:5274320/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:1404212/40‑1 (MQ=255)

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

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

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

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

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

                                          cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                                <  1:3015882/51‑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)

                                             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)

                                                     gTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAagag                 <  1:4348460/71‑1 (MQ=255)

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

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

                                                          gCAGGTAACGAACCTGCGACCGTGGATATCGACAGTATTCCCGGTTTTATTc                               >  1:2844949/1‑52 (MQ=255)

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

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

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

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

                                                              gTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTagaa                    <  1:3319912/59‑1 (MQ=255)

                                                              gTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACg         <  1:3729720/70‑1 (MQ=255)

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

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

                                                              gTAACGAACATGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGt        >  1:4708744/1‑71 (MQ=255)

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

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

                                                                aaCGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAgcgc              >  1:5675277/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)

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

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

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

TTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTCTGCTG  >  Exported/1201000‑1201138

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