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

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

Only 100 of 133 total aligned reads displayed.

AGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCG  >  NC_000913/1871980‑1872077

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

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

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

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

    cGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgc                                              >  1:692989‑M1/1‑46 (MQ=255)

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

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

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

      cAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgt                                            <  1:6161427‑M1/50‑7 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

           gAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtata                                         >  1:5684723‑M1/1‑39 (MQ=255)

           gAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtata                                         >  1:3085737‑M1/1‑39 (MQ=255)

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

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

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

              cTTTGTCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagct                                    <  1:5388789‑M1/50‑15 (MQ=255)

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

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

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

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

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

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

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

                 tATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagctt                                   <  1:1004596‑M1/48‑16 (MQ=255)

                 tATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagctt                                   <  1:5974900‑M1/48‑16 (MQ=255)

                 tATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagctt                                   <  1:39185‑M1/48‑16 (MQ=255)

                  aTCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttcc                                 >  1:2486746‑M1/1‑32 (MQ=255)

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

                  aTCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttcc                                 >  1:3685363‑M1/1‑32 (MQ=255)

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

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

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

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

                     tGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgat                              >  1:2075031‑M1/1‑29 (MQ=255)

                      ggTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgat                              >  1:7486109‑M1/1‑28 (MQ=255)

                      ggTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgat                              >  1:2185996‑M1/1‑28 (MQ=255)

                       gTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgag                           <  1:4166237‑M1/50‑24 (MQ=255)

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

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

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

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

                          ttCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcg                         <  1:324781‑M1/49‑26 (MQ=255)

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

                           tctcGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgt                        <  1:187491‑M1/49‑27 (MQ=255)

                           tctcGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgt                        <  1:2386254‑M1/49‑27 (MQ=255)

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

                             tcGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgttta                     >  1:5654482‑M1/1‑21 (MQ=255)

                             tcGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgttta                     >  1:4278323‑M1/1‑21 (MQ=255)

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

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

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

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

                               gTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgttta                     <  1:461108‑M1/48‑30 (MQ=255)

                               gTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgttta                     <  1:1924098‑M1/48‑30 (MQ=255)

                               gTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgttta                     <  1:3930351‑M1/48‑30 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

                                         ccGACGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataag         <  1:5732563‑M1/50‑42 (MQ=255)

                                         ccGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataag         >  1:4281319‑M1/1‑9 (MQ=255)

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

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

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

                                          cGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataag         >  1:4906808‑M1/1‑8 (MQ=255)

                                          cGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataag         >  1:355946‑M1/1‑8 (MQ=255)

                                          cGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataag         >  1:2534206‑M1/1‑8 (MQ=255)

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

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

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

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

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

                                              ggAGatgcgtataaagcttccgatgagcgtttatatgtcaataagcagaa    >  1:1365097‑M1/1‑4 (MQ=255)

                                              ggAGatgcgtataaagcttccgatgagcgtttatatgtcaataagcagaa    >  1:4758037‑M1/1‑4 (MQ=255)

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

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

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

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

                                                 gatgcgtataaagcttccgatgagcgtttatatgtcaataagcagaaca  >  1:5825011‑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` `< 32 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`