New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	5552330 =	42 (0.760)	3 (0.060)	3/242	NT	7.5%	coding (1559/1566 nt)	PP_4882	iron ABC transporter permease
	?	NC_002947	5552377 =	34 (0.650)					intergenic (+40/+621)	PP_4882/PP_5717	iron ABC transporter permease/hypothetical protein

AGCCCGACCGGCATGACCACCCTGGCCACCGAGGTGTGGGCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑GAGTG‑CCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAGGG  >  NC_002947/5552182‑5552475                                                                                                                                                     |                                                                                                                                                    aGCCCGACCGGCATGACCACCCTGGCCACCGAGGTGTGGGCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACaa                                                                                                                                                     <  7:232155/149‑1 (MQ=255)   tttcagatctcgagggtgtcccaaggtctaagaaaagagtgttgttaccctcttgggtggggcttgtagtatctacgaccttgttctttgctttcaactggctagtgagagttttaaagtgaacttgtgatgaagtctaggcactcAgg                                                                                                                                                   >  4:36834‑M1/147‑149 (MQ=255)    ccGACCGGCATGACCACCCTGGCCACCGAGGTGTGGGCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGc                                                                                                                                                  <  3:33979/149‑1 (MQ=255)       agatctcgagggtttcccaaggtctatgaaaagagttttgttaccctcttgggtggggcttgtagtatctacgaccttgttctttgctttcaactggctagtgagagttttaaagtgaacttgtgatgaagtctaggcactcAGGCTTg                                                                                                                                               <  5:95609‑M1/7‑1 (MQ=255)         cGGCATGACCACCCTGGCCACCGAGGTGTGGGCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGa                                                                                                                                              >  1:295153/1‑148 (MQ=255)          ggCATGACCACCCTGGCCACCGAGGTGTGGGCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGa                                                                                                                                              >  5:162406/1‑147 (MQ=255)                   accCTGGCCACCGAGGTGTGGGCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGcc                                                                                                                                     >  4:136828/1‑147 (MQ=255)                    cccTGGCCACCGAGGTGTGGGCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAg                                                                                                                                  >  8:37066/1‑149 (MQ=255)                    cccTGGCCACCGAGGTGTGGGCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAg                                                                                                                                  >  4:83888/1‑149 (MQ=255)                           cACCGAGGTGTGGGCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGcc                                                                                                                                     >  3:264713/1‑139 (MQ=255)                               gAGGTGTGGGCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCAtgtg                                                                                                                       <  7:37066/149‑1 (MQ=255)                                 tgaaaagagttttgttaccctcttgggtggggcttgtagtatctacgaccttgttctttgctttcaactggctagtgagagttttaaagtgaacttgtgatgaagtctaggcactcAGGc                                                                                                                                                  >  2:301536‑M1/117‑120 (MQ=255)                                     tgGGCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTgg                                                                                                                 <  2:160398/149‑1 (MQ=255)                                     tgGGCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTgg                                                                                                                 <  6:361530/149‑1 (MQ=255)                                     tgGGCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTg                                                                                                                  <  3:107796/148‑1 (MQ=255)                                     tgGGCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTg                                                                                                                  <  7:302438/148‑1 (MQ=255)                                       ggCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCACATGGCGGCCTGGGc                                                                                                               <  3:238999/149‑1 (MQ=255)                                        gCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGca                                                                                                              >  5:370549/1‑149 (MQ=255)                                                 gCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑Ga                                                                                                    <  3:83888/149‑1 (MQ=255)                                                 gCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑Ga                                                                                                    <  8:141611/149‑1 (MQ=255)                                                 accctcttgggtggggcttgtagtatctacgaccttgttctttgctttcaactggctagtgagagttttaaagtgaacttgtgatgaagtctaggcactcAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑Ga                                                                                                    <  3:213406‑M1/49‑1 (MQ=255)                                                  ccctcttgggtggggcttgtagtatctacgaccttgttctttgctttcaactggctagtgagagttttaaagtgaacttgtgatgaagtctaggcactcAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑Ga                                                                                                    >  4:213406‑M1/100‑148 (MQ=255)                                                   cctcttgggtggggcttgtagtatctacgaccttgttctttgctttcaactggctagtgagagttttaaactgaacttgtgatgaagtctaggcactcAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑GAGt                                                                                                  >  3:264591‑M1/99‑149 (MQ=255)                                                        tCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑GAGTG‑CCTa                                                                                            <  5:6649/149‑1 (MQ=255)                                                            gTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑GAGt                                                                                                  >  3:331213/1‑140 (MQ=255)                                                            gTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑GAGt                                                                                                  <  4:331213/140‑1 (MQ=255)                                                             ttCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑GAGTG‑CCTAGACtt                                                                                       >  7:177152/1‑149 (MQ=255)                                                                gcggcggcggcgCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑GAGTG‑CCTAGACTtcat                                                                                    >  2:254053/1‑149 (MQ=255)                                                                     ggcggcgCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑GAGTG‑CCTAGACTTCATCACaa                                                                               <  8:53046/149‑1 (MQ=255)                                                                            ccGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑GAGTG‑CCTAGACTTCATCACAAGTTCACt                                                                        >  4:240773/1‑149 (MQ=255)                                                                                    ggcgTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑GAGTG‑CCTAGACTTCATCACAAGTTCACTTTAAAAct                                                                >  6:95609/1‑149 (MQ=255)                                                                                    ggcgTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑GAGTG‑CCTAGACTTCATCACAAGTTCACTTTAAAAct                                                                <  3:136828/149‑1 (MQ=255)                                                                                    ggcgTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑GAGTG‑CCTAGACTTCATCACAAGTTCACTTTAAAAc                                                                 >  2:89529/1‑148 (MQ=255)                                                                                      cgTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑GAGTG‑CCTAGACTTCATCACAAGTTCACTTTAAAActct                                                              <  8:177152/149‑1 (MQ=255)                                                                                           ttgATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑GAGTG‑CCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTa                                                         <  3:240773/149‑1 (MQ=255)                                                                                              aTTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGGAGTG‑CCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGc                                                       <  1:254053/149‑1 (MQ=255)                                                                                                    gtgtCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑GAGTG‑CCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTg                                                 >  3:333529/1‑148 (MQ=255)                                                                                                       tCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑GAGTG‑CCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGaaa                                              <  6:370549/148‑1 (MQ=255)                                                                                                         ggggTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑GAGTG‑CCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGc                                            >  3:147824/1‑148 (MQ=255)                                                                                                           tgtttacGGTGTATTTTCTTAATATGCTGATTTTTTTGAACAAGGCTTTAAATCTTTCCCATGCCGCCATTTGTCCTTGGCACAAGCCT‑GAGTGCCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCaa                                          <  1:89529/143‑1 (MQ=255)                                                                                                           ggTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑GAGTG‑CCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCaaa                                         <  2:295153/149‑1 (MQ=255)                                                                                                                        tttGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑GAGTG‑CCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCg                              >  2:254649/1‑147 (MQ=255)                                                                                                                          tGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCACATGGCGGCCTGGGcaag                                                                                                            <  7:81472/69‑3 (MQ=21)                                                                                                                          tGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCACATGGCGGCCTGGGcaag                                                                                                            >  8:81472/1‑67 (MQ=21)                                                                                                                                    gcgGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑GAGTG‑CCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGc                >  2:299513/1‑149 (MQ=255)                                                                                                                                      ggATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑GAGTG‑CCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAg                                        <  7:195857/123‑1 (MQ=255)                                                                                                                                      ggATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑GAGTG‑CCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAg                                        >  8:195857/1‑123 (MQ=255)                                                                                                                                      ggATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑GAGTG‑CCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGccc              <  6:162406/149‑1 (MQ=255)                                                                                                                                             tttGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑GAGTG‑CCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCaaa                                         >  2:274736/1‑115 (MQ=255)                                                                                                                                             tttGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑GAGTG‑CCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCaaa                                         <  1:274736/115‑1 (MQ=255)                                                                                                                                                  aCAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑GAGTG‑CCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAggg  >  3:149898/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    AGCCCGACCGGCATGACCACCCTGGCCACCGAGGTGTGGGCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCT‑GAGTG‑CCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAGGG  >  NC_002947/5552182‑5552475

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