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

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

Only 100 of 129 total aligned reads displayed.

GGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAA  >  NC_000913/1873418‑1873516

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

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

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

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

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

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

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

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

   cgcaatcgaactttatctggttttctcgtttcactaacagaaggagATGc                                                <  1:3719289‑M2/4‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

      aatcgaactttatctggttttctcgtttcactaaccgaaggagATGCGta                                             <  1:1007621‑M2/7‑1 (MQ=255)

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

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

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

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

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

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

        tcgaactttatctggttttctcgtttcactaaccgaaggagATGCGtata                                           <  1:5621716‑M2/9‑1 (MQ=255)

        tcgaactttatctggttttctcgtttcactaaccgaaggagATGCGtata                                           <  1:883362‑M2/9‑1 (MQ=255)

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

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

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

          gaactttatctggttttctcgtttcactaaccgaaggagATGTGTATaaa                                         <  1:694716‑M2/11‑1 (MQ=255)

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

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

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

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

            actttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGc                                       >  1:4522944‑M2/38‑50 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                 atctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTcc                                   >  1:3861889‑M2/33‑49 (MQ=255)

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

                 atctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTcc                                   >  1:483219‑M2/33‑49 (MQ=255)

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

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

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

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

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

                       ttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGc                            >  1:5346211‑M2/27‑50 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                 atgCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACaa   >  1:2626885/1‑49 (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: `‑`