New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 374957240 (0.880)8 (0.180) 4/58 5.2 16.5% intergenic (‑65/+33) mutM/rpmG formamidopyrimidine‑DNA glycosylase/50S ribosomal protein L33
?REL606 = 4157631 41 (0.740)coding (278/546 nt) nusG transcription antitermination protein NusG
Rejected: Position hash score below cutoff.

AGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCG...................................  >  REL606/3749533‑3749572
..................................CCGCCGATGAAGCCCATCACACGCGGTACGCTGCGCACCAG  <  REL606/4157631‑4157597
                                                                           
                       AACAAAAACCCCCGCCGATGAAGCCCATCACACGCG                  <  1:685847/36‑1 (MQ=0)
                        ACAGAAAACCCCGCCGATGAAGCCCATCACACGCGG                 <  1:2709149/36‑1 (MQ=0)
                              acCCCCGCCGATGAAGCCCATCACACGCGGTACGCT           <  2:1925014/34‑1 (MQ=1)
                              acCCCCGCCGATGAAGCCCATCACACGCGGTACGCT           <  2:2316896/34‑1 (MQ=1)
                               ctCCCGCCGATGAAGCCCATCACACGCGGTACGCTG          <  1:2425187/34‑1 (MQ=1)
                                CCCCGCCGATGAAGCCCATCACACGCGGTACGCTGC         <  2:2501921/36‑1 (MQ=1)
                                CCCCGCCGATGAAGCCCATCACACGCGGTACGCTGC         <  2:3145682/36‑1 (MQ=1)
                                 CCCGCCGATGAAGCCCATCACACGCGGTACGCTGCG        <  1:241148/36‑1 (MQ=1)
                                                                           
AGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCG...................................  >  REL606/3749533‑3749572
..................................CCGCCGATGAAGCCCATCACACGCGGTACGCTGCGCACCAG  <  REL606/4157631‑4157597

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 12 ≤  ATCG  < 27 ≤  ATCG  < 32 ≤  ATCG