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

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

Only 100 of 188 total aligned reads displayed.

GCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAA  >  NC_000913/1873419‑1873516

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

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

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

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

    caatcgaactttatctggttttctcgtttcactaaccgaaggagATGCg                                               >  2:1390309‑M2/45‑49 (MQ=255)

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

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

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

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

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

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

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

      atcgaactttatctggttttctcgtttcactaaccgaaggagATGCGta                                             >  2:4522614‑M2/43‑49 (MQ=255)

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

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

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

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

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

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

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

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

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

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

           actttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGc                                       >  2:1866682‑M2/38‑50 (MQ=255)

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

            ctttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCt                                      >  2:4366096‑M2/37‑50 (MQ=255)

            ctttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCt                                      >  2:1137404‑M2/37‑50 (MQ=255)

               tatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTcc                                   >  2:3587608‑M2/34‑50 (MQ=255)

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

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

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

                atctggttttctcgtttcactaaccgaaggagATACGTATAAAGCTTcc                                   >  2:1411041‑M2/33‑49 (MQ=255)

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

                  ctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGAt                                >  2:380575‑M2/31‑50 (MQ=255)

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

                  ctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGAt                                >  2:839122‑M2/31‑50 (MQ=255)

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

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

                    ggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGAt                                >  2:4639462‑M2/29‑48 (MQ=255)

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

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

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

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

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

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

                     gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAg                             >  2:2265931‑M2/28‑50 (MQ=255)

                     gttttctcgtttcactaacagaaggagATGCGTATAAAGCTTCCGATg                               <  1:1540334‑M2/21‑1 (MQ=255)

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

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

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

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

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

                       tttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCg                           >  2:2238350‑M2/26‑50 (MQ=255)

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

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

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

                        ttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATCAGCGt                          >  2:364064‑M2/25‑50 (MQ=255)

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

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

                            cgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTta                       >  2:4718815‑M2/21‑49 (MQ=255)

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

                             gtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTtata                     >  2:507128‑M2/20‑50 (MQ=255)

                             gtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTtata                     >  2:924543‑M2/20‑50 (MQ=255)

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

                              tttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTatat                    >  2:1050611‑M2/19‑50 (MQ=255)

                              tttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTatat                    >  2:4143592‑M2/19‑50 (MQ=255)

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

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

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

                                  actaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTc                 >  2:1158900‑M2/15‑49 (MQ=255)

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

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

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

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

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

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

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

                                     aaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAAt              >  2:1266931‑M2/12‑49 (MQ=255)

                                     aaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAAt              >  2:714304‑M2/12‑49 (MQ=255)

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

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

                                      accgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAAt              >  2:5153085‑M2/11‑48 (MQ=255)

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

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

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

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

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

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

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

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

                                           aggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAg        >  2:325078‑M2/6‑49 (MQ=255)

                                           aggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAg        >  2:3703460‑M2/6‑49 (MQ=255)

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

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

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

                                             gagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGaa      >  2:3878121‑M2/4‑49 (MQ=255)

                                              agATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAAc     >  2:4068793‑M2/3‑49 (MQ=255)

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

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

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

GCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAA  >  NC_000913/1873419‑1873516

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