BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,201,067	G→C	13.8%	E74Q (GAA→CAA) ‡	sdaA →	L‑serine deaminase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,201,067	0	G	C	13.8%	312.4 / 32.9	123	E74Q (GAA→CAA) ‡	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand): ref base G (52/54); new base C (6/11); total (58/65)
Fisher's exact test for biased strand distribution p-value = 3.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

Only 100 of 128 total aligned reads displayed.

GTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTCTG  >  Exported/1200999‑1201135

gTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGa                                                                     <  1:912467/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:2059994/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)

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

                   tAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACg                                                                      <  1:2963199/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)

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

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

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

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

                      aGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGatt                                              <  1:1721752/71‑3 (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)

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

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

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

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

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

                             cacaCCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAc                                                                   >  1:4101675/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)

                               cacCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCg                                     >  1:4176217/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)

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

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

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

                                      atCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAt                              >  1:2619947/1‑71 (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:1733972/70‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                              tattatGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACg                      >  1:5009563/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:5295182/70‑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:4706057/1‑71 (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)

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

                                                               gTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGt     <  1:1867898/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)

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

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

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

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

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

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

GTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTCTG  >  Exported/1200999‑1201135

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