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


TCCATTTCATCGTCCAACAAAATGGGTGCAGTACAAC...................................  <  REL606/588529‑588495


.................................CAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCAT  >  REL606/1733647‑1733683


TCCATTTCATCGTCCCACAAAATGGGTGCAGTACAA                                      <  1:1319931/36‑1 (MQ=0)


TCCATTTCATCGTCCAACAAAATGGGTGCAGTACAA                                      >  1:3096993/1‑36 (MQ=1)


TCCATTTCATCGTCCAACAAAATGGGTGCAGTACAA                                      <  1:6720398/36‑1 (MQ=1)


TCCATTTCATCGTCCAACAAAATGGGTGCAGTACAA                                      >  1:1602596/1‑36 (MQ=1)


TCCATTTCATCGTCCAACAAAATGGGTGCAGTACAA                                      >  1:3233510/1‑36 (MQ=1)


 CCATTTCATCGTCCAACAAAATGGGTGCAGTACAAC                                     <  1:5979141/36‑1 (MQ=1)


 CCATTTCATCGTCCAACAAAATGGGTGCAGTACAAC                                     >  1:1965067/1‑36 (MQ=1)


 CCATTTCATCGTCCAACAAAATGGGTGCAGTACAAC                                     <  1:518720/36‑1 (MQ=1)


   ATTTCATCGTCCAACAAAATGGGTGCAGTACAACAA                                   >  1:5813488/1‑36 (MQ=1)


    TTTCATCGTCCAACAAAATGGGTGCAGTACAACAAC                                  >  1:2312313/1‑36 (MQ=1)


      atcTCGTCCAACAAAATGGGTGCAGTACAACAACTT                                <  1:2182853/33‑1 (MQ=0)


      TCATCGTCCAACAAAATGGGTGCAGTACAACAACTT                                <  1:5927855/36‑1 (MQ=1)


      TCATCGTCCAACAAAATGGGTGCAGTACAACAACTT                                <  1:5929671/36‑1 (MQ=1)


      TCATCGTCCAACAAAATGGGTGCAGTACAACAACTT                                >  1:2082809/1‑36 (MQ=1)


      TCATCGTCCAACAAAATGGGTGCAGTACAACAACTT                                >  1:1859173/1‑36 (MQ=1)


       CATCGTCCAACAAAATGGGTGCAGTACAACAACTTC                               >  1:6624203/1‑36 (MQ=1)


       CATCGTCCAACAAAATGGGTGCAGTACAACAACTTC                               <  1:3920212/36‑1 (MQ=1)


        ATCGTCCAACAAAATGGGTGCAGTACAACAACTTCG                              <  1:669229/36‑1 (MQ=1)


        ATCGTCCAACAAAATGGGTGCAGTACAACAACTTCG                              >  1:3648480/1‑36 (MQ=1)


         TCGTCCAACAAAATGGGTGCAGTACAACAACTTCGA                             <  1:703351/36‑1 (MQ=32)


          CGTCCAACAAAATGGGTGTAGTACAACAACTTCGAC                            <  1:428739/36‑1 (MQ=255)


          CGTCCAACAAAATGGGTGCAGTACAACACCTTCGAC                            >  1:6126827/1‑36 (MQ=255)


          CGTCCAACAAAATGGGTGCAGTACAACAACTTAGAt                            >  1:5921877/1‑35 (MQ=0)


           GTCCAACAAAATGGGTGCAGTACAACAACTTCGACG                           >  1:2894656/1‑36 (MQ=38)


           GTCCAACAAAATGGGTGCAGTACAACAACTTCGACG                           <  1:3376699/36‑1 (MQ=38)


           GTCCAACAAAATGGGTGCAGTACAACAACTTCGACG                           >  1:1863032/1‑36 (MQ=38)


            TCCAACAAAATGGGTGCAGTACAACAACTTCGACGA                          <  1:237377/36‑1 (MQ=40)


            TCCAACAAAATGGGTGCAGTACAACAACTTCGACGA                          >  1:5188340/1‑36 (MQ=40)


             CCAACAAAATGGGTGCAGTACAACAACTTCGACGAA                         >  1:4388468/1‑36 (MQ=255)


             CCAACAAAATGGGTGCAGTACAACAACTTCGACGAA                         <  1:3168735/36‑1 (MQ=255)


             CCAACAAAATGGGTGCAGTACAACAACTTCGACGAA                         >  1:3753760/1‑36 (MQ=255)


              CAACAAAATGGGTGCAGTACAACAACTTCGACGAAA                        <  1:2016766/36‑1 (MQ=255)


              CAACAAAATGGGTGCAGTACAACAACTTCGACGAAA                        <  1:1254894/36‑1 (MQ=255)


                 CAAAATGGGTGCAGTACAACAACTTCGACGAAATCC                     >  1:4799721/1‑36 (MQ=255)


                 CAAAATGGGTGCAGTACAACAACTTCGACGAAATCC                     <  1:5156634/36‑1 (MQ=255)


                 CAAAATGGGTGCAGTACAACAACTTCGACGAAATCC                     >  1:6681125/1‑36 (MQ=255)


                 CAAAATGGGTGCAGTACAACAACTTCGACGAAATCC                     <  1:2607963/36‑1 (MQ=255)


                  AAAGTGGGTGCAGTACAACAACTTCGACGAAATCCT                    <  1:4501355/36‑1 (MQ=12)


                  AAAATGGGTGCAGTACAACAACTTCGACGAAATCCT                    >  1:817429/1‑36 (MQ=40)


                   atATGGGTGCAGTACAACAACTTCGACGAAATCCTC                   <  1:1242169/34‑1 (MQ=17)


                   AAATGGGTGCAGTACAACAACTTCGACGAAATCCTC                   <  1:4328438/36‑1 (MQ=38)


                    AATGGGTGCAGTACAACAACTTCGACGAAATCCTCG                  >  1:2078444/1‑36 (MQ=35)


                    AATGGGTGCAGTACAACAACTTCGACGAAATCCTCG                  >  1:4707910/1‑36 (MQ=35)


                      TGGGTGCAGTACAACAACTTCGACGAAATCCTCGca                >  1:2228108/1‑34 (MQ=1)


                      TGGGTGCAGTACAACAACTTCGACGAAATCCTCGAA                <  1:1349842/36‑1 (MQ=1)


                      TGGGTGCAGTACAACAACTTCGACGAAATCCTCGAA                >  1:1951812/1‑36 (MQ=1)


                      TGGGTGCAGTACAACAACTTCGACGAAATCCTCGAA                >  1:6111275/1‑36 (MQ=1)


                      TGGGTGCAGTACAACAACTTCGACGAAATCCTCGAA                >  1:5954341/1‑36 (MQ=1)


                      TGGGTGCAGTACAACAACTTCGACGAAATCCTCGAA                <  1:4284154/36‑1 (MQ=1)


                        ggaGCAGTACAACAACTTCGACGAAATCCTCGAAGC              <  1:6210779/33‑1 (MQ=0)


                        GGTGCAGTACAACAACTTCGACGAAATCCTCGAAGC              >  1:3097903/1‑36 (MQ=1)


                        GGTGCAGTACAACAACTTCGACGAAATCCTCGAAGC              >  1:4520737/1‑36 (MQ=1)


                         GTGCAGTACAACAACTTCGACGAAATCCTCGAAGCC             >  1:4053987/1‑36 (MQ=1)


                          TGCAGTACAACAACTTCGACGAAATCCTCGAAGCCT            >  1:1373527/1‑36 (MQ=1)


                          TGCAGTACAACAACTTCGACGAAATCCTCGAAGCCT            <  1:1972381/36‑1 (MQ=1)


                            CAGTACAACAACTTCGACGAAATCCTCGAAGCCTCT          >  1:5084052/1‑36 (MQ=1)


                            CAGTACAACAACTTCGACGAAATCCTCGAAGCCTCT          <  1:2805093/36‑1 (MQ=1)


                            CAGTACAACAACTTCGACGAAATCCTCGAAGCCTCT          >  1:4023453/1‑36 (MQ=1)


                            CAGTAAAACAACTTCGACGAAATCCTCGAAGCCTCT          <  1:2151486/36‑1 (MQ=0)


                             AGTACAACAACTTCGACGAAATCCTCGCAGCCTCTG         >  1:381863/1‑36 (MQ=0)


                             AGTACAACAACTTCGACGAAATCCTCGAAGCCTCTG         <  1:5598230/36‑1 (MQ=1)


                              gaACAACAACTTCGACGAAATCCTCGAAGCCTCTGA        <  1:4296085/34‑1 (MQ=1)


                              GTACAACAACTTCGACGAAATCCTCGAAGCCTCTGA        >  1:3396807/1‑36 (MQ=1)


                              GTACAACAACTTCGACGAAATCCTCGAAGCCTCTGA        >  1:6122840/1‑36 (MQ=1)


                               taaAACAACTTCGACGAAATCCTCGAAGCCTCTGAC       <  1:3893712/33‑1 (MQ=0)


                               aACAACAACTTCGACGAAATCCTCGAAGCCTCTGAC       <  1:6419369/35‑1 (MQ=1)


                                 CAACATCTTCGACGAAATCCTCGAAGCCTCTGACGG     <  1:3230652/36‑1 (MQ=0)


                                  AACAACTTCGACGAAATCCTCGAAGCCTCTGACGGC    <  1:2591329/36‑1 (MQ=1)


                                  AACAACTTCGACGAAATCCTCGAAGCCTCTGACGGC    >  1:3834598/1‑36 (MQ=1)


                                   ACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCA   >  1:4455150/1‑36 (MQ=1)


                                    CAACTTCGACGAAATCCTCGAAGCCTCTGACGGCAT  <  1:324202/36‑1 (MQ=255)


TCCATTTCATCGTCCAACAAAATGGGTGCAGTACAAC...................................  <  REL606/588529‑588495


.................................CAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCAT  >  REL606/1733647‑1733683

Base quality scores:


ATCG

< 3 ≤


ATCG

< 11 ≤


ATCG

< 22 ≤


ATCG

< 30 ≤


ATCG

< 32 ≤


ATCG