New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 5328097	21 (0.910)	4 (0.180)	3/236	NT	18.3%	coding (290/1038 nt)	PP_4688	heme/hemin ABC transporter permease
	?	NC_002947	= 5328132	16 (0.730)					coding (255/1038 nt)	PP_4688	heme/hemin ABC transporter permease

CCCAAGCCCGCCGACAAACGCGCAGAACGACAGCAGATAAGGCGCGAAAGCATCCGGCATGCCGCCAAACCAGGCGCCGCCGACAATGGCAACCGCTGCCCCCACTGCGGCGCCGGCGGCAACCCCTACCAGCCCGGGGTCGGCCAGCGGGTTGCGGAACAGCCCCTGCATGGCCACCCCGGACAGCGCCAGCACCGCGCCCACCGCCAACCCCAACAGGGTGCGCGGCAGGCGGATCTGACCAAGGATCATCTCGGCCTGCTGCAGGCCATCTTCGGAAACGGGCAG  >  NC_002947/5327958‑5328245                                                                                                                                            |                                                                                                                                                     cccAAGCCCGCCGACAAACGCGCAGAACGACAGCAGATAAGGCGCGAAAGCATCCGGCATGCCGCCAAACCAGGCGCCGCCGACAATGGCAACCGCTTCCCCCACTGCGGCGCCGGCGGCAACCCCTACCACCCCGGCGCCGGCCAGCgg                                                                                                                                            >  2:202517/1‑150 (MQ=255)  ccAAGCCCGCCGACAAACGCGCAGAACGACAGCAGATAAGGCGAGACAGCATCAGGCATACTGCCAAACCAGGCACCGCCGACAATGGCAACCGCTGCCCCCACTGCGGCGCCGGCGGCAACCCCTACCAGCCCGGGGTCGGCCAGCggg                                                                                                                                           >  1:231240/1‑150 (MQ=255)             aacaacagcgcACAACGACAGAACATAAGGCGCGAAAGCATCCGGCATGCCGCCAAACCAGGCGCCGCCGACAATGGCAACCGCTGCCCCCACTGCGGTGCCGGCGGCAACCCCTACCAGCCCGGGGTCGGCCAGCGGGTTGCGGAACAg                                                                                                                                <  2:228598/143‑1 (MQ=255)                aaaCGCGCAGAACGACAGCAGATAAGGCGCGAAAGCATCCGGCATGCCGCCAAACCAGGCGCCGCCGACAATGGCAACCGCTGCCCCCACTGCGGCGCCGGCGGCAACCCCTACCAGCCCGGGGTCGGCCAGCGGGTTGCGGAACAGccc                                                                                                                             <  2:147577/150‑1 (MQ=255)                  aCGCGCAGAACGACAGCAGATAAGGCGCGAAAGCATCCGGCATGCCGCCAAACCAGGCGCCGCCGACAATGGCAACCGCTGCCCCCACTGCGGCGCCGGCGGCAACCCCTACCAGCCCGGGGACGGCCAGCGGGTTGCGGAACAGCCCCt                                                                                                                           >  1:529953/1‑150 (MQ=255)                       cAGAACGACAGCAGATAAGGCGCGAAAGCATCCGGCATGCCGCCAAACCAGGCGCCGCCGACAATGGCAACCGCTGCCCCCACTGCGGCGCCGGCGGCAACCCCTACCAGCCCGGGGTCGGCCAGCGGGTTGCGGAACAGCCCCTGCATg                                                                                                                      >  2:43201/1‑150 (MQ=255)                         gAACGACAGCAGATAAGGCGCGAAAGCATCCGGCATGCCGCCAAACCAGGCGCCGCCGACAATGGCAACCGCTGCCCCCACTGCGGCGCCGGCGGCAACCCCTACCAGCCCGGGGTCGGCCAGCGGGTTGCGGAACAGCCCCTGCATGGc                                                                                                                    <  1:43201/150‑1 (MQ=255)                         gAACGACAGCAGATAAGGCGCGAAAGCATCCGGCATGCCGCCAAACCAGGCGCCGCCGACAATGGCAACCGCTGCCCCCACTGCGGCGCCGGCGGCAACCCCTACCAGCCCGGGGTCGGCCAGCGGGTTGCGGAACAGCCCCTGCATGGc                                                                                                                    <  1:472709/150‑1 (MQ=255)                            cGACAGCAGATAAGGCGCGAAAGCATCCGGCATGCCGCCAAACCAGGCGCCGCCGACAATGGCAACCGCTGCCCCCACTGCGGCGCCGGCGGCAACCCCTACCAGCCCGGGGTggccatgcaggggctgttccgcaacccgctggccgac                                                                                                                 <  1:61081‑M1/150‑38 (MQ=255)                            cGACAGCAGATAAGGCGCGAAAGCATCCGGCATGCCGCCAAACCAGGCGCCGCCGACAATGGCAACCGCTGCCCCCACTGCGGCGCCGGCGGCAACCCCTACCAGCCCGGGGTggccatgcaggggctgttccgcaacccgctggccgac                                                                                                                 <  2:612739‑M1/150‑38 (MQ=255)                                                aaGCATCCGGCATGCCGCCAAACCAGGCGCCGCCGACAATGGCAACCGCTGCCCCCACTGCGGCGCCGGCGGCAACCCCTACCAGCCCGGGGTCGGCCAGCGGGTTGCGGAACAGCCCCTGCATGGCCACCCCGGACAGCGCCAGCACcg                                                                                             >  2:203728/1‑150 (MQ=255)                                                               cgccAAACCAGGCGCCGCCGACAATGGCAACCGCTGCCCCCACTGCGGCGCCGGCGGCAACCCCTACCAGCCCGGGGTCGGCCAGCGGGTTGCGGAACAGCCCCTGCATGGCCACCCCGGACAGCGCCAACACCGCGCCCACCGCCAAcc                                                                              <  1:392313/150‑1 (MQ=255)                                                                gccAAACCAGGCGCCGCCGACAATGGCAACCGCTGCCCCCACTGCGGCGCCGGCGGCAACCCCTACCAGCCCGGGGTCGGCCAGCGGGTTGCGGAACAGCGCCTGCATGGCCACCCCGGACAGCGCCAGCACCGCGCCCACCGCCAAccc                                                                             >  1:538758/1‑150 (MQ=255)                                                                 ccAAACCAGGCGCCGCCGACAATGGCAACCGCTGCCCCCACTGCGGCGCCGGCGGCAACCCCTACCAGCCCGGGGTCGGCCAGCGGGTTGCGGAACAGCCCCTGCATGGCCACCCCGGACAGCGCCAGCACCGCGCCCACCGCCAAcccc                                                                            >  2:247348/1‑150 (MQ=255)                                                                       cAGGCGCCGCCGACAATGGCAACCGCTGCCCCCACTGCGGCGCCGGCGGCAACCCCTACCAGCCCGGGGTCGGCCAGCGGGTTGCGGAACAGCCCCTGCATGGCCACCCCGGACAGCGCCAGCACCGCGCCCACCGCCAACCCCAACAgg                                                                      <  2:216129/150‑1 (MQ=255)                                                                       cAGGCGCCGCCGACAATGGCAACCGCTGCCCCCACTGCGGCGCCGGCGGCAACCCCTACCAGCCCGGGGTCGGCCAGCGGGTTGCGGAACAGCCCCTGCATGGCCACCCCGGACAGCGCCAGCACCGCGCCCACCGCCAACCCCAACAgg                                                                      >  2:437402/1‑150 (MQ=255)                                                                                  gACAATGGCAACCGCTGCCCCCACTGCGGCGCCGGCGGCAACCCCTACCAGCCCGGGGTCGGCCAGCGGGTTGCGGAACAGCCCCTGCATGGCCACCCCGGACAGCGCCGGCACCGCGCCCACCGCCAACCCCAACAGGGTGCGCGGCAg                                                           <  1:588477/150‑1 (MQ=255)                                                                                   aCAATGGCAACCGCTGCCCCCACTGCGGCGCCGGCGGCAACCCCTACCATCCCGGGGTggccatgcaggggctgttccgcaacccgctggccgaccccg                                                                                                             <  1:168851‑M1/99‑42 (MQ=255)                                                                                   aCAATGGCAACCGCTGCCCCCACTGCGGCGCCGGCGGCAACCCCTACCAGCCCGGGGTggccatgcaggggctgttccgcaacccgctggccgaccccg                                                                                                             >  2:168851‑M1/1‑58 (MQ=255)                                                                                                   cccccACTGCGGCGCCGGCGGCAACCCCTACCAGCCCGGGGCCGGCCAGCGGGTTGCGGAACAGCCCCTGCATGGCCACCCCGGACAGCGCCAGCACCGCGCCCACCGCCAACCCCAACAGGGTGCGCGGCAGGCGGATCTGACCCAGGa                                          >  2:217514/1‑150 (MQ=255)                                                                                                        aCTGCGGCGCCGGCGGCAACCCCTACCAGCCCGGGGTCGGCCAGCGGGTTGCGGAACAGCCCCTGCATGGCCACCCCGGACAGCGCCAGCACCGCGCCCACCGCCAACCCCAACAGGGTGCGCGGCAGGCGGATCTGACCAAGGatcatc                                     <  1:437402/150‑1 (MQ=255)                                                                                                                     cggcAACCCCTACCAGCCCGGGGTCGGCCAGCGGGTTGCGGAACAGCCCCTGCATGGCCACCCCGGACAGCGCCAGCACCGCGCCCACCGCCAACCCCAACAGGGTGCGCGGCAGGCGGATCTGACCAAGGATCATCTCGGCCTGCTGCa                        <  2:529953/150‑1 (MQ=255)                                                                                                                          aCCCCTACCAGCCCGGGGTCGGCCCGCGCGTTGCGGAACAGCCCCTGCATGGCCACCCCGGACAGCGCCAGCACCCCACAAACCGCCAACACCAGAAACCGATGCGCCTTTTTTTGCTTACCAAGGATCATCTCGGCCTGCTGCAGGCCa                   <  1:568748/150‑1 (MQ=255)                                                                                                                              cTACCAGCCCGGGGTCGGCCAGCGGGTTGCGGAACAGCCCCTGCATGGCCACCCCGGACAGCGCCAGCACCGCGCCCACCGCCAACCCCAACAGGGTGCGCGGCAGGCGGATCTGACCAAGGATCATCTCGGCCTGCTGCAGGCCATCtt               >  1:377401/1‑150 (MQ=255)                                                                                                                                 ccAGCCCGGGGTCGGCCAGCGGGTTGCGGAACAGCCCCTGCATGGCCACCCCAGACAGCGCCCGCACCGCGCCCACCGCCAACCCCAACAGGGTGCGCGGCAGGCGGATCTGACCAAGGATCATCTCGGCCTGCTGCAGGCCATCTTCgg            >  1:354380/1‑150 (MQ=255)                                                                                                                                       cGGGGTCGGCCAGCGGGTTGCGGAACAGCCCCTGCATGGCCaccccgggctggtaggggttgccgccggcgccgcagtgggggcagcggttgccattgt                                                         >  2:168851‑M2/1‑41 (MQ=255)                                                                                                                                       cGGGGTCGGCCAGCGGGTTGCGGAACAGCCCCTGCATGGCCaccccgggatggtaggggttgccgccggcgccgcagtgggggcagcggttgccattgt                                                         <  1:168851‑M2/99‑59 (MQ=255)                                                                                                                                           gTCGGCCAGCGGGTTGCGGAACAGCCCCTGCATGGCCaccccgggctggtaggggttgccgccggcgccgcagtgggggcagcggttgccattgtcggcggcgcctggtttggcggcatgccggatgctttcgcgccttatctgctgtcg  <  2:612739‑M2/150‑114 (MQ=255)                                                                                                                                           gTCGGCCAGCGGGTTGCGGAACAGCCCCTGCATGGCCaccccgggctggtaggggttgccgccggcgccgcagtgggggcagcggttgccattgtcggcggcgcctggtttggcggcatgccggatgctttcgcgccttatctgctgtcg  <  1:61081‑M2/150‑114 (MQ=255)                                                                                                                                            |                                                                                                                                                     CCCAAGCCCGCCGACAAACGCGCAGAACGACAGCAGATAAGGCGCGAAAGCATCCGGCATGCCGCCAAACCAGGCGCCGCCGACAATGGCAACCGCTGCCCCCACTGCGGCGCCGGCGGCAACCCCTACCAGCCCGGGGTCGGCCAGCGGGTTGCGGAACAGCCCCTGCATGGCCACCCCGGACAGCGCCAGCACCGCGCCCACCGCCAACCCCAACAGGGTGCGCGGCAGGCGGATCTGACCAAGGATCATCTCGGCCTGCTGCAGGCCATCTTCGGAAACGGGCAG  >  NC_002947/5327958‑5328245

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