BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	2,466,545	Δ9,634 bp		[intS]–[ypdJ]	15 genes [intS], yfdG, yfdH, yfdI, yfdK, yfdL, yfdM, yfdN, yfdO, yfdP, yfdQ, yfdR, yfdS, [yfdT], [ypdJ]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	2466545	2476178	9634	35 [0]	[0] 34	[intS]–[ypdJ]	[intS],yfdG,yfdH,yfdI,yfdK,yfdL,yfdM,yfdN,yfdO,yfdP,yfdQ,yfdR,yfdS,[yfdT],[ypdJ]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2466544	0 (0.000)	34 (1.240)	29/98	0.0	100%	intergenic (+161/‑1)	argW/intS	tRNA‑Arg/CPS‑53 (KpLE1) prophage; prophage CPS‑53 integrase
*	?	NC_000913	2476179 =	0 (0.000)	34 (1.240)	29/98	0.0	100%	coding (86/141 nt)	ypdJ	CPS‑53 (KpLE1) prophage; putative uncharacterized protein YpdJ

ATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAA  >  NC_000913/2466495‑2466589

aTGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTAccc                                               <  1:619971/50‑1 (MQ=255)

    cGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCta                                             >  1:1190752‑M1/1‑46 (MQ=255)

     gggTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCtagag                                          >  1:1313269‑M1/1‑45 (MQ=255)

        tAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCtagagcca                                       >  1:1839009‑M1/1‑42 (MQ=255)

        tAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCtagagcca                                       <  1:786036‑M1/50‑9 (MQ=255)

         aaaaaGTGGATAAAATAATTTTACCCACCGGATTTTTACCCtagagcca                                       >  1:1152919‑M1/1‑41 (MQ=255)

            aaGTGGATAAAATAATTTTACCCACCGGATTTTTACCCtagagccataca                                   <  1:2064637‑M1/50‑13 (MQ=255)

            aaGTGGATAAAATAATTTTACCCACCGGATTTTTACCCtagagccata                                     >  1:1703427‑M1/1‑38 (MQ=255)

             aGTGGATAAAATAATTTTACCCACCGGATTTTTACCCtagagccatacag                                  >  1:1298512‑M1/1‑37 (MQ=255)

               tGGATAAAATAATTTTACCCACCGGATTTTTACCCtagagccatacagcc                                >  1:442364‑M1/1‑35 (MQ=255)

                 gATAAAATAATTTTACCCACCGGATTTTTACCCtagagccatacagcct                               >  1:499545‑M1/1‑33 (MQ=255)

                  aTAAAATAATTTTACCCACCGGATTTTTACCCtagagccatacagcctca                             >  1:761437‑M1/1‑32 (MQ=255)

                   tAAAATAATTTTACCCACCGGATTTTTACCCtagagccatacagcctca                             >  1:160116‑M1/1‑31 (MQ=255)

                       ataatTTTACCCACCGGATTTTTACCCtagagccatacagcctcacactc                        >  1:917253‑M1/1‑27 (MQ=255)

                          atTTTACCCACCGGATTTTTACCCtagagccatacagcctcacactcgat                     >  1:149373‑M1/1‑24 (MQ=255)

                          atTTTACCCACCGGATTTTTACCCtagagccatacagcctcacactcgat                     <  1:1143267‑M1/50‑27 (MQ=255)

                           ttttACCCACCGGATTTTTACCCtagagccatacagcctcacactcgat                     >  1:666769‑M1/1‑23 (MQ=255)

                              tACCCACCGGATTTTTACCCtagagccatacagcctcacactcgatgag                  <  1:1112350‑M1/49‑30 (MQ=255)

                               aCCCACCGGATTTTTACCCtagagccatacagcctcacactcgatgaggc                >  1:1558236‑M1/1‑19 (MQ=255)

                               aCCCACCGGATTTTTACCCtagagccatacagcctcacactcgatgaggc                >  1:1722862‑M1/1‑19 (MQ=255)

                               aCCCACCGGATTTTTACCCtagagccatacagcctcacactcgatgaggc                >  1:1462281‑M1/1‑19 (MQ=255)

                                 ccACCGGATTTTTACCCtagagccatacagcctcacactcgatgaggcct              <  1:732971‑M1/50‑34 (MQ=255)

                                  cACCGGATTTTTACCCtagagccatacagcctcacactcgatgaggcctg             >  1:94452‑M1/1‑16 (MQ=255)

                                   aCCGGATTTTTACCCtagagccatacagcctcacactcgatgaggcctgt            >  1:1477760‑M1/1‑15 (MQ=255)

                                    ccGGATTTTTACCCtagagccatacagcctcacactcgatgaggcctgtc           <  1:177977‑M1/50‑37 (MQ=255)

                                      ggATTTTTACCCtagagccatacagcctcacactcgatgaggcctgtcag         >  1:1615977‑M1/1‑12 (MQ=255)

                                       gATTTTTACCCtagagccatacagcctcacactcgatgaggcctgtcagt        >  1:461398‑M1/1‑11 (MQ=255)

                                         tttttACCCtagagccatacagcctcacactcgatgaggcctgtcagttt      >  1:2158654‑M1/1‑9 (MQ=255)

                                         tttttACCCtagagccatacagcctcacactcgatgaggcctgtcagttt      <  1:2083458‑M1/50‑42 (MQ=255)

                                         tttttACCCtagagccatacagcctcacactcgatgaggcctgtcagttt      <  1:2080939‑M1/50‑42 (MQ=255)

                                          ttttACCCtagagccatacagcctcacactcgatgaggcctgtcagtttc     >  1:1630168‑M1/1‑8 (MQ=255)

                                           tttACCCtagagccatacagcctcacactcgatgaggcctgtcagtttc     <  1:1336914‑M1/49‑43 (MQ=255)

                                              aCCCtagagccatacagcctcacactcgatgaggcctgtcagtttctta  >  1:2296269‑M1/1‑4 (MQ=255)

                                              aCCCtagagccatacagcctcacactcgatgaggcctgtcagtttctta  >  1:785987‑M1/1‑4 (MQ=255)

                                              aCCCtagagccatacagcctcacactcgataaggcctgtcagtttctta  >  1:1774757‑M1/1‑4 (MQ=255)

ATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAA  >  NC_000913/2466495‑2466589

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.