New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 592057NA (NA)5 (0.100) 5/60 NT 11.4% noncoding (1343/1343 nt) IS186 repeat region
?REL606 1391825 = 39 (0.640)coding (230/1614 nt) mppA murein tripeptide (L‑ala‑gamma‑D‑glutamyl‑meso‑DAP) transporter subunit

TCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAG...................................  >  REL606/592023‑592057
................................GGGAGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAA  >  REL606/1391825‑1391861
                                                                        
          CAACCCTTAAGTTAGCGCTTATGGGAGATTGTCCCC                            <  1:6483187/36‑1 (MQ=1)
              CCTTAAGTTAGCGCTTATGGGAGATTGTCCCCGGCG                        <  1:2515403/36‑1 (MQ=1)
                             TATGGGAGATTGTCCCCGGCGTTGCGACTCAGTGGA         <  1:6746571/36‑1 (MQ=1)
                              ATGGGAGATTGTCCCCGGCGTTGCGACTCAGTGGAA        <  1:53044/36‑1 (MQ=1)
                               TGGGAGATTGTCCCCGGCGTTGCGACTCAGTGGAAA       <  1:5299704/36‑1 (MQ=1)
                                                                        
TCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAG...................................  >  REL606/592023‑592057
................................GGGAGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAA  >  REL606/1391825‑1391861

Base quality scores:  ATCG  < 3 ≤  ATCG  < 27 ≤  ATCG  < 30 ≤  ATCG  < 35 ≤  ATCG  < 37 ≤  ATCG