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] 99	[dgcJ]	[dgcJ]

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

AGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCG  >  NC_000913/1871980‑1872077

aGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAG                                                  <  1:991594‑M1/50‑1 (MQ=255)

aGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAG                                                  <  1:6609457‑M1/50‑1 (MQ=255)

  gCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGat                                                <  1:5178224‑M1/50‑3 (MQ=255)

   ccGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatg                                               <  1:5007406‑M1/50‑4 (MQ=255)

    cGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgc                                              <  1:6278950‑M1/50‑5 (MQ=255)

       aaTCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgta                                           >  1:249205‑M1/1‑43 (MQ=255)

       aaTCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgta                                           >  1:7182134‑M1/1‑43 (MQ=255)

       aaTCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgta                                           >  1:3126456‑M1/1‑43 (MQ=255)

        aTCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtat                                          <  1:3550643‑M1/50‑9 (MQ=255)

        aTCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtat                                          <  1:4287110‑M1/50‑9 (MQ=255)

        aTCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtat                                          <  1:4424348‑M1/50‑9 (MQ=255)

        aTCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgta                                           >  1:1273447‑M1/1‑42 (MQ=255)

        aTCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgta                                           >  1:2397274‑M1/1‑42 (MQ=255)

          cGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtata                                         >  1:942116‑M1/1‑40 (MQ=255)

           gAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataa                                        <  1:1049340‑M1/49‑11 (MQ=255)

           gAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataa                                        <  1:4929164‑M1/49‑11 (MQ=255)

           gAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataa                                        <  1:6626176‑M1/49‑11 (MQ=255)

            aaCTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaag                                      >  1:3563329‑M1/1‑38 (MQ=255)

            aaCTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaag                                      >  1:5093115‑M1/1‑38 (MQ=255)

            aaCTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaag                                      >  1:2593104‑M1/1‑38 (MQ=255)

            aaCTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaag                                      <  1:1541907‑M1/50‑13 (MQ=255)

              cTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagct                                    >  1:1169462‑M1/1‑36 (MQ=255)

               tttATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagctt                                   >  1:410007‑M1/1‑35 (MQ=255)

               tttATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagctt                                   >  1:6708390‑M1/1‑35 (MQ=255)

                ttATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttc                                  >  1:6947647‑M1/1‑34 (MQ=255)

                 tATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttc                                  <  1:5486490‑M1/49‑17 (MQ=255)

                 tATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttc                                  <  1:480687‑M1/49‑17 (MQ=255)

                 tATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttc                                  <  1:299776‑M1/49‑17 (MQ=255)

                 tATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttc                                  <  1:1365779‑M1/49‑17 (MQ=255)

                 tATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttc                                  <  1:1984157‑M1/49‑17 (MQ=255)

                  aTCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttcc                                 <  1:558249‑M1/49‑18 (MQ=255)

                  aTCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttcc                                 <  1:4144392‑M1/49‑18 (MQ=255)

                   tCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccg                                <  1:4718591‑M1/49‑19 (MQ=255)

                    cTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgat                              >  1:49057‑M1/1‑30 (MQ=255)

                       gTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatga                            <  1:894779‑M1/49‑23 (MQ=255)

                       gTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatga                            <  1:2195767‑M1/49‑23 (MQ=255)

                       gTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatg                             <  1:2969730‑M1/48‑22 (MQ=255)

                       gTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatg                             <  1:6000854‑M1/48‑22 (MQ=255)

                       gTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgat                              <  1:3817565‑M1/47‑21 (MQ=255)

                       gTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccga                               <  1:7155123‑M1/46‑20 (MQ=255)

                        ttttCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagc                          <  1:3396829‑M1/50‑25 (MQ=255)

                         tttCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcg                         >  1:3009939‑M1/1‑25 (MQ=255)

                         tttCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcg                         >  1:684971‑M1/1‑25 (MQ=255)

                           tctcGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtt                       >  1:2735714‑M1/1‑23 (MQ=255)

                            ctcGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtt                       <  1:6303600‑M1/49‑28 (MQ=255)

                            ctcGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtt                       >  1:2621470‑M1/1‑22 (MQ=255)

                            ctcGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtt                       <  1:4703445‑M1/49‑28 (MQ=255)

                            ctcGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtt                       >  1:7268801‑M1/1‑22 (MQ=255)

                            ctcGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtt                       <  1:3875276‑M1/49‑28 (MQ=255)

                              cGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgttta                     >  1:1836341‑M1/1‑20 (MQ=255)

                               gTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttata                   <  1:3222071‑M1/50‑32 (MQ=255)

                               gTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttata                   >  1:6933371‑M1/1‑19 (MQ=255)

                               gTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttata                   <  1:3215658‑M1/50‑32 (MQ=255)

                               gTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttata                   <  1:1878401‑M1/50‑32 (MQ=255)

                               gTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttata                   <  1:5791099‑M1/50‑32 (MQ=255)

                               gTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttat                    <  1:1817603‑M1/49‑31 (MQ=255)

                               gTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgttt                      <  1:5577952‑M1/47‑29 (MQ=255)

                               gTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgttt                      <  1:6342631‑M1/47‑29 (MQ=255)

                               gTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgttt                      <  1:3535900‑M1/47‑29 (MQ=255)

                                tttCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatat                  >  1:7607786‑M1/1‑18 (MQ=255)

                                 ttCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatg                 >  1:5104641‑M1/1‑17 (MQ=255)

                                 ttCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatg                 >  1:2145948‑M1/1‑17 (MQ=255)

                                 ttCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatg                 >  1:1297239‑M1/1‑17 (MQ=255)

                                  tCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgt                >  1:1563299‑M1/1‑16 (MQ=255)

                                  tCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgt                >  1:7579047‑M1/1‑16 (MQ=255)

                                  tCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatg                 <  1:5238223‑M1/49‑34 (MQ=255)

                                  tCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatg                 <  1:4570225‑M1/49‑34 (MQ=255)

                                   cACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgt                >  1:5427029‑M1/1‑15 (MQ=255)

                                   cACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgt                <  1:6405675‑M1/49‑35 (MQ=255)

                                    aCTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtca              <  1:2013652‑M1/50‑37 (MQ=255)

                                    aCTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtca              <  1:6875084‑M1/50‑37 (MQ=255)

                                    aCTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtc               >  1:3244650‑M1/1‑14 (MQ=255)

                                    aCTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtc               >  1:1617653‑M1/1‑14 (MQ=255)

                                    aCTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtc               <  1:5558821‑M1/49‑36 (MQ=255)

                                    aCTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtc               <  1:6752318‑M1/49‑36 (MQ=255)

                                     cTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaa             <  1:2338742‑M1/50‑38 (MQ=255)

                                     cTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaa             >  1:5083397‑M1/1‑13 (MQ=255)

                                     cTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaa             >  1:4192075‑M1/1‑13 (MQ=255)

                                      tAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaa             >  1:7432479‑M1/1‑12 (MQ=255)

                                       aaCCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaat            >  1:780349‑M1/1‑11 (MQ=255)

                                        aCCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaata           <  1:2347462‑M1/49‑40 (MQ=255)

                                        aCCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaat            <  1:2792273‑M1/48‑39 (MQ=255)

                                        aCCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaat            >  1:6943371‑M1/1‑10 (MQ=255)

                                         ccGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataag         <  1:6810925‑M1/50‑42 (MQ=255)

                                         ccGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataag         <  1:6478098‑M1/50‑42 (MQ=255)

                                         ccGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataag         <  1:870354‑M1/50‑42 (MQ=255)

                                           gAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataagca       <  1:2731390‑M1/50‑44 (MQ=255)

                                           gAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataag         >  1:2550698‑M1/1‑7 (MQ=255)

                                           gAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataag         >  1:6052463‑M1/1‑7 (MQ=255)

                                           gAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataag         >  1:5646970‑M1/1‑7 (MQ=255)

                                            aaGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataagcag      <  1:4765534‑M1/50‑45 (MQ=255)

                                            aaGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataag         >  1:3708644‑M1/1‑6 (MQ=255)

                                               gagatgcgtataaagcttccgatgagcgtttatatgtcaataagcagaac   <  1:1881358‑M1/50‑48 (MQ=255)

                                               gagatgcgtataaagcttccgatgagcgtttatatgtcaataagcagaac   <  1:5323835‑M1/50‑48 (MQ=255)

                                               gagatgcgtataaagcttccgatgagcgtttatatgtcaataagcagaa    >  1:3369039‑M1/1‑3 (MQ=255)

                                                agatgcgtataaagcttccgatgagcgtttatatgtcaataagcagaaca  <  1:2189287‑M1/50‑49 (MQ=255)

                                                agatgcgtataaagcttccgatgagcgtttatatgtcaataagcagaaca  <  1:5085086‑M1/50‑49 (MQ=255)

                                                agatgcgtataaagcttccgatgagcgtttatatgtcaataagcagaaca  <  1:7115214‑M1/50‑49 (MQ=255)

                                                agatgcgtataaagcttccgatgagcgtttatatgtcaataagcagaac   >  1:1496301‑M1/1‑2 (MQ=255)

                                                 gatgcgtataaagcttccgatgagcgtttatatgtcaataagcagaaca  >  1:5991315‑M1/1‑1 (MQ=255)

AGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCG  >  NC_000913/1871980‑1872077

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: `‑`