New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4200145 =	38 (1.280)	5 (0.180)	5/276	2.0	13.5%	intergenic (+51/‑136)	yjaG/hupA	DUF416 domain protein/HU, DNA‑binding transcriptional regulator, alpha subunit
	?	NC_000913	4200182 =	29 (1.050)					intergenic (+88/‑99)	yjaG/hupA	DUF416 domain protein/HU, DNA‑binding transcriptional regulator, alpha subunit
Rejected: Frequency below cutoff threshold.

CCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGTAATTGCAGAGAAAGCAGAAC  >  NC_000913/4200034‑4200326                                                                                                                                                     |                                                                                                                                                 ccTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaaaaa                                                                                                                                                >  2:206010/1‑151 (MQ=255)      gggCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTgg                                                                                                                                           >  2:222973/1‑151 (MQ=255)      gggCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTgg                                                                                                                                           >  2:398651/1‑151 (MQ=255)           gACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTAt                                                                                                                                       >  1:164196/1‑150 (MQ=255)                     cttacaatgtgtttatcgcttgctaagcatcgagtgcgacggaaatgccagaatagcactctcctgcatacacgcaccgatagccacttttTTTCGCCCCCCAAATGTAGACCAGACGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTa                                                                                                                            >  1:330402‑M1/92‑151 (MQ=255)                     aGGCGGGGGGGAGGAAGATTGGTATAATTTTTTAGGAATAAGACCAGGAAACGTGATTTAACGCCTGATTTTTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGGAAAAAAGTTGCTATCGGTGCgggga                                                                                                                            >  1:110084/1‑147 (MQ=255)                       gCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTa                                                                                                                            <  1:306699/149‑1 (MQ=255)                        acaatgtgtttatcgcttgctaagcatcgagtgcgacggaaatgccagaatagcactctcctgcatacacgcaccgatagccacttttTTTCGCCCCCCAAATGTAGACCAGACGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGc                                                                                                                         >  2:376010‑M1/89‑151 (MQ=255)                           atgtgtttatcgcttgctaagcatcgagtgcgacggaaatgccagaatagcactctcctgcatacacgcaccgatagccacttttTTTCGCCCCCCAAATGTAGACCAGACGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAgg                                                                                                                      >  2:135807‑M1/86‑151 (MQ=255)                            tGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAgg                                                                                                                      <  1:226255/150‑1 (MQ=255)                              tgtttatcgcttgctaagcatcgagtgcgacggaaatgccagaatagcactctcctgcatacacgcaccgatagccacttttTTTCGCCCCCCAAATGTAGACCAGACGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGgag                                                                                                                    >  1:15293‑M1/83‑150 (MQ=255)                              agagCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGgag                                                                                                                    <  2:230335/150‑1 (MQ=255)                              agagCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGgag                                                                                                                    <  2:336359/150‑1 (MQ=255)                              agagCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTATTCCCTTTCGCCCCCCGTCTGGTCTGCATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGgat                                                                                                                    <  1:184218/150‑2 (MQ=255)                               gagCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGagag                                                                                                                  >  1:115932/1‑151 (MQ=255)                                agCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGagag                                                                                                                  <  1:38609/150‑1 (MQ=255)                                    atatTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTa                                                                                                             >  2:365309/1‑151 (MQ=255)                                       ttGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTc                                                                                                          <  1:469752/151‑1 (MQ=255)                                       ttGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTc                                                                                                          <  1:354316/151‑1 (MQ=255)                                            ataATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCa                                                                                                     <  1:41338/151‑1 (MQ=255)                                                      gtgcgacggaaatgccagaatagcactctcctgcatacacgcaccgatagccactttttttcgccccccgtctggtctacatttggggggcgaaaAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCg                                                                                            >  2:381265‑M2/96‑150 (MQ=255)                                                        cAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCa                                                                                          >  2:25947/1‑150 (MQ=255)                                                         cgacggaaatgccagaatagcactctcctgcatacacgcaccgatagccacttttTTTCGCCCCCCAAATGTAGACCAGACGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCAc                                                                                         >  2:381265‑M1/56‑150 (MQ=255)                                                         aaTAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCa                                                                                          >  2:271689/1‑149 (MQ=255)                                                              gACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATg                                                                                   >  2:89963/1‑151 (MQ=255)                                                              gACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATg                                                                                   >  2:179380/1‑151 (MQ=255)                                                                    aaaaCGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAg                                                                              <  2:75426/150‑1 (MQ=255)                                                                                 ctcctgcatacacgcaccgatagccactttttttcgccccccgtctggtctacatttggggggcgaaaAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAaca                                                                 >  1:15293‑M2/69‑150 (MQ=255)                                                                                 aCGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAacac                                                                <  1:255423/151‑1 (MQ=255)                                                                                   cctgcatacacgcaccgatagccactttttttcgccccccgtctggtctacatttggggggcgaaaAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACAt                                                              >  2:135807‑M2/67‑151 (MQ=255)                                                                                      gcatacacgcaccgatagccactttttttcgccccccgtctggtctacatttggggggcgaaaAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGt                                                           >  2:376010‑M2/64‑151 (MQ=255)                                                                                       gATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTa                                                          >  2:332224/1‑151 (MQ=255)                                                                                         tacacgcaccgatagccactttttttcgccccccgtctggtctacatttggggggcgaaaAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAg                                                        >  1:330402‑M2/61‑151 (MQ=255)                                                                                           tGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGa                                                      <  1:129635/151‑1 (MQ=255)                                                                                            gTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGAt                                                     >  1:121949/1‑151 (MQ=255)                                                                                                         tCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAAc                                          <  1:12337/149‑1 (MQ=255)                                                                                                                 TTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCaa                                 <  2:207886‑M1/150‑1 (MQ=255)                                                                                                                 TTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAAc                                >  2:469170‑M1/1‑151 (MQ=255)                                                                                                                        CCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACtgattga                         <  2:426795‑M1/151‑1 (MQ=255)                                                                                                                              tGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGTAAt                   <  2:177318/151‑1 (MQ=255)                                                                                                                               ggTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGTAAtt                  >  2:91755/1‑151 (MQ=255)                                                                                                                               ggTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGTAAtt                  >  1:58464/1‑151 (MQ=255)                                                                                                                                           gggggcgaaaAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGTACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGTAATTGCAGAGAAAGc       <  2:443811‑M2/140‑1 (MQ=255)                                                                                                                                                cgaaaAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGTAATTGCAGAGAAAGCAGAAc  <  1:222969‑M2/145‑1 (MQ=255)                                                                                                                                                cgaaaAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGTAATTGCAGAGAAAGCAGAAc  <  1:239918‑M2/145‑1 (MQ=255)                                                                                                                                                cgaaaAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGTAATTGCAGAGAAAGCAGAAc  <  1:371103‑M2/145‑1 (MQ=255)                                                                                                                                                cgaaaAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGTAATTGCAGAGAAAGCAGAAc  <  2:330411‑M2/145‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 CCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGTAATTGCAGAGAAAGCAGAAC  >  NC_000913/4200034‑4200326

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