BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	4,185,993	Δ18 bp	coding (644‑661/4224 nt)	rpoC →	RNA polymerase subunit beta'

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	4185993	4186010	18	55 [0]	[0] 55	rpoC	RNA polymerase subunit beta'

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4185992	0 (0.000)	48 (1.030)	29/108	0.1	100%	coding (643/4224 nt)	rpoC	RNA polymerase subunit beta'
*	?	NC_000913	4186011 =	0 (0.000)	48 (1.030)	29/108	0.1	100%	coding (662/4224 nt)	rpoC	RNA polymerase subunit beta'

TCTGGAGCAAGAGTGCGAACAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGA  >  NC_000913/4185928‑4186053

tCTGGAGCAAGAGTGCGAACAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTA                                                               >  1:3980278‑M1/1‑65 (MQ=255)

     aGCAAGAGTGCGAACAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAtcaaa                                                          >  1:1056226‑M1/1‑60 (MQ=255)

     aGCAAGAGTGCGAACAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAtcaaa                                                          >  1:804481‑M1/1‑60 (MQ=255)

     aGCAAGAGTGCGAACAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAtcaaa                                                          >  1:1594600‑M1/1‑60 (MQ=255)

     aGCAAGAGTGCGAACAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAtcaaa                                                          <  1:3615809‑M1/65‑6 (MQ=255)

     aGCAAGAGTGCGAACAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAtcaaa                                                          >  1:2160450‑M1/1‑60 (MQ=255)

      gCAAGAGTGCGAACAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAtcaaac                                                         >  1:4387086‑M1/1‑59 (MQ=255)

          gagTGCGAACAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAtcaaactgct                                                     >  1:3304986‑M1/1‑55 (MQ=255)

          gagTGCGAACAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAtcaaactgct                                                     <  1:3790102‑M1/65‑11 (MQ=255)

                gAACAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAtcaaactgctggaagc                                               >  1:650150‑M1/1‑49 (MQ=255)

                 aaCAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAtcaaactgctggaagcg                                              >  1:2808202‑M1/1‑48 (MQ=255)

                  aCAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAtgcaactgctggaagcgt                                             >  1:16970‑M1/1‑47 (MQ=255)

                  aCAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAtcaaactgctggaagcgt                                             >  1:2610044‑M1/1‑47 (MQ=255)

                  aCAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAtcaaactgctggaagcgt                                             >  1:3152905‑M1/1‑47 (MQ=255)

                  aCAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAtcaaactgctggaagcgt                                             >  1:3645424‑M1/1‑47 (MQ=255)

                  aCAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAtcaaactgctggaagcgt                                             >  1:2557872‑M1/1‑47 (MQ=255)

                  aCAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAtcaaactgctggaagcgt                                             >  1:4146154‑M1/1‑47 (MQ=255)

                   cAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAtcaaactgctggaagcgtt                                            >  1:3665387‑M1/1‑46 (MQ=255)

                    aGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAtcaaactgctggaagcgttc                                           >  1:2811387‑M1/1‑45 (MQ=255)

                     gCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAtcaaactgctggaagcgttcg                                          >  1:1829974‑M1/1‑44 (MQ=255)

                       tGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAtcaaactgctggaagcgttcgtt                                        >  1:2410930‑M1/1‑42 (MQ=255)

                         cGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAtcaaactgctggaagcgttcgttca                                      <  1:1757016‑M1/65‑26 (MQ=255)

                           tGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAtcaaactgctggaagcgttcgttcagt                                    <  1:3568461‑M1/65‑28 (MQ=255)

                           tGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAtcaaactgctggaagcgttcgttcagt                                    <  1:569213‑M1/65‑28 (MQ=255)

                            gAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAtcaaactgctggaagcgttcgttcagtc                                   <  1:3373480‑M1/65‑29 (MQ=255)

                                 gCTGAACGAAACCAACTCCGAAACCAAGCGTAtcaaactgctggaagcgttcgttcagtctggta                              <  1:1054122‑M1/65‑34 (MQ=255)

                                 gCTGAACGAAACCAACTCCGAAACCAAGCGTAtcaaactgctggaagcgttcgttcagtctggta                              <  1:3197876‑M1/65‑34 (MQ=255)

                                 gCTGAACGAAACCAACTCCGAAACCAAGCGTAtcaaactgctggaagcgttcgttcagtctggta                              >  1:3018670‑M1/1‑32 (MQ=255)

                                 gCTGAACGAAACCAACTCCGAAACCAAGCGTAtcaaactgctggaagcgttcgttcagtctggta                              >  1:4042340‑M1/1‑32 (MQ=255)

                                      aCGAAACCAACTCCGAAACCAAGCGTAtcaaactgctggaagcgttcgttcagtctggtaacaaa                         <  1:3268344‑M1/65‑39 (MQ=255)

                                      aCGAAACCAACTCCGAAACCAAGCGTAtcaaactgctggaagcgttcgttcagtctggtaacaaa                         <  1:1368939‑M1/65‑39 (MQ=255)

                                        gAAACCAACTCCGAAACCAAGCGTAtcaaactgctggaagcgttcgttcagtctggtaacaaacc                       >  1:3074794‑M1/1‑25 (MQ=255)

                                        gAAACCAACTCCGAAACCAAGCGTAtcaaactgctggaagcgttcgttcagtctggtaacaaacc                       <  1:3631271‑M1/65‑41 (MQ=255)

                                           aCCAACTCCGAAACCAAGCGTAtcaaactgctggaagcgttcgttcagtctggtaacaaaccaga                    <  1:3180464‑M1/65‑44 (MQ=255)

                                           aCCAACTCCGAAACCAAGCGTAtcaaactgctggaagcgttcgttcagtctggtaacaaaccaga                    <  1:151284‑M1/65‑44 (MQ=255)

                                           aCCAACTCCGAAACCAAGCGTAtcaaactgctggaagcgttcgttcagtctggtaacaaaccaga                    >  1:2659859‑M1/1‑22 (MQ=255)

                                           aCCAACTCCGAAACCAAGCGTAtcaaactgctggaagcgttcgttcagtctggtaacaaaccaga                    >  1:2215345‑M1/1‑22 (MQ=255)

                                             cAACTCCGAAACCAAGCGTAtcaaactgctggaagcgttcgttcagtctggtaacaaaccagagt                  <  1:316789‑M1/65‑46 (MQ=255)

                                              aaCTCCGAAACCAAGCGTAtcaaactgctggaagcgttcgttcagtctggtaacaaaccagagtg                 <  1:1840951‑M1/65‑47 (MQ=255)

                                               aCTCCGAAACCAAGCGTAtcaaactgctggaagcgttcgttcagtctggtaacaaaccagagtgg                >  1:3217774‑M1/1‑18 (MQ=255)

                                               aCTCCGAAACCAAGCGTAtcaaactgctggaagcgttcgttcagtctggtaacaaaccagagtgg                <  1:766890‑M1/65‑48 (MQ=255)

                                               aCTCCGAAACCAAGCGTAtcaaactgctggaagcgttcgttcagtctggtaacaaaccagagtgg                <  1:1640044‑M1/65‑48 (MQ=255)

                                                 tCCGAAACCAAGCGTAtcaaactgctggaagcgttcgttcagtctggtaacaaaccagagtggat              >  1:2329329‑M1/1‑16 (MQ=255)

                                                   cGAAACCAAGCGTAtcaaactgctggaagcgttcgttcagtctggtaacaaaccagagtggatga            <  1:937030‑M1/65‑52 (MQ=255)

                                                    gAAACCAAGCGTAtcaaactgctggaagcgttcgttcagtctggtaacaaaccagagtggatgat           <  1:4126241‑M1/65‑53 (MQ=255)

                                                    gAAACCAAGCGTAtcaaactgctggaagcgttcgttcagtctggtaacaaaccagagtggatgat           <  1:2478274‑M1/65‑53 (MQ=255)

                                                    gAAACCAAGCGTAtcaaactgctggaagcgttcgttcagtctggtaacaaaccagagtggatgat           <  1:4272395‑M1/65‑53 (MQ=255)

                                                    gAAACCAAGCGTAtcaaactgctggaagcgttcgttcagtctggtaacaaaccagagtggatgat           <  1:2200037‑M1/65‑53 (MQ=255)

                                                    gAAACCAAGCGTAtcaaactgctggaagcgttcgttcagtctggtaacaaaccagagtggatgat           <  1:2809796‑M1/65‑53 (MQ=255)

                                                          aaGCGTAtcaaactgctggaagcgttcgttcagtctggtaacaaaccagagtggatgatcctgac     <  1:3833508‑M1/65‑59 (MQ=255)

                                                           aGCGTAtcaaactgctggaagcgttcgttcagtctggtaacaaaccagagtggatgatcctgacc    <  1:750224‑M1/65‑60 (MQ=255)

                                                           aGCGTAtcaaactgctggaagcgttcgttcagtctggtaacaaaccagagtggatgatcctgacc    <  1:3338874‑M1/65‑60 (MQ=255)

                                                            gCGTAtcaaactgctggaagcgttcgttcagtctggtaacaaaccagagtggatgatcctgaccg   <  1:3933485‑M1/65‑61 (MQ=255)

                                                            gCGTAtcaaactgctggaagcgttcgttcagtctggtaacaaaccagagtggatgatcctgaccg   <  1:3060599‑M1/65‑61 (MQ=255)

                                                             cGTAtcaaactgctggaagcgttcgttcagtctggtaacaaaccagagtggatgatcctgaccgt  >  1:4105074‑M1/1‑4 (MQ=255)

TCTGGAGCAAGAGTGCGAACAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGA  >  NC_000913/4185928‑4186053

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