BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_017635	4,127,945	(A)_10→9	intergenic (‑43/+92)	bglF ← / ← bglG	PTS beta‑glucoside transporter subunit IIABC/transcriptional antiterminator BglG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_017635	4,127,936	0	A	.	100.0%	157.5 / NA	36	intergenic (‑34/+101)	bglF/bglG	PTS beta‑glucoside transporter subunit IIABC/transcriptional antiterminator BglG
Reads supporting (aligned to +/- strand): ref base A (0/0); new base . (19/17); total (19/17)

GGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATT  >  NC_017635/4127888‑4127980

ggCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCa                                              >  1:592432/1‑49 (MQ=255)

   tAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCaaaa                                           >  2:721927/1‑49 (MQ=255)

   tAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCaaaa                                           >  2:696090/1‑49 (MQ=255)

     aCTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCaaaaaaa                                        <  2:543199/50‑1 (MQ=255)

        ccGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCaaaaaaaaac                                     >  2:577714/1‑49 (MQ=255)

         cGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCaaaaaaaaacc                                    >  1:1003251/1‑48 (MQ=255)

            cATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAAACCTGa                                <  2:735390/50‑1 (MQ=39)

             aTAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAAACCTGaa                               <  2:1225310/50‑1 (MQ=39)

             aTAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAAACCTGa                                >  1:57724/1‑49 (MQ=39)

             aTAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAAACCTGa                                >  2:575898/1‑49 (MQ=39)

              tAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAAACCTGaaa                              <  1:61369/50‑1 (MQ=39)

              tAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAAACCTGaa                               >  2:1063315/1‑49 (MQ=39)

              tAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAAACCTGaa                               >  2:1560134/1‑49 (MQ=39)

               aaCTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAAACCTGaaa                              >  2:1158183/1‑49 (MQ=39)

                  ttGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAAACCTGaaaaaaa                          >  2:1949946/1‑50 (MQ=39)

                    gCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAAACCTGaaaaaaa                          <  1:280915/48‑1 (MQ=39)

                      ccTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAAACCTGAAAAAAACGGc                      >  1:1805384/1‑50 (MQ=39)

                       ctctACCGCTTTGCGGCAAAACTCC‑AAAAAAAAACCTGAAAAAAACGGcc                     <  1:56146/50‑1 (MQ=39)

                       ctctACCGCTTTGCGGCAAAACTCC‑AAAAAAAAACCTGAAAAAAACGGcc                     <  1:554940/50‑1 (MQ=39)

                       ctctACCGCTTTGCGGCAAAACTCC‑AAAAAAAAACCTGAAAAAAACGGcc                     <  2:164999/50‑1 (MQ=39)

                        tctACCGCTTTGCGGCAAAACTCC‑AAAAAAAAACCTGAAAAAAACGGCCt                    >  1:241799/1‑50 (MQ=39)

                          tACCGCTTTGCGGCAAAACTCC‑AAAAAAAAACCTGAAAAAAACGGCCt                    <  1:1518042/48‑1 (MQ=39)

                            ccGCTTTGCGGCAAAACTCC‑AAAAAAAAACCTGAAAAAAACGGCCTg                   <  1:1370102/47‑1 (MQ=39)

                              gCTTTGCGGCAAAACTCC‑AAAAAAAAACCTGAAAAAAACGGCCTGACGt               >  2:1505133/1‑49 (MQ=39)

                              gCTTTGCGGCAAAACTCC‑AAAAAAAAACCTGAAAAAAACGGCCTGACGt               >  2:1322300/1‑49 (MQ=39)

                              gCTTTGCGGCAAAACTCC‑AAAAAAAAACCTGAAAAAAACGGCCTGACGt               >  1:630760/1‑49 (MQ=39)

                              gCTTTGCGGCAAAACTCC‑AAAAAAAAACCTGAAAAAAACGGCCTGACGt               >  2:210314/1‑49 (MQ=39)

                              gCTTTGCGGCAAAACTCC‑AAAAAAAAACCTGAAAAAAACGGCCTGACGTg              <  2:1479977/50‑1 (MQ=39)

                              gCTTTGCGGCAAAACTCC‑AAAAAAAAACCTGAAAAAAACGGCCTGACGTg              <  2:1767155/50‑1 (MQ=39)

                              gCTTTGCGGCAAAACTCC‑AAAAAAAAACCTGAAAAAAACGGCCTGACGTg              <  2:1556374/50‑1 (MQ=39)

                               cTTTGCGGCAAAACTCC‑AAAAAAAAACCTGAAAAAAACGGCATGACGTg              >  2:83117/1‑49 (MQ=38)

                                   gCGGCAAAACTCC‑AAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATc          >  2:1993186/1‑49 (MQ=39)

                                   gCGGCAAAACTCC‑AAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATc          <  1:2061808/49‑1 (MQ=39)

                                   gCGGCAAAACTCC‑AAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATc          <  1:1386860/49‑1 (MQ=39)

                                   gCGGCAAAACTCC‑AAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATc          >  2:979155/1‑49 (MQ=39)

                                    cGGCAAAACTCC‑AAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCa         <  2:1034875/49‑1 (MQ=39)

                                     ggCAAAACTCC‑AAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAg       >  1:660324/1‑50 (MQ=255)

                                        aaaaCTCC‑AAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCaa    >  2:1673305/1‑50 (MQ=255)

                                         aaaCTCC‑AAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAAt   >  2:346014/1‑50 (MQ=255)

                                         aaaCTCC‑AAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAAt   <  1:545080/50‑1 (MQ=255)

                                          aaCTCC‑AAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAAtt  >  2:312372/1‑50 (MQ=255)

                                          aaCTCC‑AAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAAtt  <  2:357829/50‑1 (MQ=255)

GGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATT  >  NC_017635/4127888‑4127980

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: `‑`