New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP009273	2871349 =	35 (1.700)	3 (0.160)	3/260	NT	8.6%	intergenic (+372/+579)	iap/ygbF	aminopeptidase in alkaline phosphatase isozyme conversion/CRISPR adaptation ssRNA endonuclease
	?	CP009273	2871404 =	31 (1.610)					intergenic (+427/+524)	iap/ygbF	aminopeptidase in alkaline phosphatase isozyme conversion/CRISPR adaptation ssRNA endonuclease

GCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTA  >  CP009273/2871200‑2871493                                                                                                                                                      |                                                                                                                                                 accgcgaatcgcgcataccctgcgcgtcgccgcctgcgagttccccgcgccagcggggataaaccgtcagctttataaatccggagatacggaaactagagttccccgcgccagcggggataaaccggactcaccccgaaagagattgcc                                                                                                                                                  <  2:283187‑M1/1‑1 (MQ=255)  ccgcgaaccgcgcataccctgcgcgtcgccgcctgcgagttccccgcgccagcggggataaaccgtcagctttataaatccggagctacggaaactagagttccccgcgccagcggggataaaccggactcaccccgaaagagattgcca                                                                                                                                                 >  2:158465‑M1/149‑150 (MQ=255)   cgcgactggggcgtacccggcgcgtcgccgcctgcgagtcccccgccccagccgggataaaccgtcagctttataaatccggagacacggaaactagagttccccgcgccagcggggataaaccggactcaccccgaaagagattgccag                                                                                                                                                >  1:69611‑M1/148‑150 (MQ=255)   cgcgaatcgcggctgcccggcgcggcgccgcctgcgaggtccccgcgccagcggggataaaccgtcagctttataaacccggagatacgggaactagagttcccggcgccagcggggataaaccggactcaccccgaaagagattgccag                                                                                                                                                >  1:69622‑M1/148‑150 (MQ=255)    ggtgagtCGTAGTCCATCACTCCATCTATGTTGGGACTCCCGGTTTATCCCCGCTGGCGCGGGGAGCTCCCGGGGGATAACGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGcagc                                                                                                                                               <  1:298274/144‑1 (MQ=255)     ggAACTCGTAGTCCATCATTCCACCTATGTCTGAGCTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAACGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCCTGCCTTGcagcc                                                                                                                                              <  1:185344/150‑1 (MQ=255)      gagtcgcgcatacccggcgcgtcgcgccctgcgagttccccgcgccagcggggataaaccgtcagctttataaatccggagatacggaaactagagttccccgcgccagcggggataaaccggactcaccccgaaagagattgccagcca                                                                                                                                             >  1:130918‑M1/145‑150 (MQ=255)      gAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGcagcca                                                                                                                                             >  2:158225/1‑150 (MQ=255)        aCTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGcagccagc                                                                                                                                           >  2:4089/1‑150 (MQ=255)           cGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTcc                                                                                                                                        >  1:104904/1‑150 (MQ=255)            gTAGTCCATCATTCCACCTATGTCTGACCTCCCGGTTTATCCCCGCTGGCGCGGGGAGGGTTCGGGGGATAATGTTGACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCca                                                                                                                                       <  2:83625/150‑1 (MQ=255)                taccctgcgcgtcgccgcctgcgagttccccgcgccagcggggataaaccgtcagctttataaatccggagatacggaaactagagttccccgcgccagcggggataaaccggactcaccccgaaagagattgcCagccagctgtagttc                                                                                                                                   <  1:158225‑M1/16‑8 (MQ=255)                 ccATCATTCCACCTATGTCTGAACTCCCGGTTAATCCCCGCTGGCGCGGGGCACTCCCGGGGGATAAGGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGCGGAAATCTGGGCGCCTTGCCTTGcagcccgccgcacgagc                                                                                                                                  >  1:127968/1‑141 (MQ=255)                    tcatTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCGTGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCACCGGt                                                                                                                               >  1:73462/1‑150 (MQ=255)                          cACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATcc                                                                                                                         <  2:201016/150‑1 (MQ=255)                           aCCTGTGCCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATccc                                                                                                                        <  1:208166/150‑1 (MQ=255)                                         tCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGcgccagcgggga                                                                                                          <  2:29001/150‑12 (MQ=255)                                                  ccagcggggataaaccgtcagctttataaatccggagatacggaaacgagagttccccgcgccagcggggataaaccggacgcaccccgaaagagattgcCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAgct                                                                                                 >  2:353047‑M1/101‑150 (MQ=255)                                                         ggataaaccgtcagctttataaacccggcgatacggaaactagagttccccgcgccagcggggataaaccggactcaccccgaaagagatgtcCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCtggctggc                                                                                          >  2:39492‑M1/94‑150 (MQ=255)                                                              gcgGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAAtct                                                                                     <  2:143790/150‑1 (MQ=255)                                                                      cTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCTAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCggg                                                                             <  1:177350/150‑1 (MQ=255)                                                                      cTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCggg                                                                             >  1:358002/1‑150 (MQ=255)                                                                      cTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCggg                                                                             >  2:293494/1‑150 (MQ=255)                                                                      cTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCGCTTTCggg                                                                             >  2:119331/1‑150 (MQ=255)                                                                       tCCCGGGGGATACTGTTTACGGGCCTGCGCCGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCgggg                                                                            <  1:336011/150‑1 (MQ=255)                                                                           gggggATAATGTTTACGGTCATGCGTCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTATGCGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAg                                                                        <  2:83961/150‑1 (MQ=255)                                                                           gggggATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAg                                                                        >  2:283028/1‑150 (MQ=255)                                                                                    tGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTa                                                               >  1:235645/1‑150 (MQ=255)                                                                                    gggtttcGGTCATGCGGCCCCCGGTTTGTCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCGGCCAGCTCCAGCAGCGGGTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTa                                                               <  1:85820/144‑1 (MQ=255)                                                                                               aTGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCGCGCTGgc                                                    >  1:258628/1‑150 (MQ=255)                                                                                                     gttccccgcgccagcggggataaaccggactcacccagaaagagattgcCAGCCCGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCCTTCGGGGGGAGTCCGGTTTCTCCCCGCTGGCGCgggg                                              <  1:222559‑M1/101‑1 (MQ=255)                                                                                                      ccccGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGa                                             <  2:244964/150‑1 (MQ=255)                                                                                                      ccccGGTTCATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGGAGCGGTTTATTCCCGCGGGCGCGGGGACCTCCAGCTCGGTGGCAAACTCCTACGGGGTGAGGCCGGGTTATCCCCGccgcgacggggc                                             >  1:75792/1‑139 (MQ=255)                                                                                                       cccGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGaa                                            >  1:156930/1‑150 (MQ=255)                                                                                                              tATCCCCGCGGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTAGt                                     <  1:341762/150‑1 (MQ=255)                                                                                                              tATCCCCCCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCTGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTAGt                                     <  1:329965/150‑1 (MQ=255)                                                                                                                 ccccGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTAGTTTc                                  >  1:61263/1‑150 (MQ=255)                                                                                                                 ccccGCGGGCGCGGGGAACTCTGGGCGGCTGGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTAGTTTc                                  <  2:127968/150‑1 (MQ=255)                                                                                                                    cgggggcggcgggACCTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCGCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTAGTTTCCGt                               <  2:185925/140‑1 (MQ=255)                                                                                                                    cGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTAGTTTCCGt                               <  1:283028/150‑1 (MQ=255)                                                                                                                              gggAACGCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGAt                     <  2:190228/150‑1 (MQ=255)                                                                                                                                gAACTCTGGGCGGCTTACCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTATATCCCCGATGGCGCGGGGAACCCTCGTTTCCGTATCGCCGGAttt                   >  2:214210/1‑150 (MQ=255)                                                                                                                                       gggCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGc            >  1:283187/1‑150 (MQ=255)                                                                                                                                         gcggcTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTACCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTGTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCtt          <  1:4089/150‑2 (MQ=255)                                                                                                                                             cttgccttgcAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTTAGTCCGGTTTATCCCCGCGGGCGCGGGGAACTCTAGTTTCCGTAGCTCCGGATTTATAAAGCTGACgg      <  2:104904/150‑1 (MQ=255)                                                                                                                                             cttcccttgcAGCCAGCTCCAGCAGCGCTTTATCCCCGCTGGCGCGGCGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGACTCCGTTTTATCCCCAATGGCGTGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACgg      <  2:140598/150‑1 (MQ=255)                                                                                                                                                 ccttgcAGCCAGCTCCAGCAGCGGTTTATCAACTCAGCGGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGCTTTa  >  1:346256/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                 GCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTA  >  CP009273/2871200‑2871493

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