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

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

Only 100 of 223 total aligned reads displayed.

AGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCG  >  NC_000913/1871980‑1872077

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

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

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

  gCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGat                                                <  2:1377690‑M1/50‑3 (MQ=255)

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

  gCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGat                                                <  2:4417057‑M1/50‑3 (MQ=255)

   ccGCAATCGAGCTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatg                                               >  1:3542816‑M1/1‑47 (MQ=255)

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

   ccGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatg                                               <  2:3776718‑M1/50‑4 (MQ=255)

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

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

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

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

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

     gCAATCGAACTTTATCTGGTTTTCTCGTTTCACGAACCGAAGGAGatgcg                                             <  1:1691597‑M1/50‑6 (MQ=255)

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

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

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

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

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

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

         tCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGttgcgtat                                          >  1:104780‑M1/1‑41 (MQ=255)

           gAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaa                                       <  2:1197635‑M1/50‑12 (MQ=255)

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

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

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

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

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

             aCTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagc                                     >  2:1259564‑M1/1‑37 (MQ=255)

              cTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagct                                    <  2:3856543‑M1/50‑15 (MQ=255)

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

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

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

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

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

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

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

                  aTCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttcc                                 <  2:2617901‑M1/49‑18 (MQ=255)

                   tCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccga                               >  1:5332313‑M1/1‑31 (MQ=255)

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

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

                       gTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgag                           >  1:1476909‑M1/1‑27 (MQ=255)

                       gTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgag                           <  2:5314155‑M1/50‑24 (MQ=255)

                       gTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgag                           >  2:2958647‑M1/1‑27 (MQ=255)

                       gTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgag                           >  1:5509470‑M1/1‑27 (MQ=255)

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

                       gTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatga                            <  2:3512127‑M1/49‑23 (MQ=255)

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

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

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

                       gTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgat                              <  2:367619‑M1/47‑21 (MQ=255)

                       gTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgat                              <  2:5578373‑M1/47‑21 (MQ=255)

                       gTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgat                              <  2:2832681‑M1/47‑21 (MQ=255)

                       gTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccga                               <  2:1018‑M1/46‑20 (MQ=255)

                       gTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccga                               <  2:1173382‑M1/46‑20 (MQ=255)

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

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

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

                        ttttCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagc                          <  2:1452837‑M1/50‑25 (MQ=255)

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

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

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

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

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

                          ttCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcg                         <  2:1286706‑M1/49‑26 (MQ=255)

                            ctcGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtt                       <  2:5264077‑M1/49‑28 (MQ=255)

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

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

                               gTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttata                   <  2:704152‑M1/50‑32 (MQ=255)

                               gTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttat                    <  2:1120924‑M1/49‑31 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                        aCCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataa          <  2:2112297‑M1/50‑41 (MQ=255)

                                        aCCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaata           <  2:796624‑M1/49‑40 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

                                              ggAGatgcgtataaagcttccgatgagcgtttatatgtcaataagcagaa    <  2:104844‑M1/50‑47 (MQ=255)

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

                                               gagatgcgtataaagcttccgatgagcgtttatatgtcaataagcagaac   <  2:5606771‑M1/50‑48 (MQ=255)

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

                                                agatgcgtataaagcttccgatgagcgtttatatgtcaataagcagaaca  <  1:2958571‑M1/50‑49 (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: `‑`