New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 5200727	33 (1.500)	3 (0.140)	3/246	NT	8.9%	coding (324/774 nt)	yafJ	amidotransferase
	?	NC_002947	= 5200748	30 (1.420)					coding (303/774 nt)	yafJ	amidotransferase

ACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGCGGCCCCACATTTCCCGCACGAACGGGTGGGTGTTGGACAGGCACACCCTGCCGACGTTGGCCTGACGGATATGACCGATCACCACTTCACTTTTGATCGGGTAGCGCTGCACGAGGTTGGCCACTTC  >  NC_002947/5200602‑5200876                                                                                                                                                   |                                                                                                                                 accggcaccggTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGcgg                                                                                                                               <  1:364803‑M2/150‑4 (MQ=255)              ccgggAAGGCGACGTGGATTCGGTTTAGCAAGTCGCAGAAGGGCGCTTCCCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGCGGCCCCACATTTccc                                                                                                                  <  1:280191/148‑1 (MQ=255)                            cGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGTTGTGCGCGCAGCACCAGTTGCGGCCCCACAGTTCCCGCAGGAGCgggtgg                                                                                                    <  2:364491/150‑1 (MQ=255)                            cGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGCGGCCCCACATTTCCCGCACGAACgggtgg                                                                                                    <  2:323863/150‑1 (MQ=255)                                    gTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGCGGCCCCACATTTCCCGCACGAACGGGTGGGTGTTGGa                                                                                            >  2:203788/1‑150 (MQ=255)                                     ttGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGCGGCCCCACATTTCCCGCACGAACGGGTGGGTGTTGGAc                                                                                           <  1:133651/150‑1 (MQ=255)                                        aGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGCGGCCCCACATTTCCCGCACGAACGGGTGGGTGTTGGACAgg                                                                                        <  1:203788/150‑1 (MQ=255)                                         gCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCCCCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGCGGCCCCACATTTCCCGCACGAACGGGTGGGTGTTGGACAGGc                                                                                       <  1:181347/150‑1 (MQ=255)                                         gCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGCGGCCCCACATTTCCCGCACGAACGGGTGGGTGTTGGACAGGc                                                                                       <  2:435544/150‑1 (MQ=255)                                          cAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGCGGCCCCACATTTCCCGCACGAACGGGTGGGTGTTGGACAGGca                                                                                      >  2:425550/1‑150 (MQ=255)                                          cAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTACCAACTGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGCGGCCCCACATTTCCCGCACGAACGGGTGGGTGTTGGACAGGca                                                                                      <  1:546117/150‑1 (MQ=255)                                                      ggCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGCGGCCCCACATTTCCCGCACGAACGGGTGGGTGTTGGACAGGCACACCCTGCCGAc                                                                          >  1:406699/1‑150 (MQ=255)                                                                 tGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGCGGCCCCACATTTCCCGCACGAACGGGTGGGTGTTGGACAGGCACACGCTGCCGACGTTGGCCTGAc                                                               <  2:382024/150‑1 (MQ=255)                                                                         tCGCCCACTGGCCGGTAGGAACTGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGCGGCCCCACATTTCCCGCACGAACGGGTGGGTGTTGGACAGGCACACCCTGCCGACGTTGGCCTGACGGATATGa                                                       <  2:293072/150‑1 (MQ=255)                                                                           gCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGCGGCCCCACATTTCCCGCACGAACgggtgg                                                                                                    >  2:15921/1‑103 (MQ=255)                                                                           gCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGCGGCCCCACATTTCCCGCACGAACgggtgg                                                                                                    <  1:15921/103‑1 (MQ=255)                                                                           gCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGCGGCCCCACATTTCCCGCACGAACGGGTGGGTGTTGGACAgg                                                                                        <  2:160870/115‑1 (MQ=255)                                                                           gCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGCGGCCCCACATTTCCCGCACGAACGGGTGGGTGTTGGACAgg                                                                                        >  1:160870/1‑115 (MQ=255)                                                                               aCTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGCGGCCCCACATTTCCCGCACGAACGGGTGGGTGTTGGACAGGCACACCCTGCCGACGTTGGCCTGACGGATATGACCGATc                                                 >  2:471103/1‑150 (MQ=255)                                                                                   gCCGGTAGAAACGGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGCGGCCCAACATTCCCCGCACGAACGGGTGGGTGGCCGACAGGCACACCCTGCCGACTTTGGCCTGACGGATATGACCGATcaccc                                             <  2:177594/150‑2 (MQ=255)                                                                                           tctcTGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGCGGCCCCACATTTCCCGCACGAACGGGTGGGTGTTGGACAGGCACACCCTGCCGACGTTGGCCTGACGGATATGACCGATCACcacttcactt                                     <  1:122281/147‑1 (MQ=255)                                                                                            aaCTGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGTTGTTCGCGAAGCACCAGTTGCGGCCCCACATTTCCCGCACGAACGGGTGGGTGTTGGACAGGCACACCCTGCCGACGTTg                                                                      <  1:153601/116‑1 (MQ=255)                                                                                            aaCTGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGCGGCCCCACATTTCCCGCACGAACGGGTGGGTGTTGGACAGGCACACCCTGCCGACGTTg                                                                      >  2:153601/1‑116 (MQ=255)                                                                                                 gCCAGGCCCGTGAAGTCACCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGCGGCCCCACATTTCCCGCACGAACGGGTGGGTGTTGGACAGGCACACCCTGCCGACGTTGGCCTGACGGATATGACCGATCACCACTTCACTTTTGATc                               <  1:25992/150‑1 (MQ=255)                                                                                                         cTTGAAGTCACCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGCGGCCCCACATTTCCCGCACGAACGGGTGGGTGTTGGACAGGCACACCCTGCCGACGTTGGCCTGACGGATATGACCGATCACCACTTCACTTTTGATCGGGTAgcg                       >  2:448851/1‑150 (MQ=255)                                                                                                          ttGAAGTCACCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGCGGCCCCACATTTCCCGCACGAACGGGTGGGTGTTGGACAGGCACACCCTGCCGACGTTg                                                                      >  1:445538/1‑102 (MQ=255)                                                                                                          ttGAAGTCACCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGCGGCCCCACATTTCCCGCACGAACGGGTGGGTGTTGGACAGGCACACCCTGCCGACGTTg                                                                      <  2:445538/102‑1 (MQ=255)                                                                                                           tGAAGTCACCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGCGGTCCCACATTTCCCGCACGAACGGGTGGGTGTTGGACAGGCACACCCTGCCGACGTTGGCCTGACGGATATGACCGATCACCACTTCACTTTTGATCGGGTAGCGCt                     <  2:406699/150‑1 (MQ=255)                                                                                                                tCACCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGCGGCCCCACATTTCCCGCACGAACGGGTGGGTGTTGGACAGGCACACCCTGCCGACGTTGGCCTGACGGATATGACCGATCACCACTTCACTTTTGATCGGGTAGCGCTGCACg                >  2:366943/1‑150 (MQ=255)                                                                                                                   cccAGCTGGCCGTTGTGCGCGAAGCACCAGTTGcggccagctgggtgacttcaagggcctggccagtttctaccggccagtgggcgacaccgacagcgaagcgg                                                           >  1:176104‑M2/1‑33 (MQ=255)                                                                                                                   cccAGCTGGCCGTTGTGCGCGAAGCACCAGTTGcggccagctgggtgacttcaagggcctggccagtttctaccggccagtgggcgacaccgacagcgaagcgg                                                           <  2:176104‑M2/104‑72 (MQ=255)                                                                                                                    ccAGCTGGCCGTTGTGCGCGAAGCACCAGTTGcggccagctgggtgacttcaagggcctggccagtttctaccggccagtgggcgacaccgacagcgaagcggtcttctacggactccgcaaccgcatccccagcgccttcccggaaccg            <  2:429686‑M2/150‑119 (MQ=255)                                                                                                                        cTGGCCGTTGTGCGCGAAGCACCAGTTGCGGCCCCACATTTCCCGCACGAACGGGTGGGTGTTGGACAGGCACACCCTGCCGACGTTGGCCTGACGGATATGACCGATCACCACTTCACTTTTGATCGGGTAGCGCTGCACGAGGTTGGc        <  1:425550/150‑1 (MQ=255)                                                                                                                              gTTGTGCGCGAAGCACCAGTTGCGGCCCCACATTTCCCGCACGAACGGGTGGGTGTTGGACAGGCACACCCTGCCGACGTTGGCCTGACGGATATGACCGATCACCACTTCACTTTTGATCGGGTAGCGCTGCACGAGGTTGGCCACTTc  <  1:506058/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                 ACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGTTGTGCGCGAAGCACCAGTTGCGGCCCCACATTTCCCGCACGAACGGGTGGGTGTTGGACAGGCACACCCTGCCGACGTTGGCCTGACGGATATGACCGATCACCACTTCACTTTTGATCGGGTAGCGCTGCACGAGGTTGGCCACTTC  >  NC_002947/5200602‑5200876

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.