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

GGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGC  >  NC_000913/257844‑257971

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

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

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

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

        gAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaag                                                         <  1:1536292‑M1/65‑10 (MQ=255)

         aaGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaaga                                                        <  1:2481851‑M1/65‑11 (MQ=255)

          agagCAGATTGATCAAAAAATTTACCGCCCTAGGCCCGTATATTCGTGAAGGTAagtgcaaagat                                                       >  1:2202458‑M1/1‑54 (MQ=255)

          agagCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagat                                                       >  1:891171‑M1/1‑54 (MQ=255)

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

          agagCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagat                                                       >  1:889054‑M1/1‑54 (MQ=255)

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

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

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

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

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

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

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

                    gATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattct                                             >  1:251334‑M1/1‑44 (MQ=255)

                    gATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattct                                             >  1:1954873‑M1/1‑44 (MQ=255)

                     aTCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctt                                            <  1:571736‑M1/65‑23 (MQ=255)

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

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

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

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

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

                               tgaccGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgt                                  <  1:153874‑M1/63‑33 (MQ=255)

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

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

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

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

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

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

                                    gcgctagGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggc                             <  1:2245114‑M1/62‑38 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

                                            cccGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatcag                     >  1:3073668‑M1/1‑20 (MQ=255)

                                                tataTTCGTGAAGGTAagtgcaaagataatcgattctttttcgcttgtctggctgtatgcgtcaa                 >  1:188315‑M1/1‑16 (MQ=255)

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

                                                 atatTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcttg                >  1:621462‑M1/1‑15 (MQ=255)

                                                 atatTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaac                >  1:29283‑M1/1‑15 (MQ=255)

                                                 atatTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaac                >  1:1058080‑M1/1‑15 (MQ=255)

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

                                                     tCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtga            >  1:1239445‑M1/1‑11 (MQ=255)

                                                           aGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccgg      >  1:2803486‑M1/1‑5 (MQ=255)

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

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

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

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

                                                             gTAagtgcaaagataatcgattctttttcgattgtctggccgtatgcgtcaacgtgaaaccggca    >  1:2784756‑M1/1‑3 (MQ=255)

                                                               aagtgcaaagataatcgattcttgttcgattgtctggctgtatgcgtcaacgtgaaaccggcacc  <  1:3180245‑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` `< 15 ≤` `ATCG/ATCG` `< 24 ≤` `ATCG/ATCG` `< 31 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`