BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	263,514	Δ61,897 bp		[ykfI]–ykgH	71 genes [ykfI], yafW, ykfH, ykfG, yafX, ykfF, ykfB, yafY, yafZ, ykfA, perR, insN, insI‑1, insN, eyeA, ykfC, insH‑1, mmuP, mmuM, afuC, afuB, insB‑2, insA‑2, yagB, yagA, yagE, yagF, yagG, yagH, xynR, argF, ykgS, insB‑3, insA‑3, yagJ, yagK, yagL, yagM, ykgV, yagN, intF, yagP, paoD, paoC, paoB, paoA, yagU, ykgJ, ecpE, ecpD, ecpC, ecpB, ecpA, ecpR, ykgL, ykgO, ykgM, ykgR, ykgP, eaeH, insE‑1, insF‑1, ykgA, rclC, rclB, rclA, rclR, ykgE, ykgF, ykgG, ykgH

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	263514	325410	61897	31 [0]	[0] 32	[ykfI]–ykgH	[ykfI],yafW,ykfH,ykfG,yafX,ykfF,ykfB,yafY,yafZ,ykfA,perR,insN,insI‑1,insN,eyeA,ykfC,insH‑1,mmuP,mmuM,afuC,afuB,insB‑2,insA‑2,yagB,yagA,yagE,yagF,yagG,yagH,xynR,argF,ykgS,insB‑3,insA‑3,yagJ,yagK,yagL,yagM,ykgV,yagN,intF,yagP,paoD,paoC,paoB,paoA,yagU,ykgJ,ecpE,ecpD,ecpC,ecpB,ecpA,ecpR,ykgL,ykgO,ykgM,ykgR,ykgP,eaeH,insE‑1,insF‑1,ykgA,rclC,rclB,rclA,rclR,ykgE,ykgF,ykgG,ykgH

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 263513	0 (0.000)	31 (1.060)	24/92	0.1	100%	coding (157/342 nt)	ykfI	CP4‑6 prophage; toxin of the YkfI‑YafW toxin‑antitoxin system
*	?	NC_000913	325411 =	0 (0.000)	31 (1.060)	24/92	0.1	100%	intergenic (‑47/+166)	ykgH/betA	uncharacterized protein YkgH/choline dehydrogenase

TTCTACCAGAAAATTCACGGCATCGGCTAGGGTGATACCGGCATCGAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/263466‑263513

‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCG  >  NC_000913/325411‑325455

TTCTACCAGAAAATTCACGGCATCGGCTAGGGTGATACCGGCATCGA                                                >  1:1225617/1‑47

 TCTACCAGAAAATTCACGGCATCGGCTAGGGTGATACCGGCATCGATA                                              <  1:570442/48‑1

   TACCAGAAAATTCACGGCATCGGCTAGGGTGATACCGGCATCGATAGTC                                           >  1:2023856/1‑49

     CCAGAAAATTCACGGCATCGGCTAGGGTGATACCGGCATCGATAGT                                            <  1:857484/46‑1

     CCAGAAAATTCACGGCATCGGCTAGGGTGATACCGGCATCGATAGT                                            <  1:281946/46‑1

     CCAGAAAATTCACGGCATCGGCTAGGGTGATACCGGCATCGATAGT                                            <  1:376646/46‑1

      CAGAAAATTCACGGCATCGGCTAGGGTGATACCGGCATCGATAGTCTG                                         >  1:1334359/1‑48

       AGAAAATTCACGGCATCGGCTAGGGTGATACCGGCATCGATAGTCTGTA                                       >  1:184628/1‑49

       AGAAAATTCACGGCATCGGCTAGGGTGATACCGGCATCGATAGTCTGTA                                       >  1:361080/1‑49

          AAATTCACGGCATCGGCTAGGGTGATACCGGCATCGATAGTCTGTAT                                      >  1:1955439/1‑47

          AAATTCACGGCATCGGCTAGGGTGATACCGGCATCGATAGTCTGTAT                                      >  1:782753/1‑47

           AATTCACGGCATCGGCTAGGGTGATACCGGCATCGATAGTCTGTATCAG                                   <  1:1893618/49‑1

               CACGGCATCGGCTAGGGTGATACCGGCATCGATAGTCTGTATCAGGAAT                               >  1:1800483/1‑49

                 CGGCATCGGCTAGGGTGATACCGGCATCGATAGTCTGTATCAGGAATGT                             <  1:705937/49‑1

                 CGGCATCGGCTAGGGTGATACCGGCATCGATAGTCTGTATCAGGAATGT                             <  1:397478/49‑1

                   GCATCGGCTAGGGTGATACCGGCATCGATAGTCTGTATCAGGAATGTTC                           <  1:741986/49‑1

                     ATCGGCTAGGGTGATACCGGCATCGATAGTCTGTATCAGGAATGTTCGG                         >  1:1385022/1‑49

                     ATCGGCTAGGGTGATACCGGCATCGATAGTCTGTATCAGGAATGTTCGG                         <  1:1257103/49‑1

                         GCTAGGGTGATACCGGCATCGATAGTCTGTATCAGGAATGTTCGGGTTA                     <  1:1340560/49‑1

                         GCTAGGGTGATACCGGCATCGATAGTCTGTATCAGGAATGTTCGGGTT                      >  1:2292222/1‑48

                          CTAGGGTGATACCGGCATCGATAGTCTGTATCAGGAATGTTCGGGTTAA                    <  1:1796226/49‑1

                            AGGGTGATACCGGCATCGATAGTCTGTATCAGGAATGTTCGGGTTAAA                   >  1:834330/1‑48

                             GGGTGATACCGGCATCGATAGTCTGTATCAGGAATGTTCGGGTTAAATA                 <  1:1331095/49‑1

                              GGTGATACCGGCATCGATAGTCTGTATCAGGAATGTTCGGGTTAAATAT                <  1:766916/49‑1

                               GTGATACCGGCATCGATAGTCTGTATCAGGAATGTTCGGGTTAAATAT                >  1:1861743/1‑48

                               GTGATACCGGCATCGATAGTCTGTATCAGGAATGTTCGGGTTAAATAT                >  1:954100/1‑48

                                 GATACCGGCATCGATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAG             >  1:62060/1‑49

                                     CCGGCATCGATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAA         <  1:44542/49‑1

                                     CCGGCATCGATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAA         <  1:1623370/49‑1

                                     CCGGCATCGATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAA          <  1:1988409/48‑1

                                            CGATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCG  >  1:231009/1‑49

                                             GATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCG  <  1:1453281/48‑1

TTCTACCAGAAAATTCACGGCATCGGCTAGGGTGATACCGGCATCGAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/263466‑263513

‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCG  >  NC_000913/325411‑325455

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.