BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,131,341	C→T	356 (TAG→TAA)	yghQ ←	putative inner membrane polysaccharide flippase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,131,341	0	C	T	100.0%	120.9 / NA	38	356 (TAG→TAA)	yghQ	putative inner membrane polysaccharide flippase
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (18/20); total (18/20)

CGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACAC  >  NC_000913/3131220‑3131468

cGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGc                                                                                                                >  2:341207/1‑139 (MQ=255)

cGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGc                                                                                                                >  1:670419/1‑139 (MQ=255)

 gaagaaACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCt                                                                                                               <  1:10000/139‑1 (MQ=255)

      aaCGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCgag                                                                                                          <  2:836959/139‑1 (MQ=255)

       aCGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCgaga                                                                                                         <  2:1011502/139‑1 (MQ=255)

               aaTAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAAc                                                                                                 >  1:1016454/1‑139 (MQ=255)

                     cTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCaaa                                                                                           <  1:490149/139‑1 (MQ=255)

                         ttCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACAcc                                                                                       >  1:867779/1‑139 (MQ=255)

                         ttCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACAcc                                                                                       >  2:765953/1‑139 (MQ=255)

                         ttCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACAcc                                                                                       >  1:549726/1‑139 (MQ=255)

                            aaGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGc                                                                                    <  2:670419/139‑1 (MQ=255)

                                  tttACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATg                                                                              >  2:1062313/1‑139 (MQ=255)

                                   ttACCGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCa                                                                                      >  2:906902/1‑130 (MQ=255)

                                   ttACCGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCa                                                                                      <  1:906902/130‑1 (MQ=255)

                                          agaAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGc                                                                                                                <  1:674275/97‑1 (MQ=255)

                                          agaAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGc                                                                                                                >  2:674275/1‑97 (MQ=255)

                                            aaaaaaTCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGc                                                                    >  2:164294/1‑139 (MQ=255)

                                              aaaaTCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCcg                                                                  <  1:765953/139‑1 (MQ=255)

                                                             aTGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACcgcg                                                   <  1:899107/139‑1 (MQ=255)

                                                              tGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGa                                                  >  1:174662/1‑139 (MQ=255)

                                                                  gaaatgaaatGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAgcgc                                              <  1:164294/139‑1 (MQ=255)

                                                                  gaaatgaaatGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAgcgc                                              <  1:1062313/139‑1 (MQ=255)

                                                                     atgaaatgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGccac                                           >  1:573412/1‑139 (MQ=255)

                                                                     atgaaatgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGccac                                           >  1:970989/1‑139 (MQ=255)

                                                                             aCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATAc                                   >  1:966497/1‑139 (MQ=255)

                                                                                gAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCg                                                                                                            >  1:962362/1‑63 (MQ=255)

                                                                                gAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCg                                                                                                            <  2:962362/63‑1 (MQ=255)

                                                                                      gCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCgaga                                                                                                         <  2:506730/60‑1 (MQ=255)

                                                                                      gCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCgaga                                                                                                         >  1:506730/1‑60 (MQ=255)

                                                                                       ccGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAg                                                                            >  2:489652/1‑88 (MQ=255)

                                                                                       ccGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAg                                                                            <  1:489652/88‑1 (MQ=255)

                                                                                       ccGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTccgcc                         <  1:742185/139‑1 (MQ=255)

                                                                                            aagaaTTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCa                    <  1:341207/139‑1 (MQ=255)

                                                                                                    ccccGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTGCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCACAATACCGATTCCGCCCGCCAGTAAACCg            <  2:573412/139‑1 (MQ=255)

                                                                                                              tCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCacac  <  2:867779/139‑1 (MQ=255)

                                                                                                              tCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCacac  <  2:966497/139‑1 (MQ=255)

                                                                                                                      tCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTccgcc                         >  2:29374/1‑108 (MQ=255)

                                                                                                                      tCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTccgcc                         <  1:29374/108‑1 (MQ=255)

CGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACAC  >  NC_000913/3131220‑3131468

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