New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	4561348 =	30 (1.300)	4 (0.180)	4/234	NT	13.8%	coding (245/1404 nt)	PP_4048	hypothetical protein
	?	NC_002947	4561414 =	22 (1.020)					coding (179/1404 nt)	PP_4048	hypothetical protein

CACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGCACAATCACCACCCCCGGCACCAACGACAGCAGCATGTCGGCGATGACCACGCTCAGT  >  NC_002947/4561199‑4561488                                                                                                                                                      |                                                                                                                                             gcggtggaaaaagtggccgacggctttggctggctgcgtgacctggtgcttggcgagtttgccgaggagcggccactgagcgtggtcatcgccgacatgctgctgtcgttggtgccgggggtggtgattgtgctcagtgcccgtgacctg                                                                                                                                              >  1:363123‑M1/150‑150 (MQ=255)                 tgctTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTg                                                                                                                              >  1:549351/1‑150 (MQ=255)                 tgctTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTg                                                                                                                              <  2:99514/150‑1 (MQ=255)                  gctTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTgg                                                                                                                             <  1:56681/150‑1 (MQ=255)                           ccAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTcc                                                                                                                    <  2:554900/150‑1 (MQ=255)                                         tGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCtt                                                                                                      >  1:590090/1‑150 (MQ=255)                                                 aGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGgc                                                                                              >  1:355971/1‑150 (MQ=255)                                                       gagtttgccgaggagcggccactgagcgtggtcatcgccgacatgctgctgtcgttggtgccgggggtggtgattgtgctcagtgcccgtgacctGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAg                                                                                          <  2:311014‑M1/53‑1 (MQ=255)                                                       gagtttgccgaggagcggccactgagcgtggtcatcgccgacatgctgctgtcgttggtgccgggggtggtgattgtgctcagtgcccgtgacctGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAg                                                                                          >  1:311014‑M1/96‑148 (MQ=255)                                                          tttgccgaggagcggccactgagcgtggtcatcgccgacatgctgctgtcgttggtgccgggggtggtgattgtgctcagtgcccgtgacctGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGcaca                                                                                     <  2:486662‑M1/58‑1 (MQ=255)                                                           ttgccgaggagcggccactgagcgtggtcatcgccgacatgctgctgtcgttggtgccgggggtggtgattgtgctcagtgcccgtgacctGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCacac                                                                                    >  1:486662‑M1/92‑150 (MQ=255)                                                                ccACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCACCGAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAg                                                                               >  1:15852/1‑150 (MQ=255)                                                                ccACCGCCCCCACAACAGCGCCGATGATCGGCGCGACCAGAGTCATCAGGCAGGCCCCCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGACTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACCCCGcgg                                                                               >  2:85310/1‑148 (MQ=255)                                                                    cGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCCCCCTACCAGCTGGGTCGGCTGCGTCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAg                                                                           <  1:75371/150‑1 (MQ=255)                                                                    cGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAg                                                                           <  1:274713/150‑1 (MQ=255)                                                                            ccacAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACggg                                                                   <  2:355971/150‑1 (MQ=255)                                                                                 gcgcCGATCATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCt                                                                                                                               >  1:326584/1‑85 (MQ=255)                                                                                 gcgcCGATCATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCt                                                                                                                               <  2:326584/85‑1 (MQ=255)                                                                                  cgcCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGa                                                             <  1:288565/150‑1 (MQ=255)                                                                                    ccGATGATCGGCGCGATCAGGGTCATCAGGCAGGCAACCAGCAACAACCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGc                                                           >  1:321018/1‑150 (MQ=255)                                                                                             ggCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGCACAATcacc                                                  >  1:362946/1‑150 (MQ=255)                                                                                             ggCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGCACAATcacc                                                  >  2:228208/1‑150 (MQ=255)                                                                                                      aGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTccgcc                                                                                                                 <  1:89455/78‑1 (MQ=255)                                                                                                      aGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTccgcc                                                                                                                 >  2:89455/1‑78 (MQ=255)                                                                                                       gggTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGCACAATCACCACCCCCGGCa                                        <  1:258189/150‑1 (MQ=255)                                                                                                              caggcaggCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGCACAATCACCACCCCCGGCACCAACGa                                 <  2:601639/150‑1 (MQ=255)                                                                                                               aggcaggcCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGCACAATCACCACCCCCGGCACCAACGAc                                >  2:235563/1‑150 (MQ=255)                                                                                                                     gcCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGCAAAATCACCCCCCCCGGCACCAACGACagcagc                          <  2:590090/150‑1 (MQ=255)                                                                                                                       caccaGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTa                                                                                                            >  2:403812/1‑66 (MQ=255)                                                                                                                       caccaGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTa                                                                                                            <  1:403812/66‑1 (MQ=255)                                                                                                                        accaGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGCACAATCACCACCCCCGGCACCAACGACAGCAGCATg                       >  2:608260/1‑150 (MQ=255)                                                                                                                            gCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGCACAATCACCACCCCCGGCACCAACGACAGCAGCATGTCgg                   >  2:569414/1‑150 (MQ=255)                                                                                                                                   ccACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGCACAATCACCACCCCCGGCACCAACGACAGCAGCATGTCGGCGATGAc            >  2:313494/1‑150 (MQ=255)                                                                                                                                         ttGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGCACAATCACCACCCCCGGCACCAACGACAGCAGCATGTCGGCGATGACCACGCt      >  1:96796/1‑150 (MQ=255)                                                                                                                                             cATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGCACAATCACCACCCCCGGCACCAACGACAGCAGCATGTCGGCGATGACCACGCTCAGt  <  2:321018/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                             CACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGCACAATCACCACCCCCGGCACCAACGACAGCAGCATGTCGGCGATGACCACGCTCAGT  >  NC_002947/4561199‑4561488

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