New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4160853 =	175 (2.030)	8 (0.090)	8/428	NT	4.4%	intergenic (‑63/‑271)	sthA/fabR	pyridine nucleotide transhydrogenase, soluble/transcriptional repressor of fabA and fabB
	?	NC_000913	= 4499511	NA (NA)					noncoding (1331/1331 nt)	IS2	repeat region
Rejected: Frequency below/above cutoff threshold.

ACTATACCAGAGAATGAACATACTGGATTGCTCCAGTATTGTGAATGAACGGTAACGCAATAAAACATATTTACTGCAATTCTTTGGCCGGTTCTTTTACGTACAGCGGAAACCTGCCGCTTAAACGGAGAGTATCGTCGATAAAAATCCAATAAAACGTCAGGGCAAAAGTAAGAAACAGACAAAGCAAAGGCCGCTCAGGATATAGCCAGATAAATGACGGGGATCAATTGGCt‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4161087‑4160853 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCG  <  NC_000913/4499511‑4499274                                                                                                                                                                                                                                                                                                                                                                                                                                                                                           ACTATACCAGAGAATGAACATACTGGATTGCTCCAGTATTGTGAATGAACGGTAACGCAATAAAACATATTTACTGCAATTCTTTGGCCGGTTCTTTTACGTACAGCGGAAACCTGCCGCTTAAACGGAGAGTATCGTCGATAAAAATCCAATAAAACGTCAGGGCAAAAGTAAGAACAAGACAAAGCAAAGACCGCTCAGGATATAACCAGATAAATAACGGGGAACAATAGGCTaa                                                                                                                                                                                                                                             >  2:104438/1‑236               TGAACATACTGGATTGCTCCAGTATTGTGAATGAACGGTAACGCAATAAAACATATTTACTGCAATTCTTTGGCCGGTTCTTTTACGTACAGCGGAAACCTGCCGCTTAAACGGAGAGTATCGTCGATAAAAATCCAATAAAACGTCAGGGCAAAATTAAGAAACAGACAAAGCAAAGGCCGCTCAGGATATAGCCAGATTAATAACGGGGAGCAATTGGCTGGATTTGCCCCTATATT                                                                                                                                                                                                                              >  1:690682/1‑239                                                   GGTAACGCAATAAAACATATTTACTTCAATTCCCTGGCCGGTTCTTTTACGTACAGCGGAAACCTGCCGCTTAAACGGAGAGTATCGTCGATAAAAATCCAATAAAACGTCAGGGCAAAAGTAAGAAACAGACAAAGCAAAGGCCGCTCAGGATATAGCCAGATAAATGACGGGGATCAATTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACC                                                                                                                                                                                                     <  1:504355/228‑1                                                             TAAAACATATTTACTGCAATTCTCTGGCCGGTTCTTTTACGTACAGCGGAAACCTGCCGCTTAAACGGAGAGTATCGTCGATAAAAATCCAATAAAACGTCAGGGCAAAAGTAAGAAACAGACAAAGCAAAGGCCGCTCAGGATATAGCCAGATAAATGACGGGGATCAATTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATACATTAACCCCTTAAAAcc                                                                                                                                                                                           >  2:504355/1‑226                                                                                            TTCTTTTACGTACAGCGGAAACCTGCCGCTTAAACGGAGAGTATCGTCGATAAAAATCCAATAAAACGTCAGGGCAAAAGTAAGAAACAGACAAAGCAAAGGCCGCTCAGGATATAGCCAGATAAATGACGGGGATCAATTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGC                                                                                                                                                 >  1:230428/1‑239                                                                                                    acgacaggcGAAACCTGCCGCTTAAAAGGAGAGTATCGTCGATAAAAATCCAAGAAAACGTCAGGGCAAAAGTAAGAAACAGACAAAGCAAAGGCCGCTCAGGATATAGCCAGATAAATGACGGGGATCAATTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCTCTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCG                                                                                                                                         <  2:690682/230‑1                                                                                                         AGCGGAAACCTGCCGCTTAAACGGAGAGTATCGTCGATAAAAATCCAATAAAACGTCAGGGCAAAAGTAAGAAACAGACAAAGCAAAGGCCGCTCAGGATATAGCCAGATAAATGACGGGGATCAATTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCAGCAGATATTCCCTTTGCGAGCGATAACCCAGCGCACTATGCGGAgg                                                                                                                                     >  1:398167/1‑236                                                                                                                                                                                 aacgaaaAAAGCAAAGGCCGCTCAGGATATAGCCAGAAAAATGGGGGGGATAAATTGGCTGGATTTGCCCCTTTATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATA                                                            <  2:230428/232‑1                                                                                                                                                                                                                                   TCAATTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAACGCTCTGCTATTCCGTT          >  1:768235/1‑239                                                                                                                                                                                                                                            TGGTTTTGCCCCTTTATTTCCAGACATCTTGTTTCACTTTACCCTGTAACTGCCCGCTCGCGGAGATATTCCCGTTGCGAGCGATCACCCAGCGCACTTTGCGGATGCCATTCGTTATAATGCGCGAACGCCTCTGCAAGTTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCG  <  2:859594/238‑1                                                                                                                                                                                                                                            TGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATTTAGTCACGCTTTATCGTTTTACCGAAGCTCCCTGCTATTCCGTTACTCCCCG  >  2:166565/1‑238                                                                                                                                                                                                                                            TGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCATGAAGCTCTCTGCTATTCCGTTAtcttccg  >  1:242091/1‑231                                                                                                                                                                                                                                            TGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGCCTGGTTTGGGCATGATACTGATGTGGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCG  >  2:959072/1‑238                                                                                                                                                                                                                                            TGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGG                                                                                                                                                                   <  1:556777/77‑1                                                                                                                                                                                                                                            TGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGG                                                                                                                                                                   >  2:556777/1‑77                                                                                                                                                                                                                                                                                                                                                                                                                                                                                           ACTATACCAGAGAATGAACATACTGGATTGCTCCAGTATTGTGAATGAACGGTAACGCAATAAAACATATTTACTGCAATTCTTTGGCCGGTTCTTTTACGTACAGCGGAAACCTGCCGCTTAAACGGAGAGTATCGTCGATAAAAATCCAATAAAACGTCAGGGCAAAAGTAAGAAACAGACAAAGCAAAGGCCGCTCAGGATATAGCCAGATAAATGACGGGGATCAATTGGCt‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4161087‑4160853 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCG  <  NC_000913/4499511‑4499274

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