New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 16974	0 (0.000)	40 (0.770)	29/60	0.4	100%	intergenic (‑16/‑514)	mokC/nhaA	regulatory protein for HokC, overlaps CDS of hokC/pH‑dependent sodium/proton antiporter
	?	REL606	2774435 =	NA (NA)					noncoding (1443/1443 nt)	IS150	repeat region

GCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTG................................  >  REL606/16937‑16974 ...................................GGTTGTACTGACCCCAAAAAGTTGGACAGTTAAACAC  >  REL606/2774435‑2774468                                                                               CTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGT                                 >  1:5240032/1‑36 (MQ=32)       TACATGTGTTCAGCATATAGGGGGCCTCGGGTTGTA                                <  1:1254564/36‑1 (MQ=35)       TACATGTGTTCAGCATATAGGGGGCCTCGGGTTGTA                                <  1:955183/36‑1 (MQ=35)       TACATGTGTTCAGCATATAGGGGGCCTCGGGTTGTA                                <  1:1966498/36‑1 (MQ=35)       TACATGTGTTCAGCATATAGGGGGCCNCGGGTTGTA                                <  1:2933894/36‑1 (MQ=25)        ACATGTGTTCAGCATATAGGGGGCCTCGGGTTGTAC                               >  1:214584/1‑36 (MQ=38)        ACATGTGTTCAGCATATAGGGGGCCTCGGGTTGTAC                               <  1:4023220/36‑1 (MQ=38)          ATGTGTTCAGCATATAGGGGGCCTCGGGTTGTACTG                             <  1:6563735/36‑1 (MQ=38)          ATGTGTTCAGCATATAGGGGGCCTCGGGTTGTACTG                             >  1:6474236/1‑36 (MQ=38)          ATGTGTTCAGCATATAGGGGGCCTCGGGTTGTACTG                             <  1:7059237/36‑1 (MQ=38)          ATGTGTTCAGCATATAGGGGGCCTCGGGTTGTACTG                             >  1:7261821/1‑36 (MQ=38)            GTGTTCAGCATATAGGGGGCCTCGGGTTGTACTGAC                           <  1:5334658/36‑1 (MQ=38)            GTGTTCAGCATATAGGGGGCCTCGGGTTGTACTGAC                           >  1:6503112/1‑36 (MQ=38)            GTGTTCAGCATATAGGGGGCCTCGGGTTGTAATGAC                           >  1:5304637/1‑36 (MQ=255)             TGTTCAGCATATAGGGGGCCTCGGGTTGTACTGACC                          >  1:253066/1‑36 (MQ=38)             TGTTCAGCATATAGGGGGCCTCGGGTTGTACTGACC                          <  1:1057428/36‑1 (MQ=38)              GTTCAGCATATAGGGGGCCTCGGGTTGTACTGACCC                         <  1:3991510/36‑1 (MQ=38)                  AGCATATAGGGGGCCTCGGGTTGTACTGACCCCAAA                     <  1:4979480/36‑1 (MQ=255)                   GCATATAGGGGGCCTCGGGTTGTACTGACCCCAAAA                    >  1:6088939/1‑36 (MQ=255)                   GCATATAGGGGGCCTCGGGTTGTACTGACCCCAAAA                    >  1:7027813/1‑36 (MQ=255)                    CATATAGGGGGCCTCGGGTTGTACTGACCCCAAAAA                   >  1:944497/1‑36 (MQ=255)                    CATATAGGGGGCCTCGGGTTGTACTGACCCCAAAAA                   >  1:7211107/1‑36 (MQ=255)                     ATATAGGGGGCCTCGGGTTGTACTGACCCCAAAAAG                  >  1:3583971/1‑36 (MQ=255)                       ATAGGGGGCCTCGGGTTGTACTGACCCCAAAAAGTT                >  1:2812223/1‑36 (MQ=255)                       ATAGGGGGCCTCGGGTTGTACTGACCCCAAAAAGTT                >  1:4424410/1‑36 (MQ=255)                       ATAGGGGGCCTCGGGTTGTACTGACCCCAAAAAGTT                <  1:1390360/36‑1 (MQ=255)                        TAGGGGGCCTCGGGTTGTACTGACCCCAAAAAGTTG               <  1:7793104/36‑1 (MQ=255)                        TAGGGGGCCTCGGGTTGTACTGACCCCAAAAAGTTG               >  1:1848548/1‑36 (MQ=255)                         AGGGGGCCTCGGGTTGTACTGACCCCAAAAAGTTGG              <  1:6983335/36‑1 (MQ=255)                           GGGGCCTCGGGTTGTACTGACCCCAAAAAGTTGGAC            <  1:7664972/36‑1 (MQ=255)                            GGGCCTCGGGTTGTACTGACCCCAAAAAGTTGGACA           >  1:7153219/1‑36 (MQ=255)                             GGCCTCGGGTTGTACTGACCCCAAAAAGTTGGACAG          >  1:99498/1‑36 (MQ=255)                             GGCCTCGGGTTGTACTGACCCCAAAAAGTTGGACAG          <  1:1015427/36‑1 (MQ=255)                                CTCGGGTTGTACTGACCCCAAAAAGTTGGACAGTTA       <  1:7508173/36‑1 (MQ=255)                                CTCGGGTTGTACTGACCCCAAAAAGTTGGACAGTTA       >  1:7105026/1‑36 (MQ=255)                                 TCGGGTTGTACTGACCCCAAAAAGTTGGACAGTTAA      >  1:7252821/1‑36 (MQ=255)                                 TCGGGTTGTACTGACCCCAAAAAGTTGGACAGTTAA      <  1:4621947/36‑1 (MQ=255)                                 TCGGGTTGTACTGACCCCAAAAAGTTGGACAGTTAA      >  1:4151555/1‑36 (MQ=255)                                 TCGGGTTGTACTGACCCCAAAAAGTTGGACAGTTAA      >  1:1774428/1‑36 (MQ=255)                                   GGGTTGTACTGACCCCAAAAAGTTGGACAGTTAAAC    >  1:7380887/1‑36 (MQ=38)                                                                          GCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTG................................  >  REL606/16937‑16974 ...................................GGTTGTACTGACCCCAAAAAGTTGGACAGTTAAACAC  >  REL606/2774435‑2774468 Base quality scores:  ATCG  < 3 ≤  ATCG  < 19 ≤  ATCG  < 24 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG