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

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

GCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAA  >  NC_000913/1873419‑1873516

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

               tatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTc                                    <  1:787832‑M2/16‑1 (MQ=255)

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

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

                 tctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCg                                  <  1:3141643‑M2/18‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                         tctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGtt                         >  1:940977‑M2/24‑50 (MQ=255)

                         tctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGtt                         >  1:3245663‑M2/24‑50 (MQ=255)

                         tctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGtt                         >  1:1606733‑M2/24‑50 (MQ=255)

                         tctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGtt                         >  1:873358‑M2/24‑50 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                          aaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAg           >  1:3246757‑M2/7‑47 (MQ=255)

                                          aaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAg           >  1:1953574‑M2/7‑47 (MQ=255)

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

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

                                             gagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAAc     <  1:3399624‑M2/47‑1 (MQ=255)

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

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

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

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

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

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

GCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAA  >  NC_000913/1873419‑1873516

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