BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,223,090	Δ4 bp	coding (71‑74/174 nt)	ymgI ←	uncharacterized protein YmgI

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1223090	1223093	4	71 [0]	[0] 65	ymgI	uncharacterized protein YmgI

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1223089	0 (0.000)	56 (0.800)	19/66	0.9	100%	coding (75/174 nt)	ymgI	uncharacterized protein YmgI
*	?	NC_000913	1223094 =	0 (0.000)	56 (0.800)	19/66	0.9	100%	coding (70/174 nt)	ymgI	uncharacterized protein YmgI

CTAAAGAACGCTCGCTATTATTCAGATAATTCAAAATGAGCGTGGCTGTGATGATAGGAATTATGTTTT > NC_000913/1223055‑1223123

gtctcttgatgctaagaactatCAGATAATTCaaa                                    >  1:1275416‑M2/23‑35 (MQ=255)

atctcttgatgctaagaactatCAGATAATTCaaa                                    >  1:1275368‑M2/23‑35 (MQ=255)

atctcttgatgctaagaactatCAGATAATTCaaa                                    >  1:1466304‑M2/23‑35 (MQ=255)

atctcttgatgctaagaactatCAGATAATTCaaa                                    <  1:4401343‑M2/13‑1 (MQ=255)

atctcttgatgctaagaactatCAGATAATTCaaa                                    >  1:2832837‑M2/23‑35 (MQ=255)

atctcttgatgctaagaactatCAGATAATTCaaa                                    >  1:2610237‑M2/23‑35 (MQ=255)

