BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,760,469	T→C	S261G (AGT→GGT)	selA ←	selenocysteine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,760,469	0	T	C	100.0%	201.1 / NA	57	S261G (AGT→GGT)	selA	selenocysteine synthase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (29/28); total (29/28)

ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGA  >  NC_000913/3760393‑3760546

aCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCt                                                                           <  1:670209/81‑1 (MQ=255)

  gcccgccGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTa                                                                          >  1:670269/1‑80 (MQ=255)

   cccgccGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTaaa                                                                        <  2:670209/81‑1 (MQ=255)

         gcAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGt                                                                   >  2:670312/1‑80 (MQ=255)

            aTCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACt                                                               <  1:670212/81‑1 (MQ=255)

            aTCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACt                                                               <  2:670210/81‑1 (MQ=255)

              cAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTAc                                                             >  1:670284/1‑81 (MQ=255)

               aaCTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTAc                                                             <  2:670212/80‑1 (MQ=255)

                 cTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACggg                                                          <  2:670213/81‑1 (MQ=255)

                 cTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACgg                                                           >  2:670299/1‑80 (MQ=255)

                    cTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGaa                                                        <  2:670215/80‑1 (MQ=255)

                    cTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGa                                                         <  1:670215/79‑1 (MQ=255)

                    cTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGa                                                         <  1:670217/79‑1 (MQ=255)

                        ggCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATcc                                                   <  1:670218/81‑1 (MQ=255)

                        ggCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATc                                                    >  2:670294/1‑80 (MQ=255)

                            ttGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGct                                               <  2:670218/81‑1 (MQ=255)

                             tGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGctc                                              >  2:670279/1‑81 (MQ=255)

                               gCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGctct                                             <  1:670221/80‑1 (MQ=255)

                               gCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGctct                                             <  2:670221/80‑1 (MQ=255)

                               gCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGctct                                             <  1:670220/80‑1 (MQ=255)

                               gCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCtt                                            <  1:670223/81‑1 (MQ=255)

                                ctctTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCttt                                           <  2:670223/81‑1 (MQ=255)

                                    ttCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGcc                                        >  1:670308/1‑80 (MQ=255)

                                      cGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAg                                     >  1:670282/1‑81 (MQ=255)

                                       ggCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAgc                                    >  1:670302/1‑81 (MQ=255)

                                       ggCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAgc                                    >  2:670308/1‑81 (MQ=255)

                                          aaaCCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAgcgc                                  >  2:670320/1‑80 (MQ=255)

                                          aaaCCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGcc                                 >  2:670310/1‑81 (MQ=255)

                                           aaCCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGcca                                >  2:670292/1‑81 (MQ=255)

                                             ccGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGccacc                              <  1:670225/81‑1 (MQ=255)

                                             ccGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGccac                               >  1:670292/1‑80 (MQ=255)

                                               gTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAg                            >  2:670289/1‑81 (MQ=255)

                                               gTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAg                            <  2:670225/81‑1 (MQ=255)

                                                 aCTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGtt                          >  2:670327/1‑81 (MQ=255)

                                                    ggCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCg                       >  1:670294/1‑81 (MQ=255)

                                                    ggCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCg                       >  1:670283/1‑81 (MQ=255)

                                                     gCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTcc                        <  1:670230/79‑1 (MQ=255)

                                                     gCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTcc                        <  1:670228/79‑1 (MQ=255)

                                                     gCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCg                       <  1:670227/80‑1 (MQ=255)

                                                     gCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGc                      <  2:670228/81‑1 (MQ=255)

                                                        aaGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCtt                    >  2:670325/1‑80 (MQ=255)

                                                        aaGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCtt                    >  2:670285/1‑80 (MQ=255)

                                                         aGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCtt                    <  1:670231/79‑1 (MQ=255)

                                                          gATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCAt                 <  2:670231/81‑1 (MQ=255)

                                                          gATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCAt                 >  2:670300/1‑81 (MQ=255)

                                                             cGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTa              >  2:670306/1‑81 (MQ=255)

                                                              gACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTAt             >  2:670295/1‑81 (MQ=255)

                                                                ccAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCg           <  1:670233/81‑1 (MQ=255)

                                                                  aGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCt         >  1:670303/1‑81 (MQ=255)

                                                                   gCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCtt        >  1:670336/1‑81 (MQ=255)

                                                                    cGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCttt       <  1:670234/81‑1 (MQ=255)

                                                                     gAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTg      >  1:670323/1‑81 (MQ=255)

                                                                     gAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTg      >  1:670315/1‑81 (MQ=255)

                                                                     gAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTg      >  2:670317/1‑81 (MQ=255)

                                                                     gAGCCACCGCCTAAATCAGTCACTACGGGAACATCCAGATCTTTGACGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTg      >  2:670318/1‑81 (MQ=255)

                                                                        ccACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTg   <  1:670235/81‑1 (MQ=255)

                                                                         cACCGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGa  <  2:670235/81‑1 (MQ=255)

ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGA  >  NC_000913/3760393‑3760546

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