New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_014034	= 2281196	49 (0.250)	2 (0.010)	2/184	31.5	0.02	coding (2019/2115 nt)	RCAP_rcc02121	PAS domain/exonuclease domain‑containing protein
	?	NC_014034	= 3015138	189 (0.890)					coding (955/1395 nt)	fumC	fumarate hydratase
Rejected: Position hash score below cutoff.

CGGTTCAGGCGTTCATGTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCGGCCGTGGCGACCGTGTCGCCAAGCGCGGTGTGGCGGGCTTCCTCGGGGATGGTGATGCCCAGCCGATGCGTCAGC  >  NC_014034/2281096‑2281289                                                                                                     |                                                                                              cGGTTCAGGCGTTCATGTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGgcagc                                                                                               >  2:4068231/1‑101 (MQ=37)  ggTTCAGGCGTTCATGTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGgcagca                                                                                              <  1:933479/101‑1 (MQ=34)   gTTCAGGCGTTCATGTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGcagcag                                                                                             >  2:395765/1‑101 (MQ=34)    ttCAGGCGTTCATGTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAgg                                                                                            <  2:3217256/101‑1 (MQ=32)    ttCAGGCGTTCATGTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAgg                                                                                            <  1:2931586/101‑1 (MQ=32)    ttCAGGCGTTCATGTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAgg                                                                                            >  1:2268756/1‑101 (MQ=33)     tCAGGCGTTCATGTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGgc                                                                                           <  2:2592756/101‑1 (MQ=32)     tCAGGCGTTCATGTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGgc                                                                                           >  1:2250837/1‑101 (MQ=33)      cAGGCGTTCATGTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGgcg                                                                                          <  2:928630/101‑1 (MQ=32)      cAGGCGTTCATGTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGgcg                                                                                          >  1:1027698/1‑101 (MQ=33)         gCGTTCATGTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGNGCGGGCAGCAGGCGCAg                                                                                       >  2:1019272/1‑101 (MQ=18)         gCGTTCATGTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAg                                                                                       >  2:3600718/1‑101 (MQ=33)         gCGTTCATGTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAg                                                                                       >  1:499721/1‑101 (MQ=32)           gTTCATGTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGa                                                                                     >  2:2620933/1‑101 (MQ=33)            ttCATGTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGaa                                                                                    <  1:3253017/101‑1 (MQ=32)            ttCATGTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGaa                                                                                    <  2:3378799/101‑1 (MQ=32)            ttCATGTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGaa                                                                                    <  1:3719661/101‑1 (MQ=32)            ttCATGTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGaa                                                                                    >  2:3632824/1‑101 (MQ=33)            ttCATGTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGaa                                                                                    >  2:1923154/1‑101 (MQ=33)             tCATGTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAg                                                                                   >  2:3947696/1‑101 (MQ=33)             tCATGTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCCGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAg                                                                                   <  1:3622995/101‑1 (MQ=17)              cATGTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAgg                                                                                  >  1:2792968/1‑101 (MQ=33)               aTGTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGc                                                                                 >  2:2279295/1‑101 (MQ=33)               aTGTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGc                                                                                 >  1:3671734/1‑101 (MQ=33)                tGTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGcc                                                                                >  2:1521000/1‑101 (MQ=32)                tGTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGcc                                                                                <  1:1923154/101‑1 (MQ=32)                 gTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCt                                                                               <  2:4231920/101‑1 (MQ=32)                 gTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCt                                                                               <  2:15702/101‑1 (MQ=32)                  tCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGCCGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTc                                                                              <  1:1797295/101‑1 (MQ=17)                  tCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTc                                                                              <  1:4108600/101‑1 (MQ=32)                  tCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTc                                                                              <  1:2906912/101‑1 (MQ=32)                   ctggcATGAGGCGGTCGTGTTTTCGGTCTTCGGTGAGGACGGCGCCGAGGGTGGTCGGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCGCg                                                                             <  2:2710862/98‑1 (MQ=11)                   cgcgCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCg                                                                             >  1:1893787/1‑101 (MQ=33)                    gcgcATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCgg                                                                            >  2:3061933/1‑101 (MQ=33)                    gcgcATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCgg                                                                            >  1:909720/1‑101 (MQ=33)                     cgcATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGTAGGAAGGCCTCGGc                                                                           >  1:3653299/1‑101 (MQ=17)                     cgcATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCGGc                                                                           >  2:2775168/1‑101 (MQ=33)                      gcATGAGGCGGCCGTGTTTTCGGCCTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCGGcc                                                                          <  2:3925753/101‑1 (MQ=17)                      gcATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCGGcc                                                                          >  1:1875354/1‑101 (MQ=32)                       cATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCGGCCg                                                                         >  1:2229073/1‑101 (MQ=32)                       cATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCGGCCg                                                                         <  2:2268597/101‑1 (MQ=32)                       cATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCGGCCg                                                                         <  2:1085549/101‑1 (MQ=32)                        aTGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCGGCCGt                                                                        >  2:3650955/1‑101 (MQ=32)                          gAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCGGCCGTgg                                                                      >  2:3939596/1‑101 (MQ=32)                           aGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCGGCCGTGGc                                                                     <  2:1399453/101‑1 (MQ=32)                           aGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCGGCCGTGGc                                                                     <  1:3432989/101‑1 (MQ=32)                           aGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGGGGGGAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCGGCCGCGGc                                                                     >  1:1177801/1‑101 (MQ=11)                            ggcggcCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCGGCCGTGGCg                                                                    <  2:108809/101‑1 (MQ=32)                               ggcCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCGGCCGTGGCGAcc                                                                 <  2:4029724/101‑1 (MQ=32)                               ggcCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCGGCCGTGGCGAcc                                                                 >  2:1633881/1‑101 (MQ=32)                               ggcCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCGGCCGTGGCGAcc                                                                 >  2:1214738/1‑101 (MQ=32)                               ggcCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCGGCCGTGGCGAcc                                                                 <  1:2756215/101‑1 (MQ=32)                               ggcCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCGGCCGTGGCGAcc                                                                 <  1:1990881/101‑1 (MQ=32)                               ggcCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCGGCCGTGGCGAcc                                                                 >  1:2092314/1‑101 (MQ=32)                                gcCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCGGCCGTGGCGACCg                                                                <  1:1434806/101‑1 (MQ=32)                                 ccGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCGGCCGTGGCGACCgt                                                               >  1:1845578/1‑101 (MQ=31)                                                                           aGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCGGCCGTGGCGACCGTGTCGCCAAGCGCGGTGTGGCGGGGGTTCTCCGGGGTgggggg                     >  1:501340/1‑97 (MQ=1)                                                                               cgcgGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCGGCCGTGGCGACCGTGTCGCCAAGCGCGGTGTGGCGGGCTTCCCCGGGGGGGGGGAGGccc                 >  2:386223/1‑101 (MQ=1)                                                                                gcgGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCGGCCGTGGCGACCGTGTCGCCAAGCGCGGTGTGGCGGGGGTCCTCCGGGGGGGGGGTGcccc                >  2:3600186/1‑100 (MQ=1)                                                                                gcgGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCGGCCGTGGCGACCGTGTCGCCAAGCGCGGTGTGGCGGGCGTCCTCGGGGGGGGGGGTGcccc                >  2:2234604/1‑100 (MQ=1)                                                                                  gggCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCGGCCGTGGCGACCGTGTCGCCAAGCGCGGTGTGGCGGGGGTTCTCGGGGGGGGGGGTGCCCCGc              >  2:3627623/1‑101 (MQ=1)                                                                                              cgggcAGCAGGCGCAGGAAGGCCTCGGCCGTGGCGACCGTGTCGCCAAGCGCGGTGTGGCGGGCTTCCCCGGGGGTGGGGATGCCCAGCCCGTGCGGCacg  >  2:315955/1‑99 (MQ=1)                                                                                                     |                                                                                              CGGTTCAGGCGTTCATGTCGCGCATGAGGCGGCCGTGTTTTCGGACTTCGGTGAGGACGGCGCCGAAGGTGGTCAGGCCGCGGGCCTTGAGGGCGGGCAGCAGGCGCAGGAAGGCCTCGGCCGTGGCGACCGTGTCGCCAAGCGCGGTGTGGCGGGCTTCCTCGGGGATGGTGATGCCCAGCCGATGCGTCAGC  >  NC_014034/2281096‑2281289 Base quality scores:  ATCG  < 3 ≤  ATCG  < 23 ≤  ATCG  < 30 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG