New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_014034	1640006 =	70 (0.360)	35 (0.180)	13/184	22.3	0.34	intergenic (+86/‑59)	mccB/mccA	methylcrotonoyl‑CoA carboxylase subunit beta/methylcrotonoyl‑CoA carboxylase subunit alpha
	?	NC_014034	= 2590196	64 (0.300)					intergenic (‑95/‑161)	glnE/RCAP_rcc02415	glutamate‑ammonia‑ligase adenylyltransferase/hypothetical protein
Rejected: Position hash score below cutoff.

CCCGACAGGCGATCTCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGC....................................................................................................  <  NC_014034/1640113‑1640006 ....................................................................................................GGGGGCGCTGCCCCCCGCACCCCCCGGGATATTCGCGCCAAGATGACAGGGGCGGGGCATTTGCGCGCCTTGACCCGGGGCTTGCGGCGGCTAGTCTGGCGCCATGAG  <  NC_014034/2590196‑2590097                                                                                                                                                                                                                          GCGATCTCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCT                                                                                                     <  2:4070547/101‑1 (MQ=1)         GCGATCTCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCT                                                                                                     <  1:2273980/101‑1 (MQ=1)          CGATCTCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTG                                                                                                    <  2:306223/101‑1 (MQ=1)             TCTCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCC                                                                                                 <  2:2053851/101‑1 (MQ=1)                cgcCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCC                                                                                              <  2:3565146/98‑1 (MQ=1)                CGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCC                                                                                              <  2:4049892/101‑1 (MQ=1)                CGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCC                                                                                              <  2:1130248/101‑1 (MQ=1)                 GGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCG                                                                                             <  2:1975314/101‑1 (MQ=1)                  GCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGg                                                                                            <  1:3313296/101‑2 (MQ=1)                      TGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGggccc                                                                                        <  2:2985159/101‑6 (MQ=1)                      TGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGggccc                                                                                        <  2:2196054/101‑6 (MQ=1)                       GTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCC                                                                                       <  2:853025/101‑1 (MQ=1)                        TTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCC                                                                                      <  2:3486806/101‑1 (MQ=1)                         TGGCGATCAGGATCTTTTGAAACATCGCGTTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCC                                                                                     <  2:527490/101‑1 (MQ=1)                           GCGATCAGGATCTTTTGANACATCGCGGTCNCCTTCGTCAACNGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGG                                                                                   <  1:247721/101‑1 (MQ=1)                            nGATCNGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGn                                                                                  <  1:963020/100‑2 (MQ=1)                              ATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGAT                                                                                <  1:4056427/101‑1 (MQ=1)                              ATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGAT                                                                                <  2:2341522/101‑1 (MQ=1)                               TCAGGATCTTTTGAAACATCGCGGTCTCCTTCGACAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATA                                                                               <  1:3677459/101‑1 (MQ=1)                                CAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATAT                                                                              <  2:47642/101‑1 (MQ=1)                                CAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATAT                                                                              <  1:2178667/101‑1 (MQ=1)                                 AGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATT                                                                             <  1:2303192/101‑1 (MQ=1)                                  GGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTt                                                                            <  1:939176/101‑2 (MQ=1)                                  GGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTt                                                                            <  2:3708778/101‑2 (MQ=1)                                      CTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTtgggg                                                                        <  2:1861649/101‑6 (MQ=1)                                                                                            gacccGGCGGGGGGCGCTGCCCCCCGCACCCCCCGGGATATTCGCGCCAAGATGACAGGGGCGGGGCATTTGCGCGCCTTGACCCGGGGCTTGCGGCGGCT                  <  2:1964188/96‑1 (MQ=11)                                                                                             gcccGGCGGGGGGCGCTGCCCCCCGCACCCCCCGGGATATTCGCGCCAAGATGACAGGGGCGGGGCATTTGCGCGCCTTGACCCGGGGCTTGCGGCGGCTA                 <  1:2765393/97‑1 (MQ=17)                                                                                             acccGGCGGGGGGCGCTGCCCCCCGCACCCCCCGGGATATTCGCGCCAAGATGACAGGGGCGGGGCATTTGCGCGCCTTGACCCGGGGCTTGCGGCGGCTA                 <  1:1340096/97‑1 (MQ=14)                                                                                               ccGGCGGGGGGCGCTGCCCCCCGCACCCCCCGGGATATTCGCGCCAAGATGACAGGGGCGGGGCATTTGCGCGCCTTGACCCGGGGCTTGCGGCGGCTAGT               <  1:415760/99‑1 (MQ=18)                                                                                                 GGCGGGGGGCGCTGCCCCCCGCACCCCCCGGGATATTCGCGCCAAGATGACAGGGGCGGGGCATTTGCGCGCCTTGACCCGGGGCTTGCGGCGGCTAGTCT             <  1:2915689/101‑1 (MQ=33)                                                                                                 GGCGGGGGGCGCTGCCCCCCGCACCCCCCGGGATATTCGCGCCAAGATGACAGGGGCGGGGCATTTGCGCGCCTTGACCCGGGGCTTGCGGCGGCTAGTCT             <  1:92352/101‑1 (MQ=33)                                                                                                  GCGGGGGGCGCTGCCCCCCGCACCCCCCGGGATATTCGCGCCAAGATGACAGGGGCGGGGCATTTGCGCGCCTTGACCCGGGGCTTGCGGCGGCTAGTCTG            <  1:1100577/101‑1 (MQ=32)                                                                                                   CGGGGGGCGCTGCCCCCCGCACCCCCCGGGATATTCGCGCCAAGATGACAGGGGCGGGGCATTTGCGCGCCTTGACCCGGGGCTTGCGGCGGCTAGTCTGG           <  2:2818596/101‑1 (MQ=31)                                                                                                   CGGGGGGCGCTGCCCCCCGCACCCCCCGGGATATTCGCGCCAAGATGACAGGGGCGGGGCATTTGCGCGCCTTGACCCGGGGCTTGCGGCGGCTAGTCTGG           >  1:716963/1‑101 (MQ=31)                                                                                                    GGGGGGCGCTGCCCCCCGCACCCCCCGGGATATTCGCGCCAAGATGACAGGGGCGGGGCATTTGCGCGCCTTGACCCGGGGCTTGCGGCGGCTAGTCTGGC          >  1:2189388/1‑101 (MQ=11)                                                                                                                                                                                                                  CCCGACAGGCGATCTCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGC....................................................................................................  <  NC_014034/1640113‑1640006 ....................................................................................................GGGGGCGCTGCCCCCCGCACCCCCCGGGATATTCGCGCCAAGATGACAGGGGCGGGGCATTTGCGCGCCTTGACCCGGGGCTTGCGGCGGCTAGTCTGGCGCCATGAG  <  NC_014034/2590196‑2590097 Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 33 ≤  ATCG  < 37 ≤  ATCG