New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	5910975 =	32 (1.230)	3 (0.120)	3/280	3.4	8.3%	coding (1183/1359 nt)	spuB	glutamylpolyamine synthetase
	?	NC_002947	5911007 =	35 (1.400)					coding (1151/1359 nt)	spuB	glutamylpolyamine synthetase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CACGAACTGATGACGCGCTTGAAGTTTTCATGCTCGGCGCGCTTGGTAGCGACATACCCGGTGATGAATTTCTTGCCCAGGTACTGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGGCGGCGATGGCCAGGTACGGGTTGGCGTCGGCGCCCGGCAGGCGGTTCTCGACGC  >  NC_002947/5910830‑5911118                                                                                                                                                  |                                                                                                                                                cACGAACTGATGACGCGCTTGAAGTTTTCATGCTCGGCGCGCTTGGTAGCGACGTACCCGGTGATGAATTTCTTGCCCAGGTACTGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCa                                                                                                                                               <  5:131688/148‑1 (MQ=255)               gcgcTTGAAGTTTTCATGCTCGGCGCGCTTGGTAGCGACATACCCGGTGATGAATTTCTTGCCCAGGTACTGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTc                                                                                                                                  >  2:114461/1‑147 (MQ=255)                 gcTTGAAGTTTTCATGCTCGGCGCGCTTGGTAGCGACATACCCGGTGATGAATTTCTTGCCCAGGTACTGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTa                                                                                                                               <  3:55446/148‑1 (MQ=255)                     gAAGTTTTCATGCTCGGCGCGCTTGGTAGCGACATACCCGGTGATGAATTTCTTGCCCAGGTACTGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCa                                                                                                                           >  5:95261/1‑148 (MQ=255)                          tttCATGCTCGGCGCGCTTGGTAGCGACATACCCGGTGATGAATTTCTTGCCCAGGTACTGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGcc                                                                                                                      <  4:140685/148‑1 (MQ=255)                                                  cGACATACCCGTTGATGAATTTCTTGCCCACGTACTGCACCAGCGCACGGCTGTTTTCCCTGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGAt                                                                                              <  7:27025/148‑1 (MQ=255)                                                       tACCCGGTGATGAATTTCTTGCCCAGGTACTGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCtt                                                                                         <  6:90358/148‑1 (MQ=255)                                                          ccGGTGATGAATTTCTTGCCCAGGTACTGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGa                                                                                      >  2:59793/1‑148 (MQ=255)                                                                    aTTTCTTGCCCAGGTACTGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGa                                                                             <  2:116891/147‑1 (MQ=255)                                                                      ttCTTGCCCAGGTACTGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGt                                                                          <  1:105295/148‑1 (MQ=255)                                                                           gCCCAGGTACTGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCca                                                                     >  8:45191/1‑148 (MQ=255)                                                                             ccAGGTACTGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCcac                                                                    >  2:52307/1‑147 (MQ=255)                                                                               aGGTACTGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGgcgc                                                                                                          >  4:77813/1‑107 (MQ=255)                                                                               aGGTACTGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCcacaca                                                                 >  4:65583/1‑148 (MQ=255)                                                                                ggTACTGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGgcgcg                                                                                                         <  3:77813/107‑1 (MQ=255)                                                                                   aCTGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAg                                                             <  2:27863/148‑1 (MQ=255)                                                                                     tGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGc                                                           >  6:46345/1‑148 (MQ=255)                                                                                     tGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGc                                                           >  3:115581/1‑148 (MQ=255)                                                                                              gcACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGGCGGCGa                                                  <  8:20709/148‑1 (MQ=255)                                                                                                    cTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTCGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGGCGGCGATGGCCa                                            <  2:123021/148‑1 (MQ=255)                                                                                                       ttttCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCCCACAGTAGGCTGGCGGCGATGGCCAGGt                                         >  5:79457/1‑148 (MQ=255)                                                                                                       ttttCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGGCGGCGATGGCCAGGt                                         >  3:106332/1‑148 (MQ=255)                                                                                                             aTGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGGCGGCGATGGCCAGGTACggg                                    <  8:46615/147‑1 (MQ=255)                                                                                                                  ttCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGGCGGCGATGGCCAGGTACGGGTTGGCg                              >  7:42991/1‑148 (MQ=255)                                                                                                                     cAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGGCGGCGATGGCCAGGTACGGGTTGGCGTc                            <  6:45761/147‑1 (MQ=255)                                                                                                                              cTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGGCGGCGATGGCCa                                            >  8:62072/1‑122 (MQ=255)                                                                                                                              cTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGGCGGCGATGGCCAGGTACGGGTTGGCGTCGGCGCCCg                    >  1:100275/1‑146 (MQ=255)                                                                                                                               tCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGGCGGCGATGGCCAg                                           <  7:62072/122‑1 (MQ=255)                                                                                                                                 gATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGGCGGCGATGGCCAGGTACGGGTTGGCGTCGGCGCCCGGCAgg               >  5:71018/1‑148 (MQ=255)                                                                                                                                    gggcAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGGCGGCGATGGCCAGGTACGGGTTGGCGTCGGCGCCCGGCAGGCgg            <  3:95427/145‑1 (MQ=255)                                                                                                                                         ggggCAGGCGCAGGTTGCGGCGTTCGGAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGGCGGCGATGGCCAGGTACGGGTTGGCGTCGGCGCCCGGCAGGCGGTtctc       <  3:117653/148‑1 (MQ=255)                                                                                                                                              aGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGGCGGCGATGGCCAGGTACGGGTTGGCGTCGGCGCCCGGCAGGCGGTTCTCGAcgc  <  7:125482/148‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                CACGAACTGATGACGCGCTTGAAGTTTTCATGCTCGGCGCGCTTGGTAGCGACATACCCGGTGATGAATTTCTTGCCCAGGTACTGCACCAGCGCACGGCTGTTTTCCATGCGTTCCAAGGCGTCCTCGATGGTCAGGGGCAGGCGCAGGTTGCGGCGTTCGTAGCCACGGCCCTGCACTGGCGCGCTGGCCTCGATACCTTCGACCATGCCGATGTAGCCACACAGTAGGCTGGCGGCGATGGCCAGGTACGGGTTGGCGTCGGCGCCCGGCAGGCGGTTCTCGACGC  >  NC_002947/5910830‑5911118

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