BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	2,756,160	Δ35,091 bp		[intA]–[lhgO]	33 genes [intA], yfjH, alpA, yfjI, yfjJ, abpB, abpA, yfjM, rnlA, rnlB, yfjP, yfjQ, yfjR, ypjK, yfjS, yfjT, yfjU, yfjV, yfjW, ypjI, yfjX, yfjY, ypjJ, yfjZ, ypjF, ypjA, pinH, ypjB, ypjC, ileY, ygaQ, csiD, [lhgO]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	2756160	2791250	35091	28 [0]	[0] 29	[intA]–[lhgO]	[intA],yfjH,alpA,yfjI,yfjJ,abpB,abpA,yfjM,rnlA,rnlB,yfjP,yfjQ,yfjR,ypjK,yfjS,yfjT,yfjU,yfjV,yfjW,ypjI,yfjX,yfjY,ypjJ,yfjZ,ypjF,ypjA,pinH,ypjB,ypjC,ileY,ygaQ,csiD,[lhgO]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2756159	0 (0.000)	27 (0.540)	23/96	0.6	100%	coding (1/1242 nt)	intA	CP4‑57 prophage; integrase
*	?	NC_000913	2791251 =	0 (0.000)	27 (0.540)	23/96	0.6	100%	intergenic (+1/‑22)	lhgO/gabD	L‑2‑hydroxyglutarate oxidase/NADP(+)‑dependent succinate‑semialdehyde dehydrogenase

GAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2756112‑2756159

‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTC  >  NC_000913/2791251‑2791299

GAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCACCT                                                <  1:2900219/50‑1

GAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCACCT                                                <  1:2895201/50‑1

 AAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCACCTT                                               <  1:827002/50‑1

     GCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCACCTTTGAA                                           <  1:2208699/50‑1

         TATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCACCTTTGAAAA                                         <  1:573473/48‑1

           TTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCACCTTTGAAAACAG                                      <  1:677844/49‑1

           TTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCACCTTTGAAAAC                                        <  1:1503126/47‑1

              TTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCACCTTTGAAAACAGGATG                                  <  1:2487391/50‑1

                  CGTTGGG‑CTAAACGCAGGGAGACTCCCCACCTTTGAAAACAGGATGTAGC                              <  1:2592615/50‑1

                  CGTTGGG‑CTAAACGCAGGGAGACTCCCCACCTTTGAAAACAGGATGTAG                               <  1:1667151/49‑1

                             AACGCAGGGAGACTCCCCACCTTTGAAAACAGGATGTAGCGATGAAACTT                    <  1:905192/50‑1

                             AACGCAGGGAGACTCCCCACCTTTGAAAACAGGATGTAGCGATGAAACTT                    <  1:3250342/50‑1

                             AACGCAGGGAGACTCCCCACCTTTGAAAACAGGATGTAGCGATGAAACT                     <  1:1908074/49‑1

                             AACGCAGGGAGACTCCCCACCTTTGAAAACAGGATGTAGCGATGAAAC                      <  1:2351681/48‑1

                              ACGCAGGGAGACTCCCCACCTTTGAAAACAGGATGTAGCGATGAAACTTA                   <  1:1603226/50‑1

                               CGCAGGGAGACTCCCCACCTTTGAAAACAGGATGTAGCGATGAAACTTAA                  <  1:3738266/50‑1

                                 CAGGGAGACTCCCCACCTTTGAAAACAGGATGTAGCGATGAAACTTAACG                <  1:1442834/50‑1

                                 CAGGGAGACTCCCCACCTTTGAAAACAGGATGTAGCGATGAAACTTAACG                <  1:1383992/50‑1

                                  AGGGAGACTCCCCACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGA               <  1:3506558/50‑1

                                    GGAGACTCCCCACCTTTGAAAACAGGATGTAGCGATGAAACTTATCGACA             <  1:3294611/50‑1

                                     GAGACTCCCCACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAG            <  1:3466151/50‑1

                                        ACTCCCCACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAA         <  1:3849145/50‑1

                                         CTCCCCACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAAC        >  1:3600557/1‑50

                                            CCCACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTT      <  1:2402490/49‑1

                                            CCCACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACT       <  1:272990/48‑1

                                            CCCACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACT       <  1:1832627/48‑1

                                              CACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATT   <  1:1654827/50‑1

                                               ACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTC  <  1:3006005/50‑1

GAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2756112‑2756159

‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTC  >  NC_000913/2791251‑2791299

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.