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

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

AGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCG  >  NC_000913/1871980‑1872077

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

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

 ggCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGa                                                 <  1:1447464‑M1/50‑2 (MQ=255)

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

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

   ccGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatg                                               >  1:5433378‑M1/1‑47 (MQ=255)

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

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

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

     gCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcg                                             >  1:2616007‑M1/1‑45 (MQ=255)

     gCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcg                                             <  1:4500443‑M1/50‑6 (MQ=255)

     gCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcg                                             >  1:3196944‑M1/1‑45 (MQ=255)

     gCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcg                                             >  1:2764277‑M1/1‑45 (MQ=255)

      cAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcg                                             >  1:3905799‑M1/1‑44 (MQ=255)

       aaTCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgta                                           <  1:477077‑M1/50‑8 (MQ=255)

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

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

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

         tCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtata                                         <  1:903133‑M1/50‑10 (MQ=255)

         tCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtata                                         <  1:3208463‑M1/50‑10 (MQ=255)

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

           gAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaa                                       <  1:4382909‑M1/50‑12 (MQ=255)

           gAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaa                                       <  1:1841961‑M1/50‑12 (MQ=255)

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

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

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

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

             aCTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagc                                     <  1:3278952‑M1/50‑14 (MQ=255)

             aCTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagc                                     <  1:1264079‑M1/50‑14 (MQ=255)

              cTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagct                                    <  1:3668804‑M1/50‑15 (MQ=255)

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

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

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

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

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

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

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

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

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

                        ttttCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagc                          >  1:3941859‑M1/1‑26 (MQ=255)

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

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

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

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

                          ttCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgt                        >  1:1021607‑M1/1‑24 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                                tttCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatat                  <  1:5366653‑M1/50‑33 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                        aCCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataa          <  1:3403310‑M1/50‑41 (MQ=255)

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

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

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

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

                                          cGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataagc        <  1:3541262‑M1/50‑43 (MQ=255)

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

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

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

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

                                             aGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataagcaga     <  1:1367188‑M1/50‑46 (MQ=255)

                                             aGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataagcaga     <  1:806564‑M1/50‑46 (MQ=255)

                                             aGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataagcaga     <  1:4399570‑M1/50‑46 (MQ=255)

                                             aGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataagcaga     <  1:5374033‑M1/50‑46 (MQ=255)

                                             aGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataagcag      >  1:1639989‑M1/1‑5 (MQ=255)

                                             aGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataagcag      >  1:5504943‑M1/1‑5 (MQ=255)

                                              ggAGatgcgtataaagcttccgatgagcgtttatatgtcaataagcagaa    <  1:2046934‑M1/50‑47 (MQ=255)

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

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

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

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

AGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCG  >  NC_000913/1871980‑1872077

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