New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 = 899457281 (1.150)2 (0.010) 2/198 29.6 0.01 coding (550/831 nt) RCAP_rcc00816 extradiol ring‑cleavage dioxygenase subunit B
?NC_014034 3145491 = 270 (1.100)coding (218/1029 nt) cysK2 cysteine synthase
Rejected: Position hash score below cutoff.

AGCCGCGATCTCGCCCATCAGCTGGAGATCGGCCGGACGCTGTCAGAGCTGCGCAACCGCGGCGTGCTGATCCTCGGCTCGGGCAATGTGGTGCACAACC.........................................................................  >  NC_014034/899358‑899457
....................................................................................................TTGCCCGCCGTGCCCTCGACGATGGTGCCGCCGGGTTTCAGCGCACCGCGCGCGATCGCGTCGCGGATGATGT  >  NC_014034/3145491‑3145563
                                                                                                                                                                             
AGCCGCGATCTCGCCCATCAGCTGGAGATCGGCCGGACGCTGTCAGAGCTGCGCAACCGCGGCGTGCTGATCCTCGGCTCGGGCAATGTGGTGCACAACCT                                                                          <  1:2244977/101‑1 (MQ=255)
AGCCGCGATCTCGCCCATCAGCTGGAGATCGGCCGGACGCTGTCAGAGCTGCGCAACCGCGGCGTGCTGATCCTCGGCTCGGGCAATGTGGTGCACAACCT                                                                          <  1:3410484/101‑1 (MQ=255)
AGCCGCGATCTCGCCCATCAGCTGGAGATCGGCCGGACGCTGTCAGAGCTGCGCAACCGCGGCGTGCTGATCCTCGGCTCGGGCAATGTGGTGCACAACCT                                                                          <  1:3550672/101‑1 (MQ=255)
AGCCGCGATCTCGCCCATCAGCTGGAGATCGGCCGGACGCTGTCAGAGCTGCGCAACCGCGGCGTGCTGATCCTCGGCTCGGGCAATGTGGTGCACAACCT                                                                          <  1:384245/101‑1 (MQ=255)
AGCCGCGATCTCGCCCATCAGCTGGAGATCGGCCGGACGCTGTCAGAGCTGCGCAACCGCGGCGTGCTGATCCTCGGCTCGGGCAATGTGGTGCACAACCT                                                                          <  2:3479708/101‑1 (MQ=255)
AGCCGCGATCTCGCCCATCAGCTGGAGATCGGCCGGACGCTGTCAGAGCTGCGCAACCGCGGCGTGCTGATCCTCGGCTCGGGCAATGTGGTGCACAACCT                                                                          <  2:4666805/101‑1 (MQ=255)
AGCCGCGATCTCGCCCATCAGCTGGAGATCGGCCGGACGCTGTCAGAGCTGCGCAACCGCGGCGTGCTGATCCTCGGCTCGGGCAATGTGGTGCACAACCT                                                                          <  2:661059/101‑1 (MQ=255)
                                                          GCGGCGTGCTGATCCTCGGCTCGGGCAATGTGGTGCACAACC.TGCCCGCCGTGCCCTCGACGATGGTGCCGCCGGGTTTCAGCGCACCGCGCGCGATCGCG               >  1:2685081/1‑101 (MQ=255)
                                                                       CCTCGGCTCGGGCAATGTGGTGCAC.ACCCTGCCCGCCGTGCCCTCGACGATGGTGCCGCCGGGTTTCAGCGCACCGCGCGCGATCGCGTCGCGGATGATGT  <  2:2684896/101‑1 (MQ=255)
                                                                                                                                                                             
AGCCGCGATCTCGCCCATCAGCTGGAGATCGGCCGGACGCTGTCAGAGCTGCGCAACCGCGGCGTGCTGATCCTCGGCTCGGGCAATGTGGTGCACAACC.........................................................................  >  NC_014034/899358‑899457
....................................................................................................TTGCCCGCCGTGCCCTCGACGATGGTGCCGCCGGGTTTCAGCGCACCGCGCGCGATCGCGTCGCGGATGATGT  >  NC_014034/3145491‑3145563

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG