New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_014034	= 2330175	4 (0.020)	77 (0.390)	29/186	14.7	0.95	intergenic (+150/+162)	RCAP_rcc02161/mrcA	hypothetical protein/penicillin‑binding protein 1A
	?	NC_014034	2330176 =	4 (0.020)					intergenic (+151/+161)	RCAP_rcc02161/mrcA	hypothetical protein/penicillin‑binding protein 1A
Rejected: Position hash score below cutoff.

TGGCCTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.G.......CGGCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACC  >  NC_014034/2330080‑2330244                                                                                                         |                                                                     tGGCCTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.G.......Ggggg                                                                  >  2:783732/1‑99 (MQ=255)  ggCCTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.G.......Ggggct                                                                 <  1:2023306/101‑4 (MQ=255)   gCCTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.G.......Ggggctc                                                                <  2:3030/101‑5 (MQ=255)   gCCTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.G.......Ggggctc                                                                <  2:2589968/101‑5 (MQ=255)    ccTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.G.......Ggggctcc                                                               <  1:1738679/101‑6 (MQ=255)    ccTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.G.......Ggggctcc                                                               <  1:3210951/101‑6 (MQ=255)       ccGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.G.......Ggggggcgcgc                                                            >  1:2374538/1‑93 (MQ=255)       ccGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.G.......Ggggctccgcc                                                            <  1:2267797/101‑9 (MQ=255)           cAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGCGG.......GGGCGCCGCCcccc                                                         >  2:1376707/1‑101 (MQ=38)                gcggcAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCCGGC.G.......GGGGCGCCCCccccccaccc                                                   >  1:261050/1‑96 (MQ=255)                             gcAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTgc                                             <  1:2508570/101‑1 (MQ=18)                              cAGAGCCCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTgcg                                            <  1:143234/101‑1 (MQ=18)                                 agCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCa                                         <  2:137965/101‑1 (MQ=18)                                    gccgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGa                                      <  1:791891/101‑1 (MQ=18)                                    gccgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGa                                      <  1:3550770/101‑1 (MQ=18)                                    gccgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGa                                      <  1:254394/101‑1 (MQ=18)                                     ccgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGaa                                     <  1:2771336/101‑1 (MQ=18)                                     ccgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGaa                                     <  2:3300611/101‑1 (MQ=18)                                      cgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGaaa                                    <  2:158467/101‑1 (MQ=18)                                       gcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGaaaa                                   <  2:174791/101‑1 (MQ=18)                                       gcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTCGCGCCAAGaaaa                                   <  1:2287543/101‑1 (MQ=18)                                        cACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGaaaaa                                  <  1:1565773/101‑1 (MQ=17)                                         aCGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCACGGCGTATTTGCGCCAAGAAAAAg                                 <  1:1158062/101‑1 (MQ=18)                                          cGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGc                                <  2:470202/101‑1 (MQ=18)                                          cGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGc                                <  1:4231225/101‑1 (MQ=18)                                          cGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGc                                <  1:4109795/101‑1 (MQ=18)                                          cGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGc                                <  1:2804222/101‑1 (MQ=18)                                           ggATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCTCCCCCGGCGTATTTGCGCCAAGAAAAAGcc                               <  2:3682415/101‑1 (MQ=17)                                            gATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCg                              >  2:1947682/1‑101 (MQ=17)                                             aTCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCgg                             <  2:2950817/101‑1 (MQ=18)                                             aTCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCgg                             <  2:2525313/101‑1 (MQ=18)                                               cGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCgggg                           <  2:2374362/101‑1 (MQ=17)                                                ggCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGgc                          <  2:2248284/101‑1 (MQ=17)                                                 gCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGgcg                         <  2:1650494/101‑1 (MQ=17)                                                  ccGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGgcgc                        <  1:68240/101‑1 (MQ=17)                                                   cGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGgcgcg                       <  2:2013409/101‑1 (MQ=17)                                                     aTGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCCCACCCCCGGCGGATTTGCGCCAAGAAAAAGCCGGGGGGCGCa                     >  2:633029/1‑101 (MQ=12)                                                      tGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCAt                    <  1:3119717/101‑1 (MQ=17)                                                       ggCGGAGAGTAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATg                   <  2:433213/101‑1 (MQ=17)                                                       ggCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATg                   <  2:2385734/101‑1 (MQ=17)                                                       ggCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATg                   <  2:1124815/101‑1 (MQ=17)                                                        gCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGc                  <  1:60966/101‑1 (MQ=18)                                                         cGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCg                 <  2:1152945/101‑1 (MQ=17)                                                         cGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGA.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGGGACAAGAAAAAGACGCGGCGCGGATGCg                 >  2:1102766/1‑101 (MQ=9)                                                               ggAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGccc           <  1:1376807/101‑1 (MQ=17)                                                                gAATTTCGCCGCGCAGGCGCGGCGCCCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGcccc          <  1:938465/101‑1 (MQ=17)                                                                 aaTTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCg         <  2:2083351/101‑1 (MQ=17)                                                                  aTTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCgg        <  2:2754206/101‑1 (MQ=17)                                                                    ttCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGAt      >  1:2454994/1‑101 (MQ=17)                                                                      cgcagcGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCCCCTCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCa    <  2:969230/101‑1 (MQ=12)                                                                       gccgcGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCAc   <  1:606938/101‑1 (MQ=18)                                                                        ccgcGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCAcc  <  2:3432700/101‑1 (MQ=18)                                                                                                         |                                                                     TGGCCTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.G.......CGGCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACC  >  NC_014034/2330080‑2330244 Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG