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] 35	[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 (0.910)	28/98	0.1	100%	intergenic (+161/‑1)	argW/intS	tRNA‑Arg/CPS‑53 (KpLE1) prophage; prophage CPS‑53 integrase
*	?	NC_000913	2476179 =	0 (0.000)	34 (0.910)	28/98	0.1	100%	coding (86/141 nt)	ypdJ	CPS‑53 (KpLE1) prophage; putative uncharacterized protein YpdJ

TGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCAT  >  NC_000913/2466496‑2466593

tGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTAccc                                                   >  1:390556/1‑49 (MQ=255)

tGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCt                                                  <  1:2195401‑M1/50‑2 (MQ=255)

  aCGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCtag                                                <  1:1765072‑M1/50‑4 (MQ=255)

   cGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCtaga                                               <  1:2295057‑M1/50‑5 (MQ=255)

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

     ggTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCtagagc                                             >  1:1288056‑M1/1‑44 (MQ=255)

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

         aaaaGTGGATAAAATAATTTTACCCACCGGATTTTTACCCtagagccata                                         >  1:1920613‑M1/1‑40 (MQ=255)

          aaaGTGGATAAAATAATTTTACCCACCGGATTTTTACCCtagagccata                                         >  1:2296854‑M1/1‑39 (MQ=255)

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

             gTGGATAAAATAATTTTACCCACCGGATTTTTACCCtagagccatacagc                                     >  1:1073900‑M1/1‑36 (MQ=255)

             gTGGATAAAATAATTTTACCCACCGGATTTTTACCCtagagccatacagc                                     >  1:1398437‑M1/1‑36 (MQ=255)

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

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

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

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

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

                      ataatTTTACCCACCGGATTTTTACCCtagagccatacagcctcacactc                            <  1:1252765‑M1/50‑24 (MQ=255)

                        aatTTTACCCACCGGATTTTTACCCtagagccatacagcctcacactcga                          >  1:1995254‑M1/1‑25 (MQ=255)

                           tttACCCACCGGATTTTTACCCtagagccatacagcctcacactcgat                         >  1:633988‑M1/1‑22 (MQ=255)

                           tttACCCACCGGATTTTTACCCtagagccatacagcctcacactcgat                         >  1:675569‑M1/1‑22 (MQ=255)

                             tACCCACCGGATTTTTACCCtagagccatacagcctcacactcgatgagg                     >  1:904456‑M1/1‑20 (MQ=255)

                             tACCCACCGGATTTTTACCCtagagccatacagcctcacactcgatgagg                     >  1:1260292‑M1/1‑20 (MQ=255)

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

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

                               cccACCGGATTTTTACCCtagagccatacagcctcacactcgatgaggcc                   >  1:481169‑M1/1‑18 (MQ=255)

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

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

                                            tACCCtagagccatacagcctcacactcgatgaggcctgtcagtttctt       >  1:355434‑M1/1‑5 (MQ=255)

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

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

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

                                              ccctagagccatacagcctcacactcgatgaggcctgtcagtttcttaa     <  1:2494417‑M1/49‑47 (MQ=255)

                                               cctagagccatacagcctcacactcgatgaggcctgtcagtttcttaaaa   <  1:1261202‑M1/50‑49 (MQ=255)

                                                ctagagccatacagcctcacactcgatgaggcctgtcagtttcttaaaat  <  1:1231746‑M1/50‑50 (MQ=255)

TGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCAT  >  NC_000913/2466496‑2466593

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 18 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`