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

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

Only 100 of 132 total aligned reads displayed.

GCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAA  >  NC_000913/1873419‑1873516

tccgcaatcgaactttatctggttttctcgtttcactaaccgaaggagAt                                                  <  1:5899631‑M2/2‑1 (MQ=255)

gccgcaatcgaactttatctggttttctcgtttcactaaccgaaggagAt                                                  <  1:5348140‑M2/2‑1 (MQ=255)

gccgcaatcgaactttatctggttttctcgtttcactaaccgaaggagAt                                                  <  1:7591036‑M2/2‑1 (MQ=255)

 ccgcaatcgaactttatctggttttctcgtttcactaaccgaaggagATg                                                 <  1:7283640‑M2/3‑1 (MQ=255)

 ccgcaatcgaactttatctggttttctcgtttcactaaccgaaggagATg                                                 <  1:2487149‑M2/3‑1 (MQ=255)

 ccgcaatcgaactttatctggttttctcgtttcactaaccgaaggagATg                                                 <  1:4131023‑M2/3‑1 (MQ=255)

 ccgcaatcgaactttatctggttttctcgtttcactaaccgaaggagATg                                                 >  1:3602801‑M2/48‑50 (MQ=255)

 ccgcaatcgaactttatctggttttctcgtttcactaaccgaaggagATg                                                 >  1:6112049‑M2/48‑50 (MQ=255)

  cgcaatcgaactttatctggttttctcgtttcactaaccgaaggagATGc                                                <  1:4140022‑M2/4‑1 (MQ=255)

   gcaatcgaactttatctggttttctcgtttcactaaccgaaggagATGCg                                               >  1:6557523‑M2/46‑50 (MQ=255)

   gcaatcgaactttatctggttttctcgtttcactaaccgaaggagATGCg                                               >  1:5328945‑M2/46‑50 (MQ=255)

    caatcgaactttatctggttttctcgtttcactaaccgaaggagATGCg                                               >  1:1604379‑M2/45‑49 (MQ=255)

     aatcgaactttatctggttttctcgtttcactaaccgaaggagATGCGta                                             <  1:8345696‑M2/7‑1 (MQ=255)

     aatcgaactttatctggttttctcgtttcactaaccgaaggagATGCGta                                             >  1:7611742‑M2/44‑50 (MQ=255)

     aatcgaactttatctggttttctcgtttcactaaccgaaggagATGCGta                                             >  1:726214‑M2/44‑50 (MQ=255)

     aatcgaactttatctggttttctcgtttcactaaccgaaggagATGCGta                                             <  1:6241851‑M2/7‑1 (MQ=255)

      atcgaactttatctggttttctcgtttcactaaccgaaggagATGCGtat                                            <  1:7732330‑M2/8‑1 (MQ=255)

         gaactttatctggttttctcgtttcactaaccgaaggagATGCGtata                                           >  1:7467051‑M2/40‑48 (MQ=255)

         gaactttatctggttttctcgtttcactaaccgaaggagATGCGTATaa                                          <  1:4254404‑M2/10‑1 (MQ=255)

          aactttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAg                                        >  1:4903655‑M2/39‑50 (MQ=255)

          aactttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAg                                        <  1:288884‑M2/12‑1 (MQ=255)

          aactttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAg                                        >  1:3724724‑M2/39‑50 (MQ=255)

          aactttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAg                                        >  1:1204928‑M2/39‑50 (MQ=255)

           actttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGc                                       >  1:7198941‑M2/38‑50 (MQ=255)

           actttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGc                                       >  1:390944‑M2/38‑50 (MQ=255)

           actttatctggttttctcgtttcaccaaccgaaggagATGCGTATAAAGc                                       <  1:7005637‑M2/13‑1 (MQ=255)

            ctttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCt                                      >  1:5120223‑M2/37‑50 (MQ=255)

            ctttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCt                                      >  1:5698710‑M2/37‑50 (MQ=255)

            ctttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCt                                      >  1:6522601‑M2/37‑50 (MQ=255)

              ttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTc                                    >  1:7314468‑M2/35‑50 (MQ=255)

               tatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCtt                                     <  1:6004162‑M2/15‑1 (MQ=255)

               tatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCtt                                     <  1:6511949‑M2/15‑1 (MQ=255)

                atctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTcc                                   <  1:5538505‑M2/17‑1 (MQ=255)

                  ctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGAt                                >  1:6741259‑M2/31‑50 (MQ=255)

                  ctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGAt                                >  1:2573202‑M2/31‑50 (MQ=255)

                     gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGa                                 <  1:4002250‑M2/19‑1 (MQ=255)

                     gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGa                                 <  1:1776682‑M2/19‑1 (MQ=255)

                     gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGAt                                <  1:7859733‑M2/20‑1 (MQ=255)

                     gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGa                              <  1:4193322‑M2/22‑1 (MQ=255)

                     gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAg                             <  1:3012917‑M2/23‑1 (MQ=255)

                     gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAg                             >  1:72111‑M2/28‑50 (MQ=255)

                     gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAg                             <  1:2048681‑M2/23‑1 (MQ=255)

                      ttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGc                            <  1:2122785‑M2/24‑1 (MQ=255)

                       tttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCg                           >  1:2498376‑M2/26‑50 (MQ=255)

                       tttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCg                           >  1:6006228‑M2/26‑50 (MQ=255)

                        ttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCg                           <  1:5206857‑M2/25‑1 (MQ=255)

                        ttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCg                           <  1:8123554‑M2/25‑1 (MQ=255)

                          ctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGTGCGtt                         <  1:8418427‑M2/27‑1 (MQ=255)

                          ctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGtt                         >  1:2114082‑M2/23‑49 (MQ=255)

                            cgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTta                       >  1:8371379‑M2/21‑49 (MQ=255)

                            cgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTta                       >  1:1101028‑M2/21‑49 (MQ=255)

                            cgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTta                       >  1:3287049‑M2/21‑49 (MQ=255)

                             gtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGttt                        <  1:1574180‑M2/28‑1 (MQ=255)

                             gtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTtata                     <  1:1301973‑M2/31‑1 (MQ=255)

                             gtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTtata                     >  1:7992107‑M2/20‑50 (MQ=255)

                             gtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTtat                      <  1:8239042‑M2/30‑1 (MQ=255)

                             gtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTtat                      <  1:406736‑M2/30‑1 (MQ=255)

                             gtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTta                       <  1:1070285‑M2/29‑1 (MQ=255)

                             gtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTta                       <  1:3377051‑M2/29‑1 (MQ=255)

                             gtttcactaaccgaaagagATGCGTATAAAGCTTCCGATGAGCGttt                        <  1:3050888‑M2/28‑1 (MQ=255)

                              tttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTatat                    <  1:1980664‑M2/32‑1 (MQ=255)

                              tttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTatat                    <  1:4571159‑M2/32‑1 (MQ=255)

                              tttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTatat                    <  1:7354038‑M2/32‑1 (MQ=255)

                               ttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGtttatacg                   >  1:5516579‑M2/18‑48 (MQ=255)

                                tcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGt                  >  1:7869212‑M2/17‑50 (MQ=255)

                                 cactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGt                  <  1:5599148‑M2/34‑1 (MQ=255)

                                 cactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGt                  >  1:7627348‑M2/16‑49 (MQ=255)

                                  actaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTc                 <  1:7888951‑M2/35‑1 (MQ=255)

                                  actaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTc                 <  1:528789‑M2/35‑1 (MQ=255)

                                   ctaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCaa               >  1:5530575‑M2/14‑50 (MQ=255)

                                   ctaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCaa               <  1:4121539‑M2/37‑1 (MQ=255)

                                    taaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCaa               >  1:8555324‑M2/13‑49 (MQ=255)

                                    taaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCaa               >  1:1384768‑M2/13‑49 (MQ=255)

                                    taaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCaa               >  1:6858000‑M2/13‑49 (MQ=255)

                                     aaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAAt              >  1:2178935‑M2/12‑49 (MQ=255)

                                     aaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAAt              >  1:7539381‑M2/12‑49 (MQ=255)

                                      accgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAAt              <  1:1861907‑M2/38‑1 (MQ=255)

                                      accgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAAt              <  1:2682093‑M2/38‑1 (MQ=255)

                                      accgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAAt              >  1:2693907‑M2/11‑48 (MQ=255)

                                      accgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATaa            <  1:3353849‑M2/40‑1 (MQ=255)

                                      accgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATaa            <  1:2191209‑M2/40‑1 (MQ=255)

                                      accgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATaa            <  1:3998552‑M2/40‑1 (MQ=255)

                                       ccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAg           <  1:8349617‑M2/41‑1 (MQ=255)

                                       ccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAg           <  1:7989622‑M2/41‑1 (MQ=255)

                                       ccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAg           <  1:5593260‑M2/41‑1 (MQ=255)

                                         gaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCa         <  1:6805947‑M2/43‑1 (MQ=255)

                                         gaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCa         <  1:6872956‑M2/43‑1 (MQ=255)

                                          aaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAg           >  1:181835‑M2/7‑47 (MQ=255)

                                          aaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAg        <  1:2938543‑M2/44‑1 (MQ=255)

                                           aggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGa       <  1:5951211‑M2/45‑1 (MQ=255)

                                           aggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGa       <  1:1107325‑M2/45‑1 (MQ=255)

                                           aggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGa       <  1:6025507‑M2/45‑1 (MQ=255)

                                           aggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGa       <  1:3556853‑M2/45‑1 (MQ=255)

                                              agATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAAc     >  1:5202729‑M2/3‑49 (MQ=255)

                                              agATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAAc     >  1:1368104‑M2/3‑49 (MQ=255)

                                              agATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACa    <  1:3010399‑M2/48‑1 (MQ=255)

                                              agATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACa    <  1:5871107‑M2/48‑1 (MQ=255)

                                                atgCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACaaa  <  1:4692117/50‑1 (MQ=255)

                                                atgCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACaa   >  1:529166/1‑49 (MQ=255)

                                                atgCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACaa   >  1:1703387/1‑49 (MQ=255)

GCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAA  >  NC_000913/1873419‑1873516

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