BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC JC	REL606	2,397,529	IS150 (–) Δ14 bp	coding (1124‑1137/2646 nt)	yfcTU ←	predicted outer membrane export usher protein

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	REL606	2397529	2397542	14	55 [0]	[0] 90	yfcTU	predicted outer membrane export usher protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 2397528	0 (0.000)	54 (0.690)	43/98	0.6	100%	coding (1138/2646 nt)	yfcTU	predicted outer membrane export usher protein
*	?	REL606	2774435 =	NA (NA)	54 (0.690)	43/98	0.6	100%	noncoding (1443/1443 nt)	IS150	repeat region
*	?	REL606	2397543 =	0 (0.000)	83 (1.090)	52/96	0.2	100%	coding (1123/2646 nt)	yfcTU	predicted outer membrane export usher protein
*	?	REL606	= 2775877	NA (NA)	83 (1.090)	52/96	0.2	100%	noncoding (1/1443 nt)	IS150	repeat region

GGAGGCTTCAGCACCAGAGAAAAATCCCCCCTCGACATGGTGGCCCCACTCTTGCGGACGTCCCATCATGACTTTATAACGAACCTGGCCTGGGCGA  >  REL606/2397478‑2397574

ggaggCTTCAGCACCAGAGAAAAATCCCCCCTCGACATGGTGGCCCCActc                                                >  1:3000213/1‑51 (MQ=255)

  aGGCTTCAGCACCAGAGAAAAATCCCCCCTCGACATGGTGGCCCCACTCtg                                              >  1:5088670‑M1/1‑49 (MQ=255)

  aGGCTTCAGCACCAGAGAAAAATCCCCCCTCGACATGGTGGCCCCACTCtg                                              <  1:6603582‑M1/51‑3 (MQ=255)

   ggCTTCAGCACCAGAGAAAAATCCCCCCTCGACATGGTGGCCCCACTCtgt                                             >  1:2542310‑M1/1‑48 (MQ=255)

   ggCTTCAGCACCAGAGAAAAATCCCCCCTCGACATGGTGGCCCCACTCtgt                                             <  1:1841845‑M1/51‑4 (MQ=255)

     cTTCAGCACCAGAGAAAAATCCCCCCTCGACATGGTGGCCCCACTCtgtac                                           <  1:2627770‑M1/51‑6 (MQ=255)

      ttCAGCACCAGAGAAAAATCCCCCCTCGACATGGTGGCCCCACTCtgtact                                          >  1:2485518‑M1/1‑45 (MQ=255)

      ttCAGCACCAGAGAAAAATCCCCCCTCGACATGGTGGCCCCACTCtgtact                                          <  1:6269662‑M1/51‑7 (MQ=255)

           cACCAGAGAAAAATCCCCCCTCGACATGGTGGCCCCACTCtgtactgaccc                                     >  1:1488407‑M1/1‑40 (MQ=255)

           cACCAGAGAAAAATCCCCCCTCGACATGGTGGCCCCACTCtgtactgaccc                                     <  1:6851666‑M1/51‑12 (MQ=255)

            aCCAGAGAAAAATCCCCCCTCGACATGGTGGCCCCACTCtgtactgacccc                                    >  1:6839150‑M1/1‑39 (MQ=255)

            aCCAGAGAAAAATCCCCCCTCGACATGGTGGCCCCACTCtgtactgacccc                                    >  1:5451778‑M1/1‑39 (MQ=255)

             ccAGAGAAAAATCCCCCCTCGACATGGTGGCCCCACTCtgtactgacccca                                   <  1:4161441‑M1/51‑14 (MQ=255)

              cAGAGAAAAATCCCCCCTCGACATGGTGGCCCCACTCtgtactgaccccaa                                  <  1:3665930‑M1/51‑15 (MQ=255)

               agagAAAAATCCCCCCTCGACATGGTGGCCCCACTCtgtactgaccccaaa                                 <  1:5779631‑M1/51‑16 (MQ=255)

               agagAAAAATCCCCCCTCGACATGGTGGCCCCACTCtgtactgaccccaaa                                 <  1:3916039‑M1/51‑16 (MQ=255)

               agagAAAAATCCCCCCTCGACATGGTGGCCCCACTCtgtactgaccccaaa                                 <  1:6609783‑M1/51‑16 (MQ=255)

                gagaAAAATCCCCCCTCGACATGGTGGCCCCACTCtgtactgaccccaaaa                                >  1:2764003‑M1/1‑35 (MQ=255)

                 agaAAAATCCCCCCTCGACATGGTGGCCCCACTCtgtactgaccccaaaaa                               >  1:2096668‑M1/1‑34 (MQ=255)

                 agaAAAATCCCCCCTCGACATGGTGGCCCCACTCtgtactgaccccaaaaa                               <  1:6632142‑M1/51‑18 (MQ=255)

                   aaaaaTCCCCCCTCGACATGGTGGCCCCACTCtgtactgaccccaaaaagt                             >  1:1382050‑M1/1‑32 (MQ=255)

                   aaaaaTCCCCCCTCGACATGGTGGCCCCACTCtgtactgaccccaaaaagt                             >  1:396857‑M1/1‑32 (MQ=255)

                    aaaaTCCCCCCTCGACATGGTGGCCCCACTCtgtactgaccccaaaaagtt                            >  1:1565489‑M1/1‑31 (MQ=255)

                    aaaaTCCCCCCTCGACATGGTGGCCCCACTCtgtactgaccccaaaaagtt                            >  1:3800926‑M1/1‑31 (MQ=255)

                      aaTCCCCCCTCGACATGGTGGCCCCACTCtgtactgaccccaaaaagttgg                          <  1:2700337‑M1/51‑23 (MQ=255)

                       aTCCCCCCTCGACATGGTGGCCCCACTCtgtactgaccccaaaaagttgga                         >  1:710027‑M1/1‑28 (MQ=255)

                         ccccccTCGACATGGTGGCCCCACTCtgtactgaccccaaaaagttggaca                       >  1:5772451‑M1/1‑26 (MQ=255)

                         ccccccTCGACATGGTGGCCCCACTCtgtactgaccccaaaaagttggaca                       >  1:4710494‑M1/1‑26 (MQ=255)

                         ccccccTCGACATGGTGGCCCCACTCtgtactgaccccaaaaagttggac                        <  1:811240‑M1/50‑25 (MQ=255)

                            cccTCGACATGGTGGCCCCACTCtgtactgaccccaaaaagttggacagtt                    >  1:6583840‑M1/1‑23 (MQ=255)

                            cccTCGACATGGTGGCCCCACTCtgtactgaccccaaaaagtt                            <  1:519540‑M1/43‑21 (MQ=255)

                             ccTCGACATGGTGGCCCCACTCtgtactgaccccaaaagttggacagttaa                   <  1:1263470‑M1/51‑30 (MQ=255)

                             ccTCGACATGGTGGCCCCACTCtgtactgaccccaaaaagttggacagtta                   >  1:3262745‑M1/1‑22 (MQ=255)

                             ccTCGACATGGTGGCCCCACTCtgtactgaccccaaaaagttggacagtta                   <  1:2609976‑M1/51‑30 (MQ=255)

                             ccTCGACATGGTGGCCCCACTCtgtactgaccccaaaaagttggacagtta                   <  1:4455531‑M1/51‑30 (MQ=255)

                              cTCGACATGGTGGCCCCACTCtgtactgaccccaaaaagttggacagttaa                  <  1:3011246‑M1/51‑31 (MQ=255)

                              cTCGACATGGTGGCCCCACTCtgtactgaccccaaaaagttggacagttaa                  <  1:1885329‑M1/51‑31 (MQ=255)

                              cTCGACATGGTGGCCCCACTCtgtactgaccccaaaaagttggacagttaa                  <  1:3468992‑M1/51‑31 (MQ=255)

                               tCGACATGGTGGCCCCACTCtgtactgaccccaaaaagttggacagttaaa                 >  1:4730869‑M1/1‑20 (MQ=255)

                               tCGACATGGTGGCCCCACTCtgtactgaccccaaaaagttggacagttaaa                 <  1:5261948‑M1/51‑32 (MQ=255)

                                cGACATGGTGGCCCCACTCtgtactgaccccaaaaagttggacagttaaac                <  1:6697409‑M1/51‑33 (MQ=255)

                                cGACATGGTGGCCCCACTCtgtactgaccccaaaaagttggacagttaaac                <  1:7227021‑M1/51‑33 (MQ=255)

                                cGACATGGTGGCCCCACTCtgtactgaccccaaaaagttggacagttaaa                 >  1:3813920‑M1/1‑19 (MQ=255)

                                 gACATGGTGGCCCCACTCtgtactgaccccaaaaagttggacagttaaaca               >  1:1084294‑M1/1‑18 (MQ=255)

                                 gACATGGTGGCCCCACTCtgtactgaccccaaaaagttggacagttaaaca               <  1:2296283‑M1/51‑34 (MQ=255)

                                  aCATGGTGGCCCCACTCtgtactgaccccaaaaagttggacagttaaacac              <  1:2852399‑M1/51‑35 (MQ=255)

                                    aTGGTGGCCCCACTCtgtactgaccccaaaaagttggacagttaaacacga            >  1:967155‑M1/1‑15 (MQ=255)

                                    aTGGTGGCCCCACTCtgtactgaccccaaaaagttggacagttaaacacga            <  1:6421881‑M1/51‑37 (MQ=255)

                                      ggtggCCCCACTCtgtactgaccccaaaaagttggacagttaaacacgagg          <  1:5877883‑M1/51‑39 (MQ=255)

                                       gtggCCCCACTCtgtactgaccccaaaaagttggacagttaaacacgaggc         <  1:3720266‑M1/51‑40 (MQ=255)

                                          gCCCCACTCtgtactgaccccaaaaagttggacagttaaacacgaggcata      <  1:3464693‑M1/51‑43 (MQ=255)

                                            cccACTCtgtactgaccccaaaaagttggacagttaaacacgaggcatata    <  1:3611424‑M1/51‑45 (MQ=255)

                                            cccACTCtgtactgaccccaaaaagttggacagttaaacacgaggcatata    >  1:4887339‑M1/1‑7 (MQ=255)

                                             ccACTCtgtactgaccccaaaaagttggacagttaaacacgaggcatatag   >  1:1775377‑M1/1‑6 (MQ=255)

                                              cACTCtgtactgaccccaaaaagttggacagttaaacacgaggcatatagg  <  1:2877587‑M1/51‑47 (MQ=255)

GGAGGCTTCAGCACCAGAGAAAAATCCCCCCTCGACATGGTGGCCCCACTCTTGCGGACGTCCCATCATGACTTTATAACGAACCTGGCCTGGGCGA  >  REL606/2397478‑2397574

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