BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,720,634	Δ1,442 bp	IS150‑mediated	insJ–insK	insJ, insK

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3720634	3722075	1442	41 [0]	[0] 37	insJ–insK	insJ,insK

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3720633	0 (0.000)	37 (1.000)	29/96	0.1	100%	noncoding (1/1443 nt)	IS150	repeat region
*	?	NC_000913	3722076 =	0 (0.000)	37 (1.000)	29/96	0.1	100%	intergenic (+27/+252)	insK/glyS	IS150 conserved protein InsB/glycine‑‑tRNA ligase subunit beta

CCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGG  >  NC_000913/3722030‑3722122

tatttctgcgcattcccagcctcctgtgttgatttccaacgagtataa                                               <  1:2359070‑M2/2‑1 (MQ=255)

tatttctgcgcattcccagcctcctgtgttgatttccaacgagtataa                                               <  1:270745‑M2/2‑1 (MQ=255)

tatttctgcgcattcccagcctcctgtgttgatttccaacgagtata                                                <  1:829440‑M2/1‑1 (MQ=255)

tatttctgcgcattcccagcctcctgtgttgatttccaacgagtata                                                <  1:2435992‑M2/1‑1 (MQ=255)

tatttctgcgcattcccagcctcctgtgttgatttccaacgagtatAAc                                              <  1:2651105‑M2/3‑1 (MQ=255)

 atttctgcgcattcccagcctcctgtgttgatttccaacgagtatAACt                                             >  1:2320739‑M2/46‑49 (MQ=255)

  tttctgcgcattcccagcctcctgtgttgatttccaacgagtatAACtt                                            <  1:479162‑M2/5‑1 (MQ=255)

  tttctgcgcattcccagcctcctgtgttgatttccaacgagtatAACtt                                            >  1:820723‑M2/45‑49 (MQ=255)

   ttctgcgcattcccagcctcctgtgttgatttccaacgagtatAACttt                                           >  1:1481047‑M2/44‑49 (MQ=255)

   ttctgcgcattcccagcctcctgtgttgatttccaacgagtatAACttt                                           >  1:1671213‑M2/44‑49 (MQ=255)

         gcattcccagcctcctgtgttgatttccaacgagtatAACTTTGATTta                                     <  1:1086052‑M2/12‑1 (MQ=255)

          cattcccagcctcctgtgttgatttccaacgagtatAACTTTGATTta                                     >  1:647105‑M2/37‑48 (MQ=255)

            ttcccagcctcctgtgttgatttccaacgagtatAACTTTGATTTATAg                                  <  1:1447405‑M2/15‑1 (MQ=255)

                cagcctcctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAg                              <  1:2349573‑M2/19‑1 (MQ=255)

                 agcctcctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAgg                             >  1:1562722‑M2/30‑49 (MQ=255)

                 agcctcctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAgg                             <  1:2709218‑M2/20‑1 (MQ=255)

                 agcctcctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAgg                             >  1:203046‑M2/30‑49 (MQ=255)

                  gcctcctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAGGt                            <  1:2467829‑M2/21‑1 (MQ=255)

                     tcctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAGGTggg                         >  1:2702950‑M2/26‑49 (MQ=255)

                          tgttgatttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCttt                    <  1:86270‑M2/29‑1 (MQ=255)

                          tgttgatttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCtt                     <  1:1377655‑M2/28‑1 (MQ=255)

                          tgttgatttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCt                      <  1:592992‑M2/27‑1 (MQ=255)

                          tgttgatttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCt                      <  1:512242‑M2/27‑1 (MQ=255)

                              gatttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTtct                 <  1:1758883‑M2/32‑1 (MQ=255)

                              gatttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTtct                 <  1:640644‑M2/32‑1 (MQ=255)

                               atttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTtctc                <  1:1369343‑M2/33‑1 (MQ=255)

                               atttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTTctct               <  1:2073262‑M2/34‑1 (MQ=255)

                                  tccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCtctgtc            >  1:1046981‑M2/13‑49 (MQ=255)

                                   ccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCtctgtct           >  1:2108278‑M2/12‑49 (MQ=255)

                                    caacgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCtctgtct           <  1:495921‑M2/38‑1 (MQ=255)

                                      acgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGcc        >  1:1071808‑M2/9‑49 (MQ=255)

                                         agtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCttt     <  1:2840947‑M2/44‑1 (MQ=255)

                                         agtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCttt     <  1:1998499‑M2/44‑1 (MQ=255)

                                          gtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTc    <  1:2292675‑M2/45‑1 (MQ=255)

                                          gtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTc    <  1:2517514‑M2/45‑1 (MQ=255)

                                           tatAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCg   >  1:627014‑M2/4‑49 (MQ=255)

                                             tAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCgg  <  1:1722587‑M2/47‑1 (MQ=255)

CCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGG  >  NC_000913/3722030‑3722122

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