New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 381306239 (0.850)8 (0.180) 4/58 5.2 16.7% intergenic (‑143/‑791) ivbL/emrD ilvB operon leader peptide/multidrug efflux system protein
?REL606 = 4157631 41 (0.740)coding (278/546 nt) nusG transcription antitermination protein NusG
Rejected: Position hash score below cutoff.

CACGTTTTGTACTGAATTGCAGATAACAAAAAACCCCGCCG...................................  >  REL606/3813022‑3813062
...................................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)
                                                                            
CACGTTTTGTACTGAATTGCAGATAACAAAAAACCCCGCCG...................................  >  REL606/3813022‑3813062
...................................CCGCCGATGAAGCCCATCACACGCGGTACGCTGCGCACCAG  <  REL606/4157631‑4157597

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