New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_014034	2330179 =	64 (0.230)	34 (0.120)	7/154	30.5	0.22	intergenic (+154/+158)	RCAP_rcc02161/mrcA	hypothetical protein/penicillin‑binding protein 1A
	?	NC_014034	= 2577171	183 (0.510)					intergenic (+93/‑336)	fdx/purU	ferredoxin VI/formyltetrahydrofolate deformylase
Rejected: Position hash score below cutoff.

CCCCGGGACCCCCGGGATCTTTGGCCCCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGG...................................................................................................  <  NC_014034/2330295‑2330179 ..............................................................................................GCAAATACGCCGGGGGTGCGGGGGCG..GCCCCCGCCGGGTCGCCTTGGCGAAGCCAAGGCGAAACGCTTTCGCGCCCAGGCAAAAGGCCGCCCGAAGGCGGCCGTTTCATCGCGATGTCAG  <  NC_014034/2577171‑2577075                                                                                                                                                                                                                                           TCTTTGGCCCCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGG                                                                                                    >  1:5421380/1‑101 (MQ=1)                   CTTTGGCCCCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGC                                                                                                   >  1:2605348/1‑101 (MQ=1)                   CTTTGGCCCCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGC                                                                                                   >  2:6820494/1‑101 (MQ=1)                   CTTTGGCCCCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGCCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGC                                                                                                   <  1:6435370/101‑1 (MQ=1)                    TTTGGCCCCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGCG                                                                                                  >  1:2781970/1‑101 (MQ=1)                    TTTGGCCCCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGCG                                                                                                  >  1:3955461/1‑101 (MQ=1)                    TTTGGCCCCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGCG                                                                                                  <  2:535469/101‑1 (MQ=1)                    TTTGGCCCCAGGTGAAAACCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGCG                                                                                                  >  2:3757785/1‑101 (MQ=1)                     TTGGCCCCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGNCGGGGGTGCGGGGGCG..G                                                                                               >  2:1824999/1‑101 (MQ=1)                     TTGGCCCCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGCG..G                                                                                               >  2:3209959/1‑101 (MQ=1)                     TTGGCCCCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGCG..G                                                                                               <  1:317458/101‑1 (MQ=1)                      gGGCCCCAGGTAAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGCG..Ga                                                                                              <  1:3577230/100‑2 (MQ=1)                      TGGCGCCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGCG..Ga                                                                                              <  2:3088632/101‑2 (MQ=1)                      TGGCCCCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGCG..Ga                                                                                              <  1:5169302/101‑2 (MQ=1)                      TGGCCCCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGCG..Ga                                                                                              >  2:5328451/1‑100 (MQ=1)                      TGGCCCCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGCG..Ga                                                                                              <  2:6665029/101‑2 (MQ=1)                      TGGCCCCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGCG..Ga                                                                                              >  1:869739/1‑100 (MQ=1)                      TGGCCCCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGCG..Ga                                                                                              <  2:3577990/101‑2 (MQ=1)                      TGGCCCCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGCG..Ga                                                                                              >  1:2405457/1‑100 (MQ=1)                       GGCCCCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGCG..Gag                                                                                             <  1:6053257/101‑3 (MQ=1)                       GGCCCCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGCGCGGGGGCG..Gag                                                                                             >  1:6558912/1‑99 (MQ=1)                        GCCCCAGGTGAAGGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGCG..Gagc                                                                                            <  1:2311046/101‑4 (MQ=1)                        GCCCCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGCG..Gagc                                                                                            <  2:497259/101‑4 (MQ=1)                        GCCCCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGCG..Gagc                                                                                            <  2:6986451/101‑4 (MQ=1)                        GCCCCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGCG..Gagc                                                                                            >  1:6913031/1‑98 (MQ=1)                        GCCCCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGCG..Gagc                                                                                            <  1:2660593/101‑4 (MQ=1)                         CCCCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCGCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGCG..Gagcc                                                                                           <  2:1232489/101‑5 (MQ=1)                         CCCCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGCG..Gagcc                                                                                           <  2:5646517/101‑5 (MQ=1)                         CCCCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGCG..Gagcc                                                                                           <  1:2633850/101‑5 (MQ=1)                           CCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGCG..Gagcccc                                                                                         <  2:1465193/101‑7 (MQ=1)                           CCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGCGCGGGGGCG..Gagcccc                                                                                         >  1:723002/1‑95 (MQ=1)                             AGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGCG..Gagcccccg                                                                                       <  1:6598021/101‑9 (MQ=1)                             AGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGCG..Gagcccccg                                                                                       <  1:3194736/101‑9 (MQ=1)                               GTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGGGCGGAGCCCCCGCC                                                                                       >  1:6015718/1‑101 (MQ=31)                                                                                                                                                                                                                          CCCCGGGACCCCCGGGATCTTTGGCCCCAGGTGAAAGCCGCCCGGGCCGCCGGTGATCCGGGGCCCCGCATGCGCGCCCCGGCTTTTTCTTGGCGCAAATACGCCGGGGGTGCGGGG...................................................................................................  <  NC_014034/2330295‑2330179 ..............................................................................................GCAAATACGCCGGGGGTGCGGGGGCG..GCCCCCGCCGGGTCGCCTTGGCGAAGCCAAGGCGAAACGCTTTCGCGCCCAGGCAAAAGGCCGCCCGAAGGCGGCCGTTTCATCGCGATGTCAG  <  NC_014034/2577171‑2577075 Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 40 ≤  ATCG