BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,816,579	G→A	intergenic (‑30/‑97)	rph ← / → yicC	truncated RNase PH/UPF0701 family protein YicC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,816,579	0	G	A	100.0%	200.1 / NA	51	intergenic (‑30/‑97)	rph/yicC	truncated RNase PH/UPF0701 family protein YicC
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (18/33); total (18/33)

CGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAG  >  NC_000913/3816533‑3816624

cGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACt                                             <  1:2719415/49‑1 (MQ=255)

cGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTg                                            <  1:2252842/50‑1 (MQ=255)

 ggCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgc                                           >  1:3602199/1‑50 (MQ=255)

 ggCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgc                                           >  1:3092993/1‑50 (MQ=255)

 ggCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgc                                           >  1:278760/1‑50 (MQ=255)

 ggCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgc                                           <  1:2018153/50‑1 (MQ=255)

   ccTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgcg                                          >  1:1323158/1‑49 (MQ=255)

     tgctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAt                                       >  1:3630496/1‑50 (MQ=255)

       ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAt                                       <  1:262803/48‑1 (MQ=255)

       ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATta                                     <  1:831542/50‑1 (MQ=255)

        tgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATta                                     >  1:2806522/1‑49 (MQ=255)

        tgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATta                                     >  1:2822132/1‑49 (MQ=255)

        tgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATta                                     >  1:473670/1‑49 (MQ=255)

         ggACGCATATTGAAATCTCCGGCTTGAAACTAATGTGACTGCGCATtata                                   >  1:2841512/1‑50 (MQ=39)

          gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATtata                                   <  1:2530397/49‑1 (MQ=255)

          gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATtata                                   <  1:292974/49‑1 (MQ=255)

          gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATtat                                    <  1:672918/48‑1 (MQ=255)

           aCGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACg                                 <  1:3769803/50‑1 (MQ=255)

            cGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACgg                                <  1:2492073/50‑1 (MQ=255)

            cGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACg                                 >  1:3472828/1‑49 (MQ=255)

              cATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGAc                              <  1:3742960/50‑1 (MQ=255)

               atatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACt                             <  1:3808089/50‑1 (MQ=255)

                  ttGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGTCTTcc                          <  1:3249605/50‑1 (MQ=255)

                  ttGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTcc                          <  1:1443767/50‑1 (MQ=255)

                  ttGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTcc                          <  1:1927327/50‑1 (MQ=255)

                  ttGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTcc                          <  1:1638202/50‑1 (MQ=255)

                     aaaTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcgg                        <  1:959855/49‑1 (MQ=255)

                      aaTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcggcg                      <  1:270649/50‑1 (MQ=255)

                         ctcCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGt                    >  1:2761995/1‑49 (MQ=255)

                            cGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTa                  <  1:2810711/48‑1 (MQ=255)

                            cGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTAtt                <  1:2564723/50‑1 (MQ=255)

                             ggCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTc               >  1:146400/1‑50 (MQ=255)

                             ggCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTc               <  1:2151633/50‑1 (MQ=255)

                              gCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTcc              >  1:2175564/1‑50 (MQ=255)

                               cTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCt             <  1:1989248/50‑1 (MQ=255)

                                  gAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTa            <  1:1041815/48‑1 (MQ=255)

                                  gAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTAt           <  1:2217512/49‑1 (MQ=255)

                                   aaaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATc          <  1:169509/49‑1 (MQ=255)

                                    aaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCt        <  1:2547881/50‑1 (MQ=255)

                                     aCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTg       <  1:3266161/50‑1 (MQ=255)

                                     aCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTg       <  1:2346064/50‑1 (MQ=255)

                                     aCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTg       <  1:629692/50‑1 (MQ=255)

                                      cAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGa      >  1:1617448/1‑50 (MQ=255)

                                       aaaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGa      >  1:3727511/1‑49 (MQ=255)

                                        aaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGAc     >  1:3410687/1‑49 (MQ=255)

                                         aTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGAc     >  1:1532606/1‑48 (MQ=255)

                                          tgtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAg  >  1:385451/1‑50 (MQ=255)

                                          tgtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAg  <  1:3150289/50‑1 (MQ=255)

                                          tgtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAg  <  1:773422/50‑1 (MQ=255)

                                          tgtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAg  <  1:2411200/50‑1 (MQ=255)

                                          tgtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAg  <  1:3237627/50‑1 (MQ=255)

CGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAG  >  NC_000913/3816533‑3816624

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