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	57 [0]	[0] 57	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)	55 (0.770)	40/96	0.1	100%	intergenic (+2/‑201)	nikR/rhsB	DNA‑binding transcriptional repressor NikR/rhs element protein RhsB
*	?	NC_000913	3625680 =	NA (NA)	55 (0.770)	40/96	0.1	100%	coding (1124/1125 nt)	yhhJ	ABC transporter family protein YhhJ

GCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACAC  >  NC_000913/3618943‑3619040

gcggcgTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGT                                                   >  1:1235698‑M1/1‑49 (MQ=255)

 cggcgTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTta                                                 <  1:1087644‑M1/50‑3 (MQ=255)

 cggcgTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTta                                                 >  1:379223‑M1/1‑48 (MQ=255)

  ggcgTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtac                                                <  1:5358629‑M1/50‑4 (MQ=255)

    cgTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgc                                              <  1:1409919‑M1/50‑6 (MQ=255)

    cgTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgc                                              >  1:3632848‑M1/1‑45 (MQ=255)

     gTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgc                                              >  1:1092105‑M1/1‑44 (MQ=255)

     gTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgc                                              >  1:4108138‑M1/1‑44 (MQ=255)

     gTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgc                                              >  1:2296528‑M1/1‑44 (MQ=255)

      tGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgcca                                            >  1:936627‑M1/1‑43 (MQ=255)

      tGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgcca                                            <  1:2002734‑M1/50‑8 (MQ=255)

        cGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccatt                                          >  1:775023‑M1/1‑41 (MQ=255)

         ggCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattg                                         <  1:1203556‑M1/50‑11 (MQ=255)

         ggCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattg                                         >  1:3689684‑M1/1‑40 (MQ=255)

          gCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgt                                        <  1:500360‑M1/50‑12 (MQ=255)

           cACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtc                                       >  1:2957627‑M1/1‑38 (MQ=255)

              gggCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtccca                                    <  1:4132259‑M1/50‑16 (MQ=255)

              gggCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtccca                                    <  1:2569130‑M1/50‑16 (MQ=255)

              gggCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtccc                                     >  1:5229039‑M1/1‑35 (MQ=255)

                gCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaat                                  <  1:1305306‑M1/50‑18 (MQ=255)

                 cATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatt                                 >  1:3710041‑M1/1‑32 (MQ=255)

                 cATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatc                                 <  1:224349‑M1/50‑19 (MQ=255)

                  aTTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcg                                >  1:3539403‑M1/1‑31 (MQ=255)

                   tttGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgt                               >  1:3640274‑M1/1‑30 (MQ=255)

                   tttGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgt                               >  1:703711‑M1/1‑30 (MQ=255)

                      gCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtctt                            >  1:4973831‑M1/1‑27 (MQ=255)

                       cAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttc                           <  1:4309627‑M1/50‑25 (MQ=255)

                       cAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtctt                            <  1:2048600‑M1/49‑24 (MQ=255)

                            cTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcctgaa                      >  1:2567640‑M1/1‑21 (MQ=255)

                               gCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcctgaatcg                   >  1:4497060‑M1/1‑18 (MQ=255)

                               gCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcctgaatcg                   >  1:3124853‑M1/1‑18 (MQ=255)

                               gCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcctgaatc                    <  1:2037861‑M1/49‑32 (MQ=255)

                               gCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcctgaatc                    <  1:1906342‑M1/49‑32 (MQ=255)

                               gCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcctgaat                     <  1:5435731‑M1/48‑31 (MQ=255)

                                ccGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcctgaatcgc                  <  1:4473250‑M1/50‑34 (MQ=255)

                                ccGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcctgaatcgc                  <  1:280188‑M1/50‑34 (MQ=255)

                                ccGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcctgaatcg                   <  1:2650607‑M1/49‑33 (MQ=255)

                                 cGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcctgaatcg                   >  1:2289812‑M1/1‑16 (MQ=255)

                                 cGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcctgaatcg                   >  1:3439418‑M1/1‑16 (MQ=255)

                                 cGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcctgaatcg                   >  1:3130392‑M1/1‑16 (MQ=255)

                                  gaaggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcag                >  1:4214966‑M1/1‑15 (MQ=255)

                                  gaaggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcag                >  1:3350057‑M1/1‑15 (MQ=255)

                                   aaggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagc               <  1:366834‑M1/50‑37 (MQ=255)

                                   aaggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcag                >  1:5026549‑M1/1‑14 (MQ=255)

                                    aggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcag                >  1:5057772‑M1/1‑13 (MQ=255)

                                       aagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcg           >  1:5255947‑M1/1‑10 (MQ=255)

                                       aagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcg           >  1:5138656‑M1/1‑10 (MQ=255)

                                        agATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcg           >  1:1360906‑M1/1‑9 (MQ=255)

                                         gATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgc          >  1:3859131‑M1/1‑8 (MQ=255)

                                          aTTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaa        >  1:4631034‑M1/1‑7 (MQ=255)

                                             gAGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaatg      <  1:5520342‑M1/49‑46 (MQ=255)

                                              aGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaatggt    <  1:4691‑M1/50‑48 (MQ=255)

                                              aGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaatggt    >  1:4084121‑M1/1‑3 (MQ=255)

                                              aGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaatggt    <  1:354602‑M1/50‑48 (MQ=255)

                                              aGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaatggt    <  1:2998861‑M1/50‑48 (MQ=255)

                                              aGTtacgccattgtcccaatcgtcttcctgaatcgcagcagcgcaatggt    <  1:2164944‑M1/50‑48 (MQ=255)

                                                ttacgccattgtcccaatcttcttcctgaatcgcagcagcgcaatggtga  <  1:4668930‑M1/50‑50 (MQ=255)

GCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACAC  >  NC_000913/3618943‑3619040

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