BRESEQ :: Evidence

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	annotation	gene	product
*	?	REL606	= 339475	52 (1.510)	5 (0.150)	4/44	4.0	coding (592/1083 nt)	lacI	lac repressor
*	?	REL606	= 2022999	50 (0.910)	5 (0.150)	4/44	4.0	coding (600/1071 nt)	hisC	histidinol‑phosphate aminotransferase
Rejected: Position hash score below cutoff.


TTGATTGCGAGTGAGATATTTATGCCAGCCAGCCAGACGCAGACGCGCCGAGACAGAACTTAATGGGCCCG  >  REL606/339440‑339510


ttGATTGCGAGTGAGATATTTATGccagccagccag                                     >  1:5224842/1‑36 (MQ=255)


 tGATTGCGAGTGAGATATTTATGCCAGCTAGCcaga                                    >  1:4349453/1‑36 (MQ=255)


 tGATTGCGAGTGAGATATTTATGCCAGCCAGCcaga                                    >  1:4166280/1‑36 (MQ=255)


 tGATTGCGAGTGAGATATTTATGACAGCCAGCcaga                                    >  1:1346047/1‑36 (MQ=255)


  gATTGCGAGTGAGATATTTATGCCAGCCAGCcagac                                   <  1:4906291/36‑1 (MQ=255)


   aTTGCTAGTAAGATATTTATGCCAGCCAGCcagacg                                  <  1:2878502/36‑1 (MQ=25)


   aTTGCGAGTGAGATATTTATGCCAGCCAGCcagacg                                  >  1:3000090/1‑36 (MQ=255)


    ttGCGAGTGAGATATTTATGCCAGCCAGCcagacgc                                 >  1:3817585/1‑36 (MQ=255)


     tGCGAGTGAGATATTTATGCCAGCCAGCcagacgca                                <  1:831420/36‑1 (MQ=255)


      gCGAGTGAGATATTTATGCCAGCCAGCcagacgcag                               <  1:7309260/36‑1 (MQ=255)


        gAGTGAGATATTTATGCCAGCCAGCcagacgcagac                             <  1:2416018/36‑1 (MQ=255)


         aGTGAGATATTTATGCCAGCCAGCcagacacagacg                            >  1:1893361/1‑36 (MQ=255)


         aGTGAGATATTTATGCCAGCCAGCcacacgcagacg                            >  1:1521392/1‑36 (MQ=255)


          gTGAAATATTTATGCCAGCCAGCCagacgcagacgc                           <  1:7452457/36‑1 (MQ=255)


           tGAGATATTTATGCCAGCCAGCCAGACGCAGAcgcg                          >  1:3270343/1‑36 (MQ=255)


             agaTGTTTATGCCAGCCAGCCAGACGCAGACGCGcc                        <  1:1631573/36‑1 (MQ=255)


             agaTATTTATGCCAGCCAGCCAGCCGCAGACGCGcc                        >  1:3417803/1‑36 (MQ=255)


             agaTATTTATGCCAGCCAGCCAGACGCAGACGCGcc                        >  1:2372048/1‑36 (MQ=255)


             agaTATTTATGCCAGCCAGCCAGACGCAGACGCGcc                        <  1:246630/36‑1 (MQ=255)


             agaTATTTATGCCAGCCAGCCAGACGCAGACGCGcc                        >  1:5581393/1‑36 (MQ=255)


             agaTATTTATGCCAGCCAGCCAGACGCAGACGCGcc                        >  1:2924686/1‑36 (MQ=255)


              gctatTTATGCCATCCAGCCAGACGCACCCGCGCCg                       >  1:1668367/3‑36 (MQ=18)


              gaTATTTATGCCAGCCAGCCAGACGCAGACGCGCCg                       >  1:1657251/1‑36 (MQ=255)


              gaTATTTATGCCAGCCAGCCAGACGCAGACGCGCCg                       >  1:1073584/1‑36 (MQ=255)


                tatTTATGCCAGACACCAAGACGCAGACGCGCCgaa                     >  1:3781998/1‑35 (MQ=18)


                 atTTATGCCAGCCAGCCAGACGCAGACGCGCCgaga                    <  1:949227/36‑1 (MQ=255)


                   ttATGCCAGCCAGCCAGACGCAGACGCGCCGAGACa                  <  1:4486592/36‑1 (MQ=255)


                    tATGCCAGCCAGCCAGACGCAGACGCGCCGAGACAg                 >  1:7194864/1‑36 (MQ=255)


                    tATGCCAGCCAGCCAGACGCAGACGCGCCGAGACAg                 >  1:1450263/1‑36 (MQ=255)


                     aTGCCAGCCAGCCAGACGCAGACGCGCCTTGACAGa                >  1:1370891/1‑36 (MQ=25)


                     aTGCCAGCCAGCCAGACGCAGACGCGCCGAGACAGa                >  1:501476/1‑36 (MQ=255)


                        ccagccagccagACGCAGACGCGCCGAGCCAGAACt             >  1:2928062/1‑36 (MQ=255)


                        ccagccagccagACGCAGACGCGCCGAGACAGAACt             >  1:120683/1‑36 (MQ=255)


                        caggccagccagTCGCAGACGCGCCGAGACAGAACt             <  1:3392161/33‑1 (MQ=25)


                         cagccagccagACGCAGACGCGCCGAGACAGAACtt            >  1:7304737/1‑36 (MQ=255)


                         cagccagccagACGCAGACGCGCCGAGACAGAACtt            >  1:4946347/1‑36 (MQ=255)


                          agccagccagACGCAGACGCGCCGAGACAGAACTTa           >  1:1191572/1‑36 (MQ=255)


                          agccagccagACGCAGACGCGCCGAGACAGAACTTa           >  1:1821134/1‑36 (MQ=255)


                          agccagccagACGCAGACGCGCCGAGACAGAACTTa           >  1:3636796/1‑36 (MQ=255)


                          agccagcaagACGCCGACGCGCCGAGACAGTACCTa           >  1:6581604/1‑36 (MQ=18)


                           gccagccagACGCAGACGCGCCGAGACAGAACTTaa          <  1:5904347/36‑1 (MQ=255)


                           gccagccagACGCAGACGCGCCGAGACAGAACTTaa          >  1:3693055/1‑36 (MQ=255)


                            ccagccggACGCAGACGCGCCGAGACAGAACTTAAt         <  1:1626414/36‑1 (MQ=255)


                              agccagACGCAGACGCGCCGAGACAGAACTTAATgg       >  1:4566813/1‑36 (MQ=255)


                              agccagACGCAGACGCGCCGAGACAGAACTTAATgg       >  1:2502391/1‑36 (MQ=255)


                               gccagACGCAGACGCGCGGAGACAGAACTTAATggg      <  1:5093142/36‑1 (MQ=255)


                               gccagACGCAGACGCGCCGAGACAGAACTTAATggg      >  1:5573579/1‑36 (MQ=255)


                                ccagACGCAGACGCGCTGAGACAGAACTTAATGGGc     >  1:1386824/1‑36 (MQ=255)


                                 cagacgcagacgCGCCGAGACAGAACTTAATGGGcc    >  1:17157/1‑36 (MQ=255)


                                 cagacgcagacgCGCCGAGACAGAACTTAATGGGcc    <  1:1787464/36‑1 (MQ=255)


                                 cagacgcagacgCGCCGAGACAGAACTTAATGGGcc    <  1:2060529/36‑1 (MQ=255)


                                   gacgcagacgcGCCGAGACAGAACTTAATGGGCCCg  >  1:4962414/1‑36 (MQ=255)


                                   gacgcagacgcGCCGAGACAGAACTTAATGGGCCCg  <  1:4349366/36‑1 (MQ=255)


TTGATTGCGAGTGAGATATTTATGCCAGCCAGCCAGACGCAGACGCGCCGAGACAGAACTTAATGGGCCCG  >  REL606/339440‑339510

Base quality scores:


ATCG

< 3 ≤


ATCG

< 8 ≤


ATCG

< 16 ≤


ATCG

< 27 ≤


ATCG

< 32 ≤


ATCG