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

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

Only 100 of 101 total aligned reads displayed.

CTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAA  >  NC_000913/1873420‑1873516

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                             tttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGAGTTTatat                    >  1:3999532‑M2/19‑50 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                      ccgaaggagATGCGTATAAAGCTTCCGTTGAGCGTTTATATGTCAATAAg           >  1:6425499‑M2/10‑50 (MQ=255)

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

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

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

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

                                        gaaggagATGTGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAg           >  1:5629426‑M2/8‑48 (MQ=255)

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

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

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

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

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

                                           ggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGTAGaa      >  1:2490784‑M2/5‑50 (MQ=255)

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

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

                                            gagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGaa      >  1:3215673‑M2/4‑49 (MQ=255)

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

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

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

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

CTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAA  >  NC_000913/1873420‑1873516

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