New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4160853 =99 (1.410)5 (0.070) 5/436 NT 4.8% 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.

TGCTCCAGTATTGTGAATGAACGGTAACGCAATAAAACATATTTACTGCAATTCTTTGGCCGGTTCTTTTACGTACAGCGGAAACCTGCCGCTTAAACGGAGAGTATCGTCGATAAAAATCCAATAAAACGTCAGGGCAAAAGTAAGAAACAGACAAAGCAAAGGCCGCTCAGGATATAGCCAGATAAATGACGGGGATCAATTGGCt‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4161059‑4160853
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCG  <  NC_000913/4499511‑4499274
                                                                                                                                                                                                                                                                                                                                                                                                                                                              
TGCTCCAGTATTGTGAATGAACGGTAACGCAATAAAACATATTTACTGCAATTCTTTGGCCGGTTCTTTTACGTACAGCGGAAACCTGCCGCTTAAACGGAGAGTATCGTCGATAAAAATCCAATAAAACGTCAGGGCAAAAGTAAGCAACAGACAAAGCAAAGGCCGCTCAGGAGATAGCCAGATAAATTATGGGGAACAATTGGCTGGATTTGCCCCTATATTTCCAGACCTCTGTT                                                                                                                                                                                                                >  1:857890/1‑239
TGCTCCAGTATTGTGAATGAACGGTAACGCAATAAAACATATTTACTGCAATTCTTTGGCCGGTTCTTTTACGTACAGCGGAAACCTGCCGCTTAAACGGAGAGTATCGTCGATAAAAATCCAATAAAACGTCAGGGCAAAAGTAAGAAACAGACAAAGCAAAGGCCGCTCAGGATATAGCCAGATAAATGACGGGGATCAATTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTT                                                                                                                                                                                                                <  2:305036/239‑1
                                                        TGGCCGGTTCTTTTACGTACAGCGGAAACCTGCCGCTTAAACGGAGAGTATCGTCGATAAAAATCCAATAAAACGTCAGGGCAAAAGTAAGAAACAGACAAAGCAAAGGCCGCTCAGGATATAGCCAGATAAATGACGGGGATCAATTGGCTGGATTTGCCCCTATATTTCCAGACATCTGGTATCACTTAACCCATTACAAGCCCGCGGCCCCAAATATTCCCTTGGCGAGCGAAAA                                                                                                                                                         >  2:77294/1‑238
                                                                ttTTTTAAGTACCAGCGAAACCTGTCCTTTAAACGTAAGATTTCCTTGATAAAAATCCAATAAAACGTCAGGGCAAAAGTAAAAAACAGACAAAGCAAAGGCCGCTCAGGATATAGCCAGATAAATGACGGGAATAAATTGGCTGGATTTGCCCCTATATTTCCATACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCA                                                                                                                                                <  2:857890/237‑1
                                                                                                 ACGGAGAGTATCGTCGATAAAAATCCAATAAAACGTCAGGGCAAAAGTAAGAAACAGACAAAGCAAAGGCCGCTCAGGATATAGCCAGATAAATGACGGGGATCAATTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAA                                                                                                                <  1:77294/239‑1
                                                                                                                                                                                                aatccccctccTTGGTTGGTTTTGCCCCTATATTCCCAGACCTCTTTTATCACTTAACCCATTACAAGCCCGCTGCCGAAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTA                 <  2:137099/228‑1
                                                                                                                                                                                                            TGGCTGGATTTGCCCCTATATTTCCAGACCTCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTAGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCT     <  1:201041/239‑1
                                                                                                                                                                                                               CTGGATTTTCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCC   <  1:468091/238‑1
                                                                                                                                                                                                                TGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAACGTCTTTGCTATTCCGTTACCCTCCG  >  1:404121/1‑238
                                                                                                                                                                                                                                                                                                                                                                                                                                                              
TGCTCCAGTATTGTGAATGAACGGTAACGCAATAAAACATATTTACTGCAATTCTTTGGCCGGTTCTTTTACGTACAGCGGAAACCTGCCGCTTAAACGGAGAGTATCGTCGATAAAAATCCAATAAAACGTCAGGGCAAAAGTAAGAAACAGACAAAGCAAAGGCCGCTCAGGATATAGCCAGATAAATGACGGGGATCAATTGGCt‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4161059‑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 < 27 ≤ 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.