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 (1.010)	23/96	0.1	100%	coding (1/1242 nt)	intA	CP4‑57 prophage; integrase
*	?	NC_000913	2791251 =	0 (0.000)	27 (1.010)	23/96	0.1	100%	intergenic (+1/‑22)	lhgO/gabD	L‑2‑hydroxyglutarate oxidase/NADP(+)‑dependent succinate‑semialdehyde dehydrogenase

AACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2756109‑2756159

‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTT  >  NC_000913/2791251‑2791295

AACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCA                                               <  1:1745788/50‑1

  CGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCACC                                             >  1:538304/1‑50

   GAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCACCT                                            <  1:2129800/50‑1

    AAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCACCTT                                           >  1:754986/1‑50

     AAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCACCTTT                                          <  1:1821415/50‑1

      AAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCACCTTT                                          >  1:828661/1‑49

        GCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCACCTTTGA                                        >  1:1877224/1‑49

         CTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCACCTTTGAA                                       >  1:1714579/1‑49

             ATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCACCTTTGAAAACAG                                  >  1:845882/1‑50

              TTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCACCTTTGAAAACAGG                                 <  1:1070461/50‑1

                     CGTTGGG‑CTAAACGCAGGGAGACTCCCCACCTTTGAAAACAGGATGTAGC                          >  1:1281510/1‑50

                     CGTTGGG‑CTAAACGCAGGGAGACTCCCCACCTTTGAAAACAGGATGTAG                           <  1:586313/49‑1

                      GTTGGG‑CTAAACGCAGGGAGACTCCCCACCTTTGAAAACAGGATGTAGCG                         <  1:1318174/50‑1

                      GTTGGG‑CTAAACGCAGGGAGACTCCCCACCTTTGAAAACAGGAGGTAGCG                         <  1:89542/50‑1

                           ggCTAAACGCAGGGAGACTCCCCACCTTTGAAAACAGGATGTAGCGATG                      <  1:1270533/47‑1

                           ggCTAAACGCAGGGAGACTCCCCACCTTTGAAAACAGGATGTAGCGATG                      <  1:1334403/47‑1

                              TAAACGCAGGGAGACTCCCCACCTTTGAAAACAGGATGTAGCGATGAAAC                  >  1:1982750/1‑50

                                   GCAGGGAGACTCCCCACCTTTGAAAACAGGATGTAGCGATGAAACTTAAC             <  1:1852320/50‑1

                                    CAGGGAGACTCCCCACCTTTGAAAACAGGATGTAGCGATGAAACTTAACG            <  1:36632/50‑1

                                    CAGGGAGACTCCCCACCTTTGAAAACAGGATGTAGCGATGAAACTTAAC             >  1:424384/1‑49

                                     AGGGAGACTCCCCACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGA           >  1:364427/1‑50

                                       GGAGACTCCCCACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACA         <  1:1776818/50‑1

                                        GAGACTCCCCACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAG        <  1:901472/50‑1

                                         AGACTCCCCACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAG        >  1:275955/1‑49

                                          GACTCCCCACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTA      <  1:1851782/50‑1

                                          GACTCCCCACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTA      <  1:1442157/50‑1

                                          GACTCCCCACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGT       >  1:2162910/1‑49

                                               CCCACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTT  <  1:451054/49‑1

AACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2756109‑2756159

‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTT  >  NC_000913/2791251‑2791295

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.