BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,197,137	25956 bp→CAGATAATTCAA		[ymfD]–[ymgI]	42 genes [ymfD], ymfE, lit, intE, xisE, ymfH, ymfI, ymfJ, ymfK, ymfT, ymfL, ymfM, oweE, ymfN, aaaE, ymfR, beeE, jayE, ymfQ, ycfK, tfaP, tfaE, stfE, pinE, mcrA, icdC, iraM, ymgK, ymgL, ycgX, bluR, bluF, ycgZ, ymgA, ariR, ymgC, pdeG, ymgF, ycgH, ymgD, ymgG, [ymgI]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1197137	1223092	25956	43 [0]	[0] 45	[ymfD]–[ymgI]	[ymfD],ymfE,lit,intE,xisE,ymfH,ymfI,ymfJ,ymfK,ymfT,ymfL,ymfM,oweE,ymfN,aaaE,ymfR,beeE,jayE,ymfQ,ycfK,tfaP,tfaE,stfE,pinE,mcrA,icdC,iraM,ymgK,ymgL,ycgX,bluR,bluF,ycgZ,ymgA,ariR,ymgC,pdeG,ymgF,ycgH,ymgD,ymgG,[ymgI]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1197136	0 (0.000)	32 (1.160) +CAGATAATTCAA	23/72	0.0	100%	coding (397/666 nt)	ymfD	e14 prophage; putative SAM‑dependent methyltransferase
*	?	NC_000913	1223093 =	0 (0.000)	32 (1.160) +CAGATAATTCAA	23/72	0.0	100%	coding (71/174 nt)	ymgI	uncharacterized protein YmgI

CTGCTTAACTAAAGAACGCTCGCTATTATTCAGATAATTCAAAATGAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGAATA  >  NC_000913/1223047‑1223140

ctgattttagtaatctcttgatgctaagaactatcagataattcaaa                                                 <  1:1388785‑M2/1‑1 (MQ=255)

   attttagtaatctcttgatgctaagaactatcagataattcaaAGCGt                                             <  1:687104‑M2/5‑1 (MQ=255)

   attttagtaatctcttgatgctaagaactatcagataattcaaAGCGTg                                            <  1:1447437‑M2/6‑1 (MQ=255)

    ttttagtaatctcttgatgctaagaactatcagataattcaaAGCGTgg                                           >  1:2708773‑M2/43‑49 (MQ=255)

      ttagtaatctcttgatgctaagaactatcagataattcaaAGCGTGGc                                          <  1:2769671‑M2/8‑1 (MQ=255)

      ttagtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCt                                         >  1:386525‑M2/41‑49 (MQ=255)

      ttagtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCt                                         >  1:293813‑M2/41‑49 (MQ=255)

       tagtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCtg                                        >  1:1953105‑M2/40‑49 (MQ=255)

         gtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCtgtg                                      <  1:1277104‑M2/12‑1 (MQ=255)

         gtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCtgtg                                      <  1:2112812‑M2/12‑1 (MQ=255)

          taatctcttgatgctaagaactatcagataattcaaAGCGTGGCTGtga                                     <  1:2245688‑M2/13‑1 (MQ=255)

           aatctcttgatgctaagaactatcagataattcaaAGCGTGGCTGtgat                                    >  1:913397‑M2/36‑49 (MQ=255)

           aatctcttgatgctaagaactatcagataattcaaAGCGTGGCTGtgat                                    >  1:564084‑M2/36‑49 (MQ=255)

               tcttgatgctaagaactatcagataattcaaAGCGTGGCTGtgatga                                  >  1:391730‑M2/32‑47 (MQ=255)

                 ttgatgctaagaactatcagataattcaaAGCGTGGCTGtgatga                                  <  1:2226102‑M2/16‑1 (MQ=255)

                 ttgatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATa                                <  1:2524517‑M2/18‑1 (MQ=255)

                 ttgatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATa                                <  1:2516098‑M2/18‑1 (MQ=255)

                  tgatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAgg                              >  1:1719237‑M2/29‑48 (MQ=255)

                  tgatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAgg                              >  1:2732987‑M2/29‑48 (MQ=255)

                   gatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGa                             <  1:1074994‑M2/21‑1 (MQ=255)

                   gatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGa                             <  1:2822269‑M2/21‑1 (MQ=255)

                    atgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAAt                           <  1:2502896‑M2/23‑1 (MQ=255)

                    atgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAAt                           <  1:205939‑M2/23‑1 (MQ=255)

                     tgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAAtt                          <  1:2634855‑M2/24‑1 (MQ=255)

                       ctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTAt                        >  1:2205239‑M2/24‑49 (MQ=255)

                       ctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTAt                        >  1:901816‑M2/24‑49 (MQ=255)

                       ctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTAt                        >  1:314863‑M2/24‑49 (MQ=255)

                          agaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATg                       <  1:252166‑M2/27‑1 (MQ=255)

                          agaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGtt                     >  1:815007‑M2/21‑49 (MQ=255)

                            aactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGtttt                   <  1:2585750‑M2/31‑1 (MQ=255)

                              ctatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGttttt                  >  1:866601‑M2/17‑48 (MQ=255)

                                 tcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACg              <  1:2445574‑M2/36‑1 (MQ=255)

                                   agataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGt             >  1:216550‑M2/12‑48 (MQ=255)

                                     ataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGaa          >  1:1116585‑M2/10‑49 (MQ=255)

                                     ataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGaa          <  1:1439259‑M2/40‑1 (MQ=255)

                                        attcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATga       <  1:770166‑M2/43‑1 (MQ=255)

                                         ttcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAAt         >  1:252586‑M2/6‑46 (MQ=255)

                                         ttcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAAt         >  1:777084‑M2/6‑46 (MQ=255)

                                           caaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATgaga     <  1:963352‑M2/45‑1 (MQ=255)

                                           caaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGaa    >  1:2543792‑M2/4‑49 (MQ=255)

                                            aaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGaa    >  1:2142621‑M2/3‑48 (MQ=255)

                                            aaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGAAt   <  1:891397‑M2/47‑1 (MQ=255)

                                             aAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGAAt   >  1:1804551‑M2/2‑48 (MQ=255)

                                             aAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGAATa  <  1:334216‑M2/48‑1 (MQ=255)

                                             aAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGAATa  <  1:258930‑M2/48‑1 (MQ=255)

CTGCTTAACTAAAGAACGCTCGCTATTATTCAGATAATTCAAAATGAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGAATA  >  NC_000913/1223047‑1223140

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 34 ≤` `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.