BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,720,634	Δ1,442 bp	IS150‑mediated	insJ–insK	insJ, insK

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3720634	3722075	1442	88 [0]	[0] 87	insJ–insK	insJ,insK

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3720633	0 (0.000)	87 (0.890)	54/98	0.0	100%	noncoding (1/1443 nt)	IS150	repeat region
*	?	NC_000913	3722076 =	0 (0.000)	87 (0.890)	54/98	0.0	100%	intergenic (+27/+252)	insK/glyS	IS150 conserved protein InsB/glycine‑‑tRNA ligase subunit beta

GACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTG  >  NC_000913/3722028‑3722124

tatatttctgcgcattcccagcctcctgtgttgatttccaacgagtataa                                                 >  1:5390924‑M2/49‑50 (MQ=255)

  tatttctgcgcattcccagcctcctgtgttgatttccaacgagtatAAc                                                <  1:6325180‑M2/3‑1 (MQ=255)

  tatttctgcgcattcccagcctcctgtgttgatttccaacgagtatAAc                                                <  1:1940447‑M2/3‑1 (MQ=255)

  tatttctgcgcattcccagcctcctgtgttgatttccaacgagtatAACt                                               >  1:127983‑M2/47‑50 (MQ=255)

  tatttctgcgcattcccagcctcctgtgttgatttccaacgagtatAACt                                               <  1:2100711‑M2/4‑1 (MQ=255)

   atttctgcgcattcccagcctcctgtgttgatttccaacgagtatAACtt                                              >  1:7406155‑M2/46‑50 (MQ=255)

     ttctgcgcattcccagcctcctgtgttgatttccaacgagtatAACttt                                             >  1:6692748‑M2/44‑49 (MQ=255)

     ttctgcgcattcccagcctcctgtgttgatttccaacgagtatAACttt                                             >  1:4588673‑M2/44‑49 (MQ=255)

     ttctgcgcattcccagcctcctgtgttgatttccaacgagtatAACttt                                             >  1:6767880‑M2/44‑49 (MQ=255)

      tctgcgcattcccagcctcctgtgttgatttccaacgagtatAACTTTg                                            <  1:4079106‑M2/7‑1 (MQ=255)

      tctgcgcattcccagcctcctgtgttgatttccaacgagtatAACTTTGa                                           >  1:5420523‑M2/43‑50 (MQ=255)

      tctgcgcattcccagcctcctgtgttgatttccaacgagtatAACTTTGa                                           >  1:6311609‑M2/43‑50 (MQ=255)

          cgcattcccagcctcctgtgttgatttccaacgagtatAACTTTGAtt                                         <  1:2224843‑M2/10‑1 (MQ=255)

          cgcattcccagcctcctgtgttgatttccaacgagtatAACTTTGATTta                                       <  1:2833557‑M2/12‑1 (MQ=255)

          cgcattcccagcctcctgtgttgatttccaacgagtatAACTTTGATTta                                       >  1:3035188‑M2/39‑50 (MQ=255)

          cgcattcccagcctcctgtgttgatttccaacgagtatAACTTTGATTta                                       <  1:4235289‑M2/12‑1 (MQ=255)

          cgcattcccagcctcctgtgttgatttccaacgagtatAACTTTGATTta                                       <  1:7443578‑M2/12‑1 (MQ=255)

           gcattcccagcctcctgtgttgatttccaacgagtatAACTTTGATTtat                                      <  1:5885046‑M2/13‑1 (MQ=255)

           gcattcccagcctcctgtgttgatttccaacgagtatAACTTTGATTtat                                      <  1:2158982‑M2/13‑1 (MQ=255)

           gcattcccagcctcctgtgttgatttccaacgagtatAACTTTGATTtat                                      <  1:2900451‑M2/13‑1 (MQ=255)

           gcattcccagcctcctgtgttgatttccaacgagtatAACTTTGATTta                                       >  1:979530‑M2/38‑49 (MQ=255)

            cattcccagcctcctgtgttgatttccaacgagtatAACTTTGATttatc                                     >  1:5733412‑M2/37‑49 (MQ=255)

             attcccagcctcctgtgttgatttccaacgagtatAACTTTGATTTATAg                                    <  1:2035428‑M2/15‑1 (MQ=255)

             attcccagcctcctgtgttgatttccaacgagtatAACTTTGATTTATAg                                    <  1:4634451‑M2/15‑1 (MQ=255)

             attcccagcctcctgtgttgatttccaacgagtatAACTTTGATTTATAg                                    >  1:2013996‑M2/36‑50 (MQ=255)

             attcccagcctcctgtgttgatttccaacgagtatAACTTTGATTTATAg                                    >  1:3086825‑M2/36‑50 (MQ=255)

             attcccagcctcctgtgttgatttccaacgagtatAACTTTGATTTATAg                                    >  1:4521040‑M2/36‑50 (MQ=255)

             attcccagcctcctgtgttgatttccaacgagtatAACTTTGATTTATAg                                    <  1:381358‑M2/15‑1 (MQ=255)

              ttcccagcctcctgtgttgatttccaacgagtatAACTTTGATTTATAGt                                   >  1:5835659‑M2/35‑50 (MQ=255)

               tcccagcctcctgtgttgatttccaacgagtatAACTTTGATTTATAGt                                   >  1:1503604‑M2/34‑49 (MQ=255)

                cccagcctcctgtgttgatttccaacgagtatAACTTTGATTTATAGTc                                  <  1:7313774‑M2/17‑1 (MQ=255)

                 ccagcctcctgtgttgatttccaacgagtatAACTTTGATTTCTAgtcct                                >  1:21864‑M2/32‑48 (MQ=255)

                 ccagcctcctgtgttgatttccaacgagtatAACTTTGATTTATAGTCa                                 <  1:253428‑M2/18‑1 (MQ=255)

                 ccagcctcctgtgttgatttccaacgagtatAACTTTGATTTATAGTCa                                 <  1:1667206‑M2/18‑1 (MQ=255)

                  cagcctcctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAgg                               <  1:6173754‑M2/20‑1 (MQ=255)

                  cagcctcctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAgg                               <  1:4466821‑M2/20‑1 (MQ=255)

                  cagcctcctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAgg                               <  1:2759746‑M2/20‑1 (MQ=255)

                   agcctcctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAGGt                              >  1:6491572‑M2/30‑50 (MQ=255)

                    gcctcctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAGGTg                             >  1:407235‑M2/29‑50 (MQ=255)

                     cctcctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAGGTgg                            >  1:6604364‑M2/28‑50 (MQ=255)

                      ctcctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAGGTggg                           >  1:1709661‑M2/27‑50 (MQ=255)

                       tcctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAGGTgggg                          >  1:1080681‑M2/26‑50 (MQ=255)

                       tcctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAGGTgggg                          >  1:2968446‑M2/26‑50 (MQ=255)

                       tcctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAGGTgggg                          >  1:279421‑M2/26‑50 (MQ=255)

                        cctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGc                         >  1:6509575‑M2/25‑50 (MQ=255)

                         ctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAGGTgggg                          <  1:4638492‑M2/25‑1 (MQ=255)

                         ctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCt                        <  1:2633105‑M2/27‑1 (MQ=255)

                         ctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCt                        <  1:2578820‑M2/27‑1 (MQ=255)

                         ctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCt                        <  1:288901‑M2/27‑1 (MQ=255)

                           gtgttgatttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCttt                      >  1:2100016‑M2/22‑50 (MQ=255)

                            tgttgatttccaacgagtatAACTTTGTTTTATAGTCAGGTGGGGCtttt                     >  1:7400718‑M2/21‑50 (MQ=255)

                            tgttgatttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCttt                      <  1:440833‑M2/29‑1 (MQ=255)

                            tgttgatttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCttt                      <  1:2671035‑M2/29‑1 (MQ=255)

                             gttgatttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTtc                    >  1:6779010‑M2/20‑50 (MQ=255)

                              ttgatttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTtct                   >  1:552941‑M2/19‑50 (MQ=255)

                              ttgatttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTtct                   <  1:3144905‑M2/32‑1 (MQ=255)

                              ttgatttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTtc                    <  1:7365849‑M2/31‑1 (MQ=255)

                              ttgatttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTtc                    <  1:5111379‑M2/31‑1 (MQ=255)

                               tgatttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTtctc                  >  1:1091647‑M2/18‑50 (MQ=255)

                                gatttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTtctc                  <  1:2135797‑M2/33‑1 (MQ=255)

                                 atttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTTctct                 <  1:2781185‑M2/34‑1 (MQ=255)

                                 atttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTTctct                 <  1:5174908‑M2/34‑1 (MQ=255)

                                 atttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCtctg                >  1:5464434‑M2/16‑50 (MQ=255)

                                 atttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCtctg                <  1:5904960‑M2/35‑1 (MQ=255)

                                  tttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCtctgt               >  1:7005001‑M2/15‑50 (MQ=255)

                                  tttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCtctgt               >  1:3672293‑M2/15‑50 (MQ=255)

                                    tccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCtctgtc              <  1:2431463‑M2/37‑1 (MQ=255)

                                       aacgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGcc          >  1:4765998‑M2/10‑50 (MQ=255)

                                        acgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCt         >  1:6649629‑M2/9‑50 (MQ=255)

                                        acgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCt         <  1:6755557‑M2/42‑1 (MQ=255)

                                         cgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCtt        >  1:60184‑M2/8‑50 (MQ=255)

                                         cgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCtt        <  1:3691467‑M2/43‑1 (MQ=255)

                                         cgagtatAACGTTGAGTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCtt        <  1:5097468‑M2/43‑1 (MQ=255)

                                          gagtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCttt       <  1:790012‑M2/44‑1 (MQ=255)

                                          gagtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCttt       <  1:6619547‑M2/44‑1 (MQ=255)

                                          gagtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCttt       >  1:6627732‑M2/7‑50 (MQ=255)

                                          gagtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCttt       >  1:3134666‑M2/7‑50 (MQ=255)

                                           agtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTc      <  1:1004728‑M2/45‑1 (MQ=255)

                                           agtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTc      <  1:3949874‑M2/45‑1 (MQ=255)

                                            gtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCg     >  1:1936241‑M2/5‑50 (MQ=255)

                                            gtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCg     <  1:433271‑M2/46‑1 (MQ=255)

                                            gtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCg     <  1:6177896‑M2/46‑1 (MQ=255)

                                             tatAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCgg    <  1:1739797‑M2/47‑1 (MQ=255)

                                             tatAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCgg    >  1:4819978‑M2/4‑50 (MQ=255)

                                              atAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGt   >  1:5336038‑M2/3‑50 (MQ=255)

                                               tAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGt   <  1:4627357‑M2/48‑1 (MQ=255)

                                               tAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTg  <  1:6131354‑M2/49‑1 (MQ=255)

GACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTG  >  NC_000913/3722028‑3722124

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