New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	BS168_reference_genome	= 2653333	22 (0.700)	3 (0.100)	3/224	NT	17.6%	pseudogene (341/471 nt)	sigKn	RNA polymerase sporulation‑specific sigma‑K factor precursor (Sigma‑27) (N‑terminal half)
	?	BS168_reference_genome	2701367 =	7 (0.230)					pseudogene (30/417 nt)	sigKc	RNA polymerase sporulation‑specific sigma‑K factor precursor (Sigma‑27) (C‑terminal fragment)

ATGTAGCTTGGTTATTTAAGGTATCCCCTTCATTCGCTAAGAAACGTTTAAGCATTTATTACAGGAAACTGACTCAACATTTATTTAGTCAAACTGTTACGGGAAATTTATGCCCCCCCTTTTTGTAATGAATAATGAAATCCTCATGCATTTGCGCGCATTGAAAAAAACAAAAAAAGACGTCTCCCTTCACGATCCGATCGGGCAGGATAAGGAAGGGAATGAAATAAGCCTGATTGACGTTTTGAAATCAGAAAACGAAGACGTGATAGATA  >  BS168_reference_genome/2701229‑2701503                                                                                                                                           |                                                                                                                                         aTGTAGCTTGGTTATTTAAGGTATCCCCTTCATTCGCTAAGAAACGTTTAAGCATTTATTACAGGAAACTGACTCAACATTTATTTAGTCAAACTGTTACGGGAAATTTATGCCCCCCCTTTTTGtaatgaataatgaa                                                                                                                                          <  8:205892/139‑1 (MQ=255)                      tATCCCCTTCATTCGCTAAGAAACGTTTAAGCATTTATTACAGGAAACTGACTCAACATTTATTTAGTCAAACTGTTACGGGAAATTTATGCCCCCCCTTTTTGTAATGAATAATGAAATCCTCATGCATTTGCGCGCa                                                                                                                     <  1:17068/139‑1 (MQ=255)                                     cTAAGAAACGTTTAAGCATTTTTTACAGGAAACTGACTCAACATTTATTGAGTCAAACTGTTACGGGAATTTTATGCCCCCCCTTTTTGTAATGAATTATGAAATCCTCATGCATTTGCGAGCATTGATAAAAACatt                                                                                                       >  4:264276/1‑136 (MQ=255)                                     cTAAGAAACGTTTAAGCATTTATTACAGGAAACTGACTCAACATTTATTTAGTCAAACTGTTACGGGAAATTTATGCCCCCCCTTTTTGTAATGAATAATGAAATCCTCATGCATTTGCGCGCATTGAAAAAAACaaa                                                                                                       >  1:231976/1‑138 (MQ=255)                                           catcggaacgatcgggcttattaaaggaattgaaagctattccgctggaaaagggacaaagctggcgacgtatgcagcgaggtgtattgaaaatgaAATCCTCATGCATTTGCGCGCATTGAAAAAAACAAAAAAAGAc                                                                                                >  2:53444‑M2/97‑139 (MQ=255)                                             tcggaacgatcgggcttattaaaggaattgaaagctattccgctggaaaagggacaaagctggcgacgtatgcagcgaggtgtattgaaaatgaAATCCTCATGCATTTGCGCGCATTGAAAAAAACAAAAAAAGACGt                                                                                              >  5:65179‑M2/95‑139 (MQ=255)                                                     atcgggcttattaaaggaattgaaagctattccgctggaaaagggacaaagctggcgacgtatgcagcgaggtgtattgaaaatgaAATCCTCATGCATTTGCGCGCATTGAAAAAAACAAAAAAAGACGTCTCCCTTc                                                                                      >  8:175954‑M2/87‑139 (MQ=255)                                                        ttattaCAGGAAACTGACTCAACATTTATTTAGTCAAACTGTTACGGGAAATTTATGCCCCCCCTTTTTGTAATGAATAATGAAATCCTCATGCATTTGCGCGCATTGAAAAAAACAAAAAAAGACGTCTCCCTTCAc                                                                                    >  3:254233/1‑138 (MQ=255)                                                                                                         aaTTTATGCCCCCCCTTTTTGTAATGAATAATGAAATCCTCATGCATTTGCGCGCATTGAAAAAAACAAAAAAAGACGTCTCCCTTCACGATCCGCTCGGGCAGGATAAGTAAGGGAATGAAATAAGCCTGATTGACGt                                  >  1:244042/1‑139 (MQ=255)                                                                                                         aaTTTATCCCCCCCCTTTTTGTAATGAATAATGAAATCCTCATGCATTTGCGCGCAATGAAAAAAACAAAAAAAGACTTCTCCCTTCACGATCTGATCGGGCAGGATAAGTAAGGGATTGAAATAAGCCTGATGCACGt                                  >  3:170695/1‑139 (MQ=255)                                                                                                                   cccccTTTTTGTAATGAAAAATGAAATCCTCATGCATTTGCGCGCATTGAAAAAAACAAAAAAAGACGTCTCCCTTCCCGATCCGATCGTGCAGGATAAGGAAGGGAATGAAATACGCCTGATTGACGTTTTGAAAt                          >  4:7710/1‑137 (MQ=255)                                                                                                                                        tgaAATCCTCATGCATTTGCGCGCATTGAAAAAAACAAAAAAAGACGTCTCCCTTCACGATCCGATCGGGCAGGATAAGGAAGGGAATGAAATAAGCCTGATTGACGTTTTGAAATCAGAAAACGAAGACGTgataga    >  5:75899‑M2/4‑138 (MQ=255)                                                                                                                                         taAAACCCCATCCATTGGCGCCAATACAAAAAAAAAAAAACCCACGGCTCCCTTCACGATCCGATCGGGCAGGATAAGGAAGGGAATGAAATAAGCCTGATTGACGTTTTGAAATCAGAAAACGAAGACGTgatagata  <  4:254233‑M2/137‑1 (MQ=255)                                                                                                                                           |                                                                                                                                         ATGTAGCTTGGTTATTTAAGGTATCCCCTTCATTCGCTAAGAAACGTTTAAGCATTTATTACAGGAAACTGACTCAACATTTATTTAGTCAAACTGTTACGGGAAATTTATGCCCCCCCTTTTTGTAATGAATAATGAAATCCTCATGCATTTGCGCGCATTGAAAAAAACAAAAAAAGACGTCTCCCTTCACGATCCGATCGGGCAGGATAAGGAAGGGAATGAAATAAGCCTGATTGACGTTTTGAAATCAGAAAACGAAGACGTGATAGATA  >  BS168_reference_genome/2701229‑2701503

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