New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 72407 =27 (0.610)4 (0.090) 4/56 5.1 9.1% coding (446/1701 nt) araB ribulokinase
?REL606 3650765 = 53 (0.960)intergenic (+46/+141) cspA/hokA major cold shock protein/toxic polypeptide, small
Rejected: Position hash score below cutoff.

GATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGC  >  REL606/72373‑72442
                                  |                                   
gATTTTTGCCCAGAACCATTCGCTGGAATAAATAcc                                    >  1:1315257/1‑36 (MQ=255)
gATATTTCCCCAGAACCATTCGCTGGAATAAATAcc                                    <  2:922772/36‑1 (MQ=37)
 aTTTTTGACCAGAACCATTCGCTGGAATAAATACCg                                   <  1:2825887/36‑1 (MQ=255)
  tttttGCCCAGAACCATTCGCTGGAATAAATACCGc                                  >  2:1977564/1‑36 (MQ=255)
  tttttGCCCAGAACCATTCGCTGGAATAAATACCGc                                  <  1:2430167/36‑1 (MQ=255)
   ttttGCCCAGAACCATTCGCTGGAATAAATACCGcc                                 <  1:1095852/36‑1 (MQ=255)
     ttGCCCAAAACCATTCGCTGGAATAAATACCGCCaa                               <  2:2983616/36‑1 (MQ=255)
       gCCCAGAACCATTCGATGGAATAAATACCGCCAata                             >  1:949519/1‑36 (MQ=255)
          cAGAACCATTCGCTGGAATAAATACCGCCAATATAg                          >  2:3274435/1‑36 (MQ=255)
          cAGAACCATTCGCTGGAATAAATACCGCCAATATAg                          >  1:2873232/1‑36 (MQ=255)
          aaaaaCCATTCGCAGGAAAAAACCCCGCCAATATAg                          <  2:1628282/33‑1 (MQ=16)
          aaaaaCAATCCGCGGAAAAAAATCCCGCCAATATAg                          <  2:2173324/33‑1 (MQ=2)
           aGAACCATTCGCTGGAATAAATACCGCCAATATAGc                         >  1:3110073/1‑36 (MQ=255)
           aGAACCATTCGCTGGAATAAATACCGCCAATATAGc                         >  1:3374924/1‑36 (MQ=255)
            gAACCATTCGCTGGAATAAATACCGCCAATATAGCg                        >  2:1094817/1‑36 (MQ=255)
            gAACCATTCGCTGGAATAAATACCGCCAATATAGCg                        >  1:2558788/1‑36 (MQ=255)
              aCCATTCGCTGGCATAAATACCGCCAATATCGCggg                      >  2:2776919/1‑36 (MQ=25)
              aCCATTCGCTGGAATAAATACCGCCAATATAGCggg                      >  1:3355617/1‑36 (MQ=255)
               ccATTCGCTGGAATAAATACCGCCAATATAGCGGGa                     >  2:1392790/1‑36 (MQ=255)
                 aTTCGCTGGAATAAATACCGCCAATATAGCGGGggt                   >  2:444240/1‑36 (MQ=255)
                 aTTCGCTGGAATAAATACCGCCAATATAGCGGGagt                   >  2:2672852/1‑36 (MQ=255)
                   tCGCTGGAATAAATACCGCCAATATAGCGGGagtag                 >  2:2725761/1‑36 (MQ=255)
                    cGCTGGAATAAATACTGCCAATATAGCGGGggtagt                >  2:2331602/1‑36 (MQ=37)
                           ataaataCCGCCAATATAGCGGGAGTAGTTAACAtt         >  1:1338343/1‑36 (MQ=255)
                           ataaaccCCGCCAATATAGCGGGAGTAGTCAACAtt         <  1:1311005/36‑1 (MQ=37)
                                taCCGCCAATATAGCGGGAGTAGTCAACATTGcctg    >  1:2881839/1‑34 (MQ=255)
                                taCCGCCAATATAGCGGGAGTAGTCAACATTGCCCg    >  1:3216195/1‑36 (MQ=255)
                                ccccGCCAATATAGCGGGAGTAGTCAACATTGCCCg    <  1:444221/34‑1 (MQ=255)
                                  ccGCCAATATAGCGGGAGTAGTCAACATTGCCCGgc  >  2:103482/1‑36 (MQ=255)
                                  ccGCCAATATAGCGGGAGTAGTCAACATTGCCCGgc  <  2:811690/36‑1 (MQ=255)
                                  |                                   
GATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGC  >  REL606/72373‑72442

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 16 ≤  ATCG  < 29 ≤  ATCG  < 32 ≤  ATCG