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.

AGCTGGACTGGGAAAACCTGGCGGCAGAAAAACAGACGCTGGTGTTCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGCTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAGCAAATGAACAGTCCATCGCTAATTATTATTGGTCGGGTTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGG  >  NC_000913_3_vch_pgi/3498921‑3499211                                                                                                                                               |                                                                                                                                                     aggtGGACTGGGAGAGCCAGGCGGGAGAAAGACAGACGCTGGTGTTCTATAGGGGGTTGAATCAGGCCGCGACTATTCGGCAAAAGCTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTc                                                                                                                                                <  2:195609/146‑1 (MQ=255)   cTGGACTGGGAAAACCTGGCGGCAGAAAAACAGACGCTGGTGTTCTATATGGGGATGAATCAGGCCGCGACTATTCAGCAAAAGCTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCAc                                                                                                                                              >  2:188561/1‑149 (MQ=255)        cTGGGAAAACCTGGCGGCAGAAAAACAGACGCTGGTGTTCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGCTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGtcaatcacgcgct                                                                                                                                         <  1:301759‑M1/149‑14 (MQ=255)             aaaaCCTGGCGGCAGAAAAACAGACGCTGGCGTTCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGCTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTg                                                                                                                                    >  2:207946/1‑149 (MQ=255)                  cTGGCGGCAGAAAAACAGACGGTGGTGTTCGATATGGGTGTGAATGAGGCCGCCACGATTCAGCAATATCTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGgcaatcacgcgctgcgtgactgc                                                                                                                               <  2:42936‑M1/149‑24 (MQ=255)                   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)                                                            aagcAGGCCGCGACTATTCAGCAAGAGCTGATTGAACACGGAATGACAGGAAAAATGCCGGTGGCAATTGTCGAAAACGGGACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAGCAAATGaa                                                                                     <  1:188751/146‑1 (MQ=255)                                                                       gACTATTCAGCAAAAGCTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGtcaatcacgcgctgcgtgactgccgtaccgttttcgacaattgccaccggcatttcgc                                                                                            >  1:159053‑M1/1‑73 (MQ=255)                                                                       gACTATTCAGCAAAAGCTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGtcaatcacgcgctgcgtgactgccgtaccgttttcgacaattgccaccggcatttcgc                                                                                            <  2:159053‑M1/131‑59 (MQ=255)                                                                           aTTCAGCAAAAGCTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAGCAAATGAACAGTCCATCGCTGat                                                                      >  2:389320/1‑149 (MQ=255)                                                                               aGCAAACGCTGATTGACCACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAGCAAATGAACAGTCCATCGCTAattatt                                                                  <  1:136307/149‑1 (MQ=255)                                                                                gCAAAAGCTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGga                                                                                                         >  1:352662/1‑109 (MQ=255)                                                                                gCAAAAGCTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGga                                                                                                         <  2:352662/109‑1 (MQ=255)                                                                                    aaGCTGATTGAACACCGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAGCAAATGAACAGTCCATCGCTAATTATTATTgg                                                             <  1:320952/149‑1 (MQ=255)                                                                                       cTGATTGAACACGGAATGCCAGGCGAAATGCAGGTGCCAAGTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAGCAAATGAACAGTCCATCGCTAATTATTATTGGTCg                                                          >  2:281130/1‑149 (MQ=255)                                                                                         ggtgggacacGGAAGGCCAGGCGAAATGGCGGTGGAAATTGTCGAAAAGGGGTCGGCAGTCACGCAGCGCCTGATTGACGGTACGCTCACACAGCTGGGAGTACTGGCGCAGCAAATGAGCAGTCCATCGCTAATTCTTAGTGGTCggg                                                        <  2:167508/143‑1 (MQ=255)                                                                                             gAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTAACACAGCTGGGAGAACTGGCGCAGCAAATGAACAGTCCATCGCTAATTATTATTGGTCGGgttgt                                                    >  1:190510/1‑149 (MQ=255)                                                                                                gacGGAATGACAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGGCACGCAGCGGGTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGTAGCTAATGAACAGTCCATCGCTAATTATTATTGGTCGGGTTGTTgg                                                 <  2:116115/148‑1 (MQ=255)                                                                                                 acGGAAGGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAGCAAATGAACAGTCCATCGCTAATTATTATTGGTCGGGTTGTTGGc                                                <  1:186875/149‑1 (MQ=255)                                                                                                             gCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGAcgtaccgttttcgacaattgccaccggcatttcgcctggcattccgtgttcaatcagcttttgctgaatagtc                                                      <  2:159053‑M2/131‑74 (MQ=255)                                                                                                             gCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGAcgtaccgttttcgacaattgccaccggcatttcgcctggcattccgtgttcaatcagcttttgctgaatagtc                                                      >  1:159053‑M2/1‑58 (MQ=255)                                                                                                                     ccGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACCGTACGCTCACACAGCTGGGAGAACTGGCGCAGCAAATGAACAGTCCATCGCTAATTATTATTGGTCGGGTCGTTGGCCTGCGCGATAaactgaactg                            >  2:380767/1‑149 (MQ=255)                                                                                                                                              cGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAGCAAATGAACAGTCCATCGCTAATTATTATTGGTCGGGTTGTTGGCCTGCGCGATCATCCTGACTGGTTCTCCAACCATTAATTTAACCCg   <  1:291715/149‑1 (MQ=255)                                                                                                                                               gtcaatcacgcgctgcgtgactgccgtaccgttttcgacaattgccaccggcatttcgcctggcattccgtgttcaatcagcttttgctgaatagtcgcggccggactcaaccccacatagaacatcagcgtctgtttttctgccgcca  <  1:310553‑M1/149‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                     AGCTGGACTGGGAAAACCTGGCGGCAGAAAAACAGACGCTGGTGTTCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGCTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAGCAAATGAACAGTCCATCGCTAATTATTATTGGTCGGGTTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGG  >  NC_000913_3_vch_pgi/3498921‑3499211

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