New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	3828223 =	10 (0.110)	29 (0.350)	21/252	4.1	69.6%	intergenic (‑204/+155)	PP_3381/PP_3382	transposase/gluconate 2‑dehydrogenase cytochrome c subunit
	?	NC_002947	3828259 =	16 (0.190)					intergenic (‑240/+119)	PP_3381/PP_3382	transposase/gluconate 2‑dehydrogenase cytochrome c subunit
Rejected: Coverage evenness skew score above cutoff.

CCCTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAGGGCGCTGAAAAACACGTCAGCGAAGCGACTTTACATCACCG  >  NC_002947/3828211‑3828402                                                 |                                                                                                                                                cccTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGa                                            <  1:617828/150‑1 (MQ=32) cccTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGa                                            >  2:1178248/1‑150 (MQ=32)  ccTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAg                                           >  1:1985135/1‑150 (MQ=32)   cTAACAAAGCAAAAGGCGGCTGGCCCTCCCGGGATAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAgg                                          <  2:1985135/150‑1 (MQ=16)    tAACAAATCAAAAGTCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAggg                                         <  2:754180/150‑1 (MQ=21)       cAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAGGgcgc                                      >  2:1416112/1‑150 (MQ=37)        aaaGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAGGGCGCt                                     >  2:301770/1‑150 (MQ=38)         aaGCAAAGGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCCCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAGGGCGCTg                                    >  2:1609401/1‑150 (MQ=255)            cAAAGGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAGGGCGCTGaaa                                 <  1:1609401/150‑1 (MQ=255)            cAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAGGGCGCTGaaa                                 <  1:1232126/150‑1 (MQ=40)              aaaGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAGGGCGCTGaaaaa                               <  1:1483761/150‑1 (MQ=255)                aGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAGGGCGCTGAAAAaca                             <  2:116504/150‑1 (MQ=255)                 ggcggcTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAGGGCGCTGAAAAacac                            <  1:2061123/150‑1 (MQ=255)                           ccggtagggCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAGGGCGCTGAAAAACACGTCagcgaag                  <  2:906671/144‑1 (MQ=255)                            cTACAGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAGGGCGCTGAAAAACACGTCagcgaaac                 >  1:1085335/1‑148 (MQ=255)                              aCCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAGGGCGCTGAAAAACACGTCagcgaagcga               >  2:142825/1‑150 (MQ=255)                               ccGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGGCGCCGGACACCTCATGGTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAGGGCGCTGAAAAACACGTCAGCGAAGCGAc              <  1:255383/150‑1 (MQ=255)                                         cGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAGGGCGCTGAAAAACACGTCAGCGAAGCGACTTTACATCAc    >  2:1289760/1‑150 (MQ=255)                                           ccTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAGGGCGCTGAAAAACACGTCAGCGAAGCGACTTTACATCACcg  >  2:2117528/1‑150 (MQ=255)                                                 |                                                                                                                                                CCCTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGCGAGGGCGCTGAAAAACACGTCAGCGAAGCGACTTTACATCACCG  >  NC_002947/3828211‑3828402

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.