New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1649883	35 (1.040)	5 (0.150)	5/260	3.4	11.8%	coding (90/192 nt)	ydfD	Qin prophage; uncharacterized protein
	?	NC_000913	= 1649893	41 (1.260)					coding (100/192 nt)	ydfD	Qin prophage; uncharacterized protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCA  >  NC_000913/1649737‑1650027                                                                                                                                                   |                                                                                                                                                 aataTGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGa                                                                                                                                                <  5:78680/149‑1 (MQ=255)       cATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTa                                                                                                                                           <  5:102922/148‑1 (MQ=255)                  tGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCaa                                                                                                                               >  5:152317/1‑149 (MQ=255)                    cTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAAc                                                                                                                              <  1:40547/148‑1 (MQ=255)                         tCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGaa                                                                                                                          >  3:100639/1‑147 (MQ=255)                            tttACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGa                                                                                                                      >  6:72737/1‑148 (MQ=255)                               aCGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGaaaa                                                                                                                   <  4:100639/148‑1 (MQ=255)                                   cTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTccc                                                                                                              <  8:123861/149‑1 (MQ=255)                                      aTGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCg                                                                                                             <  5:72737/147‑1 (MQ=255)                                         ccAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAAc                                                                                                        >  5:24774/1‑149 (MQ=255)                                           aaaaGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTg                                                                                                      >  8:12333/1‑149 (MQ=255)                                               ggATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTAc                                                                                                  >  4:36038/1‑149 (MQ=255)                                               ggATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTAc                                                                                                  >  4:33957/1‑149 (MQ=255)                                                        caATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGaa                                                                                                                          >  4:96976/1‑116 (MQ=255)                                                        caATGAATTCAGCAGTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGaa                                                                                                                          <  3:96976/116‑1 (MQ=255)                                                               ttCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTa                                                                                  >  5:125180/1‑149 (MQ=255)                                                                    cATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCAc                                                                             <  6:152317/149‑1 (MQ=255)                                                                            cTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGataa                                                                     <  6:24774/149‑1 (MQ=255)                                                                                  cTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGATAATATCg                                                                >  1:20343/1‑148 (MQ=255)                                                                                     aCAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGa                                                                                          >  3:54967/1‑119 (MQ=255)                                                                                     aCAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGa                                                                                          <  4:54967/119‑1 (MQ=255)                                                                                                 gTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGt                                                  <  6:125180/147‑1 (MQ=255)                                                                                                 gTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTc                                                 <  3:36038/148‑1 (MQ=255)                                                                                                 gTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTc                                                 <  3:33957/148‑1 (MQ=255)                                                                                                              cAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTcc                                   >  1:97115/1‑149 (MQ=255)                                                                                                                 ttGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCAc                                                                             <  1:65263/104‑1 (MQ=255)                                                                                                                 ttGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCAc                                                                             >  2:65263/1‑104 (MQ=255)                                                                                                                             tcagtGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCgg                    <  2:20343/149‑1 (MQ=255)                                                                                                                               agtGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGATAATATCGaa                                                              >  3:57839/1‑105 (MQ=255)                                                                                                                               agtGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGATAATATCGaa                                                              <  4:57839/105‑1 (MQ=255)                                                                                                                                          gACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTtata       >  1:152130/1‑149 (MQ=255)                                                                                                                                             aaTGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGc    >  6:108430/1‑149 (MQ=255)                                                                                                                                               tGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCa  >  6:85219/1‑149 (MQ=255)                                                                                                                                                   ggAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCt           <  2:1750/136‑1 (MQ=255)                                                                                                                                                   ggAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCt           >  1:1750/1‑136 (MQ=255)                                                                                                                                                   |                                                                                                                                                 AATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAATATCCGGTTTCATTCTTATATGCCA  >  NC_000913/1649737‑1650027

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