BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	257,908	Δ776 bp		insB9–[crl]	insB9, insA9, [crl]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	257908–258676	258683	8–776	52 [0]	[0] 52	insB9–[crl]	insB9,insA9,[crl]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 257907	0 (0.000)	47 (1.090)	19/112	0.3	100%	intergenic (+8/+16)	crl/insB9	RNA polymerase holoenzyme assembly factor Crl/IS1 transposase B
*	?	NC_000913	258684 =	0 (0.000)	47 (1.090)	19/112	0.3	100%	pseudogene (9/331 nt)	crl	RNA polymerase holoenzyme assembly factor Crl

CTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACC  >  NC_000913/258620‑258747

ggacacccgaagagcagattgatcagaaaatttaccgcactaggcccgtatattcgtgaaggtaa                                                                 >  1:2364752‑M2/65‑65 (MQ=255)

ggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaa                                                                 >  1:590247‑M2/65‑65 (MQ=255)

ggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaa                                                                 >  1:1867735‑M2/65‑65 (MQ=255)

  acacccgaagagcagattgatcgaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGt                                                               <  1:1035691‑M2/3‑1 (MQ=255)

  acacccgaagagcagattgatcaaaagatttaccgcactaggcccgtatattcgtgaaggtaAGt                                                               <  1:2335363‑M2/3‑1 (MQ=255)

  acacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGt                                                               <  1:1588829‑M2/3‑1 (MQ=255)

  acacccgaagagcagattgaccaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGt                                                               <  1:671209‑M2/3‑1 (MQ=255)

     cccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAgtgtc                                                            >  1:2712499‑M2/60‑63 (MQ=255)

     cccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaATTGCa                                                            >  1:1674971‑M2/60‑65 (MQ=255)

     cccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCa                                                            >  1:1815988‑M2/60‑65 (MQ=255)

         aagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGa                                                        <  1:1983864‑M2/10‑1 (MQ=255)

          agagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGAt                                                       <  1:2935574‑M2/11‑1 (MQ=255)

          agagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggaaAGTGCAAAGAt                                                       >  1:2524707‑M2/55‑65 (MQ=255)

            agcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAaagactg                                                     >  1:3216848‑M2/53‑62 (MQ=255)

            agcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATaa                                                     <  1:2649240‑M2/13‑1 (MQ=255)

            agcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATaa                                                     >  1:2603745‑M2/53‑65 (MQ=255)

            agcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATaa                                                     >  1:3160227‑M2/53‑65 (MQ=255)

            agcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATaa                                                     >  1:3207743‑M2/53‑65 (MQ=255)

            agcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATaa                                                     >  1:512448‑M2/53‑65 (MQ=255)

            agcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATaa                                                     >  1:2476348‑M2/53‑65 (MQ=255)

             gcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAAt                                                    >  1:2467733‑M2/52‑65 (MQ=255)

                 attgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGAt                                                <  1:1758950‑M2/18‑1 (MQ=255)

                 attgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGAt                                                >  1:2611928‑M2/48‑65 (MQ=255)

                               ttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGt                                  <  1:2966975‑M2/32‑1 (MQ=255)

                               ttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGt                                  <  1:1611093‑M2/32‑1 (MQ=255)

                               ttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGt                                  <  1:2720302‑M2/32‑1 (MQ=255)

                               ttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGt                                  <  1:1129164‑M2/32‑1 (MQ=255)

                               ttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGt                                  <  1:760369‑M2/32‑1 (MQ=255)

                               ttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGt                                  <  1:2102957‑M2/32‑1 (MQ=255)

                               ttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGt                                  <  1:3786605‑M2/32‑1 (MQ=255)

                                   cgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTgg                              <  1:3690240‑M2/36‑1 (MQ=255)

                                   cgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTgg                              <  1:1914151‑M2/36‑1 (MQ=255)

                                   cgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTgg                              <  1:1614536‑M2/36‑1 (MQ=255)

                                   cgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCGTTTTCGATTGTCTgg                              <  1:1577087‑M2/36‑1 (MQ=255)

                                       ctaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGt                          <  1:925237‑M2/40‑1 (MQ=255)

                                       ctaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGt                          <  1:3863286‑M2/40‑1 (MQ=255)

                                       ctaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGt                          <  1:1608974‑M2/40‑1 (MQ=255)

                                       ctaggcccgtatatgcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGt                          <  1:1386573‑M2/40‑1 (MQ=255)

                                            cccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATgct                     >  1:592346‑M2/21‑64 (MQ=255)

                                            cccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCg                     >  1:256303‑M2/21‑65 (MQ=255)

                                            cccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCg                     >  1:3622284‑M2/21‑65 (MQ=255)

                                            cccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCg                     >  1:1322322‑M2/21‑65 (MQ=255)

                                                 atattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTcaaa                >  1:2023673‑M2/16‑64 (MQ=255)

                                                 atattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAAc                <  1:1992297‑M2/50‑1 (MQ=255)

                                                  tattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACg               <  1:1348992‑M2/51‑1 (MQ=255)

                                                  tattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACg               <  1:822767‑M2/51‑1 (MQ=255)

                                                      cgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGaa           >  1:2235089‑M2/11‑65 (MQ=255)

                                                           aggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCgg      <  1:3692232‑M2/60‑1 (MQ=255)

                                                             gtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCa    <  1:1535630‑M2/62‑1 (MQ=255)

                                                             gtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGCCTGTATGCGTCAACGGGAACCCGGCa    >  1:1066658‑M2/4‑65 (MQ=255)

                                                               aAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCAcc  >  1:3351257‑M2/2‑65 (MQ=255)

                                                               aAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCAcc  <  1:1225756‑M2/64‑1 (MQ=255)

CTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACC  >  NC_000913/258620‑258747

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