New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_vch_pgi	= 3499063	19 (0.860)	4 (0.180)	4/280	2.6	17.1%	coding (1236/1374 nt)	cysG	fused siroheme synthase 1,3‑dimethyluroporphyriongen III dehydrogenase and siroheme ferrochelatase/uroporphyrinogen methyltransferase
	?	NC_000913_3_vch_pgi	= 3499086	20 (0.910)					coding (1259/1374 nt)	cysG	fused siroheme synthase 1,3‑dimethyluroporphyriongen III dehydrogenase and siroheme ferrochelatase/uroporphyrinogen methyltransferase
Rejected: Frequency below/above cutoff threshold.

TGGCGGCAGAAAAACAGACGCTGGTGTTCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGCTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAG‑CAAATGAACAGTCCATCGCTAATTATTATTGGTCGGGTTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTC  >  NC_000913_3_vch_pgi/3498939‑3499231                                                                                                                                                    |                                                                                                                                                   tGGCGGCAGAAAAACAGACGCTGATGTTCTATGTGGGGTTGAGTCCGGCCGCGACTATTCAGCAAAAGCTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGAc                                                                                                                                                   <  1:310553‑M2/149‑2 (MQ=255)           aaaaCAGTCGCTGGTGTTCTATAGGGGGTTGAATCAGGCCGCGACTATTCAGCCAAGGGGCTGTTAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTca                                                                                                                                         <  1:91497/149‑1 (MQ=255)                        gtgtTCTATATGGGGTAGAATCAGGCCGCGACTATTCAGCAAAAGCTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACGGCTGGgaga                                                                                                                            >  1:341889/1‑149 (MQ=255)                            tCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGCTGATTGAACACGGAATGCCAGGCGAGATGCGGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGGAGGACTg                                                                                                                        >  2:74340/1‑149 (MQ=255)                                          aagcAGGCCGCGACTATTCAGCAAGAGCTGATTGAACACGGAATGACAGGAAAAATGCCGGTGGCAATTGTCGAAAACGGGACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAG‑CAAATGaa                                                                                                         <  1:188751/146‑1 (MQ=255)                                                         aTTCAGCAAAAGCTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAG‑CAAATGAACAGTCCATCGCTGat                                                                                          >  2:389320/1‑149 (MQ=255)                                                             aGCAAACGCTGATTGACCACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAG‑CAAATGAACAGTCCATCGCTAattatt                                                                                      <  1:136307/149‑1 (MQ=255)                                                              gCAAAAGCTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGga                                                                                                                              >  1:352662/1‑109 (MQ=255)                                                              gCAAAAGCTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGga                                                                                                                              <  2:352662/109‑1 (MQ=255)                                                                  aaGCTGATTGAACACCGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAG‑CAAATGAACAGTCCATCGCTAATTATTATTgg                                                                                 <  1:320952/149‑1 (MQ=255)                                                                     cTGATTGAACACGGAATGCCAGGCGAAATGCAGGTGCCAAGTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAG‑CAAATGAACAGTCCATCGCTAATTATTATTGGTCg                                                                              >  2:281130/1‑149 (MQ=255)                                                                       ggtgggacacGGAAGGCCAGGCGAAATGGCGGTGGAAATTGTCGAAAAGGGGTCGGCAGTCACGCAGCGCCTGATTGACGGTACGCTCACACAGCTGGGAGTACTGGCGCAG‑CAAATGAGCAGTCCATCGCTAATTCTTAGTGGTCggg                                                                            <  2:167508/143‑1 (MQ=255)                                                                           gAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTAACACAGCTGGGAGAACTGGCGCAG‑CAAATGAACAGTCCATCGCTAATTATTATTGGTCGGgttgt                                                                        >  1:190510/1‑149 (MQ=255)                                                                              gacGGAATGACAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGGCACGCAGCGGGTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGTAG‑CTAATGAACAGTCCATCGCTAATTATTATTGGTCGGGTTGTTgg                                                                     <  2:116115/148‑1 (MQ=255)                                                                               acGGAAGGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAG‑CAAATGAACAGTCCATCGCTAATTATTATTGGTCGGGTTGTTGGc                                                                    <  1:186875/149‑1 (MQ=255)                                                                                           gCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGAcgtaccgttttcgacaattgccaccggcatttcgcctggcattccgtgttcaatcagcttttgctgaatagtc                                                                           >  1:159053‑M2/1‑58 (MQ=255)                                                                                           gCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGAcgtaccgttttcgacaattgccaccggcatttcgcctggcattccgtgttcaatcagcttttgctgaatagtc                                                                           <  2:159053‑M2/131‑74 (MQ=255)                                                                                                   ccGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACCGTACGCTCACACAGCTGGGAGAACTGGCGCAG‑CAAATGAACAGTCCATCGCTAATTATTATTGGTCGGGTCGTTGGCCTGCGCGATAaactgaactg                                                >  2:380767/1‑149 (MQ=255)                                                                                                                            cGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAG‑CAAATGAACAGTCCATCGCTAATTATTATTGGTCGGGTTGTTGGCCTGCGCGATCATCCTGACTGGTTCTCCAACCATTAATTTAACCCg                       <  1:291715/149‑1 (MQ=255)                                                                                                                             tgCAGTCACGCAGCGTGTGACTGAAGGTGCGGTCACACAGATGGTAGCGTTGGCGCAA‑CAAATGAACAGTCCATCGCTAATTATTATTGGGCGGGTTGTTGGTGGTCGCGATAAACGGAACTAGTTCTACAACCATTAATTTAACCCgg                      <  2:356478/148‑1 (MQ=255)                                                                                                                              gCAGTCACGCAGCGCGTGATTGCcgtaccgttttcgacaattgccaccggcatttcgcctggcattccgtgttcaatcagatattgctgaatcgtggcggcctcattcacacccatatcgaacaccaccgtctgtttttctgccgccag                      <  2:42936‑M2/149‑127 (MQ=255)                                                                                                                                        aGCGCGTGATTGAcgtaccgttttcgacaattgccaccggcatttcgcctggcattccgtgttcaatcagcttttgctgaatagtcgcggcctgattcaaccccatatagaacaccagcgtctgtttttctgccgccaggttttcccag            <  1:301759‑M2/149‑137 (MQ=255)                                                                                                                                           gcgTGAATGACGGTACGCTCCCACAGCTGGGAGAACTGGCGGAG‑CAAAGGAACAGTCCATCGCTAATTATTATTGGTCGGGTTGTTGCGCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTg        <  2:341889/149‑1 (MQ=255)                                                                                                                                           acgTAATCGACGGTACGCTCACCAAGCTGGGAGAACTGGCGCAGACAACGCAACCG‑CCGTCGCTAATTATCATTAGTCAGGTTCTTGGCCAGCGCGATAAACTGAACGGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTg        <  1:74340/148‑1 (MQ=255)                                                                                                                                                 ttGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAA‑CAAATGAACAGTCCATCGCTAATTATTATTGGTCGGGTTGTTGGCCTTCGCGATAAACTCAACTGGTTCTCCAAACATCAGTTTCAACCGGTCCAGACCATACTCGCcttc  >  2:98791/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                   TGGCGGCAGAAAAACAGACGCTGGTGTTCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGCTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAG‑CAAATGAACAGTCCATCGCTAATTATTATTGGTCGGGTTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTC  >  NC_000913_3_vch_pgi/3498939‑3499231

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ 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.