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

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

Only 100 of 168 total aligned reads displayed.

GGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAA  >  NC_000913/1873418‑1873516

ggccgcaatcgaactttatctggttttctcgtttcactaaccgaaggaga                                                   <  1:3698982‑M2/1‑1 (MQ=255)

 gccgcaatcgaactttatctggttttctcgtttcactcaccgaaggagAt                                                  >  2:2266992‑M2/49‑50 (MQ=255)

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

 gccgcaatcgaactttatctggttttctcgtttcactaaccgaaggagAt                                                  >  1:1524744‑M2/49‑50 (MQ=255)

 gccgcaatcgaactttatctggttttctcgtttcactaaccgaagcagAt                                                  >  1:2788261‑M2/49‑50 (MQ=255)

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

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

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

     caatcgaactttatctggttttctcgtttcactaaccgaaggagATGCGt                                              <  2:2489884‑M2/6‑1 (MQ=255)

      aatcgaactttatctggttttctcgtttcactaaccgaaggagATGCGta                                             >  2:3224910‑M2/44‑50 (MQ=255)

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

        tcgaactttatctggttttctcgtttcactaaccgaaggagATGCGtata                                           <  1:2332509‑M2/9‑1 (MQ=255)

        tcgaactttatctggttttctcgtttcactaaccgaaggagATGCGtata                                           <  1:878960‑M2/9‑1 (MQ=255)

         cgaactttatctggttttctcgtttcactaaccgaaggagATGCGtata                                           >  1:737846‑M2/41‑49 (MQ=255)

          gaactttatctggttttctcgtttcactaaccgaaggagATGCGtata                                           >  2:42835‑M2/40‑48 (MQ=255)

          gaactttatctggttttctcgtttcactaaccgaaggagATGCGtata                                           >  2:956103‑M2/40‑48 (MQ=255)

          gaactttatctggttttctcgtttcactaaccgaaggagATGCGTATaaa                                         <  1:1371623‑M2/11‑1 (MQ=255)

           aactttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAg                                        >  2:3518762‑M2/39‑50 (MQ=255)

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

           aactttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAg                                        >  2:814740‑M2/39‑50 (MQ=255)

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

            actttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGc                                       <  2:1849425‑M2/13‑1 (MQ=255)

             ctttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCt                                      <  2:2473370‑M2/14‑1 (MQ=255)

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

                tatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTcc                                   >  1:1920557‑M2/34‑50 (MQ=255)

                tatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTc                                    <  2:876681‑M2/16‑1 (MQ=255)

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

                 atctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTcc                                   >  2:3103809‑M2/33‑49 (MQ=255)

                    tggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGAt                                >  2:741106‑M2/30‑49 (MQ=255)

                     ggttttctcgtttcactaaccgaaggagATGAGTATAAAGCTTCCGAt                                >  1:495629‑M2/29‑48 (MQ=255)

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

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

                      gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATg                               <  2:946324‑M2/21‑1 (MQ=255)

                      gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATg                               <  1:1554567‑M2/21‑1 (MQ=255)

                      gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATg                               <  2:2905831‑M2/21‑1 (MQ=255)

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

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

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

                       ttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGc                            >  2:2749585‑M2/27‑50 (MQ=255)

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

                       ttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGc                            >  2:2246950‑M2/27‑50 (MQ=255)

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

                       ttttcgcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGc                            <  1:3406650‑M2/24‑1 (MQ=255)

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

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

                         ttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGt                          >  2:3373800‑M2/25‑50 (MQ=255)

                          tctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGtt                         >  2:1689391‑M2/24‑50 (MQ=255)

                          tctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGtt                         >  2:1824378‑M2/24‑50 (MQ=255)

                          tctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGt                          <  1:3466723‑M2/26‑1 (MQ=255)

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

                           ctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGtt                         <  2:1162208‑M2/27‑1 (MQ=255)

                            tcgtttcactacccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTta                       >  1:3967755‑M2/22‑50 (MQ=255)

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

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

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

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

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

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

                               tttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTatat                    >  1:2117374‑M2/19‑50 (MQ=255)

                               tttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTatat                    >  1:2687638‑M2/19‑50 (MQ=255)

                                ttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATg                   >  2:3092547‑M2/18‑50 (MQ=255)

                                ttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATg                   >  1:4509733‑M2/18‑50 (MQ=255)

                                ttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATg                   >  1:1090231‑M2/18‑50 (MQ=255)

                                 tcactaaccgaaggagATGCGTATAAAGCTGCCGATGAGCGTTTATATg                   <  1:1235694‑M2/33‑1 (MQ=255)

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

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

                                   actaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTc                 >  1:2969463‑M2/15‑49 (MQ=255)

                                   actaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCa                <  2:538000‑M2/36‑1 (MQ=255)

                                   actaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCa                <  1:2023997‑M2/36‑1 (MQ=255)

                                    ctaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCaa               >  2:850056‑M2/14‑50 (MQ=255)

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

                                     taaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCaa               >  2:1210056‑M2/13‑49 (MQ=255)

                                      aaccgaaggagATGCGCATAAAGCTTCCGCTGCGCGTTTATATGTCAAt              >  2:1554719‑M2/12‑49 (MQ=255)

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

                                       accgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATa             <  2:3788616‑M2/39‑1 (MQ=255)

                                         tgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAg           >  2:1678733‑M2/9‑49 (MQ=255)

                                         cgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAg           >  2:2365535‑M2/9‑49 (MQ=255)

                                         cgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGc          <  2:136161‑M2/42‑1 (MQ=255)

                                          gaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAg           >  1:3668892‑M2/8‑48 (MQ=255)

                                          gaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAg           >  1:2813810‑M2/8‑48 (MQ=255)

                                          gaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAg           >  1:3077829‑M2/8‑48 (MQ=255)

                                           aaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAg           >  2:3692220‑M2/7‑47 (MQ=255)

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

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

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

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

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

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

                                            aggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAg        >  2:166690‑M2/6‑49 (MQ=255)

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

                                             ggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGaa      >  1:1062264‑M2/5‑50 (MQ=255)

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

                                                gATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACaa   <  2:1455113‑M2/49‑1 (MQ=255)

                                                gATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACaa   <  1:680626‑M2/49‑1 (MQ=255)

                                                gATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACa    >  2:330875‑M2/2‑49 (MQ=255)

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

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

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

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

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

GGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAA  >  NC_000913/1873418‑1873516

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