New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3948434	26 (1.030)	4 (0.170)	4/238	2.4	14.0%	intergenic (+10/+15)	yifE/yifB	UPF0438 family protein/magnesium chelatase family protein and putative transcriptional regulator
	?	NC_000913	= 3948439	25 (1.060)					intergenic (+15/+10)	yifE/yifB	UPF0438 family protein/magnesium chelatase family protein and putative transcriptional regulator
Rejected: Frequency below/above cutoff threshold.

AGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAGATGGATGAGCAAACGGTCAATCGCTCGATAGCTAACTGCCTCCTGCAAATGCTGACGTGTGATAATGTCAGACTGATCAATATCAGCAATGGTTCGTGCAACTTTC  >  NC_000913/3948291‑3948575                                                                                                                                                     |                                                                                                                                         aGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGATGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAAt                                                                                                                                          >  2:19604/1‑149 (MQ=255)             aaGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTa                                                                                                                              <  5:2070/149‑1 (MQ=255)                        gtgtGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGtt                                                                                                                   >  6:89814/1‑149 (MQ=255)                         tgtgGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGttt                                                                                                                  <  5:89814/149‑1 (MQ=255)                                        tGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAGATGGATGa                                                                                                    <  6:80074/148‑1 (MQ=255)                                          aCGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGATGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAGATGGATGAGc                                                                                                  <  1:19604/148‑1 (MQ=255)                                                 ttAAGCGTCCAAAACGTTTTCATACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTAt                                                                                                                             <  1:32161/114‑1 (MQ=255)                                                 ttAAGCGTCCAAAACGTTTTCATACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTAt                                                                                                                             >  2:32161/1‑114 (MQ=255)                                                     gCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTa                                                                                                                              >  4:81683/1‑109 (MQ=255)                                                     gCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTa                                                                                                                              <  3:81683/109‑1 (MQ=255)                                                            aaaCGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAGATGGATGAGCAAACGGTCAATCGCTCGa                                                                                >  1:102318/1‑148 (MQ=255)                                                                        acCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAGATGGATGAGCAAACGGTCAATCGCTCGATAGCTAACTGcct                                                                   >  8:43703/1‑149 (MQ=255)                                                                        acCCATTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACCCCCATTCTGACCATAAAAAAAAAGGGCGAAATGCCCTTTTTTTAAGTCAGTAGTTTCTGCAGATGGATGAGCAAACGGTCAATCGCTCCATAGCTAACTGgct                                                                   >  6:35859/1‑149 (MQ=255)                                                                                    ggTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAGATGGATGAGCAAACGGTCAATCGCTCGATAGCTAACTGCCTCCTGCAAATGCt                                                       <  7:43703/149‑1 (MQ=255)                                                                                     gTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAGATGGATGAGCAAACGGTCAATCTCTCGATAGCTAACTGCCTCCTGCAAATGCt                                                       >  3:78016/1‑148 (MQ=255)                                                                                           cGCAGGTTGAAGGTGCTGAAGACTACTCCGATTCTGACGATAAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAGTTGGATGAGCAAACGGTCAATCGCTCGATAGCTAACTGCCTCCTGCAAATGCTGACgtgt                                                <  5:35859/149‑1 (MQ=255)                                                                                                             aaGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAGATGGATGAGCAAACGGTCAATCGCTCGATAGCTAACTGCCTCCTGCAAATGCTGACGTGTGATAATGTCAGACTGa                                <  2:102318/147‑1 (MQ=255)                                                                                                             aaGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAGATGGATGAGCAAACGGTCAATCGCTCGATAGCTAACTGCCTCCTGCAAATGCTGACGTGTGATAATGTCAGACTGAt                               <  1:134627/148‑1 (MQ=255)                                                                                                                     acCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAGATGGATGAGCAAACGGTCAATCGCTCGATAGCTAACTGCCTCCTGCAAATGCTGAc                                                    >  4:104050/1‑119 (MQ=255)                                                                                                                     acCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAGATGGATGAGCAAACGGTCAATCGCTCGATAGCTAACTGCCTCCTGCAAATGCTGAc                                                    <  3:104050/119‑1 (MQ=255)                                                                                                                             tGACGATTAATAAAAAGGGCGAAATgccctttttattaatcgtcagaatcggtgtagtcttcagcaccttcaacctgcggtttaccgccggaaagg                                                                   <  5:30694‑M2/96‑72 (MQ=255)                                                                                                                             tGACGATTAATAAAAAGGGCGAAATgccctttttattaatcgtcagaatcggtgtagtcttcagcaccttcaacctgcggtttaccgccggaaagg                                                                   >  6:30694‑M2/1‑25 (MQ=255)                                                                                                                               aCGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAGATGGATGAGCAAACGGTCAATCGCTCGATAGCTAACTGCCTCCTGCaaa                                                           >  4:62954/1‑102 (MQ=255)                                                                                                                               aCGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAGATGGATGAGCAAACGGTCAATCGCTCGATAGCTAACTGCCTCCTGCaaa                                                           <  3:62954/102‑1 (MQ=255)                                                                                                                                      ataaAAAGGGCGAAATgccctttttattaatcgtcagaatcggtgtagtcttcagcaccttcaacctgcggtttaccgccggaaagggtgtgaaaacgtttt                                                    >  4:60841‑M2/1‑16 (MQ=255)                                                                                                                                      ataaAAAGGGCGAAATgccctttttattaatcgtcagaatcggtgtagtcttcagcaccttcaacctgcggtttaccgccggaaagggtgtgaaaacgtttt                                                    <  3:60841‑M2/102‑87 (MQ=255)                                                                                                                                       taaAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAGATGGATGAGCAAACGGTCAATCGCTCGATAGCTAACTGcctc                                                                  >  7:70608/1‑87 (MQ=255)                                                                                                                                       taaAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAGATGGATGAGCAAACGGTCAATCGCTCGATAGCTAACTGcctc                                                                  <  8:70608/87‑1 (MQ=255)                                                                                                                                         aaaaGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAGATGGATGAGCAAACGGTCAATCGCTCGATAGCTAACTGCCTCCTGCAAATGCTGACGTGTGATAATGTCAGACTGATCAATATCAGCAATGGTTCGTGCAACttt   >  7:71973/1‑148 (MQ=255)                                                                                                                                         aaaaGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAGATGGATGAGCAAACGGTCAATCGCTCGATAGCTAACTGCCTCCTGCAAATGCTGACGTGTGATAATGTCAGACTGATCAATATCAGCAATGGTTCGTGCAACTTTc  <  8:34303/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                         AGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAGATGGATGAGCAAACGGTCAATCGCTCGATAGCTAACTGCCTCCTGCAAATGCTGACGTGTGATAATGTCAGACTGATCAATATCAGCAATGGTTCGTGCAACTTTC  >  NC_000913/3948291‑3948575

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