New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	6121904 =	29 (1.240)	5 (0.230)	5/242	NT	16.0%	coding (302/1386 nt)	PP_5371	rubredoxin/rubredoxin reductase
	?	NC_002947	6121950 =	25 (1.130)					coding (256/1386 nt)	PP_5371	rubredoxin/rubredoxin reductase

TGCGGCGGTGCGCCCCTGAGCCAGGGCGATGGATAGCGCGGGCTTGGAATAAAGATGGCCCGACTCCTGGGTCAGTACGCGAATCTCCAACTGCGGGGCGGCCTTGCGCAGGGCCTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAG  >  NC_002947/6121807‑6122094                                                                                                                                                |                                                                                                                                                 ctcagcaccccgatcgcagcagtggctaccgcgccggttatccaggcgcaaccactggttgccgctcgtcccgcccaaccggtggtcatcatcggtaGCGGCCTTGCGCAGGGCCTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCgat                                                                                                                                            >  2:509695‑M1/98‑150 (MQ=255)    ggcggTGCGCCCCTGAGCCAGGGCGATGGATAGCGCGGGCTTGGAATAAAGATGGCCCGACTCCTGGGTCAGTACGCGAATCTCCAACTGCGGGGCGGCCTTGCGCAGGGCCTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCgatgat                                                                                                                                         <  1:528046/150‑1 (MQ=255)              cccTGAGCCAGGGCGATGGATAGCGCGGGCTTGGAATACAGAAAGCCGAATTCTTGGGTTAGTACGCTAAAAGCGAACGGCGGGGGTGCCGTGCTCAGGGCCTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGt                                                                                                                               >  2:70474/1‑150 (MQ=255)               ccTGAGCCAGGGCGATGGATAGCGCGGGCTTGGAATAAAGATGGCCCGACTCCTGGGTCAGTACGCGAATCTCCAACTGCGGGGCGGCCTTGCGCAGGGCCTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGtt                                                                                                                              <  1:210573/150‑1 (MQ=255)                                 aTAGCGCGGGCTTGGAATAAAGATGGCCCGACTCCTGGGTCAGTACGCGAATCTCCAACTGCGGGGCGGCCTTGCGCCGAGCCTGTGCCAGGCCGTAGCCGGCATGGCCGCTCCCGCTGATGGACACCGGTTGTGTGCGATGACCGGCaa                                                                                                            >  1:167145/1‑150 (MQ=255)                                     cgcgGGCTTGGAATAAAGATGGCCCGACTCCTGGGTCAGTACGCGAATCTCCAACTGCGGGGCGGCCTTGCGCAGGGCCTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAg                                                                                                        <  2:306657/150‑1 (MQ=255)                                      gcgGGCTTGGAATAAAGATGGCCCGACTCCTGGGTCAGTACGCGAATCTCCAACTGCGGGGCGGCCTTGCGCAGGGCCTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGt                                                                                                       <  1:238481/150‑1 (MQ=255)                                       cgGGCTTGGAATAAAGATGGCCCGACTCCTGGGTCAGTACGCGAATCTCCAACTGCGGGGCGGCCTTGCGCAGGGCCTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTg                                                                                                      <  1:280175/150‑1 (MQ=255)                                              gaaaTAAAGATGGCCCGACTCCTGGGTCAGTACGCGAATCTCCAACTGCGGGGCGGCCTTGCGCAGGGCCTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTgcgc                                                                                               >  2:547414/3‑150 (MQ=255)                                                aaTAAAGATGGCCCGACTCCTGGGTCAGTACGCGAATCTCCAACTGCGGGGCGGCCTTGCGCAGGGCCTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCt                                                                                             >  1:479084/1‑150 (MQ=255)                                                 aTAAAGATGGCCCGACTCCTGGGTCAGTACGCGAATCTCCAACTGCGGGGCGGCCTTGCGCAGGGCCTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGc                                                                                                                          >  1:42549/1‑120 (MQ=255)                                                 aTAAAGATGGCCCGACTCCTGGGTCAGTACGCGAATCTCCAACTGCGGGGCGGCCTTGCGCAGGGCCTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGc                                                                                                                          <  2:42549/120‑1 (MQ=255)                                                            ccGACTCCTGGGTCAGTACGCGAATCTCCAACTGCGGGGCGGCCTTGCGCAGGGCCTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCgg                                                                                                                                <  2:497569/103‑1 (MQ=255)                                                            ccGACTCCTGGGTCAGTACGCGAATCTCCAACTGCGGGGCGGCCTTGCGCAGGGCCTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCgg                                                                                                                                >  1:497569/1‑103 (MQ=255)                                                             cGACTCCTGGGTCAGTACGCGAATCTCCAACTGCGGGGCGGCCTTGCGCAGGGCCTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGgcgc                                                                                <  2:531504/150‑1 (MQ=255)                                                                   cTGGGTCAGTACGCGAATCTCCAACTGCGGGGCGGCCTTGCGCAGGGCCTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGc                                                                          >  2:75418/1‑150 (MQ=255)                                                                                 gAATCTCCAACTGCGGGGCGGCCTTGCGCAGGGCCTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGAACACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCg                                                            <  1:75418/150‑1 (MQ=255)                                                                                    tctcCAACTGCGGGGCGGCCTTGCGCAGGGCCTGTGCCAGGCCGTAGCCGGACTGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCgggg                                                         >  2:151450/1‑150 (MQ=255)                                                                                       ccAACTGCGGGGCGGCCTTGCGCAGGGCCTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGc                                                      >  2:448978/1‑150 (MQ=255)                                                                                        cAACTGCGGGGCGGCCTTGCGCAGGGCCTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGAATACCGGTTGGGAGGGACGAGCGGAAACCAGTGGTTGCGCCTGGATAACCGTCGCGGTAGACACTGCTTCGATCCGGGTGCt                                                     >  2:38129/1‑150 (MQ=255)                                                                                         aaCTGCGGGGCGGCCTTGCGCAGGGCCTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCGatgatg                                                                                                                                        <  2:235332/66‑1 (MQ=255)                                                                                         aaCTGCGGGGCGGCCTTGCGCAGGGCCTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCGatgatg                                                                                                                                        >  1:235332/1‑66 (MQ=255)                                                                                           atcggtagcggccatgccggctacggcctggcacaggccctgcgcaaggccgcTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAg                                                  >  2:509695‑M2/54‑150 (MQ=255)                                                                                              cGGGGCGGCCTTGCGCAGGGCCTGTGCCAGGCCGTAGCCGGCATGGCCGATGCCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTc                                               <  1:38129/150‑1 (MQ=255)                                                                                                 ccatgccggctacggcctggcacaggccctgcgcaaggcatggccgcTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTg                                                                                                      >  2:349079‑M2/48‑92 (MQ=255)                                                                                                 ccatgccggctacggcctggcacaggccctgcgcaaggcatggccgcTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTg                                                                                                      <  1:349079‑M2/45‑1 (MQ=255)                                                                                                  gcggccatgccggctacggcctggcacaggccctgcgcaaggccgcTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATc                                           <  1:509695‑M2/104‑1 (MQ=255)                                                                                                   cggccatgccggctacggcctggcacaggccctgcgcaaggccgcTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCa                                          >  2:552936‑M2/46‑150 (MQ=255)                                                                                                          gcgcAGGGCCTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTaaaa                                   <  1:547414/150‑1 (MQ=255)                                                                                                                          gggCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGAcc                   <  2:167145/149‑1 (MQ=255)                                                                                                                                    cGGCCTGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAg         <  2:129516/150‑1 (MQ=255)                                                                                                                                          ggccgcTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATa   >  1:186735‑M2/7‑150 (MQ=255)                                                                                                                                           gccgcTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAg  >  1:345823‑M2/6‑150 (MQ=255)                                                                                                                                           gccgcTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGACGACCTTGCACTCAGGGCATAg  >  2:38963‑M2/6‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                 TGCGGCGGTGCGCCCCTGAGCCAGGGCGATGGATAGCGCGGGCTTGGAATAAAGATGGCCCGACTCCTGGGTCAGTACGCGAATCTCCAACTGCGGGGCGGCCTTGCGCAGGGCCTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAG  >  NC_002947/6121807‑6122094

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