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

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

AGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCG  >  NC_000913/1871980‑1872077

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

  gCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGat                                                >  1:1104752‑M1/1‑48 (MQ=255)

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

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

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

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

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

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

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

         tCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtat                                          >  1:1714677‑M1/1‑41 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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