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	91 [0]	[0] 93	[dgcJ]	[dgcJ]

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

GCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAA  >  NC_000913/1873419‑1873516

gccgcaatcgaactttatctggttttctcgtttcactaaccgaaggagAt                                                  >  1:689415‑M2/49‑50 (MQ=255)

gccgcaatcgaactttatctggttttctcgtttcactaaccgaaggagAt                                                  <  1:3275130‑M2/2‑1 (MQ=255)

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

 ccgcaatcgaactttatctggttttctcgtttcactaaccgaaggagATg                                                 >  1:6043998‑M2/48‑50 (MQ=255)

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

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

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

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

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

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

      atcgaactttatctggttttctcgtttcactaaccgaaggagATGCGta                                             >  1:584928‑M2/43‑49 (MQ=255)

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

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

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

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

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

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

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

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

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

           actttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGc                                       <  1:2462285‑M2/13‑1 (MQ=255)

           actttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGc                                       <  1:4762334‑M2/13‑1 (MQ=255)

           actttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGc                                       <  1:5813098‑M2/13‑1 (MQ=255)

            ctttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCt                                      <  1:3687179‑M2/14‑1 (MQ=255)

             tttatctggttttctcgtttcactaaccgcaggagATGCGTATAAAGCtt                                     >  1:3720672‑M2/36‑50 (MQ=255)

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

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

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

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

                 tctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGa                                 >  1:269014‑M2/32‑50 (MQ=255)

                 tctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGa                                 >  1:5891986‑M2/32‑50 (MQ=255)

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

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

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

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

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

                     gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAg                             <  1:3045022‑M2/23‑1 (MQ=255)

                     gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAg                             <  1:5179618‑M2/23‑1 (MQ=255)

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

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

                       tttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCg                           >  1:6331455‑M2/26‑50 (MQ=255)

                       tttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCg                           >  1:5333650‑M2/26‑50 (MQ=255)

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

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

                        ttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGt                          >  1:1477168‑M2/25‑50 (MQ=255)

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

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

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

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

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

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

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

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

                             gtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTtata                     <  1:5523941‑M2/31‑1 (MQ=255)

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

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

                             gtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTta                       <  1:654479‑M2/29‑1 (MQ=255)

                              tttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTatat                    >  1:2585313‑M2/19‑50 (MQ=255)

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

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

                                tcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGt                  >  1:2330345‑M2/17‑50 (MQ=255)

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

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

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

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

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

                                  actaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTc                 <  1:451589‑M2/35‑1 (MQ=255)

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

                                  actaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCa                <  1:3625999‑M2/36‑1 (MQ=255)

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

                                    taaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCaa               >  1:1137765‑M2/13‑49 (MQ=255)

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

                                      accgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAAt              >  1:4198497‑M2/11‑48 (MQ=255)

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

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

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

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

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

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

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

                                        cgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGc          <  1:590051‑M2/42‑1 (MQ=255)

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

                                           aggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAg        >  1:3668994‑M2/6‑49 (MQ=255)

                                           aggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAg        >  1:1702372‑M2/6‑49 (MQ=255)

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

                                            ggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGaa      >  1:40941‑M2/5‑50 (MQ=255)

                                            ggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGaa      <  1:1964321‑M2/46‑1 (MQ=255)

                                            ggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGaa      <  1:4586360‑M2/46‑1 (MQ=255)

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

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

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

                                                atgCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACaa   >  1:3047374/1‑49 (MQ=255)

                                                atgCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACaa   >  1:4243145/1‑49 (MQ=255)

GCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAA  >  NC_000913/1873419‑1873516

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