atctcttgatgctaagaactatCAGATAATTCaaa                                    >  1:4636635‑M2/23‑35 (MQ=255)

 tctcttgatgctaagaactatCAGATAATTcaaag                                   >  1:6217493‑M2/22‑34 (MQ=255)

 tctcttgatgctaagaactatCAGATAATTcaaag                                   <  1:3468191‑M2/14‑2 (MQ=255)

 tctcttgatgctaagaactatCAGATAATTcaaag                                   >  1:3422620‑M2/22‑34 (MQ=255)

                    atcAGATAATTCAAAgcgtggctgtgatgatagga                >  1:1322301‑M1/2‑15 (MQ=255)

                    atcAGATAATTCAAAgcgtggctgtgatgatagga                >  1:6098993‑M1/2‑15 (MQ=255)

                    atcAGATAATTCAAAgcgtggctgggatgatagga                >  1:1757663‑M1/2‑15 (MQ=255)

                     tCAGATAATTCAAAgcgtggctgtgatgataggaa               >  1:3499093‑M1/1‑14 (MQ=255)

                     tCAGATAATTCAAAgcgtggctgtgatgataggaa               >  1:4228867‑M1/1‑14 (MQ=255)

                     tCAGATAATTCAAAgcgtggctgtgatgataggaa               >  1:4431885‑M1/1‑14 (MQ=255)

                     tCAGATAATTCAAAgcgtggctgtgatgataggaa               >  1:4724031‑M1/1‑14 (MQ=255)

                     tCAGATAATTCAAAgcgtggctgtgatgataggaa               >  1:2109429‑M1/1‑14 (MQ=255)

                     tCAGATAATTCAAAgcgtggctgtgatgataggaa               >  1:7293419‑M1/1‑14 (MQ=255)

                     tCAGATAATTCAAAgcgtggctgtgatgataggaa               <  1:535161‑M1/35‑22 (MQ=255)

                     tCAGATAATTCAAAgcgtggctgtgatgataggaa               <  1:1163172‑M1/35‑22 (MQ=255)

                      cAGATAATTCAAAgcgtggctgtgatgataggaat              <  1:499972‑M1/35‑23 (MQ=255)

                      cAGATAATTCAAAgcgtggctgtgatgataggaat              <  1:684722‑M1/35‑23 (MQ=255)

                      cAGATAATTCAAAgcgtggctgtgatgataggaat              <  1:3789439‑M1/35‑23 (MQ=255)

                       aGATAATTCAAAgcgtggctgtgatgataggaatt             <  1:2655770‑M1/35‑24 (MQ=255)

                       aGATAATTCAAAgcgtggctgtgatgataggaatt             >  1:2296607‑M1/1‑12 (MQ=255)

                       aGATAATTCAAAgcgtggctgtgatgataggaatt             <  1:1027277‑M1/35‑24 (MQ=255)

                        gATAATTCAAAgcgtggctgtgatgataggaatta            <  1:5232949‑M1/35‑25 (MQ=255)

                        gATAATTCAAAgcgtggctgtgatgataggaatta            <  1:6194133‑M1/35‑25 (MQ=255)

                        gATAATTCAAAgcgtggctgtgatgataggaatta            <  1:543449‑M1/35‑25 (MQ=255)

                        gATAATTCAAAgcgtggctgtgatgataggaatta            <  1:3472222‑M1/35‑25 (MQ=255)

                        gATAATTCAAAgcgtggctgtgatgataggaatta            <  1:3472221‑M1/35‑25 (MQ=255)

                        gATAATTCAAAgcgtggctgtgatgataggaatta            <  1:3472155‑M1/35‑25 (MQ=255)

                        gATAATTCAAAgcgtggctgtgatgataggaatta            <  1:3071650‑M1/35‑25 (MQ=255)

                         aTAATTCAAAgcgtggctgtgatgataggaattat           <  1:7444659‑M1/35‑26 (MQ=255)

                          tAATTCAAAgcgtggctgtgatgataggaattatg          >  1:4011341‑M1/1‑9 (MQ=255)

                          tAATTCAAAgcgtggctgtgatgataggaattatg          >  1:1842236‑M1/1‑9 (MQ=255)

                          tAATTCAAAgcgtggctgtgatgataggaattatg          >  1:720952‑M1/1‑9 (MQ=255)

                          tAATTCAAAgcgtggctgtgatgataggaattatg          >  1:2525828‑M1/1‑9 (MQ=255)

                           aaTTCAAAgcgtggctgtgatgataggaattatgt         <  1:4401517‑M1/35‑28 (MQ=255)

                           aaTTCAAAgcgtggctgtgatgataggaattatgt         <  1:4670213‑M1/35‑28 (MQ=255)

                            aTTCAAAgcgtggctgtgatgataggaattatgtt        >  1:2873781‑M1/1‑7 (MQ=255)

                             ttCAAAgcgtggctgtgatgataggaattatgttt       <  1:2872815‑M1/35‑30 (MQ=255)

                             ttCAAAgcgtggctgtgatgataggaattatgttt       <  1:484903‑M1/35‑30 (MQ=255)

                             ttCAAAgcgtggctgtgatgataggaattatgttt       <  1:4179727‑M1/35‑30 (MQ=255)

                             ttCAAAgcgtggctgtgatgataggaattatgttt       <  1:4170254‑M1/35‑30 (MQ=255)

                             ttCAAAgcgtggctgtgatgataggaattatgttt       <  1:3970626‑M1/35‑30 (MQ=255)

                              tCAAAgcgtggctgtgatgataggaattatgtttt      >  1:6913812‑M1/1‑5 (MQ=255)

                              tCAAAgcgtggctgtgatgataggaattatgtttt      <  1:4563477‑M1/35‑31 (MQ=255)

                              tCAAAgcgtggctgtgatgataggaattatgtttt      <  1:6086735‑M1/35‑31 (MQ=255)

                               cAAAgcgtgggtgtgatgataggaattatgttttt     <  1:6413468‑M1/35‑32 (MQ=255)

                               cAAAgcgtggctgtgatgataggaattatgttttt     >  1:7062281‑M1/1‑4 (MQ=255)

                               cAAAgcgtggctgtgatgataggaattatgttttt     >  1:972044‑M1/1‑4 (MQ=255)

                               cAAAgcgtggctgtgatgataggaattatgttttt     <  1:1390664‑M1/35‑32 (MQ=255)

                               cAAAgcgtggctgtgatgataggaattatgttttt     <  1:2362540‑M1/35‑32 (MQ=255)

                               cAAAgcgtggctgtgatgataggaattatgttttt     <  1:3377590‑M1/35‑32 (MQ=255)

                                aaagcgtggctgtgatgataggaattatgtttttt    >  1:5377645‑M1/1‑3 (MQ=255)

                                aaagcgtggctgtgatgataggaattatgtttttt    >  1:5629380‑M1/1‑3 (MQ=255)

                                aaagcgtggctgtgatgataggaattatgtttttt    <  1:436814‑M1/35‑33 (MQ=255)

                                aaagcgtggctgtgatgataggaattatgtttttt    >  1:3521376‑M1/1‑3 (MQ=255)

                                aaagcgtggctgtgatgataggaattatgtttttt    <  1:2498296‑M1/35‑33 (MQ=255)

                                 aagcgtggctgtgatgataggaattatgtttttta   <  1:5486904‑M1/35‑34 (MQ=255)

                                 aagcgtggctgtgatgataggaattatgtttttta   <  1:7664555‑M1/35‑34 (MQ=255)

                                 aagcgtggctgtgatgataggaattatgtttttta   >  1:1328035‑M1/1‑2 (MQ=255)

                                 aagcgtggctgtgatgataggaattatgtttttta   <  1:698134‑M1/35‑34 (MQ=255)

                                 aagcgtggctgtgatgataggaattatgtttttta   <  1:2796061‑M1/35‑34 (MQ=255)

                                  agcgtggctgtgatgataggaattatgttttttac  <  1:3803411‑M1/35‑35 (MQ=255)

                                  agcgtggctgtgatgataggaattatgttttttac  >  1:3982856‑M1/1‑1 (MQ=255)

                                  agcgtggctgtgatgataggaattatgttttttac  >  1:6169481‑M1/1‑1 (MQ=255)

                                  agcgtggctgtgatgataggaattatgttttttac  <  1:4769941‑M1/35‑35 (MQ=255)

                                  agcgtggctgtgatgataggaattatgttttttac  >  1:1088334‑M1/1‑1 (MQ=255)

CTAAAGAACGCTCGCTATTATTCAGATAATTCAAAATGAGCGTGGCTGTGATGATAGGAATTATGTTTT > NC_000913/1223055‑1223123

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