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%	452.8 / NA	114	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 (54/60); total (54/60)

Only 100 of 114 total aligned reads displayed.

ACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAGG  >  NC_000913/3816532‑3816625

aCGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACt                                              >  1:4799577/1‑50 (MQ=255)

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

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

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

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

 cGGCCTGCTGGAAGCATATTGAAATCTCCGGCTTGAAACAAATGTGACt                                              <  1:7602676/49‑1 (MQ=38)

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

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

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

    ccTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgcgc                                          <  1:4279081/50‑1 (MQ=255)

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

     ctgctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCa                                         >  1:2077103/1‑50 (MQ=255)

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

        ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAtt                                       <  1:4189527/49‑1 (MQ=255)

        ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAtt                                       <  1:6795930/49‑1 (MQ=255)

        ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAtt                                       <  1:7213070/49‑1 (MQ=255)

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

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

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

          ggACGCATATTGAAATCTCCGGCTTGAAACAAATGCGACTGCGCATtata                                    >  1:8476608/1‑50 (MQ=39)

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

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

           gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATAc                                   <  1:4599761/50‑1 (MQ=255)

           gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATAc                                   <  1:3748686/50‑1 (MQ=255)

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

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

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

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

              gCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGa                                >  1:5178014/1‑50 (MQ=255)

              gCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGa                                >  1:2091095/1‑50 (MQ=255)

              gCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGa                                >  1:5220771/1‑50 (MQ=255)

              gCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGa                                <  1:5972218/50‑1 (MQ=255)

               cATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGAc                               >  1:7173331/1‑50 (MQ=255)

               cATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGAc                               >  1:7161632/1‑50 (MQ=255)

               cATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGAc                               >  1:929457/1‑50 (MQ=255)

                atatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACt                              >  1:7916933/1‑50 (MQ=255)

                atatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACt                              >  1:2719545/1‑50 (MQ=255)

                 tatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACtt                             >  1:5685324/1‑50 (MQ=255)

                 tatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACtt                             >  1:5845798/1‑50 (MQ=255)

                 tatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACtt                             <  1:5011707/50‑1 (MQ=255)

                 tatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACtt                             >  1:2245780/1‑50 (MQ=255)

                  atTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTc                            >  1:2472184/1‑50 (MQ=255)

                   ttGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTcc                           >  1:7463715/1‑50 (MQ=255)

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

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

                    tGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcg                          >  1:6390867/1‑50 (MQ=255)

                     gAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcgg                         >  1:49685/1‑50 (MQ=255)

                     gAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcgg                         >  1:3741205/1‑50 (MQ=255)

                     gAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcg                          <  1:1739352/49‑1 (MQ=255)

                     gAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcg                          <  1:257712/49‑1 (MQ=255)

                     gAAATCTCCGGCTTGAAACAAATGTGACAGCGCATTATACGGACTTCcg                          <  1:1692822/49‑1 (MQ=255)

                      aaaTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcggc                        >  1:6205829/1‑50 (MQ=255)

                      aaaTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcggc                        <  1:249012/50‑1 (MQ=255)

                      aaaTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcggc                        <  1:1965023/50‑1 (MQ=255)

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

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

                        aTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcggcgg                      >  1:689365/1‑50 (MQ=255)

                         tctcCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGt                     >  1:6489261/1‑50 (MQ=255)

                           tcCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTa                   >  1:467526/1‑50 (MQ=255)

                           tcCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTa                   >  1:7305955/1‑50 (MQ=255)

                            ccGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTAt                  >  1:4360732/1‑50 (MQ=255)

                            ccGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTAt                  >  1:1829488/1‑50 (MQ=255)

                            ccGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTAt                  >  1:8368634/1‑50 (MQ=255)

                             cGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGtt                    <  1:1849872/47‑1 (MQ=255)

                             cGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGt                     <  1:8175462/46‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                 ttGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTa             >  1:6653022/1‑50 (MQ=255)

                                 ttGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTa             >  1:4449452/1‑50 (MQ=255)

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

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

                                    aaaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATcc          <  1:1281354/50‑1 (MQ=255)

                                    aaaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATcc          >  1:3797805/1‑50 (MQ=255)

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

                                     aaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCt         >  1:6238485/1‑50 (MQ=255)

                                     aaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCt         >  1:2241578/1‑50 (MQ=255)

                                      aCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCt         >  1:2660509/1‑49 (MQ=255)

                                      aCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCt         >  1:232675/1‑49 (MQ=255)

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

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

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

                                        aaaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGAc      <  1:2465977/50‑1 (MQ=255)

                                        aaaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGAc      <  1:6700711/50‑1 (MQ=255)

                                        aaaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGAc      <  1:4954836/50‑1 (MQ=255)

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

                                         aaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACa     <  1:5200560/50‑1 (MQ=255)

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

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

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

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

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

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

                                            gtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAgg  >  1:7523574/1‑50 (MQ=255)

ACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAGG  >  NC_000913/3816532‑3816625

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