BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,872,030	Δ1,437 bp		[dgcJ]	[dgcJ]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1872030	1873466	1437	87 [0]	[0] 88	[dgcJ]	[dgcJ]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1872029	0 (0.000)	81 (1.150)	42/66	0.0	100%	intergenic (+169/‑12)	cdgI/dgcJ	putative c‑di‑GMP binding protein CdgI/putative diguanylate cyclase DgcJ
*	?	NC_000913	1873467 =	0 (0.000)	81 (1.150)	42/66	0.0	100%	coding (1427/1491 nt)	dgcJ	putative diguanylate cyclase DgcJ

TTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTGCCATTTATCATGAAATTGCACCATAGAATGCT > NC_000913/1871995‑1872063

tttATCTGGTTTTCTCGTTTCACTAACCGAAGGAG                                    <  1:1255792‑M1/35‑1 (MQ=255)

 ttATCTGGTTTTCTCGTTTCACTAACCGAAGGAGa                                   >  1:6186550‑M1/1‑34 (MQ=255)

  tATCTGGTTTTCTCGTTTCACTAACCGAAGGAGat                                  <  1:271423‑M1/35‑3 (MQ=255)

  tATCTGGTTTTCTCGTTTCACTAACCGAAGGAGat                                  >  1:1322166‑M1/1‑33 (MQ=255)

  tATCTGGTTTTCTCGTTTCACTAACCGAAGGAGat                                  >  1:839715‑M1/1‑33 (MQ=255)

  tATCTGGTTTTCTCGTTTCACTAACCGAAGGAGat                                  <  1:858987‑M1/35‑3 (MQ=255)

   aTCTGGTTTTCTCGTTTCACTAACCGAAGGAGatg                                 >  1:3928767‑M1/1‑32 (MQ=255)

    tCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgc                                >  1:687632‑M1/1‑31 (MQ=255)

    tCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgc                                <  1:3036626‑M1/35‑5 (MQ=255)

    tCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgc                                >  1:1557289‑M1/1‑31 (MQ=255)

     cTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcg                               <  1:5956589‑M1/35‑6 (MQ=255)

      tGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgt                              <  1:4083701‑M1/35‑7 (MQ=255)

       ggTTTTCTCGTTTCACTAACCGAAGGAGatgcgta                             <  1:7027036‑M1/35‑8 (MQ=255)

       ggTTTTCTCGTTTCACTAACCGAAGGAGatgcgta                             <  1:5283044‑M1/35‑8 (MQ=255)

        gTTTTCTCGTTTCACTAACCGAAGGAGatgcgtat                            <  1:4361525‑M1/35‑9 (MQ=255)

        gTTTTCTCGTTTCACTAACCGAAGGAGatgcgtat                            >  1:1443727‑M1/1‑27 (MQ=255)

         ttttCTCGTTTCACTAACCGAAGGAGatgcgtata                           <  1:7287664‑M1/35‑10 (MQ=255)

         ttttCTCGTTTCACTAACCGAAGGAGatgcgtata                           >  1:519659‑M1/1‑26 (MQ=255)

          tttCTCGTTTCACTAACCGAAGGAGatgcgtataa                          <  1:7013282‑M1/35‑11 (MQ=255)

          tttCTCGTTTCACTAACCGAAGGAGatgcgtataa                          <  1:6804709‑M1/35‑11 (MQ=255)

           ttCTCGTTTCACTAACCGAAGGAGatgcgtataaa                         >  1:6236986‑M1/1‑24 (MQ=255)

           ttCTCGTTTCACTAACCGAAGGAGatgcgtataaa                         >  1:1732977‑M1/1‑24 (MQ=255)

           ttCTCGTTTCACTAACCGAAGGAGatgcgtataaa                         >  1:890861‑M1/1‑24 (MQ=255)

           ttCTCGTTTCACTAACCGAAGGAGatgcgtataaa                         >  1:7774179‑M1/1‑24 (MQ=255)

            tctcGTTTCACTAACCGAAGGAGatgcgtataaag                        >  1:5584671‑M1/1‑23 (MQ=255)

            tctcGTTTCACTAACCGAAGGAGatgcgtataaag                        >  1:3256628‑M1/1‑23 (MQ=255)

            tctcGTTTCACTAACCGAAGGAGatgcgtataaag                        <  1:821988‑M1/35‑13 (MQ=255)

              tcGTTTCACTAACCGAAGGAGatgcgtataaagct                      >  1:7729914‑M1/1‑21 (MQ=255)

               cGTTTCACTAACCGAAGGAGatgcgtataaagctt                     <  1:177084‑M1/35‑16 (MQ=255)

                gTTTCACTAACCGAAGGAGatgcgtataaagcttc                    >  1:3521433‑M1/1‑19 (MQ=255)

                gTTTCACTAACCGAAGGAGatgcgtataaagcttc                    <  1:2249449‑M1/35‑17 (MQ=255)

                gTTTCACTAACCGAAGGAGatgcgtataaagcttc                    <  1:2298441‑M1/35‑17 (MQ=255)

                gTTTCACTAACCGAAGGAGatgcgtataaagcttc                    >  1:2105258‑M1/1‑19 (MQ=255)

                gTTTCACTAACCGAAGGAGatgcgtataaagcttc                    >  1:3455636‑M1/1‑19 (MQ=255)

                gTTTCACTAACCGAAGGAGatgcgtataaagcttc                    <  1:2849941‑M1/35‑17 (MQ=255)

                 tttCACTAACCGAAGGAGatgcgtataaagcttcc                   >  1:5556703‑M1/1‑18 (MQ=255)

                 tttCACTAACCGAAGGAGatgcgtataaagcttcc                   >  1:6276274‑M1/1‑18 (MQ=255)

                 tttCACTAACCGAAGGAGatgcgtataaagcttcc                   <  1:2144296‑M1/35‑18 (MQ=255)

                 tttCACTAACCGAAGGAGatgcgtataaagcttcc                   >  1:1696094‑M1/1‑18 (MQ=255)

                  ttCACTAACCGAAGGAGatgcgtataaagcttccg                  >  1:5603015‑M1/1‑17 (MQ=255)

                  ttCACTAACCGAAGGAGatgcgtataaagcttccg                  >  1:4296437‑M1/1‑17 (MQ=255)

                  ttCACTAACCGAAGGAGatgcgtataaagcttccg                  >  1:1329596‑M1/1‑17 (MQ=255)

                  ttCACTAACCGAAGGAGatgcgtataaagcttccg                  >  1:1950818‑M1/1‑17 (MQ=255)

                   tCACTAACCGAAGGAGatgcgtataaagcttccga                 >  1:6986247‑M1/1‑16 (MQ=255)

                   tCACTAACCGAAGGAGatgcgtataaagcttccga                 >  1:6256128‑M1/1‑16 (MQ=255)

                   tCACTAACCGAAGGAGatgcgtataaagcttccga                 <  1:7435904‑M1/35‑20 (MQ=255)

                    cACTAACCGAAGGAGatgcgtataaagcttccgat                <  1:2789181‑M1/35‑21 (MQ=255)

                    cACTAACCGAAGGAGatgcgtataaagcttccgat                >  1:2380463‑M1/1‑15 (MQ=255)

                    cACTAACCGAAGGAGatgcgtataaagcttccgat                <  1:2311753‑M1/35‑21 (MQ=255)

                    cACTAACCGAAGGAGatgcgtataaagcttccgat                <  1:2680496‑M1/35‑21 (MQ=255)

                     aCTAACCGAAGGAGatgcgtataaagcttccgatg               >  1:377558‑M1/1‑14 (MQ=255)

                       tAACCGAAGGAGatgcgtataaagcttccgatgag             <  1:5601011‑M1/35‑24 (MQ=255)

                         aCCGAAGGAGatgcgtataaagcttccgatgagcg           >  1:3531127‑M1/1‑10 (MQ=255)

                         aCCGAAGGAGatgcgtataaagcttccgatgagcg           >  1:4111820‑M1/1‑10 (MQ=255)

                         aCCGAAGGAGatgcgtataaagcttccgatgagcg           >  1:2155840‑M1/1‑10 (MQ=255)

                         aCCGAAGGAGatgcgtataaagcttccgatgagcg           >  1:2930316‑M1/1‑10 (MQ=255)

                          ccGAAGGAGatgcgtataaagcttccgatgagcgt          <  1:446016‑M1/35‑27 (MQ=255)

                          ccGAAGGAGatgcgtataaagcttccgatgagcgt          <  1:5796969‑M1/35‑27 (MQ=255)

                           cGAAGGAGatgcgtataaagcttccgatgagcgtt         <  1:7676638‑M1/35‑28 (MQ=255)

                           cGAAGGAGatgcgtataaagcttccgatgagcgtt         <  1:4618676‑M1/35‑28 (MQ=255)

                           cGAAGGAGatgcgtataaagcttccgatgagcgtt         <  1:4607234‑M1/35‑28 (MQ=255)

                           cGAAGGAGatgcgtataaagcttccgatgagcgtt         <  1:775640‑M1/35‑28 (MQ=255)

                           cGAAGGAGatgcgtataaagcttccgatgagcgtt         >  1:3078051‑M1/1‑8 (MQ=255)

                            gAAGGAGatgcgtataaagcttccgatgagcgttt        <  1:5637506‑M1/35‑29 (MQ=255)

                            gAAGGAGatgcgtataaagcttccgatgagcgttt        <  1:6921942‑M1/35‑29 (MQ=255)

                             aaGGAGatgcgtataaagcttccgatgagcgttta       >  1:7600142‑M1/1‑6 (MQ=255)

                             aaGGAGatgcgtataaagcttccgatgagcgttta       <  1:6973692‑M1/35‑30 (MQ=255)

                             aaGGAGatgcgtataaagcttccgatgagcgttta       <  1:942373‑M1/35‑30 (MQ=255)

                             aaGGAGatgcgtataaagcttccgatgagcgttta       >  1:6504559‑M1/1‑6 (MQ=255)

                             aaGGAGatgcgtataaagcttccgatgagcgttta       >  1:4272765‑M1/1‑6 (MQ=255)

                             aaGGAGatgcgtataaagcttccgatgagcgttta       <  1:4516283‑M1/35‑30 (MQ=255)

                             aaGGAGatgcgtataaagcttccgatgagcgttta       <  1:4884888‑M1/35‑30 (MQ=255)

                              aGGAGatgcgtataaagcttccgatgagcgtttat      >  1:2264011‑M1/1‑5 (MQ=255)

                              aGGAGatgcgtataaagcttccgatgagcgtttat      >  1:604884‑M1/1‑5 (MQ=255)

                              aGGAGatgcgtataaagcttccgatgagcgtttat      <  1:5997526‑M1/35‑31 (MQ=255)

                               ggAGatgcgtataaagcttccgatgagcgtttata     <  1:1049015‑M1/35‑32 (MQ=255)

                               ggAGatgcgtataaagcttccgatgagcgtttata     <  1:5377362‑M1/35‑32 (MQ=255)

                               ggAGatgcgtataaagcttccgatgagcgtttata     <  1:7133191‑M1/35‑32 (MQ=255)

                                gagatgcgtataaagcttccgatgagcgtttatat    >  1:4994127‑M1/1‑3 (MQ=255)

                                gagatgcgtataaagcttccgatgagcgtttatat    >  1:2947531‑M1/1‑3 (MQ=255)

                                 agatgcgtataaagcttccgatgagcgtttatatg   >  1:6478097‑M1/1‑2 (MQ=255)

                                 agatgcgtataaagcttccgatgagcgtttatatg   >  1:1082509‑M1/1‑2 (MQ=255)

                                  gatgcgtataaagcttccgatgagcgtttatatgt  >  1:5122929‑M1/1‑1 (MQ=255)

                                  gatgcgtataaagcttccgatgagcgtttatatgt  <  1:4365636‑M1/35‑35 (MQ=255)

                                  gatgcgtataaagcttccgatgagcgtttatatgt  <  1:5571300‑M1/35‑35 (MQ=255)

                                  gatgcgtataaagcttccgatgagcgtttatatgt  <  1:6724747‑M1/35‑35 (MQ=255)

                                  gatgcgtataaaccttccgatgagcgtttatatgc  >  1:3779798‑M1/1‑1 (MQ=255)

TTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTGCCATTTATCATGAAATTGCACCATAGAATGCT > NC_000913/1871995‑1872063

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 29 ≤` `ATCG/ATCG` `< 35 ≤` `ATCG/ATCG` `< 38 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`