BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,201,068	A→T	14.2%	E74V (GAA→GTA) ‡	sdaA →	L‑serine deaminase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,201,068	0	A	T	14.2%	305.8 / 35.5	119	E74V (GAA→GTA) ‡	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand): ref base A (52/50); new base T (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 = 9.04e-01

Only 100 of 125 total aligned reads displayed.

GTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTCTGCTG  >  minE/1200999‑1201138

gTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGa                                                                        <  1:912467/70‑1 (MQ=255)

 ttCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAc                                                                      <  1:481076/71‑1 (MQ=255)

   cACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCt                                                                    >  1:1876963/1‑71 (MQ=255)

    aCTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAcc                                                                     <  1:1358070/69‑1 (MQ=255)

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

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

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

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

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

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

       gTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAc                                                                      >  1:2793122/1‑65 (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)

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

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

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

                     aaGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGa                                                   >  1:797924/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:822517/71‑4 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                 ccGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGt                                      >  1:1125878/1‑71 (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:2652223/71‑1 (MQ=255)

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

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

                                       tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAt                                 <  1:349867/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)

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

                                            attattATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                                <  1:5573431/50‑1 (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)

                                                  atGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTag                      <  1:5295182/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)

                                                           gCAGGTAACGAACCTGCGACCGTGGATATCGACAGTATTCCCGGTTTTATTc                               >  1:2844949/1‑52 (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:3789326/70‑1 (MQ=255)

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

                                                               gTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACg         <  1:3729720/70‑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:2158501/1‑63 (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:5413966/1‑71 (MQ=255)

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

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

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

GTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTCTGCTG  >  minE/1200999‑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: `‑`