New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 370838 =31 (0.550)4 (0.070) 3/64 NT 11.1% intergenic (+30/‑89) phoA/psiF bacterial alkaline phosphatase/hypothetical protein
?REL606 3051517 = 33 (0.530)coding (835/1224 nt) glcF glycolate oxidase iron‑sulfur subunit

AAATAAAACCGCGCCCGGCGGTGAATTTTCGCTGTCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAAT  >  REL606/370803‑370873
                                   |                                   
caaTAAAACCGCGCCCGGCGGTGAATTTTCGCTGTc                                     <  1:1793455/35‑1 (MQ=255)
atatgccgataaagcgcgtcaggtcagtgaactggc                                     >  1:6682266‑M1/36‑36 (MQ=255)
aaaTAAAACCGCGCCCGGCGGTGAATTTTCGCTGTc                                     <  1:4226388/36‑1 (MQ=255)
aaaTAAAACCGCGCCCGGCGGTGAATTTTCGCTGTc                                     <  1:279916/36‑1 (MQ=255)
aaaTAAAACCGCGCCCGGCGGTGAATTTTCGCTGTc                                     <  1:3244785/36‑1 (MQ=255)
 tatgccgataaagcgcgtcaggtcagtgaactggCg                                    >  1:1128895‑M1/35‑36 (MQ=255)
 tatgccgataaagcgcgtcaggtcagtgaactggCg                                    >  1:7808986‑M1/35‑36 (MQ=255)
 tatgccgataaagcgcgtcaggtcagtgaactggCg                                    <  1:2392775‑M1/2‑1 (MQ=255)
  atgccgataaagcgcgtcaggtcagtgaacttgCgg                                   <  1:4504355‑M1/3‑1 (MQ=255)
  aTAAAACCGCGCCCGGCGGTGAATTTTCGCTGTCgg                                   >  1:5611666/1‑36 (MQ=255)
    aaaaCCGCGCCCGGCGGTGAATTTTCGCTTTCgggg                                 >  1:5708465/1‑35 (MQ=255)
    aaaaCCGCGCCCGGCGGCGAATTTTCGCTGTCGggt                                 >  1:6731716/1‑36 (MQ=255)
    aaaaACGCGCCCGGCGGGGACTTTTAGTTGTCgggg                                 >  1:4344222/1‑35 (MQ=9)
     aaaCCGCGCCCGGCGGTGAATTTTCGCTGTCGggtg                                >  1:587035/1‑36 (MQ=255)
     aaaCCGCGCCCGGCGGTGAATTTTCGCTGTCGggtg                                >  1:7014228/1‑36 (MQ=255)
      cgataaagcgcgtcaggtcagtgaactggcgggggg                               <  1:4997079‑M1/7‑4 (MQ=255)
      acccGCGCCCGGCGGTGAATTTTCGCTGTCGggtgg                               <  1:3095552/34‑1 (MQ=255)
       aCCGCGCCCGGCGGTGAATTTTCGCTGTCGggtggt                              <  1:1662131/36‑1 (MQ=255)
       aCCGCGCCCGGCGGCGAATTTTCGCTGTCGggtggt                              >  1:7402176/1‑36 (MQ=255)
        ccGCGCCCGGCGGTGAATTTTCGCTGTCGGGTGGtt                             <  1:4422469/36‑1 (MQ=255)
        ataaagcgcgtcaggtcagtgaactggCGGTTGGtt                             <  1:8073706‑M1/9‑1 (MQ=255)
         cgcgCCCGGCGGTGAATTTTCGCTGTCGGGTGGttt                            <  1:3379044/36‑1 (MQ=255)
         cgcgCCCGGCGGTGAATTTTCGCTGTCGGGTGGttt                            >  1:8136379/1‑36 (MQ=255)
          gcgcCCGGCGGTGAATTTTCGCTGTCGGGTGGtttt                           >  1:7400848/1‑36 (MQ=255)
           aagcgcgtcagggcagtgaactggCGGGTGGttttt                          <  1:660530‑M1/12‑1 (MQ=255)
            gcCCGGCGGTGAATTTTCGCTGTCGGGTGGtttttt                         >  1:6845270/1‑36 (MQ=255)
              cgcgtcaggtcagtgaactggCGGGGGGTTTTTTTg                       <  1:2595051‑M1/15‑1 (MQ=255)
                    gTGAATTTTCGCTGTCGGGTGGTTTTTTTGCTGTTa                 <  1:7534613/36‑1 (MQ=255)
                     tGAATTTTCGCTGTCGGGTGGTTTTTTTGCTGTTAg                >  1:1026385/1‑36 (MQ=255)
                     tGAATTTTCGCTGTCGGGTGGTTTTTTTGCTGTTAg                >  1:3427460/1‑36 (MQ=255)
                     tGAATTTTCGCTGTCGGGTGGTTTTTTTGCTGTTAg                >  1:537266/1‑36 (MQ=255)
                      gAATTTTCGCTGTCGGGTGGTTTTTTTGCTGTTAGc               >  1:2803950/1‑36 (MQ=255)
                      gAATTTTCGCTGTCGGGTGGTTTTTTTGCTGTTAGc               >  1:7596271/1‑36 (MQ=255)
                       aaTTTTCGCTGTCGGGTGGTTTTTTTGTTGTTAGCa              >  1:4845758/1‑36 (MQ=255)
                        aTTTTCGCTGTCTGGTGGTTTTTTTGCTGTTAGCaa             <  1:2568243/36‑1 (MQ=255)
                             cGCTGTCGGGTGGTTTTTTTGCTGTTAGCAACCAGa        <  1:1135280/36‑1 (MQ=255)
                              gCTGTCGGGTGGTTTTTTTGCTGTTAGCAACCAGAc       <  1:692166/36‑1 (MQ=255)
                               cTGTCGGGTGGTTTTTTTGCTGTTAGCAACCAGACt      <  1:6727102/36‑1 (MQ=255)
                               cTGTCGGGTGGTTTTTTTGCTGTTAGCAACCAGACt      <  1:1669755/36‑1 (MQ=255)
                                 gTCGGGTGGTTTTTTTTCTGTTAGCAACCAGCCTTa    <  1:1651946/36‑1 (MQ=37)
                                 gTCGGGTGGTTTTTTTGTTGTTAGCAACCAGACTTa    >  1:6417463/1‑36 (MQ=255)
                                 gTCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTa    >  1:367805/1‑36 (MQ=255)
                                 gTCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTa    >  1:6928335/1‑36 (MQ=255)
                                 gTCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTa    >  1:2797061/1‑36 (MQ=255)
                                 gTCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTa    >  1:2369637/1‑36 (MQ=255)
                                   CGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAAt  >  1:510808‑M1/1‑36 (MQ=255)
                                   |                                   
AAATAAAACCGCGCCCGGCGGTGAATTTTCGCTGTCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAAT  >  REL606/370803‑370873

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 21 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG