New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3521815 =	34 (0.940)	4 (0.120)	4/236	3.5	11.0%	coding (158/540 nt)	hofN	DNA catabolic putative fimbrial assembly protein
	?	NC_000913	3521857 =	33 (0.990)					coding (116/540 nt)	hofN	DNA catabolic putative fimbrial assembly protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GTTGCGAACGCTGTTCGCGTAACTGCTGTTGCTCCAGCAAACGTGGCTTCGTTATCTGTAAGCTGCGGGCGAGTTGTTGTTCCGCCTGAAGCAAAACTGCGTCTATGCGCGCTTCGGCGCTGCCTGTCAGACGCAGTATTAGCGTTATCCCGACGGCCAGCAGCAGAGGCGCAACGAACATCAGCAACCAGAAACGCAGAAAAGCGGTCCGGCGTTGCTGTCGCCAGGGCAAAAAATTAATTGGCGGGTTCATCAGTACGCCTTCCCAAGCGC  >  NC_000913/3521720‑3521992                                                                                                                                          |                                                                                                                                        gTTGCGAACGCTGTTCGCGTAACTGCTGTTGCTCCAGCAAACGTGGCTTCGTTATCTGTAAGCTGCGGGCGAGTTGTTGTTCCGCCTGAAGCAAAACTGCGTCTATGCGCGCTTCGGCGCTGCCTGTCAGACGCAGTATTAGCGTTATc                                                                                                                              >  6:116486/1‑149 (MQ=255)       aaCGCTGTTCGCGTAACTGCTGTTGCTCCAGCAAACGTGGCTTCGTTATCTGTAAGCTGCGGGCGAGTTGTTGTTCCGCCTGAAGCAAAACTGCGTCTATGCGCGCTTCGGCTCTGCCTGTCAGACGCAGTATTAGCGTTATCCCGACg                                                                                                                        <  5:116486/149‑1 (MQ=255)            tGTTCGCGTAACTGCTGTTGCTCCAGCAAACGTGGCTTCGTTATCTGTAAGCTGCGGGCGAGTTGTTGTTCCGCCTGAAGCAAAACTGCGTCTATGCGCGCTTCGGCGCTGCCTGTCAGACGCAGTATTAGCGTTATCCCGACGGc                                                                                                                      >  7:57982/1‑146 (MQ=255)              ttCGCGTAACTGCTGTTGCTCCAGCAAACGTGGCTTCGTTATCTGTAAGCTGCGGGCGAGTTGTTGTTCCGCCTGAAGCAAAACTGCGTCTATGCGCGCTTCGGCGCTGCCTGTCAGACGCAGTATTAGCGTTATCCCGACGGCca                                                                                                                    <  8:57982/146‑1 (MQ=255)                             ttGCTCCAGCAAACGTGGCTTCGTTATCTGTAAGCTGCGGGCGAGTTGTTGTTCCGCCTGAAGCAAAACTGCGTCTATGCGCGCTTCGGCGCTGCCTGTCAGACGCAGTATTAGCGTTATCCCGACGGCca                                                                                                                    >  1:149483/1‑131 (MQ=255)                             ttGCTCCAGCAAACGTGGCTTCGTTATCTGTAAGCTGCGGGCGAGTTGTTGTTCCGCCTGAAGCAAAACTGCGTCTATGCGCGCTTCGGCGCTGCCTGTCAGACGCAGTATTAGCGTTATCCCGACGGCca                                                                                                                    <  2:149483/131‑1 (MQ=255)                                       aaaCGTGGCTTCGTTATCTGTAAGCTGCGGGCGAGTTGTTGTTCCGCCTGAAGCAAAACTGCGTCTATGCGCGCTTCGGCGCTGCCTGTCAGACGCAGTATTAGCGTTATCCCGAc                                                                                                                         <  5:111098/116‑1 (MQ=255)                                       aaaCGTGGCTTCGTTATCTGTAAGCTGCGGGCGAGTTGTTGTTCCGCCTGAAGCAAAACTGCGTCTATGCGCGCTTCGGCGCTGCCTGTCAGACGCAGTATTAGCGTTATCCCGAc                                                                                                                         >  6:111098/1‑116 (MQ=255)                                                         tGTAAGCTGCGGGCGAGTTGTTGTTCCGCCTGAAGCAAAACTGCGTCTATGCGCGCTTCGGCGCTGCCTGTCAGACGCAGTATTAGCGTTATCCCGACGGCCAGCAGCAGAGGCGCAACGAACATCAGCAACCAGAAACGCAGAAAAGc                                                                      >  2:98039/1‑149 (MQ=255)                                                          gTAAGCTGCGGGCGAGTTGTTGTTCCGCCTGAAGCAAAACTGCGTCTATGCGCGCTTCGGCGCTGCCTGTCAGACGCAGTATTAGCGTTATCCCGACGGCCAGCAGCAGAGgc                                                                                                         <  2:186568/113‑1 (MQ=255)                                                          gTAAGCTGCGGGCGAGTTGTTGTTCCGCCTGAAGCAAAACTGCGTCTATGCGCGCTTCGGCGCTGCCTGTCAGACGCAGTATTAGCGTTATCCCGACGGCCAGCAGCAGAGgc                                                                                                         >  1:186568/1‑113 (MQ=255)                                                          gTAAGCTGCGGGCGAGTTGTTGTTCCGCCTGAAGCAAAACTGCGTCTATGCGCGCTTCGGCGCTGCCTGTCAGACGCAGTATTAGCGTTATCCCGACGGCCAGCAGCAGAGGCGCAACGAACATCAGCAACCAGAAACGCAGAAAAGCg                                                                     >  2:69405/1‑149 (MQ=255)                                                           tAAGCTGCGGGCGAGTTGTTGTTCCGCCTGAAGCAAAACTGCGTCTATGCGCGCTTCGGCGCTGCCTGTCAGACGCAGTATTAGCGTTATCCCGACGGCCAGCAGCAGAGGCGCAACGAACATCAGCAACCAGAAACGCAGAAAAGc                                                                      >  2:77706/1‑147 (MQ=255)                                                           tAAGCTGCGGGCGAGTTGTTGTTCCGCCTGAAGCAAAACTGCGTCTATGCGCGCTTCGGCGCTGCCTGTCAGACGCAGTATTAGCGTTATCCCGACGGCCAGCAGCAGAGGCGCAACGAACATCAGCAACCAGAAACGCAGAAAAGc                                                                      <  1:77706/147‑1 (MQ=255)                                                           tAAGCTGCGGGCGAGTTGTTGTTCCGCCTGAAGCAAAACTGCGTCTATGCGCGCTTCGGCGCTGCCTGTCAGACGCAGTATTAGCGTTATCCCGACGGCCAGCAGCAGAGGCGCAACGAACATCAGCAACCAGAAACGCAGAAAAGCgg                                                                    >  7:76485/1‑149 (MQ=255)                                                                tGCGGGCGAGTTGTTGTTCCGCCTGAAGCAAAACTGCGTCTATGCGCGCTTCGGCGCTGCCTGTCAGACGCAGTATTAGCGTTATCCCGACGGCCAGCAGCAGAGGCGCAACGAACATCAGCAAc                                                                                       >  8:37410/1‑125 (MQ=255)                                                                tGCGGGCGAGTTGTTGTTCCGCCTGAAGCAAAACTGCGTCTATGCGCGCTTCGGCGCTGCCTGTCAGACGCAGTATTAGCGTTATCCCGACGGCCAGCAGCAGAGGCGCAACGAACATCAGCAAc                                                                                       <  7:37410/125‑1 (MQ=255)                                                                    ggCGAGTTGTTGTTCCGCCTGAAGCAAAACTGCGTCTATGCGCGCTTCGGCGCTGCCTGTCAGACGCAGTATTAGCGTTATCCCGACGGCCAGCAGCAGAGGCGCAACGAACATCAGCAACCAGAAACGCAGAAAAGCGGTCCGGCGtt                                                           <  1:98039/149‑1 (MQ=255)                                                                                      cTGAAGCAAAACTGCGTCTATGCGCGCTTCGGCGCTGCCTGTCAGACGCAGTATTAGCGTTATCCCGACGGCCAGCAGCAGAGGCGCAACGAACATCAGCAACCAGAAACGCAGAAAAGCGGTcc                                                                 <  8:20/125‑1 (MQ=255)                                                                                      cTGAAGCAAAACTGCGTCTATGCGCGCTTCGGCGCTGCCTGTCAGACGCAGTATTAGCGTTATCCCGACGGCCAGCAGCAGAGGCGCAACGAACATCAGCAACCAGAAACGCAGAAAAGCGGTcc                                                                 >  7:20/1‑125 (MQ=255)                                                                                       tGAAGCAAAACTGCGTCTATGCGCGCTTCGGCGCTGCCTGTCAGACGCAGTATTAGCGTTATCCCGACGGCCAGCAGCAGAGGCGCAACGAACAt                                                                                              <  5:47298/95‑1 (MQ=255)                                                                                       tGAAGCAAAACTGCGTCTATGCGCGCTTCGGCGCTGCCTGTCAGACGCAGTATTAGCGTTATCCCGACGGCCAGCAGCAGAGGCGCAACGAACAt                                                                                              >  6:47298/1‑95 (MQ=255)                                                                                           gCAAAACTGCGTCTATGCGCGCTTCGGCGCTGCCTGTCAGACGCAGTATTAGCGTTATc                                                                                                                              <  2:181528/59‑1 (MQ=255)                                                                                           gCAAAACTGCGTCTATGCGCGCTTCGGCGCTGCCTGTCAGACGCAGTATTAGCGTTATc                                                                                                                              >  1:181528/1‑59 (MQ=255)                                                                                                 cTGCGTCTATGCGCGCTTCGGCGCTGCCTGTCAGACGCAGTATTAGCGTTATCCCGACGGCCAGCAGCAGAGGCGCAACGAACATCAGCAACCAGAAACGCAGAAAAGCGGTCCGGCGTTGCTGTCGCCAGGGCAAAAAATTAATTGGc                              <  1:69405/149‑1 (MQ=255)                                                                                                         aTGCGCGCTTCGGCGCTGCCTGTCAGACGCAGTATTAGCGTTATCCCGACGGCCAGCAGCaga                                                                                                            >  6:116811/1‑63 (MQ=255)                                                                                                         aTGCGCGCTTCGGCGCTGCCTGTCAGACGCAGTATTAGCGTTATCCCGACGGCCAGCAGCaga                                                                                                            <  5:116811/63‑1 (MQ=255)                                                                                                                           ccTGTCAGACGCAGTATTAGCGTTATCCCGACGGCCAGCAGCAGAGGCGCAACGAACATCa                                                                                            >  4:139694/1‑61 (MQ=255)                                                                                                                           ccTGTCAGACGCAGTATTAGCGTTATCCCGACGGCCAGCAGCAGAGGCGCAACGAACATCa                                                                                            <  3:139694/61‑1 (MQ=255)                                                                                                                             tGTCAGACGCAGTATTAGCGTTATCCCGACGGCCAGCAGCAGAGGCGCAACGAACATCAGCAACCAGAAACGCAGAAAAGCGGTCCGGCGTTGCTGTCGCCAGGGCAAAAAATTAATTGGCGGGTTCATCAGTa                 >  7:44910/1‑134 (MQ=255)                                                                                                                             tGTCAGACGCAGTATTAGCGTTATCCCGACGGCCAGCAGCAGAGGCGCAACGAACATCAGCAACCAGAAACGCAGAAAAGCGGTCCGGCGTTGCTGTCGCCAGGGCAAAAAATTAATTGGCGGGTTCATCAGTa                 <  8:44910/134‑1 (MQ=255)                                                                                                                             tGTCAGACGCAGTATTAGCGTTATCCCGACGGCCAGCAGCAGAGGCGCAACGAACATCAGCAACCAGAAACGCAGAAAAGCGGTCCGGCGTTGCTGTCGCCAGGGCAAAAAATTAATTGGCGGGTTCATCAGTACGCCTTCCCAAgcgc  >  1:40774/1‑149 (MQ=255)                                                                                                                              gTCAGACGCAGTATTAGCGTTATCCCGACGGCCAGCAGCAGAGGCGCAACGAACATCAGCAACCAGAAACGCAGAAAAGCGGTCCGGCGTTGCTGTCGCCAGGGCAAAAAATTAATTGGCGGGTTCATCAGTACGCCTTCCCAAgcgc  >  1:70892/1‑148 (MQ=255)                                                                                                                                          |                                                                                                                                        GTTGCGAACGCTGTTCGCGTAACTGCTGTTGCTCCAGCAAACGTGGCTTCGTTATCTGTAAGCTGCGGGCGAGTTGTTGTTCCGCCTGAAGCAAAACTGCGTCTATGCGCGCTTCGGCGCTGCCTGTCAGACGCAGTATTAGCGTTATCCCGACGGCCAGCAGCAGAGGCGCAACGAACATCAGCAACCAGAAACGCAGAAAAGCGGTCCGGCGTTGCTGTCGCCAGGGCAAAAAATTAATTGGCGGGTTCATCAGTACGCCTTCCCAAGCGC  >  NC_000913/3521720‑3521992

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.