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–258675	258683	9–776	65 [0]	[0] 65	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)	58 (1.590)	19/112	0.2	100%	intergenic (+8/+16)	crl/insB9	RNA polymerase holoenzyme assembly factor Crl/IS1 transposase B
*	?	NC_000913	258684 =	0 (0.000)	58 (1.590)	19/112	0.2	100%	pseudogene (9/331 nt)	crl	RNA polymerase holoenzyme assembly factor Crl

GGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGC  >  NC_000913/257844‑257971

ggACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAa                                                                 >  1:116078‑M1/1‑64 (MQ=255)

  acacCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagt                                                               <  1:514729‑M1/65‑4 (MQ=255)

  acacCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagt                                                               <  1:1788877‑M1/65‑4 (MQ=255)

  acacCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagt                                                               <  1:271396‑M1/65‑4 (MQ=255)

  acacCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagt                                                               >  1:3048506‑M1/1‑62 (MQ=255)

        gAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaag                                                         >  1:954245‑M1/1‑56 (MQ=255)

          agagCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagat                                                       <  1:224577‑M1/65‑12 (MQ=255)

            agCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGGAGGTAagtgcaaagactg                                                     >  1:1636675‑M1/1‑52 (MQ=255)

            agCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataa                                                     >  1:1399997‑M1/1‑52 (MQ=255)

            agCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataa                                                     >  1:243551‑M1/1‑52 (MQ=255)

            agCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataa                                                     >  1:2160204‑M1/1‑52 (MQ=255)

            agCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataa                                                     >  1:462790‑M1/1‑52 (MQ=255)

            agCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataa                                                     >  1:735513‑M1/1‑52 (MQ=255)

            agCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagactg                                                     >  1:2359324‑M1/1‑52 (MQ=255)

               aGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcg                                                  >  1:2059631‑M1/1‑49 (MQ=255)

               aGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcg                                                  >  1:1731918‑M1/1‑49 (MQ=255)

                  ttGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgatt                                               >  1:486951‑M1/1‑46 (MQ=255)

                  ttGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgatt                                               >  1:1876909‑M1/1‑46 (MQ=255)

                  ttGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgatt                                               >  1:1097889‑M1/1‑46 (MQ=255)

                           aaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcga                                      >  1:782885‑M1/1‑37 (MQ=255)

                               ttACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgt                                  <  1:1704321‑M1/65‑33 (MQ=255)

                               ttACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgt                                  <  1:342680‑M1/65‑33 (MQ=255)

                               ttACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgt                                  <  1:1525427‑M1/65‑33 (MQ=255)

                               ttACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgt                                  <  1:24643‑M1/65‑33 (MQ=255)

                               ttACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgt                                  <  1:13759‑M1/65‑33 (MQ=255)

                               ttACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgt                                  <  1:101268‑M1/65‑33 (MQ=255)

                                   cGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctgg                              <  1:216538‑M1/65‑37 (MQ=255)

                                   cGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctgg                              <  1:1838137‑M1/65‑37 (MQ=255)

                                   cGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctgg                              <  1:852118‑M1/65‑37 (MQ=255)

                                    gCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggc                             <  1:207894‑M1/65‑38 (MQ=255)

                                       cTAGGCGCGTATATTCGTGAAGGTAagtgcgaagataatcgattctttttcgattgtctggctgt                          <  1:1192427‑M1/65‑41 (MQ=255)

                                       cTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgt                          <  1:1926875‑M1/65‑41 (MQ=255)

                                       cTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgt                          <  1:2816219‑M1/65‑41 (MQ=255)

                                       cTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgt                          <  1:1323901‑M1/65‑41 (MQ=255)

                                       cTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgt                          <  1:312097‑M1/65‑41 (MQ=255)

                                       cTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgt                          <  1:747193‑M1/65‑41 (MQ=255)

                                       cTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgt                          <  1:1729134‑M1/65‑41 (MQ=255)

                                       cTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgt                          <  1:177077‑M1/65‑41 (MQ=255)

                                       cTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgt                          <  1:1997060‑M1/65‑41 (MQ=255)

                                       cTAGGCCCGTATATTCGTGAAGGTAagtacaaagataatcgattctttttcgattgtctggctgt                          <  1:142740‑M1/65‑41 (MQ=255)

                                           gCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgc                      >  1:1538952‑M1/1‑21 (MQ=255)

                                            cccGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcg                     >  1:406281‑M1/1‑20 (MQ=255)

                                            cccGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcg                     >  1:816930‑M1/1‑20 (MQ=255)

                                            cccGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcg                     >  1:3108340‑M1/1‑20 (MQ=255)

                                            cccGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcg                     >  1:1019024‑M1/1‑20 (MQ=255)

                                            cccGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcg                     >  1:676902‑M1/1‑20 (MQ=255)

                                            cccGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcg                     >  1:274606‑M1/1‑20 (MQ=255)

                                            cccGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcg                     >  1:1435339‑M1/1‑20 (MQ=255)

                                            cccGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcg                     <  1:2268471‑M1/65‑46 (MQ=255)

                                            cccGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcg                     >  1:1619989‑M1/1‑20 (MQ=255)

                                                tataTTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaa                 >  1:2814692‑M1/1‑16 (MQ=255)

                                                  tatTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacg               <  1:2510924‑M1/65‑52 (MQ=255)

                                                  tatTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacg               <  1:928470‑M1/65‑52 (MQ=255)

                                                  tatTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacg               <  1:1570183‑M1/65‑52 (MQ=255)

                                                  tatTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacg               <  1:2139670‑M1/65‑52 (MQ=255)

                                                  tatTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacg               <  1:1698522‑M1/65‑52 (MQ=255)

                                                  tatTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacg               <  1:1628535‑M1/65‑52 (MQ=255)

                                                      cGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaa           >  1:1899228‑M1/1‑10 (MQ=255)

                                                           aGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccgg      <  1:96851‑M1/65‑61 (MQ=255)

                                                           aGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccgg      <  1:2419362‑M1/65‑61 (MQ=255)

                                                             gTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggca    <  1:1997487‑M1/65‑63 (MQ=255)

                                                             gTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggca    <  1:2349119‑M1/65‑63 (MQ=255)

                                                             gTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggca    >  1:1026959‑M1/1‑3 (MQ=255)

                                                             gTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggca    <  1:2618019‑M1/65‑63 (MQ=255)

                                                               aagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcacc  <  1:2473834‑M1/65‑65 (MQ=255)

GGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGC  >  NC_000913/257844‑257971

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