New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2541529	110 (0.800)	6 (0.050)	5/216	NT	5.3%	coding (151/876 nt)	cysW	sulfate/thiosulfate ABC transporter permease
	?	NC_000913	= 2541543	109 (0.830)					coding (137/876 nt)	cysW	sulfate/thiosulfate ABC transporter permease

AACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2541378‑2541529 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcatTCAGCAAGGGGCTGATGCCGGTTTTACAGAATCTGGCCGATCCGGACATGCTGCACGCCATCTGGCTGACGGTGATGATCGCGCTGATTGCCGTACCGGTAAACCTGGTGTTCGGCATTCTGCTGGCCTGGCTGGTGACGCGCTTTAA  <  NC_000913/2541543‑2541396                                                                                                                                                                                                                                                                                                              AACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGG                                                                                                                                                         >  4:823106/1‑149  ACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGC                                                                                                                                                        <  3:619648/149‑1    TGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCAT                                                                                                                                                      >  4:650272/1‑149                                                            CAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATTCAGCAAGGGGCTGATGCCGG                                                                                                                                 <  4:455199/114‑1                                                            CAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATTCAGCAAGGGGCTGATGCCGG                                                                                                                                 >  3:455199/1‑114                                                                                         TCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATTCAGCAAGGGGCTGATGCCGG                                                                                                                                 <  4:1050167/85‑1                                                                                         TCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATTCAGCAAGGGGCTGATGCCGG                                                                                                                                 >  3:1050167/1‑85                                                                                               CCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATTCAGCAAGGGGCTGA                                                                                                                                       <  2:1007163/73‑1                                                                                               CCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATTCAGCAAGGGGCTGA                                                                                                                                       >  1:1007163/1‑73                                                                                                                                                    GGCATTCAGCAAGGGGCTGATGCCGGTTTTACAGAATCTGGCCGATCCGGACATGCTGCACGCCATCTGGCTGACGGTGATGATCGCGCTGATTGCCGTACCGGTAAACCTGGTGTTCGGCATTCTGCTGGCCTGGCTGGTGAC           <  4:733745/144‑1                                                                                                                                                    GGCATTCAGCAAGGGGCTGATGCCGGTTTTACAGAATCTGGCCGATCCGGACATGCTGCACGCCATCTGGCTGACGGTGATGATCGCGCTGATTGCCGTACCGGTAAACCTGGTGTTCGGCATTCTGCTGGCCTGGCTGGTGAC           >  3:733745/1‑144                                                                                                                                                    GGCATTCAGCAAGGGGCTGATGCCGGTTTTACAGAATCTGGCCGATCCGGACATGCTGCACGCCATCTGGCTGACGGTGATGATCGCGCTGATTGCCGTACCGGTAAACCTGG                                          <  1:718430/113‑1                                                                                                                                                    GGCATTCAGCAAGGGGCTGATGCCGGTTTTACAGAATCTGGCCGATCCGGACATGCTGCACGCCATCTGGCTGACGGTGATGATCGCGCTGATTGCCGTACCGGTAAACCTGG                                          >  2:718430/1‑113                                                                                                                                                     GCATTCAGCAAGGGGCTGATGCCGGTTTTACAGAATCTGGCCGATCCGGACATGCTGCACGCCATCTGGCTGACGGTGATGATCGCGCTGATTGCCGTACCGGTAAACCTGGTGTTCGGCATTCTGCTGGCCTGGCTGGTGACGCGCT      <  4:552274/148‑1                                                                                                                                                     GCATTCAGCAAGGGGCTGATGCCGGTTTTACAGAATCTGGCCGATCCGGACATGCTGCACGCCATCTGGCTGACGGTGATGATCGCGCTGATTGCCGTACCGGTAAACCTGGTGTTCGGCATTCTGCTGGCCTGGCTGGTGACGCGCT      >  3:552274/1‑148                                                                                                                                                     GCATTCAGCAAGGGGCTGATGCCGGTTTTACAGAATCTGGCCGATCCGGACATGCTGCACGCCATCTGGCTGACGGTGATGATCGCGCTGATTGCCGTACCGGTAAACCTGGTGTTCGG                                   <  4:396871/119‑1                                                                                                                                                     GCATTCAGCAAGGGGCTGATGCCGGTTTTACAGAATCTGGCCGATCCGGACATGCTGCACGCCATCTGGCTGACGGTGATGATCGCGCTGATTGCCGTACCGGTAAACCTGGTGTTCGG                                   >  3:396871/1‑119                                                                                                                                                     GCATTCAGCAAGGGGCTGATGCCGGTTTTACAGAATCTGGCCGATCCGGACATGCTGCACGCCATCTGGCTGACGGTGATGATCGCGCTGATTGCCGTACCGGTAAACCTGGTGTTC                                     <  3:39709/117‑1                                                                                                                                                     GCATTCAGCAAGGGGCTGATGCCGGTTTTACAGAATCTGGCCGATCCGGACATGCTGCACGCCATCTGGCTGACGGTGATGATCGCGCTGATTGCCGTACCGGTAAACCTGGTGTTC                                     >  4:39709/1‑117                                                                                                                                                      CATTCAGCAAGGGGCTGATGCCGGTTTTACAGAATCTGGCCGATCCGGACATGCTGCACGCCATCTGGCTGACGGTGATGATCGCGCTGATTGCCTTACCGGTAAACCTGGTGTTCGGCATTCTGCTGGCCTGGCTGGTGACGCGCTTT    >  4:899605/1‑149                                                                                                                                                       ATTCAGCAAGGGGCTGATGCCGGTTTTACAGAATCTGGCCGATCCGGACATGCTGCACGCCATCTGGCTGACGGTGATGATCGCGCTGATTGCCGTACCGGTAAACCTGGTGTTCGGCATTCTGCTGGCCTGGCTGGTGACGCGC       >  3:555397/1‑145                                                                                                                                                       ATTCAGCAAGGGGCTGATGCCGGTTTTACAGAATCTGGCCGATCCGGACATGCTGCACGCCATCTGGCTGACGGTGATGATCGCGCTGATTGCCGTACCGGTAAACCTGGTGTTCGGCATTCTGCTGGCCTGGCTGGTGAAGCGCTTTA   >  1:1299598/1‑149                                                                                                                                                       ATTCAGCAAGGGGCTGATGCCGGTTTTACAGAATCTGGCCGATCCGGACATGCTGCACGCCATCTGGCTGACGGTGATGATCGCGCTGATTGCCGTACCGGTAAACCTGGTGTTCGGCATTCTGCTGGCCTGGCGGGTGACGCGC       <  4:555397/145‑1                                                                                                                                                       ATTCAGCAAGGGGCTGATGCCGGTTTTACAGAATCTGGCCGATCCGGACATGCTGCACGCCATCTGGCTGACGGTGATGATCGCGCTGATTGCCGTACCGGTAAACCTGGTGTTCGGCATTC                              <  4:1458626/122‑1                                                                                                                                                       ATTCAGCAAGGGGCTGATGCCGGTTTTACAGAATCTGGCCGATCCGGACATGCTGCACGCCATCTGGCTGACGGTGATGATCGCGCTGATTGCCGTACCGGTAAACCTGGTGTTCGGCATTC                              >  3:1458626/1‑122                                                                                                                                                       ATTCAGCAAGGGGCTGATGCCGGTTTTACAGAATCTGGCCGATCCGGACATGCTGCACGCCATCTGGCTGACGGTGATGATCGCGCTGATTGCCGTACC                                                     <  2:644564/99‑1                                                                                                                                                       ATTCAGCAAGGGGCTGATGCCGGTTTTACAGAATCTGGCCGATCCGGACATGCTGCACGCCATCTGGCTGACGGTGATGATCGCGCTGATTGCCGTACC                                                     >  1:644564/1‑99                                                                                                                                                        TTCAGCAAGGGGCTGATGCCGGTTTTACAGAATCTGGCCGATCCGGACATGCTGCACGCCATCTGGCTGACGGTGATGATCGCGCTGATTGCCGTACCGGTAAACCTGGTGTTCGGCATTCTGCTGGCCTGGCTGGTGACGCGCTTTAA  >  4:111762/1‑149                                                                                                                                                        TTCAGCAAGGGGCTGATGCCGGTTTTACAGAATCTGGCCGATCCGGACATGCTGCACGCCATCTGGCTGACGGTGATGATCGCGCTGATTGCCGTACCGGTAAACCTGGTGTTCGGCATTCTGCTGGCCTGGCTGGTGACGCGCTTTAA  >  2:629047/1‑149                                                                                                                                                        TTCAGCAAGGGGCTGATGCCGGTTTTACAGAATCTGGCCGATCCGGACATGCTGCACGCCATCTGGCTGACGGTGATGATCGCGCTGATTGCCGTACCGGTAAACCTGGTGTTCGGCATTCTGCTGGCCTGGCTGGTGACGCGCTTTAA  <  1:629047/149‑1                                                                                                                                                                                                                                                                                                              AACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2541378‑2541529 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcatTCAGCAAGGGGCTGATGCCGGTTTTACAGAATCTGGCCGATCCGGACATGCTGCACGCCATCTGGCTGACGGTGATGATCGCGCTGATTGCCGTACCGGTAAACCTGGTGTTCGGCATTCTGCTGGCCTGGCTGGTGACGCGCTTTAA  <  NC_000913/2541543‑2541396

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.