New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1647054 =	22 (0.840)	3 (0.120)	3/276	2.8	14.2%	intergenic (+317/+120)	flxA/intK	Qin prophage; uncharacterized protein/pseudogene, integrase fragment, Qin prophage;Phage or Prophage Related
	?	NC_000913	1647120 =	16 (0.660)					intergenic (+383/+54)	flxA/intK	Qin prophage; uncharacterized protein/pseudogene, integrase fragment, Qin prophage;Phage or Prophage Related
Rejected: Frequency below cutoff threshold.

CAATTATGAATTGCAACTATTTCTGTAGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGG  >  NC_000913/1646909‑1647197                                                                                                                                                  |                                                                                                                                                tccctggcttacgcgccgaacgcaaagtcatcagcaagttaagcccccgaagttacctaaaaagaagcctgatcccgataaaaaagattgaaaccaagatcgattcggttgagtgcatatccattcatagggtagattcttaagtCGCGtt                                                                                                                                            <  1:142869‑M1/6‑1 (MQ=255) tccctggcttacgcgccgaacgcaaagtcatcagcaagttaagcccccgaagttacctaaaaagaagcctgatcccgataaaaaagattgaaaccaagatcgattcggttgagtgcatatccattcatagggtagattcttaagtCGCGtt                                                                                                                                            <  2:139403‑M1/6‑1 (MQ=255)   cctggcttacgcgccgaacgcaaagtcatcagcaagttaagcccccgaagttacctaaaaagaagcctgatcccgataaaaaagattgaaaccaagatcgattcggttgagtgcatatccattcatagggtagattcttaagtCGCGTTTc                                                                                                                                          >  1:139370‑M1/144‑151 (MQ=255)           ttGCAACTATTTCTGTAGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCatat                                                                                                                                  >  1:410192/1‑151 (MQ=255)                              cACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAAc                                                                                                               <  2:337846/151‑1 (MQ=255)                                cTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACaa                                                                                                             >  1:13951/1‑151 (MQ=255)                                   ttGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAacc                                                                                                          <  2:394490/151‑1 (MQ=255)                                    tgtgGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAACCAAGTAACAAacca                                                                                                         >  2:357977/1‑151 (MQ=255)                                    tgtgGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAacca                                                                                                         >  1:361319/1‑151 (MQ=255)                                        tcgACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACTCGCACCGCCCAGACTGACTCAGCGCGTTTCTTGCATATCCCCGGTAAAACAAGTAACAAACCAcgcg                                                                                                     >  1:127713/3‑151 (MQ=255)                                         ggACAGTCCACAAAACTGCCAACTTCCGTTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGa                                                                                                    >  2:67260/1‑151 (MQ=255)                                              gTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATg                                                                                               >  2:392587/1‑151 (MQ=255)                                                         tGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGaaa                                                                                    <  1:312579/151‑1 (MQ=255)                                                            cAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAAcgc                                                                                 <  2:99175/151‑1 (MQ=255)                                                               cTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGa                                                                               >  1:28982/1‑150 (MQ=255)                                                                                tAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTAt                                                             >  2:174597/1‑151 (MQ=255)                                                                                       cATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGa                                                      <  1:114712/151‑1 (MQ=255)                                                                                           aGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATg                                                  >  1:412492/1‑151 (MQ=255)                                                                                           aGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGAAATg                                                  >  2:177387/1‑151 (MQ=255)                                                                                             cATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATg                                                  <  1:83742/149‑1 (MQ=255)                                                                                              aTAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCCCCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCAc                                               <  2:318263/151‑1 (MQ=255)                                                                                               tAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACt                                              <  2:61767/151‑1 (MQ=255)                                                                                                   cTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAAc                                          <  2:410295/151‑1 (MQ=255)                                                                                                       attcggttgagtgcatatccattcatagggtagattcttaagtCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGaa                                      <  1:317220‑M1/108‑1 (MQ=255)                                                                                                            aTAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGAt                                 >  2:255077/1‑151 (MQ=255)                                                                                                                 gtgcatatccattcatagggtagattcttaagtCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATTAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTgg                            <  1:157974‑M1/118‑1 (MQ=255)                                                                                                                       atccattcatagggtagattcttaagtCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCaa                      <  1:410852‑M1/124‑1 (MQ=255)                                                                                                                                           gactCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAgg  >  1:72152/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                CAATTATGAATTGCAACTATTTCTGTAGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGG  >  NC_000913/1646909‑1647197

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