BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	REL606	1,733,647	IS150 (–) +3 bp	coding (683‑685/1413 nt)	pykF	pyruvate kinase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	annotation	gene	product
*	?	REL606	588495 =	NA (NA)	57 (1.170)	35/62	0.1	noncoding (1/1977 nt)	IS150	repeat region
*	?	REL606	1733647 =	37 (0.680)	57 (1.170)	35/62	0.1	coding (683/1413 nt)	pykF	pyruvate kinase
*	?	REL606	= 666130	NA (NA)	37 (0.690)	29/68	0.5	noncoding (1443/1443 nt)	IS150	repeat region
*	?	REL606	= 1733649	65 (1.190)	37 (0.690)	29/68	0.5	coding (685/1413 nt)	pykF	pyruvate kinase


CGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAT...................................  >  REL606/666096‑666130


...................................TGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGATGA  <  REL606/1733649‑1733614


ctTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAT                                     <  1:2474577/34‑1 (MQ=1)


CGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAT                                     <  1:4372063/36‑1 (MQ=1)


CGTGTTTAACTGTCCAACTTTTTGGGGTCAGCACAT                                     >  1:5004036/1‑36 (MQ=0)


CGTGTATAACTGTCCAACTTTTTGGGGTCAGTACAT                                     <  1:6587015/36‑1 (MQ=0)


 GTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATa                                    <  1:2026751/36‑2 (MQ=1)


 GTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATa                                    <  1:3440988/36‑2 (MQ=1)


 GTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATG                                    <  1:7054208/36‑1 (MQ=32)


  TGTTTAACTGTCCAACTTTTTGGGGTCAGTACATac                                   <  1:4330236/36‑3 (MQ=1)


  TGTTTAACTGTCCAACTTTTTGGGGTCAGTACATac                                   >  1:6331722/1‑34 (MQ=1)


  TGTTTAACTGTCCAACTTTTTGGGGTCAGTACATac                                   >  1:635345/1‑34 (MQ=1)


  TGTTTAACTGTCCAACTTTTTGGGGTCAGTACATGT                                   <  1:1265789/36‑1 (MQ=35)


   GTTTAACTGTCCAACTTTTTGGGGTCAGTACATacg                                  <  1:2886302/36‑4 (MQ=1)


   GTTTAACTGTCCAACTTTTTGGGGTCAGTACATGTT                                  <  1:2608350/36‑1 (MQ=38)


   GTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTT                                  <  1:6725608/36‑1 (MQ=38)


   GTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTT                                  >  1:4421246/1‑36 (MQ=38)


    TTTAACTGTCCAACTTTTTGGGGTCAGTAGATGTTG                                 >  1:6501866/1‑36 (MQ=255)


    TTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTTG                                 <  1:40355/36‑1 (MQ=38)


    TTTAACTGTCCAACTTTTTGCGGTCAGTACATGTTG                                 >  1:5997928/1‑36 (MQ=255)


     TTAACTGTCCAACTTTTTGGGGTCAGTACATGTTGA                                >  1:6493408/1‑36 (MQ=40)


     TTAACTGTCCAACTTTTTGGGGTCAGTACATGTTGA                                <  1:6260572/36‑1 (MQ=40)


      TAACTGTCCAACTTTTTGGGGTCAGTACATGTTGAG                               <  1:2387214/36‑1 (MQ=255)


       AACTGTCCAACTTTTTGGGGTCAGTACATGTTGAGG                              <  1:5786085/36‑1 (MQ=255)


       AACTGTCCAACTTTTTGGGGTCAGTACATGTTGAGG                              >  1:4915720/1‑36 (MQ=255)


            TCCAACTTTTTGGGGTCAGTACATGTTGAGGCCTTC                         <  1:3644252/36‑1 (MQ=255)


             CCAACTTTTTGGGGTCAGTACATGTTGAGGCCTTCC                        >  1:651316/1‑36 (MQ=255)


              CAACTTTTTGGGGTCAGTACATGTTGAGGCCTTCCT                       >  1:6960958/1‑36 (MQ=255)


              CAACTTTTTGGGGTCAGTACATGTTGAGGCCTTCCT                       <  1:3398128/36‑1 (MQ=255)


              CAACTTTTTGGGGTCAGTACATGTTGAGGCCCTCCT                       >  1:589780/1‑36 (MQ=255)


               AACTTTTTGGGGTCAGTACATGTTGAGGCCTTCCTG                      >  1:3356538/1‑36 (MQ=255)


                 CTTTTTGGGGTCAGTACATGTTGAGGCCTTCCTGGT                    <  1:4778567/36‑1 (MQ=255)


                  TTTTTGGGGTCAGTACATGTTGAGGCCTTCCTGGTT                   >  1:2471673/1‑36 (MQ=255)


                    TTTGGGGTCAGTACATGTTGAGGCCTTCCTGGTTTT                 <  1:5499920/36‑1 (MQ=255)


                    TTTGGGGTCAGTACATGTTGAGGCCTTCCTGGTTTT                 >  1:4062409/1‑36 (MQ=255)


                    TTTGGGGTCAGTACATGTTGAGGCCTTCCTGGTTTT                 >  1:3726680/1‑36 (MQ=255)


                     TTGGGGTCAGTACATGTTGAGGCCTTCCTGGTTTTC                <  1:7079733/36‑1 (MQ=255)


                       GGGGTCAGTACATGTTGAGGCCTTCCTGGTTTTCGA              >  1:2723664/1‑36 (MQ=255)


                        GGGTCAGTACATGTTGAGGCCTTCCTGGTTTTCGAT             <  1:5746124/36‑1 (MQ=255)


                        GGGTCAGTACATGTTGAGGCCTTCCTGGTTTTCGAT             >  1:691462/1‑36 (MQ=255)


                         GGTCAGTACATGTTGAGGCCTTCCTGGTTTTCGATT            <  1:4168245/36‑1 (MQ=255)


                         GGTCAGTACATGTTGAGGCCTTCCTGGTTTTCGATT            <  1:1371748/36‑1 (MQ=255)


                           TCAGTACATGTTGAGGCCTTCCTGGTTTTCGATTTT          >  1:6480009/1‑36 (MQ=255)


                             AGTACATGTTGAGGCCTTCCTGGTTTTCGATTTTGG        >  1:5974115/1‑36 (MQ=255)


                             AGTACATGTTGAGGCCTTCCTGGTTTTCGATTTTGG        <  1:4727219/36‑1 (MQ=255)


                             AGTACATGTTGAGGCCTTCCTGGTTTTCGATCTTGG        >  1:1588985/1‑36 (MQ=255)


                              GTACATGTTGAGGCCTTCCTGGTTTTCGATTTTGGA       <  1:802586/36‑1 (MQ=255)


                               TACATGTTGAGGCCTTCCTGGTTTTCGATTTTGGAG      <  1:2616775/36‑1 (MQ=255)


                                 CATGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGAT    <  1:6789813/36‑1 (MQ=255)


                                   TGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGATGA  <  1:7238822/36‑1 (MQ=255)


                                   TGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGATGA  >  1:7416103/1‑36 (MQ=255)


CGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAT...................................  >  REL606/666096‑666130


...................................TGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGATGA  <  REL606/1733649‑1733614

Base quality scores:


ATCG

< 3 ≤


ATCG

< 13 ≤


ATCG

< 22 ≤


ATCG

< 30 ≤


ATCG

< 32 ≤


ATCG