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	45 [0]	[0] 45	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)	41 (1.050)	23/118	0.2	100%	intergenic (+8/+16)	crl/insB9	RNA polymerase holoenzyme assembly factor Crl/IS1 transposase B
*	?	NC_000913	258684 =	0 (0.000)	41 (1.050)	23/118	0.2	100%	pseudogene (9/331 nt)	crl	RNA polymerase holoenzyme assembly factor Crl

GAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGC  >  NC_000913/257840‑257971

gAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTA                                                                  >  1:2374037‑M1/1‑68 (MQ=255)

   tGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagt                                                               <  1:1246646‑M1/68‑4 (MQ=255)

   tGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagt                                                               <  1:1272070‑M1/68‑4 (MQ=255)

   tGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagt                                                               <  1:507521‑M1/68‑4 (MQ=255)

   tGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagt                                                               <  1:3349249‑M1/68‑4 (MQ=255)

    ggACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtg                                                              >  1:747626‑M1/1‑64 (MQ=255)

    ggACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtg                                                              >  1:3155512‑M1/1‑64 (MQ=255)

      acacCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgca                                                            >  1:3278985‑M1/1‑62 (MQ=255)

      acacCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagggca                                                            >  1:2695752‑M1/1‑62 (MQ=255)

        acCCGAAGAGCACATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaa                                                          >  1:2124368‑M1/1‑60 (MQ=255)

         cccGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaag                                                         >  1:2516985‑M1/1‑59 (MQ=255)

             aaGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataa                                                     <  1:1687821‑M1/68‑14 (MQ=255)

              agagCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataat                                                    >  1:2974908‑M1/1‑54 (MQ=255)

              agagCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataat                                                    >  1:799689‑M1/1‑54 (MQ=255)

                agCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcg                                                  >  1:1966276‑M1/1‑52 (MQ=255)

                agCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcg                                                  >  1:3480659‑M1/1‑52 (MQ=255)

                agCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcg                                                  >  1:74973‑M1/1‑52 (MQ=255)

                    gATTGATCAAAAAATTTACCGCACTAGGCCCGTATGTTCGTGAAGGTAagtgcaaagataatcgattc                                              <  1:647207‑M1/68‑21 (MQ=255)

                        gATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttt                                          >  1:587272‑M1/1‑44 (MQ=255)

                               aaaTTTCCCGCCCTAGCCCGGTATATTCGTGAAGGTgtgcaaagataatcgattctttttcgattgt                                    <  1:3275572‑M1/68‑32 (MQ=255)

                                aaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgt                                  <  1:1462803‑M1/68‑33 (MQ=255)

                                aaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgt                                  <  1:1965949‑M1/68‑33 (MQ=255)

                                aaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgt                                  <  1:3265106‑M1/68‑33 (MQ=255)

                                    tACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctgg                              <  1:3424274‑M1/68‑37 (MQ=255)

                                        gCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgt                          <  1:2505172‑M1/68‑41 (MQ=255)

                                        gCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgt                          <  1:2263047‑M1/68‑41 (MQ=255)

                                          cctAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtat                        <  1:1221942‑M1/67‑43 (MQ=255)

                                            tAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgc                      <  1:2531786‑M1/68‑45 (MQ=255)

                                            tAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgc                      <  1:725495‑M1/68‑45 (MQ=255)

                                             aGGCGCGTATATTCGTGAAGGTGagtgcaaagataatcgattctttttcgattgtttggctgtatgcg                     <  1:2121147‑M1/68‑46 (MQ=255)

                                                cccGTATATTCGTGCAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtca                  >  1:976435‑M1/1‑20 (MQ=255)

                                                cccGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtca                  >  1:2254962‑M1/1‑20 (MQ=255)

                                                cccGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtca                  >  1:2418428‑M1/1‑20 (MQ=255)

                                                cccGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtca                  >  1:90145‑M1/1‑20 (MQ=255)

                                                  cGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaac                <  1:1309610‑M1/68‑51 (MQ=255)

                                                  cGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaac                <  1:2166550‑M1/68‑51 (MQ=255)

                                                   gTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacg               <  1:2005392‑M1/68‑52 (MQ=255)

                                                     atatTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtg             >  1:3307005‑M1/1‑15 (MQ=255)

                                                      tatTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtga            >  1:2769936‑M1/1‑14 (MQ=255)

                                                        ttCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaa          >  1:265450‑M1/1‑12 (MQ=255)

                                                         tCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaac         >  1:1554187‑M1/1‑11 (MQ=255)

                                                         tCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaac         >  1:2749233‑M1/1‑11 (MQ=255)

                                                              aaGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggca    <  1:989571‑M1/68‑63 (MQ=255)

                                                               aGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcac   >  1:3080829‑M1/1‑5 (MQ=255)

                                                                ggTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcacc  <  1:1653568‑M1/68‑65 (MQ=255)

GAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGC  >  NC_000913/257840‑257971

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