New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1340206 =27 (0.470)6 (0.100) 3/70 NT 15.4% intergenic (+107/‑86) yciM/pyrF hypothetical protein/orotidine 5'‑phosphate decarboxylase
?REL606 3170213 = 39 (0.680)coding (104/1005 nt) ygjR predicted NAD(P)‑binding dehydrogenase

CGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGC  >  REL606/3170178‑3170248
                                   |                                   
cgAGGCCGCCCATGAGAGCGGTAAATACAAGTTAAc                                     <  1:2722482/36‑1 (MQ=255)
 gAGGCCGCCCATGAGAGCGGTAAATACAAGTTAAcc                                    <  1:3710607/36‑1 (MQ=255)
 gAGGCCGCCCATGAGAGCGGTAAATACAAGTTAAcc                                    >  1:3084235/1‑36 (MQ=255)
 gAGGCCGCCCATGAGAGCGGTAAATACAAGTTAAcc                                    <  1:2087541/36‑1 (MQ=255)
 gAGGCCGCCCATGAGAGCGGTAAATACAAGTTAAcc                                    >  1:5848955/1‑36 (MQ=255)
 gAGGCCGCCCATGAGAGCGGTAAATACAAGTTAAcc                                    <  1:1282887/36‑1 (MQ=255)
 gAGGCCGCCCATGAGAGCGGTAAATACAAGTTAAcc                                    >  1:1651855/1‑36 (MQ=255)
   ggCCGCCCATGAGAGCGGTAAATACAAGTTAAccgc                                  <  1:2592336/36‑1 (MQ=255)
   ggCCGCCCATGAGAGCGGTAAATACAAGTTAAccgc                                  >  1:302368/1‑36 (MQ=255)
    gccgccCATGAGAGCGGTAAATACAAGTTAAccgcc                                 >  1:6311712/1‑36 (MQ=255)
    gccgccCATGAGAGCGGTAAATACAAGTTAAccgcc                                 >  1:6799625/1‑36 (MQ=255)
    gccgccCAGGAGAGCGGTAAATACAAGTTAAccgcc                                 <  1:3083712/36‑1 (MQ=255)
       gccCATGAGAGCGGTAAATACAAGTTAACCGCCGta                              >  1:510168/1‑36 (MQ=255)
       gccCATGAGAGCGGTAAATACAAGTTAACCGCCGta                              >  1:1581885/1‑36 (MQ=255)
          cATGAGAGCGGTAAATACAAGTTAACCGCCGTatat                           >  1:3858012/1‑36 (MQ=255)
          aaggagaGCGATAAAAAGAAATAACCCGCCGTatat                           <  1:6917194/33‑1 (MQ=2)
          aaagagaGAGGAAAAAAAAAGTAAACCGCCGTatat                           <  1:413978/33‑1 (MQ=9)
            tGAGAGCGGTAAATACAAGTTAACCGCCGTATATTc                         >  1:6174072/1‑36 (MQ=255)
              agagCGGTAAATACAAGTTAACCGCCGTATATTccc                       >  1:488643/1‑36 (MQ=255)
               gagCGGTAAATACAAGTTAACCGCTGTATATTCCCg                      >  1:2876422/1‑36 (MQ=255)
                agCGGTAAATACAAGTTAACCGCCGTATATTCCCGc                     >  1:5167618/1‑36 (MQ=255)
                  ggaaaaaaaaaattaacCCGCCGTATATTCCCGCAg                   <  1:3623525‑M2/19‑1 (MQ=255)
                   ggTAAATACAAGTTAACCGCCGTATATTCCCGCAGc                  >  1:219492/1‑36 (MQ=255)
                   ggTAAATACAAGTTAACCGCCGTATATTCCCGCAGc                  >  1:854727/1‑36 (MQ=255)
                   ggTAAATACAAGTTAACCGCCGTATANTCCCGCAGc                  >  1:4545563/1‑36 (MQ=255)
                   ggTAAAGACAAGTNAACCGACGTATATTCCCGCAGc                  >  1:2666450/1‑36 (MQ=37)
                   gaaaaaaaaaattaacCCGCCGTATATTCCCGCAGc                  <  1:3138902‑M2/20‑1 (MQ=255)
                     tAAATACAAGTTAACCGCCGTATATTCCCGCAGCCt                >  1:3512548/1‑36 (MQ=255)
                     aaaaTAAAATTAACCCGCCGTATATTCCCGCAGCCt                <  1:5496035/35‑1 (MQ=16)
                        aTACAAGTTAACCGCCGTATATTCCCGCAGCCTTGa             >  1:6521066/1‑36 (MQ=255)
                        aTACAAGTTAACCGCCGTATATTCCCGCAGCCTTGa             >  1:1983696/1‑36 (MQ=255)
                        aTACAAGTTAACCGCCGTATATTCCCGCAGCCTTGa             >  1:1688598/1‑36 (MQ=255)
                         tACAAGTTAACCGCCGTATATTCCCGCAGCCTTGaa            >  1:4185868/1‑36 (MQ=255)
                          aaaattaaaCCGCCGTATATTCCCGCAGCCTTGAAc           <  1:6781589‑M2/27‑1 (MQ=255)
                          aCAAGTTAACCGCCGTATATTCCCTCAGCCTTGAAc           >  1:3915522/1‑36 (MQ=255)
                          aCAAGTTAACCGCCGTATATTCCCGCAGCCTTGAAc           >  1:5732717/1‑36 (MQ=255)
                          aCAAGTTAACCGCCGTATACTCCCGCAGCCGTGAAc           >  1:31411/1‑36 (MQ=38)
                           cAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACa          >  1:5861964/1‑36 (MQ=255)
                            aaGTTAACCGCCGTATATTCCCGCAGCCTTGAACAg         >  1:5110249/1‑36 (MQ=255)
                            aaGTTAACCGCCGTATATTCCCGCAGCCTTGAACAg         >  1:1143583/1‑36 (MQ=255)
                              ttaaaCCGCCGTATATTCCCGCAGCCTTGAACAGGc       <  1:158885‑M2/31‑1 (MQ=255)
                               taccCCGCCGTATATTCCCGCAGCCTTGAACAGGcc      <  1:6006724‑M2/32‑1 (MQ=255)
                                 gcCCGCCGTATATTCCCGCAGCCTTGAACAGGCCCa    <  1:3028033‑M2/34‑1 (MQ=255)
                                  aCCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAg   >  1:6607745/1‑36 (MQ=255)
                                  aCCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAg   >  1:7079792/1‑36 (MQ=255)
                                   ccgccgTATATTCCCGCAGCCTTGAACAGGCCCAGc  >  1:5637673/1‑36 (MQ=255)
                                   ccgccgTATATTCCCGCAGCCTTGAACAGGCCCAGc  <  1:6871034/36‑1 (MQ=255)
                                   |                                   
CGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGC  >  REL606/3170178‑3170248

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG