BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,618,992	Δ6,688 bp		rhsB–[yhhJ]	rhsB, yhhH, yhhI, [yhhJ]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3618992	3625679	6688	77 [0]	[0] 79	rhsB–[yhhJ]	rhsB,yhhH,yhhI,[yhhJ]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3618991	0 (0.000)	77 (0.890)	48/96	0.0	100%	intergenic (+2/‑201)	nikR/rhsB	DNA‑binding transcriptional repressor NikR/rhs element protein RhsB
*	?	NC_000913	3625680 =	NA (NA)	77 (0.890)	48/96	0.0	100%	coding (1124/1125 nt)	yhhJ	ABC transporter family protein YhhJ

TTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGA  >  NC_000913/3625631‑3625728

gcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtt                                                  <  1:1138573‑M2/1‑1 (MQ=255)

 cggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTa                                                 <  1:4372099‑M2/2‑1 (MQ=255)

  ggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTAc                                                >  1:2896786‑M2/48‑50 (MQ=255)

  ggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTAc                                                >  1:1236018‑M2/48‑50 (MQ=255)

  ggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTAc                                                >  1:5270671‑M2/48‑50 (MQ=255)

  ggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTAc                                                <  1:1639141‑M2/3‑1 (MQ=255)

  ggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTAc                                                <  1:1785211‑M2/3‑1 (MQ=255)

  ggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTAc                                                >  1:4901925‑M2/48‑50 (MQ=255)

     gtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGcc                                             <  1:6489925‑M2/6‑1 (MQ=255)

     gtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGcc                                             <  1:1335807‑M2/6‑1 (MQ=255)

        cggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGCCAtt                                          <  1:4257184‑M2/9‑1 (MQ=255)

         ggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTg                                         <  1:4508878‑M2/10‑1 (MQ=255)

         ggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTg                                         <  1:4240169‑M2/10‑1 (MQ=255)

         ggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTg                                         <  1:518572‑M2/10‑1 (MQ=255)

          gcacgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGt                                        <  1:4406775‑M2/11‑1 (MQ=255)

          gcacgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGt                                        <  1:2279473‑M2/11‑1 (MQ=255)

          gcacgggcatttgcagtgctggccgaaggaagattgagtTACGCCATTGt                                        <  1:6742983‑M2/11‑1 (MQ=255)

           cacgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTc                                       >  1:1465454‑M2/39‑50 (MQ=255)

           cacgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTc                                       >  1:5558462‑M2/39‑50 (MQ=255)

           cacgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTc                                       >  1:3986126‑M2/39‑50 (MQ=255)

            acgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTcc                                      <  1:4222882‑M2/13‑1 (MQ=255)

            acgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTcc                                      <  1:4180770‑M2/13‑1 (MQ=255)

            acgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTcc                                      <  1:2850213‑M2/13‑1 (MQ=255)

             cgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTccc                                     <  1:2545927‑M2/14‑1 (MQ=255)

             cgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTcc                                      >  1:1138793‑M2/37‑49 (MQ=255)

             cgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTcc                                      >  1:3992232‑M2/37‑49 (MQ=255)

               ggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCaa                                   <  1:6619854‑M2/16‑1 (MQ=255)

                gcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAAt                                  <  1:2134031‑M2/17‑1 (MQ=255)

                gcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAAt                                  >  1:1205940‑M2/34‑50 (MQ=255)

                gcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAAt                                  >  1:6046958‑M2/34‑50 (MQ=255)

                 catttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATc                                 <  1:848964‑M2/18‑1 (MQ=255)

                 catttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATc                                 <  1:11781‑M2/18‑1 (MQ=255)

                  atttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCg                                >  1:961969‑M2/32‑50 (MQ=255)

                   tttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGt                               >  1:1522536‑M2/31‑50 (MQ=255)

                   tttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGt                               <  1:4596902‑M2/20‑1 (MQ=255)

                      gcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCt                             <  1:3606769‑M2/22‑1 (MQ=255)

                      gcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCt                             <  1:4245703‑M2/22‑1 (MQ=255)

                       cagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTc                           >  1:5967752‑M2/27‑50 (MQ=255)

                        agtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTcc                          >  1:1863460‑M2/26‑50 (MQ=255)

                          tgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCt                         >  1:1441084‑M2/24‑49 (MQ=255)

                          tgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTg                        <  1:4857346‑M2/27‑1 (MQ=255)

                          tgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTg                        <  1:2463842‑M2/27‑1 (MQ=255)

                            cttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGaa                      >  1:5351326‑M2/22‑50 (MQ=255)

                            cttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGaa                      <  1:5341692‑M2/29‑1 (MQ=255)

                            cttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGaa                      <  1:6030958‑M2/29‑1 (MQ=255)

                            cttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGa                       <  1:4319510‑M2/28‑1 (MQ=255)

                               gccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAAt                     <  1:4886578‑M2/30‑1 (MQ=255)

                               gccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATc                    <  1:594463‑M2/31‑1 (MQ=255)

                               gccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATc                    <  1:4995692‑M2/31‑1 (MQ=255)

                               gccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATc                    <  1:1808979‑M2/31‑1 (MQ=255)

                               gccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCg                   >  1:6109532‑M2/19‑50 (MQ=255)

                                ccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCg                   >  1:5150925‑M2/18‑49 (MQ=255)

                                ccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCg                   <  1:253134‑M2/32‑1 (MQ=255)

                                ccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCg                   <  1:1308286‑M2/32‑1 (MQ=255)

                                 cgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgca                 <  1:5372533‑M2/34‑1 (MQ=255)

                                 cgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCg                   >  1:5948614‑M2/17‑48 (MQ=255)

                                 cgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCg                   >  1:4651066‑M2/17‑48 (MQ=255)

                                  gaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgcag                <  1:4509477‑M2/35‑1 (MQ=255)

                                    aggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgcagca              <  1:5981581‑M2/37‑1 (MQ=255)

                                    aggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgcagca              <  1:934324‑M2/37‑1 (MQ=255)

                                    aggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgcagca              <  1:2089204‑M2/37‑1 (MQ=255)

                                    aggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgcag                >  1:3944777‑M2/14‑48 (MQ=255)

                                     ggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGcagcag             >  1:2493643‑M2/13‑50 (MQ=255)

                                     ggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGcagcag             >  1:462302‑M2/13‑50 (MQ=255)

                                      gaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCagcagc            >  1:3630564‑M2/12‑50 (MQ=255)

                                       aagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAgcg           >  1:4840283‑M2/11‑50 (MQ=255)

                                         gattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAgcgc          >  1:1985644‑M2/9‑49 (MQ=255)

                                         gattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCa         <  1:4494365‑M2/42‑1 (MQ=255)

                                          attgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCaa        >  1:5107396‑M2/8‑50 (MQ=255)

                                           ttgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAAt       >  1:4367914‑M2/7‑50 (MQ=255)

                                           ttgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAAt       >  1:5768201‑M2/7‑50 (MQ=255)

                                             gagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATgg     >  1:4260638‑M2/5‑50 (MQ=255)

                                             gagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATgg     >  1:6256883‑M2/5‑50 (MQ=255)

                                             gagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATg      <  1:2080641‑M2/45‑1 (MQ=255)

                                              agtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATgg     <  1:5754013‑M2/46‑1 (MQ=255)

                                              agtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGt    >  1:4407191‑M2/4‑50 (MQ=255)

                                              agtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGt    >  1:3685262‑M2/4‑50 (MQ=255)

                                                 tACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGa  >  1:486059/1‑49 (MQ=255)

                                                 tACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGa  >  1:6340017/1‑49 (MQ=255)

TTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGA  >  NC_000913/3625631‑3625728

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