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	164 [0]	[0] 166	[dgcJ]	[dgcJ]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1872029	0 (0.000)	157 (1.480)	73/96	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)	157 (1.480)	73/96	0.0	100%	coding (1427/1491 nt)	dgcJ	putative diguanylate cyclase DgcJ

Only 100 of 166 total aligned reads displayed.

GGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAA  >  NC_000913/1873418‑1873516

ggccgcaatcgaactttatctggttttctcgtttcactaaccgaaggaga                                                   <  1:350305‑M2/1‑1 (MQ=255)

  ccgcaatcgaactttatctggttttctcgtttcactaaccgaaggagATg                                                 <  1:6522901‑M2/3‑1 (MQ=255)

   cgcaatcgaactttatctggttttctcgtttcactaaccgaaggagATGc                                                >  1:7780757‑M2/47‑50 (MQ=255)

   cgcaatcgaactttatctggttttctcgtttcactaaccgaaggagATGc                                                >  1:1275353‑M2/47‑50 (MQ=255)

   cgcaatcgaactttatctggttttctcgtttcactaaccgaaggagATGc                                                <  1:4495043‑M2/4‑1 (MQ=255)

    gcaatcgaactttatctggttttctcgtttcactaaccgaaggagATGCg                                               >  1:3378781‑M2/46‑50 (MQ=255)

    gcaatcgaactttatctggttttctcgtttcactaaccgaaggagATGCg                                               <  1:512383‑M2/5‑1 (MQ=255)

    gcaatcgaactttatctggttttctcgtttcactaaccgaaggagATGCg                                               <  1:4926382‑M2/5‑1 (MQ=255)

     caatcgaactttatctggttttctcgtttcactaaccgaaggagATGCg                                               >  1:4558200‑M2/45‑49 (MQ=255)

     caatcgaactttatctggttttctcgtttcactaaccgaaggagATGCGt                                              <  1:4210524‑M2/6‑1 (MQ=255)

      aatcgaactttatctggttttctcgtttcactaaccgaaggagATGCGta                                             >  1:4549249‑M2/44‑50 (MQ=255)

       atcgaactttatctggttttctcgtttcactaaccgaaggagATGCGtat                                            <  1:8100528‑M2/8‑1 (MQ=255)

       atcgaactttatctggttttctcgtttcactaaccgaaggagATGCGtat                                            <  1:4039889‑M2/8‑1 (MQ=255)

        tcgaactttatctggttttctcgtttcactaaccgaaggagATGCGtat                                            >  1:7034701‑M2/42‑49 (MQ=255)

        tcgaactttatctggttttctcgtttcactaaccgaaggagATGCGtat                                            >  1:1895871‑M2/42‑49 (MQ=255)

         cgaactttatctggttttctcgtttcactaaccgaaggagATGCGtata                                           >  1:8100731‑M2/41‑49 (MQ=255)

          gaactttatctggttttctcgtttcactaaccgaaggagATGCGtata                                           >  1:2881300‑M2/40‑48 (MQ=255)

          gaactttatctggttttctcgtttcactaaccgaaggagATGCGTATaaa                                         <  1:4295632‑M2/11‑1 (MQ=255)

          gaactttatctggttttctcgtttcactaaccgaaggagATGCGTATaaa                                         <  1:6547937‑M2/11‑1 (MQ=255)

          gaactttatctggttttctcgtttcactaaccgaaggagATGCGTATaaa                                         <  1:6670647‑M2/11‑1 (MQ=255)

          gaactttatctggttttctcgtttcactaaccgaaggagATGCGTATaa                                          <  1:1500286‑M2/10‑1 (MQ=255)

          gaactttatctggttttctcgtttcactaaccgaaggagATGCGTATaa                                          <  1:1082719‑M2/10‑1 (MQ=255)

           aactttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAg                                        <  1:961344‑M2/12‑1 (MQ=255)

           aactttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAg                                        >  1:2330395‑M2/39‑50 (MQ=255)

             ctttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCt                                      >  1:1926056‑M2/37‑50 (MQ=255)

             ctttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCt                                      >  1:548590‑M2/37‑50 (MQ=255)

              tttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCtt                                     >  1:3482507‑M2/36‑50 (MQ=255)

               ttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTc                                    >  1:1888098‑M2/35‑50 (MQ=255)

               ttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTc                                    >  1:2927550‑M2/35‑50 (MQ=255)

                tatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCtt                                     <  1:5424003‑M2/15‑1 (MQ=255)

                tatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTcc                                   >  1:7624937‑M2/34‑50 (MQ=255)

                 atctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTcc                                   <  1:6654740‑M2/17‑1 (MQ=255)

                 atctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTcc                                   <  1:4834641‑M2/17‑1 (MQ=255)

                 atctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTcc                                   <  1:5380377‑M2/17‑1 (MQ=255)

                   ctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGAt                                >  1:1965008‑M2/31‑50 (MQ=255)

                    tggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGAt                                >  1:6245218‑M2/30‑49 (MQ=255)

                    tggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGAt                                >  1:733534‑M2/30‑49 (MQ=255)

                     ggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGAt                                >  1:1751006‑M2/29‑48 (MQ=255)

                      gttttctcgtttctctaagcgacggagATGCGTATAAAGCTTCCGATGa                              <  1:3800042‑M2/22‑1 (MQ=255)

                      gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGa                                 <  1:8274917‑M2/19‑1 (MQ=255)

                      gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGAt                                <  1:474111‑M2/20‑1 (MQ=255)

                      gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGAt                                <  1:2024797‑M2/20‑1 (MQ=255)

                      gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGAt                                <  1:2307743‑M2/20‑1 (MQ=255)

                      gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGAt                                <  1:3767971‑M2/20‑1 (MQ=255)

                      gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATg                               <  1:5863402‑M2/21‑1 (MQ=255)

                      gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGa                              <  1:1657412‑M2/22‑1 (MQ=255)

                      gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGa                              <  1:845245‑M2/22‑1 (MQ=255)

                      gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAg                             >  1:612998‑M2/28‑50 (MQ=255)

                       ttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGc                            <  1:6994519‑M2/24‑1 (MQ=255)

                         ttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCg                           <  1:1598899‑M2/25‑1 (MQ=255)

                         ttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCg                           <  1:1002532‑M2/25‑1 (MQ=255)

                          tctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGt                          <  1:5388727‑M2/26‑1 (MQ=255)

                           ctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGtt                         <  1:5302644‑M2/27‑1 (MQ=255)

                           ctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGtt                         >  1:7683878‑M2/23‑49 (MQ=255)

                           ctcgcttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGtt                         >  1:2850759‑M2/23‑49 (MQ=255)

                            tcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTta                       >  1:1072595‑M2/22‑50 (MQ=255)

                             cgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTta                       >  1:5879236‑M2/21‑49 (MQ=255)

                             cgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTta                       >  1:1430411‑M2/21‑49 (MQ=255)

                              gtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGttt                        <  1:20075‑M2/28‑1 (MQ=255)

                              gtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTtata                     >  1:1869792‑M2/20‑50 (MQ=255)

                              gtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTtat                      <  1:1190575‑M2/30‑1 (MQ=255)

                               tttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTatat                    <  1:3639847‑M2/32‑1 (MQ=255)

                                ttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATg                   >  1:796145‑M2/18‑50 (MQ=255)

                                 tcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATg                   <  1:7167688‑M2/33‑1 (MQ=255)

                                  cactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGt                  >  1:3808455‑M2/16‑49 (MQ=255)

                                  cactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGt                  <  1:6160770‑M2/34‑1 (MQ=255)

                                   actaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTc                 >  1:4210314‑M2/15‑49 (MQ=255)

                                    ctaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCaa               >  1:6330690‑M2/14‑50 (MQ=255)

                                    ctaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCaa               >  1:2783681‑M2/14‑50 (MQ=255)

                                    ctaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCaa               <  1:6111345‑M2/37‑1 (MQ=255)

                                    ctaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCaa               <  1:3485624‑M2/37‑1 (MQ=255)

                                      aaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAAt              >  1:3079147‑M2/12‑49 (MQ=255)

                                       accgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAAt              <  1:7752011‑M2/38‑1 (MQ=255)

                                       accgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATaa            <  1:250372‑M2/40‑1 (MQ=255)

                                       accgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATa             <  1:4010864‑M2/39‑1 (MQ=255)

                                       accgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATa             <  1:1496353‑M2/39‑1 (MQ=255)

                                       accgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATa             <  1:247538‑M2/39‑1 (MQ=255)

                                       accgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATa             <  1:5665865‑M2/39‑1 (MQ=255)

                                        ccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAg           >  1:3096885‑M2/10‑50 (MQ=255)

                                        ccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAg           >  1:88936‑M2/10‑50 (MQ=255)

                                        ccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAg           >  1:1133471‑M2/10‑50 (MQ=255)

                                        ccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAg           <  1:5694906‑M2/41‑1 (MQ=255)

                                        ccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAg           <  1:6738000‑M2/41‑1 (MQ=255)

                                         cgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAg           >  1:2596849‑M2/9‑49 (MQ=255)

                                         cgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAg           >  1:4122354‑M2/9‑49 (MQ=255)

                                          gaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAg           >  1:5273312‑M2/8‑48 (MQ=255)

                                          gaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCa         <  1:3926580‑M2/43‑1 (MQ=255)

                                           aaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAg        <  1:2657213‑M2/44‑1 (MQ=255)

                                           aaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAg        <  1:4507483‑M2/44‑1 (MQ=255)

                                            aggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGa       <  1:7201432‑M2/45‑1 (MQ=255)

                                               agATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAAc     >  1:5776150‑M2/3‑49 (MQ=255)

                                               agATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACa    <  1:6772563‑M2/48‑1 (MQ=255)

                                               agATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACa    <  1:6318412‑M2/48‑1 (MQ=255)

                                                gATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACaa   <  1:5980030‑M2/49‑1 (MQ=255)

                                                gATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACaa   <  1:6642614‑M2/49‑1 (MQ=255)

                                                gATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACa    >  1:3568087‑M2/2‑49 (MQ=255)

                                                gATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACa    >  1:2460205‑M2/2‑49 (MQ=255)

                                                gATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACa    >  1:7176827‑M2/2‑49 (MQ=255)

                                                 atgCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACaaa  <  1:2718421/50‑1 (MQ=255)

                                                 atgCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACaaa  <  1:1256957/50‑1 (MQ=255)

GGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAA  >  NC_000913/1873418‑1873516

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: `‑`