New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_014034	158231 =	64 (0.310)	25 (0.120)	4/198	30.4	0.34	intergenic (+34/‑196)	parC/tuf1	DNA topoisomerase IV subunit A/translation elongation factor Tu
	?	NC_014034	1253896 =	32 (0.150)					intergenic (+54/+37)	RCAP_rcc01181/hisH	hypothetical protein/imidazole glycerol phosphate synthase subunit HisH
Rejected: Position hash score below cutoff.

GGGCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTGG  >  NC_014034/158131‑158321                                                                                                     |                                                                                           gggCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCgg                                                                                            <  1:2754677/101‑1 (MQ=255)  gccAAGCGGGCGAGCGACGGTCGGGTGGCGGCGCGCGGTTTGCCGCGCTGCAACCGGTTCACCTGACGGGGTTTCGGGCGCTGCCGCACCCGACACGCggg                                                                                           <  2:2717809/99‑1 (MQ=255)    cAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCCGCCGCGCCCGACCCGCggggg                                                                                         <  1:2400178/101‑1 (MQ=255)      aGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCt                                                                                       <  1:1443798/101‑1 (MQ=255)       gCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTc                                                                                      <  1:3720198/101‑1 (MQ=255)       gCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTc                                                                                      <  1:4138733/101‑1 (MQ=255)        cGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAGCCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTcc                                                                                     <  1:2718556/101‑1 (MQ=255)        cGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTcc                                                                                     <  2:531793/101‑1 (MQ=255)          ggCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGc                                                                                   <  1:875704/101‑1 (MQ=255)           gCGCGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCACGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGTCCTTCCCGCGTTCGCTCCGcc                                                                                  <  2:3511186/101‑1 (MQ=255)            cGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGccc                                                                                 <  1:1590976/101‑1 (MQ=255)             gAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGcccc                                                                                <  2:1454663/101‑1 (MQ=255)                  ccGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCAc                                                                           <  1:3735030/101‑1 (MQ=255)                    gggCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCAccc                                                                         <  2:516418/101‑1 (MQ=255)                       cGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCAcccccc                                                                      <  1:2873971/101‑2 (MQ=255)                        tgATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCAGCACCCCCgg                                                                     <  1:3970716/100‑1 (MQ=255)                         gATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGc                                                                    <  1:437954/101‑1 (MQ=255)                            ggCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTa                                                                 <  2:2261613/101‑1 (MQ=255)                            ggCGCCGCGCGGTTCCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTa                                                                 <  2:1242703/101‑1 (MQ=255)                            ggCCCCGCGCGGGTTCCCGCGCCACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCCCCCGACCCGCGGGGGCTCCCGCCCCGGACCCCCGGCgga                                                                 <  2:4075815/101‑3 (MQ=255)                             gcgcCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAt                                                                <  1:2292356/101‑1 (MQ=255)                             gcgcCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAt                                                                <  1:133170/101‑1 (MQ=255)                              cgccgcGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAtt                                                               <  1:1132999/101‑1 (MQ=255)                               gccgcGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAttt                                                              <  2:2516726/101‑1 (MQ=255)                                ccgcGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAtttt                                                             <  2:3562982/101‑1 (MQ=255)                                ccgcGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAtttt                                                             <  2:228368/101‑1 (MQ=255)                                 cgcgcgGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTg                                                            <  1:3009398/101‑1 (MQ=255)                                 cgcgcgGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTg                                                            <  1:2031226/101‑1 (MQ=255)                                      ggTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAg                                                       <  2:702820/101‑1 (MQ=255)                                          tCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGaaa                                                   <  1:2208524/101‑1 (MQ=255)                                           cccGCGCGACAACCGTTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGaaaa                                                  <  2:1332295/101‑1 (MQ=255)                                           cccGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGaaaa                                                  <  1:2908053/101‑1 (MQ=255)                                           cccGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGaaaa                                                  <  1:2619926/101‑1 (MQ=255)                                               cgcgACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCaa                                              <  1:2708767/101‑1 (MQ=255)                                               cgcgACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCaa                                              <  2:2034292/101‑1 (MQ=255)                                                  gACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGt                                           <  1:3477323/101‑1 (MQ=255)                                                    cacccGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGCCCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGc                                         <  2:3753061/98‑1 (MQ=255)                                                     aaCCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGcc                                        <  2:3767569/101‑1 (MQ=255)                                                      aggggTTCCCCTGACGGGGGTGCGGGCGCTGCCGCGCCCGCCCCGCGGGGGCTCCGCCCCCGCACGCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGccc                                       <  2:1569950/98‑1 (MQ=37)                                                       ccGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGcccc                                      <  2:3233003/101‑1 (MQ=255)                                                         gtttCCCCTGACGGGGTTTCGGGCGGTGCCGCCCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGcccccc                                    <  1:3964507/99‑1 (MQ=255)                                                         ggTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGcccccc                                    <  2:3986989/101‑1 (MQ=255)                                                          gTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCg                                   <  2:691190/101‑1 (MQ=255)                                                           ttCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCgg                                  <  2:2535233/101‑1 (MQ=255)                                                           ttCCCCGGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCgg                                  <  1:755045/101‑1 (MQ=255)                                                             ccccTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGAc                                <  2:3954211/101‑1 (MQ=255)                                                               ccTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACgt                              <  1:2412699/101‑1 (MQ=255)                                                                cTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACgtg                             <  2:3994755/101‑1 (MQ=255)                                                                cTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACgtg                             <  2:918737/101‑1 (MQ=255)                                                                cTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACgtg                             <  2:812729/101‑1 (MQ=255)                                                                 tGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACgtgt                            <  2:2209994/101‑1 (MQ=255)                                                                    cGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGt                         <  2:1883179/101‑1 (MQ=255)                                                                           tCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCTCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCgggg                  <  1:3571874/101‑1 (MQ=255)                                                                           tCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCgggg                  <  1:3666088/101‑1 (MQ=255)                                                                              ggCGCTGCCGCGCCCGGCCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGAcg               <  2:3841580/101‑1 (MQ=255)                                                                              ggCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGAcg               <  2:3231975/101‑1 (MQ=255)                                                                               gcgcGGGCGCGCCCGGCCCGCGGGGGCCCCGCCCCCGGACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGAcgc              <  2:1837796/101‑1 (MQ=38)                                                                                 gcTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGAcgcgc            <  2:3410770/101‑1 (MQ=255)                                                                                 gcTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGAcgcgc            <  1:823773/101‑1 (MQ=255)                                                                                  cTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGAcgcgcg           <  2:3027202/101‑1 (MQ=255)                                                                                     ccGCGGCTGCCCCGCGGGGGCTCCGCCCCCGCTCCCCCGGCGTATTTTGCCAAGGAGAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTg        <  2:2512590/101‑1 (MQ=38)                                                                                     ccGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTg        <  2:2177597/101‑1 (MQ=255)                                                                                      cgcgCCCGACCCGCGGGGGCTCCGCCCCCGCGCCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGc       <  2:2117500/101‑1 (MQ=255)                                                                                           ccGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTgg  <  2:2242647/101‑1 (MQ=255)                                                                                           ccGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTgg  <  2:535078/101‑1 (MQ=255)                                                                                                     |                                                                                           GGGCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTGG  >  NC_014034/158131‑158321 Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 23 ≤  ATCG  < 32 ≤  ATCG  < 38 ≤  ATCG