New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP009273	2534351 =	38 (0.990)	4 (0.110)	4/302	NT	11.5%	coding (688/876 nt)	cysW	sulfate/thiosulfate ABC transporter subunit
	?	CP009273	2534553 =	25 (0.680)					coding (486/876 nt)	cysW	sulfate/thiosulfate ABC transporter subunit

GTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAGCAAAATCGCCGCTTCGTCTTCCTGGCTGCCCTGGCTTAACATCACCGGCACCAGTTCGCGCACCACAAACGGACACGTCACGAAGATGGTGACCAGCACCATTCCCGGCCAGGAGAACATAATTTGCAGGTTATGCTCGTCGAGCCAACCGCCGAGCGGGCCGTTAGAGCCGTAGAACAGCAAATACACCAGACCGGCAACCACCGGCGATACGGCAAACGGAATGTCCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTAC  >  CP009273/2534281‑2534810                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                  gtaagtgAGCCGACGTTGTTGTAGTCCGGCTCCAGCAATTAAATCGGTAACGGCGGCGACAGGGTTTTGCCGCAAATCGACCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAGCAAAATCGCCGCTTCGTCTTCCTGGCTGCCCTGGCTTAACATCACCGGCACCAGTTCGCGcac                                                                                                                                                                                                                                                                 <  1:187927/269‑1 (MQ=255)              cGGGGTTGTANGCCTGCTCCAGCCAGCCAAACTGAAACGGCGTAAACAAGGTTATGCCGCGACTCTGGCCGAAAACCACAAAAACCGGGCAAAACGCGCCAATTGCGCGCGCGTTGGCAAACAACACGCCATAAAGCTCCGCCCAGCGGATCTTCGGTCAGGTAACGCGACGGAACACTTGCCACCCGCACGCGCCACGCCAAATCGCCGCTTCGTCTTCCTGGCTGCCACGGCTTAAAAGCCCCGGCACCAGCTCGCGCACCACAAACGGaca                                                                                                                                                                                                                                                     <  2:57649/274‑1 (MQ=255)                 tgttgtAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCAAAGGAACATCTGCCAGCCGGACGCGCCAAGCAAAATCGCCGCTTCGTCTTCCTGGCTGCCCTGGCTTAACATCACCGGCACCAGTTAGCGCACCACACACCGACACGt                                                                                                                                                                                                                                                  >  1:446347/1‑274 (MQ=255)                                 actcAAGTCCATCCGTACCGGCAGCGACAGGGGTTCGACGCGACGCGAGCCGGAAACCACTGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCCCGCCATAACGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAGCAAAATCGCCGCTTCGTCTTCCTGGCTGCCCTGGCTTAACATCACCGGCACCAGTTCGCGCACCACAAACGGACACGTCACGAAGATGGTGAcc                                                                                                                                                                                                                                  <  2:287067/271‑1 (MQ=255)                                                                    cgccgcGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAGCAAAATCGCCGCTTCGTCTTCCTGGCTGCCCTGGCTTAACATCACCGGCACCAGTTCGCGCACCACAAACGGACACGCCACGAAGAGGGTGACCAGCACCATACCCGGCCAGAAGAACATAATTTGCAg                                                                                                                                                                                               >  1:237616/1‑274 (MQ=255)                                                                                     gAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAGCAAAATCGCCGCTTCGTCTTCCTGGCTGCCCTGGCTTAACATCACCGGCACCAGTTCGCGCACCACAAACGGACACGTCACGaa                                                                                                                                                                                                                                            >  2:250006/1‑212 (MQ=255)                                                                                         ccaccGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAGCAAAATCGCCGCTTCGTCTTCCTGGCTGCCCTGGCTTAACATCACCGGCACCAGTTCGCGCACCACAAACGGACACGTCACGAAGATg                                                                                                                                                                                                                                        <  1:250011/212‑1 (MQ=255)                                                                                                gccagggcagccaggaagacgaagcggcgattttgcttggcgcgtccggctggcagatgttccgtcgcgtcacattaccgaacatccgctgggcgctgctttatggcgtggtgttgaccaacgcccgcgcaattggcgagtttggcgcggtgtcggtggtttccggctcgattcgcgCACCACAAACGGACACGTCACGAAGATGGTGACCAGCACCATTCCCgg                                                                                                                                                                                                                    >  2:153965‑M2/178‑225 (MQ=255)                                                                                                    gggcagccaggaagacgaagcggcgattttgcttggcgcgtccggctggcagatgttccgtcgcgtcacattaccgaacatccgctgggcgctgctttatggcgtggtgttgaccaacgcccgcgcaattggcgagtttggcgcggtgtcggtggtttccggctcgattcgcgCACCACAAACGGACACGTCACGAAGATGGTGACCAGCACCATTCCCGGCCAg                                                                                                                                                                                                                <  1:153968‑M2/52‑1 (MQ=255)                                                                                                               gCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAACAGCGCACAGCGGGTGATCGGTAATGAGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAGCAAAATCGCCGCTTCGTCTTCCTGGCTGCCCTGGCTCAAAATAACCGGCACCAGTTCGCGCACCAAAAACGGACACGTCACGAAGATGGTGACCAGCACCATTCCCGACAAGGAGAACATCATATGCAGGTTATGCTCGTAGAACCAACCTCACAGAGGGCCgggagatgcgt                                                                                                                                                   >  2:181511/1‑265 (MQ=255)                                                                                                                       cgcggggtctcgTCAACACCCCGACATAAAGCAGCGCCCAGCGGATGTGCGGTAATGAGACGCGGCGGAACATCTGCCAGCCGGACGCGCCAAGCAAAATCGCCGCTTCGTCTTCCTGGCTGCCCTGGCTTAACATCACCGGCACCAGTTCGCGCACCACAAACGGACACGTCACGAAGATGGTGACCAGCACCATTCCCGGCCAGGAGAACATAATTTGCAGGTTATGCTCGTCGAGCCAACCGCCGAGCGGGCCGTTAGAGCCGTAGAACAGc                                                                                                                                           <  1:81809/264‑1 (MQ=255)                                                                                                                                            cGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAGCAAAATCGCCGCTTCGTCTTCCTGGCTGCCCTGGCTTAACATCACCGACACCAGTTCGCGCACCACAAACGGACACGTCACGAAGATGGTGACCAGCACCATTCCCGGCCAGGAGAAAATAATTTGCAGGTTATGCTCGTCGAGCCAACCGCACAGCGGGCCGATAGAGCCGTAGAAAAGCACATACACCAGACCGGCAac                                                                                                                       >  2:709952/1‑274 (MQ=255)                                                                                                                                                     gcagggcccgtgggATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAGCAAAATCGCCGCTTCGTCTTCCTGGCTGCCCTGGCTTAACATCACCGGCACCAGTTCGCGCACCACAAACGGACACGTCACGAAGATGGTGACCAGCACCATTCCCGGCCAGGAGAACATAATTTGCAGGTTATGCTCGTCGAGCCAACCGCCGAGCGGGCCGTTAGAGCCGTAGAACAGCAAATACACCAGACCGGCAAGCACCGGCGAt                                                                                                             <  1:709962/263‑1 (MQ=255)                                                                                                                                                          ggccAGGGGATGTTCGGTCATGTGACGCGACGGAACATCTGCCAGCCGGACTCCCCAAGCAAAATCGCCGCGTCGTCTTCCTGGCTGCCCTGGCTTAACATCACCGGCACCAGTTCGCGCACCACAAACGGACACGTCACGAAGATGGTGACCAGCACCATTCCCGGCCAGGAGAACATAATTTGCAGGTTATGCTCGTCGAGCCAACCGCCGAGCGGGCCGTTAGAGCCGTAGAACAGCAAATACACCAGACCGGCAACCACCGGCGATACGGc                                                                                                        <  1:181514/273‑1 (MQ=255)                                                                                                                                                              aGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAGCAAAATCGCCGCTTCGTCTTCCTGGCTGCCCTGGCTTAACATCACCGGCACCAGTTCGCGCACCACAAACGGACACGTCACGAAGATGGTGACCAGCACCa                                                                                                                                                                                                                           >  1:16196/1‑156 (MQ=255)                                                                                                                                                                  gATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAGCAAAATCGCCGCTTCGTCTTCCTGGCTGCCCTGGCTTAACATCACCGGCACCAGTTCGCGCACCACAAACGGACACGTCACGAAGATGGTGACCAGCACCATTcc                                                                                                                                                                                                                       <  2:16196/156‑1 (MQ=255)                                                                                                                                                                      ttCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAGCAAAATCGCCGCTTCGTCTTCCTGGCTGCCCTGGCTTAACATCACCGGCACCAGTTCGCGCACCACAAACGGACACGTCACGAAGATGGTGACCAGCACCATTCCCGGCCAGGAGAACa                                                                                                                                                                                                         >  1:324268/1‑166 (MQ=255)                                                                                                                                                                          gTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAGCAAAATCGCCGCTTCGTCTTCCTGGCTGCCCTGGCTTAACATCACCGGCACCAGTTCGCGCACCACAAACGGACACGTCACGAAGATGGTGACCAGCACCATTCCCGGCCAGGAGAACATAAt                                                                                                                                                                                                     <  2:324262/166‑1 (MQ=255)                                                                                                                                                                                                        cGGACGCGCCAAGCAAAATCGCCGCTTCGTCTTCCTGGCTGCCCTGGCTTAACATCACCGGCACCAGTTCGCGCACCACAAACGGACACGTCACGAAGATGGTGACCAGCACCATTCCCGGCCAGGAGAACATAATTTGCAGGTTATGCTCGNCGAGCCAACCGCCGAGCGGGCCGTTAGAGCCGTAGAACAGCAAATACACCAGACCGGCAACCACCGGCGATACGGCAAACGGAATGTCCAGTAGCGTCAGCAGTAACTGGCGTCCGGGGAAg                                                          >  1:615770/1‑275 (MQ=255)                                                                                                                                                                                                          gACGCGCCAAGCAAAATCGCCGCTTCGTCTTCCTGGCTGCCCTGGCTTAACATCACCGGCACCAGTTCGCGCACCACAAAc                                                                                                                                                                                                                                                          >  1:607372/1‑81 (MQ=255)                                                                                                                                                                                                             gcgcCAAGCAAAATCGCCGCTTCGTCTTCCTGGCTGCCCTGGCTTAACATCACCGGCACCAGTTCGCGCACCACAAACg                                                                                                                                                                                                                                                         >  1:560863/1‑79 (MQ=255)                                                                                                                                                                                                              cgcCAAGCAAAATCGCCGCTTCGTCTTCCTGGCTGCCCTGGCTTAACATCACCGGCACCAGTTCGCGCACCACAAACGGac                                                                                                                                                                                                                                                      <  2:607363/81‑1 (MQ=255)                                                                                                                                                                                                                 cAAGCAAAATCGCCGCTTCGTCTTCCTGGCTGCCCTGGCTTAACATCACCGGCACCAGTTCGCGCACCACAAACGGaca                                                                                                                                                                                                                                                     <  2:560854/79‑1 (MQ=255)                                                                                                                                                                                                                        acgcccgcgcaattggcgagtttggcgcggtgtcggtggtttccggctcgattcgcgCACCACAAACGGACACGTCACGAAGATGGTGACCAGCACCATTCCCGGCCAGGAGAACATa                                                                                                                                                                                                       >  1:220131‑M2/58‑118 (MQ=255)                                                                                                                                                                                                                            ccgcgcaattggcgagtttggcgcggtgtcggtggtttccggctcgattcgcgCACCACAAACGGACACGTCACGAAGATGGTGACCAGCACCATTCCCGGCCAGGAGAACATAAttt                                                                                                                                                                                                   <  2:220126‑M2/65‑1 (MQ=255)                                                                                                                                                                                                                             cctcgtCGTCTGCCTGGCTGCCCTGGCTTAACATCACCGGCACCCGTTCGCGCACCACAAACGGACACGTCACGAAGATGGTGACCAGCACCATTCCCGGCCAGGAGAACATAATTTGCNGGTTATGCTCGTCGAGCCAACCGCCGAGCGGGCCGTTAGAGCCGTAGAACAGCAAATACACCAGACCGGCAACCACCGGCGATACGGCAAACGGAATGTCCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCAccagc                                        <  1:613880/267‑1 (MQ=255)                                                                                                                                                                                                                                             gCTGCCCTGGCTTAACATCACCGGCACCAGTTCGCGCACCACAAACGGACACGTCACGAAGATGGTGACCAGCACCATTCCCGGCCAGGAGAACATAATTTGCAGGTTATGCTCGTCGAGCCAAccgc                                                                                                                                                                        >  2:721092/1‑128 (MQ=255)                                                                                                                                                                                                                                                 cccTGGCTTAACATCACCGGCACCAGTTCGCGCACCACAAACGGACACGTCACGAAGATGGTGACCAGCACCATTCCCGGCCAGGAGAACATAATTTGCAGGTTATGCTCGTCGAGCCAACCGCCGAg                                                                                                                                                                    <  1:721102/128‑1 (MQ=255)                                                                                                                                                                                                                                                                aCCGGCACCAGTTCGCGCACCACAAACGGACACGTCACGAAGATGGTGACCAGCACCATTCCCGGCCAGGAGAACATAATTTGCAGGTTATGCTCGTCGAGCCAACCGCCGAGCGGGCCGTTAGAGCCGTAGAACAGCAAATACACCAGACCGGCAACCACCGGCGATACGGCAAACGGAATGTCCAGTa                                                                                       >  1:335750/1‑190 (MQ=255)                                                                                                                                                                                                                                                                aCCGGCACCAGTTCGCGCACCACAAACGGACACGTCACGAAGATGGTGACCAGCACCATTCCCGGCCAGGAGAACATAATTTGCAGGTTATGCTCGTCGAGCCAACCGCCGAGCGGGCCGTTAGAGCCGTAGAACAGCAAATACACCAGACCGGCAACCACCGGCGATACGGCAAACGGAATGTCCAGTAGCGCCAGAAGTAACGGGCGACCAGGGAAGATAAAGCGGGTAACCAGCCAGGCCAGCAGAGTGACGAACACCAGGTTGAACgggtc  >  2:207020/1‑272 (MQ=255)                                                                                                                                                                                                                                                                              gcgCACCACAAACGGACACGTCACGAAGATGGTGACCAGCACCATTCCCGGCCAGGAGAACATAATTTGCAGGTTATGCTCGTCGAGCCAACCGCCGAGCGGGCCGTTAGAGCCGTAGAACAGCAAATACACCAGACCGGCAACCACCGGCGATACg                                                                                                          >  2:394474‑M2/4‑157 (MQ=255)                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                  GTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAGCAAAATCGCCGCTTCGTCTTCCTGGCTGCCCTGGCTTAACATCACCGGCACCAGTTCGCGCACCACAAACGGACACGTCACGAAGATGGTGACCAGCACCATTCCCGGCCAGGAGAACATAATTTGCAGGTTATGCTCGTCGAGCCAACCGCCGAGCGGGCCGTTAGAGCCGTAGAACAGCAAATACACCAGACCGGCAACCACCGGCGATACGGCAAACGGAATGTCCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTAC  >  CP009273/2534281‑2534810

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ 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.