New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	3369704 =	21 (0.710)	4 (0.140)	3/236	NT	18.5%	coding (2967/3015 nt)	PP_2964	transposase
	?	NC_002947	3369735 =	15 (0.530)					coding (2998/3015 nt)	PP_2964	transposase
Rejected: Coverage evenness skew score above cutoff.

CAGCAACGTCAAGCGCTTCGGTGAGTTCAACCTCAAGCTCAATAGGCCACCGGAGCCCTGGATCAAGGACTCGGTATTCCAACAGGCCGCCGGCTCGCTTCGAGCCAACACGGCCAGCAAGACCGGTGCCTAGGAGGCAACATGATCGAGATCCCTCCCGCGACTTTTCGAGTCACACCCTACGGAGAAGTGGATGGTGTAGCACTGGACAAGCTGCGAGCCAGCTTCGACACCTCCCAGCTTCTGCGACTGGTTGACCAATTCG  >  NC_002947/3369620‑3369884                                                                                                                    |                                                                                                                                                      cAGCAACGTCAAGCGCTTCGGTGAGTTCAACCTCAAGCTCAATAGGCCACCGGAGCCCTGGATCAAGGACTCGGTATTCCAACAGGCCGCCGGCTCGCTTCGAGCCAACACGGCCAGCAAGACCGGTGCCTAGGAGGCAACATGAt                                                                                                                         >  1:627318/1‑146 (MQ=255) cAGCAACGTCAAGCGCTTCGGTGAGTTCAACCTCAAGCTCAATAGGCCACCGGAGCCCTGGATCAAGGACTCGGTATTCCAACAGGCCGCCGGCTCGCTTCGAGCCAACACGGCCAGCAAGACCGGTGCCTAGGAGGCAACATGAt                                                                                                                         <  2:627318/146‑1 (MQ=255)        gTCAAGCGCTTCGGTGAGTTCAACCTCAAGCTCAATAGGCCACCGGAGCCCTGGATCAAGGACTCGGTATTCCAACAGGCCGCCGGCTCGCTTCGAGCCAACACGGCCAGCAAGACCGGTGCCTAGGAGGCAACATGATCGAGAtccctc                                                                                                              >  1:717530/1‑150 (MQ=255)                      tGAGTTCAACCTCAAGCTCAATAGGCCACCGGAGCCCTGGATCAAGGACTCGGTATTCCAACAGGCCGCCGGCTCGCTTCGAGCCAACACGGCCAGCAAGACCGGTGCCTAGGAGGCAACATGATCGAGATCCCTCCCGCGACTTTTCGa                                                                                                <  1:512992/150‑1 (MQ=255)                          ttCAACCTCAAGCTCAATAGGCCACCGGAGCCCTGGATCAAGGACTCGGTATTCCAACAGGCCGCCGGCTCGCTTCGAGCCAACACGGCCAGCAAGACCGGTGCCTAGGAGGCAACATGATCGAGATCCCTCCCGCGACTTTTCGAGTca                                                                                            >  1:523357/1‑150 (MQ=255)                          ttCAACCTCAAGCTCAATAGGCCACCGGAGCCCTGGATCAAGGACTCGGTATTCCAACAGGCCGCCGGCTCGCTTCGAGCCAACACGGCCAGCAAGACCGGTGCCTAGGAGGCAACATGATCGAGATCCCTCCCGCGACTTTTCGAGTca                                                                                            >  1:310290/1‑150 (MQ=255)                                       tCAATAGGCCACCGGAGCCCTGGATCAAGGACTCGGTATTCCAACAGGCCGCCGGCTCGCTTCGAGCCAACACGGCCAGCAAGACCGGTGCCTAGGAGGCAACATGATCGAGATCCCTCCCGCGACTTTTCGAGTCACACCCTACGgaga                                                                               >  2:463758/1‑150 (MQ=255)                                          aTAGGCCACCGGAGCCCTGGATCAAGGACTCGGTATTCCAACAGGCCGCCGGCTCGCTTCGAGCCAACACGGCCAGCAAGACCGGTGCCTAGGAGGCAACATGATCGAGATCCCTCCCGCGACTTTTCGAGTCACACCCTACGGAGAAGt                                                                            <  1:463758/150‑1 (MQ=255)                                               ccACCGGAGCCCTGGATCAAGGACTCGGTATTCCAACAGGCCGCCGGCTCGCTTCGAGCCAACACGGCCAGCAAGACCGGTGCCTAGGAGGCAACATGATCGAGATCCCTCCCGCGACTTTTCGAGTCACACCCTACGGAGAAGTGGATg                                                                       >  1:657392/1‑150 (MQ=255)                                                 atcGGAGCCCTGGATCAAGGACTCGGTATTCCAACAGGCCGCCGGCTCGCTTCGAGCCAACACGGCCAGCAAGACCGGTGCCTAGGAGGCAACATGATCGAGATCCCTCCCGCGACTTTTCGAGTCACACCCTACGGAGAAGTGGATGgt                                                                     >  1:34189/3‑150 (MQ=255)                                                           tGGATCAAGGACTCGGTATTCCAACAGGCCGCCGGCTCGCTTCGAGCCAACACGGCCAGCAAGACCGGTGCCTAGGAGGCAACATGATCGAGATCCCTCCCGCGACTTTTCGAGTCACACCCTACGGAGAAGTGGATGGTGTAGCACTgg                                                           >  2:473347/1‑150 (MQ=255)                                                            ggATCAAGGACTCGGTATTCCAACAGGCCGCCGGCTCGCTTCGAGCCAACACGGCCAGCAAGACCGGTGCCTAGGAGGCAACATGATCGAGATCCCTCCCGCGACTTTTCGAGTCACACCCTACGGAGAAGTGGATGGTGTAGCACTGGa                                                          <  2:34189/150‑1 (MQ=255)                                                                       tCGGTATTCCAACAGGCCGCCGGCTCGCTTCGAGCCAACACGGCCAGCAAGACCGGTGCCTAGGAGGc                                                                                                                                 >  1:754482/1‑68 (MQ=255)                                                                       tCGGTATTCCAACAGGCCGCCGGCTCGCTTCGAGCCAACACGGCCAGCAAGACCGGTGCCTAGGAGGc                                                                                                                                 <  2:754482/68‑1 (MQ=255)                                                                                  aCAGGCCGCCGGCTCGCTTCGAGCCAACACGGCCAGCAAGACCGGTGCCTAGGAGGCAACATGATCGAGATCCCTCCCGCGACTTTTCGAGTCACACCCTACGGAGAAGTGGATGGTGTAGCACTGGACAAGCTGCGAGCCAGCTTCGac                                    >  2:775469/1‑150 (MQ=255)                                                                                                           aaCACGGCCAGCAAGACCGGTGCCTAGGAGGCAACATGATCGAGATCCCTCCCGCGACTTTTCGAGTCACACCCTACGGAGAAGTGGATGGTGGAGCACTGGACCAGCTGCGAGCCAGCTGCGACACCTCGCAGCTTCTGCGACTGGTTg           >  2:92350/1‑150 (MQ=255)                                                                                                             cacGGCCAGCAAGACCGGTGCCTAGGAGGCAACATGATCGAGATCCCTCCCGCGACTTTTCGAGTCACACCCTACGGAGAAGTGGATGGTGTAGCACTGGACAAGCTGCGAGCCAGCTTCGACACCTCCCAGCTTCTGCGACTGGTTGAc         <  2:310290/150‑1 (MQ=255)                                                                                                                    aGCAAGACCGGTGCCTAGGAGGCAACATGATCGAGATCCCTCCCGCGACTTTTCGAGTCACACCCTACGGAGAAGTGGATGGTGTAGCACTGGACAAGCTGCGAGCCAGCTTCGACACCTCCCAGCTTCTGCGACTGGTTGACCAATTCg  >  1:382180/1‑150 (MQ=255)                                                                                                                    |                                                                                                                                                      CAGCAACGTCAAGCGCTTCGGTGAGTTCAACCTCAAGCTCAATAGGCCACCGGAGCCCTGGATCAAGGACTCGGTATTCCAACAGGCCGCCGGCTCGCTTCGAGCCAACACGGCCAGCAAGACCGGTGCCTAGGAGGCAACATGATCGAGATCCCTCCCGCGACTTTTCGAGTCACACCCTACGGAGAAGTGGATGGTGTAGCACTGGACAAGCTGCGAGCCAGCTTCGACACCTCCCAGCTTCTGCGACTGGTTGACCAATTCG  >  NC_002947/3369620‑3369884

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