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	38 [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.220)	25/98	0.1	100%	noncoding (1/1443 nt)	IS150	repeat region
*	?	NC_000913	3722076 =	0 (0.000)	37 (1.220)	25/98	0.1	100%	intergenic (+27/+252)	insK/glyS	IS150 conserved protein InsB/glycine‑‑tRNA ligase subunit beta

ATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATAT  >  NC_000913/3720586‑3720681

atatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac                                               >  1:2700254‑M1/1‑47 (MQ=255)

atatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac                                               >  1:414739‑M1/1‑47 (MQ=255)

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

 tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact                                              <  1:966256‑M1/50‑5 (MQ=255)

 tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact                                              >  1:1867749‑M1/1‑46 (MQ=255)

 tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact                                              <  1:1299270‑M1/50‑5 (MQ=255)

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

  atTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt                                             >  1:1615825‑M1/1‑45 (MQ=255)

   tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt                                            >  1:2010788‑M1/1‑44 (MQ=255)

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

     tCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttga                                          >  1:342087‑M1/1‑42 (MQ=255)

     tCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttg                                           <  1:2200624‑M1/49‑8 (MQ=255)

     tCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttg                                           <  1:2695457‑M1/49‑8 (MQ=255)

          cgcATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttgattta                                      >  1:179839‑M1/2‑38 (MQ=255)

          cgcATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttgattt                                       <  1:2211026‑M1/48‑12 (MQ=255)

          cgcATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttgat                                         <  1:698119‑M1/46‑10 (MQ=255)

           gCATTCCCAGCCTCCTGTATTGATTTCCAACGAGTATaactttgatttat                                     <  1:2284289‑M1/50‑14 (MQ=255)

             aTTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttgatttatag                                   >  1:1444915‑M1/1‑35 (MQ=255)

             aTTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttgatttatag                                   >  1:962591‑M1/1‑35 (MQ=255)

                  cAGCCTCCTGTGTTGATTTCCAACGAGTATaactttgatttatagtcagg                              >  1:1190511‑M1/1‑30 (MQ=255)

                   aGCCTCCTGTGTTGATTTCCAACGAGTATaactttgatttatagtcaggt                             >  1:1535691‑M1/1‑29 (MQ=255)

                   aGCCTCCTGTGTTGATTTCCAACGAGTATaactttgatttatagtcaggt                             <  1:53301‑M1/50‑22 (MQ=255)

                     cctcctGTGTTGATTTCCAACGAGTATaactttgatttatagtcaggtgg                           <  1:2189481‑M1/50‑24 (MQ=255)

                      ctcctGTGTTGATTTCCAACGAGTATaactttgatttatagtcaggtggg                          >  1:458475‑M1/1‑26 (MQ=255)

                              ttGATTTCCAACGAGTATaactttgatttatagtcaggtggggcttttct                  <  1:2600200‑M1/50‑33 (MQ=255)

                              ttGATTTCCAACGAGTATaactttgatttatagtcaggtggggcttttct                  <  1:1969014‑M1/50‑33 (MQ=255)

                                gATTTCCAACGAGTATaactttgatttatagtcaggtggggcttttctct                >  1:41919‑M1/1‑16 (MQ=255)

                                gATTTCCAACGAGTATaactttgatttatagtcaggtggggcttttctc                 <  1:1173111‑M1/49‑34 (MQ=255)

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

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

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

                                     ccAACGAGTATaactttgatttatagtcaggtggggcttttctctgtctg           >  1:1366970‑M1/1‑11 (MQ=255)

                                      cAACGAGTATaactttgatttatagtcaggtggggcttttctctgtctgc          >  1:905850‑M1/1‑10 (MQ=255)

                                        aCGAGTATaactttgatttatagtcaggtggggcttttctctgtctgcct        <  1:1799204‑M1/50‑43 (MQ=255)

                                        aCGAGTATaactttgatttatagtcaggtggggcttttctctgtctgcct        <  1:1377806‑M1/50‑43 (MQ=255)

                                             tataactttgatttatagtcaggtggggcttttctctgtctgcctttcgg   >  1:2186465‑M1/1‑3 (MQ=255)

                                             tataactttgatttatagccaggtggggcttttctctgtctgcctttcgg   >  1:112377‑M1/1‑3 (MQ=255)

                                              ataactttgatttatagtcaggtggggcttttctctgtctgcctttcggt  >  1:1113623‑M1/1‑2 (MQ=255)

ATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATAT  >  NC_000913/3720586‑3720681

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