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

TATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATA  >  NC_000913/3720585‑3720680

tataTTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt                                                 >  1:1895003/1‑48 (MQ=255)

  tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt                                                 <  1:2500732/46‑1 (MQ=255)

  tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt                                                 <  1:1779726/46‑1 (MQ=255)

  tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt                                                 <  1:1621347/46‑1 (MQ=255)

  tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac                                              <  1:2651105‑M1/49‑4 (MQ=255)

  tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa                                               <  1:2359070‑M1/48‑3 (MQ=255)

  tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa                                               <  1:270745‑M1/48‑3 (MQ=255)

  tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATa                                                <  1:829440‑M1/47‑2 (MQ=255)

  tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATa                                                <  1:2435992‑M1/47‑2 (MQ=255)

   atTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact                                             >  1:2320739‑M1/1‑45 (MQ=255)

    tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt                                            <  1:479162‑M1/49‑6 (MQ=255)

    tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt                                            >  1:820723‑M1/1‑44 (MQ=255)

     ttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt                                           >  1:1671213‑M1/1‑43 (MQ=255)

     ttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt                                           >  1:1481047‑M1/1‑43 (MQ=255)

            gCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttgattta                                     <  1:1086052‑M1/49‑13 (MQ=255)

             cATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttgattta                                     >  1:647105‑M1/1‑36 (MQ=255)

               ttCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttgatttatag                                  <  1:1447405‑M1/49‑16 (MQ=255)

                   cAGCCTCCTGTGTTGATTTCCAACGAGTATaactttgatttatagtcag                              <  1:2349573‑M1/49‑20 (MQ=255)

                    aGCCTCCTGTGTTGATTTCCAACGAGTATaactttgatttatagtcagg                             <  1:2709218‑M1/49‑21 (MQ=255)

                    aGCCTCCTGTGTTGATTTCCAACGAGTATaactttgatttatagtcagg                             >  1:1562722‑M1/1‑29 (MQ=255)

                    aGCCTCCTGTGTTGATTTCCAACGAGTATaactttgatttatagtcagg                             >  1:203046‑M1/1‑29 (MQ=255)

                     gCCTCCTGTGTTGATTTCCAACGAGTATaactttgatttatagtcaggt                            <  1:2467829‑M1/49‑22 (MQ=255)

                        tcctGTGTTGATTTCCAACGAGTATaactttgatttatagtcaggtggg                         >  1:2702950‑M1/1‑25 (MQ=255)

                             tgtTGATTTCCAACGAGTATaactttgatttatagtcaggtggggcttt                    <  1:86270‑M1/49‑30 (MQ=255)

                             tgtTGATTTCCAACGAGTATaactttgatttatagtcaggtggggctt                     <  1:1377655‑M1/48‑29 (MQ=255)

                             tgtTGATTTCCAACGAGTATaactttgatttatagtcaggtggggct                      <  1:592992‑M1/47‑28 (MQ=255)

                             tgtTGATTTCCAACGAGTATaactttgatttatagtcaggtggggct                      <  1:512242‑M1/47‑28 (MQ=255)

                                 gATTTCCAACGAGTATaactttgatttatagtcaggtggggcttttct                 <  1:1758883‑M1/48‑33 (MQ=255)

                                 gATTTCCAACGAGTATaactttgatttatagtcaggtggggcttttct                 <  1:640644‑M1/48‑33 (MQ=255)

                                  aTTTCCAACGAGTATaactttgatttatagtcaggtggggcttttctct               <  1:2073262‑M1/49‑35 (MQ=255)

                                  aTTTCCAACGAGTATaactttgatttatagtcaggtggggcttttctc                <  1:1369343‑M1/48‑34 (MQ=255)

                                     tCCAACGAGTATaactttgatttatagtcaggtggggcttttctctgtc            >  1:1046981‑M1/1‑12 (MQ=255)

                                      ccAACGAGTATaactttgatttatagtcaggtggggcttttctctgtct           >  1:2108278‑M1/1‑11 (MQ=255)

                                       cAACGAGTATaactttgatttatagtcaggtggggcttttctctgtct           <  1:495921‑M1/48‑39 (MQ=255)

                                         aCGAGTATaactttgatttatagtcaggtggggcttttctctgtctgcc        >  1:1071808‑M1/1‑8 (MQ=255)

                                            aGTATaactttgatttatagtcaggtggggcttttctctgtctgccttt     <  1:2840947‑M1/49‑45 (MQ=255)

                                            aGTATaactttgatttatagtcaggtggggcttttctctgtctgccttt     <  1:1998499‑M1/49‑45 (MQ=255)

                                             gTATaactttgatttatagtcaggtggggcttttctctgtctgcctttc    <  1:2292675‑M1/49‑46 (MQ=255)

                                             gTATaactttgatttatagtcaggtggggcttttctctgtctgcctttc    <  1:2517514‑M1/49‑46 (MQ=255)

                                              tataactttgatttatagtcaggtggggcttttctctgtctgcctttcg   >  1:627014‑M1/1‑3 (MQ=255)

                                                taactttgatttatagtcaggtggggcttttctctgtctgcctttcgg  <  1:1722587‑M1/48‑48 (MQ=255)

TATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATA  >  NC_000913/3720585‑3720680

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