New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 2330179 =41 (0.190)25 (0.110) 19/194 22.4 0.16 intergenic (+154/+158) RCAP_rcc02161/mrcA hypothetical protein/penicillin‑binding protein 1A
?NC_014034 2555842 = 223 (0.990)intergenic (‑705/+41) mraZ/RCAP_rcc02384 protein MraZ/Mrp/NBP35 family protein
Rejected: Position hash score below cutoff.

CCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.G.......CGGCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAT  >  NC_014034/2330089‑2330279
                                                                                                  |                                                                                                    
ccAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGCTG.......GGGCTCCGCCccc                                                                                             <  2:2661408/101‑1 (MQ=38)
                    cAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTgcg                                                                               <  1:4183025/101‑1 (MQ=18)
                     agagCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTgcgc                                                                              <  1:614274/101‑1 (MQ=18)
                     agagCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTgcgc                                                                              <  1:4269809/101‑1 (MQ=18)
                       agCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCa                                                                            <  1:2923123/101‑1 (MQ=18)
                       agCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCa                                                                            <  2:3173422/101‑1 (MQ=18)
                        gccccgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCaa                                                                           <  1:3155748/98‑1 (MQ=18)
                         cgccgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGTCCCCGCACCCCCGGCGTATTTGCGCCAAg                                                                          <  2:3902959/101‑1 (MQ=21)
                         cgccgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAg                                                                          <  2:3723166/101‑1 (MQ=18)
                         cgccgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAg                                                                          <  1:3367409/101‑1 (MQ=18)
                           ccgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGaa                                                                        <  1:1377464/101‑1 (MQ=18)
                            cgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGaaa                                                                       <  1:3363685/101‑1 (MQ=18)
                              cACGGATCGGCCGATTGTGGAGAGGAATTTCGCCGCGCAGGCGCGGCGTCCCGGGC.GGGGGCTCCGCCCCCGCTCCCCCGGCGTATTTGCGCCAAGaaaaa                                                                     <  1:3690095/101‑1 (MQ=12)
                              cACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCCCCCCCGGCGTATTTGCGCCAAGaaaaa                                                                     <  2:500791/101‑1 (MQ=17)
                              cACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGaaaaa                                                                     <  2:3567426/101‑1 (MQ=17)
                                cGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGc                                                                   <  1:3710638/101‑1 (MQ=18)
                                 ggttgggCCTATGGCGGAGAGGAATTTCGCCTCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGcc                                                                  <  1:3142471/96‑1 (MQ=9)
                                 ggATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGcc                                                                  <  1:2692365/101‑1 (MQ=18)
                                 ggATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGcc                                                                  <  2:2661258/101‑1 (MQ=18)
                                       gCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGgcg                                                            <  2:4341996/101‑1 (MQ=17)
                                       gCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGgcg                                                            <  1:1572619/101‑1 (MQ=17)
                                        ccGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGgcgc                                                           <  1:2169334/101‑1 (MQ=17)
                                        ccGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCCCACCCCCGGCGTATTTGCGCCAAGAAAAAAACGGGgcgc                                                           >  1:191330/1‑101 (MQ=9)
                                         cGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGgcgcg                                                          <  2:3020424/101‑1 (MQ=17)
                                         cGATGGCGCAGAGGCATTTCGCCGCGCAGGCGCAGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGgcgcg                                                          <  1:3848669/101‑1 (MQ=9)
                                            tGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCAt                                                       <  1:2744976/101‑1 (MQ=17)
                                            tGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCAt                                                       <  1:4005371/101‑1 (MQ=17)
                                             ggCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATg                                                      <  1:684923/101‑1 (MQ=17)
                                             ggCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATg                                                      >  1:3254901/1‑101 (MQ=17)
                                                ggAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCGCCCCCGGCGTATTTGCGCTAAGAAAAAGCCGGGGCGCGCATGCgg                                                   <  1:888535/101‑1 (MQ=12)
                                                ggAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCgg                                                   <  2:340920/101‑1 (MQ=17)
                                                ggAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCga                                                   <  2:3456097/101‑2 (MQ=17)
                                                 gagaGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCggg                                                  >  2:4233911/1‑101 (MQ=17)
                                                 gagaGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCggg                                                  <  2:2559754/101‑1 (MQ=17)
                                                      gAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGcccc                                             >  2:1288734/1‑101 (MQ=17)
                                                      gAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGcccc                                             <  2:1961365/101‑1 (MQ=17)
                                                      gAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGCGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGcccc                                             <  2:2521389/101‑1 (MQ=18)
                                                        aTTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCgg                                           <  1:1619891/101‑1 (MQ=17)
                                                           tCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATc                                        <  2:2522010/101‑1 (MQ=18)
                                                            cgccgcGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCa                                       <  2:928724/101‑1 (MQ=18)
                                                               cgcgCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACcg                                    <  2:1168758/101‑1 (MQ=18)
                                                                                         gggctccgcccccGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTggg           <  2:3067909/92‑1 (MQ=255)
                                                                                           gctccgcccccGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGc         >  1:3462125/8‑101 (MQ=255)
                                                                                            ctccgcccccGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGcc        <  1:3304813/95‑1 (MQ=255)
                                                                                             tccgcccccGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCa       <  2:3411448/96‑1 (MQ=255)
                                                                                             tccgcccccGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCa       >  1:2376141/6‑101 (MQ=255)
                                                                                             tccgcccccGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCa       >  1:2244058/6‑101 (MQ=255)
                                                                                             tccgcccccGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCa       >  1:1478313/6‑101 (MQ=255)
                                                                                              ccgcccccGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCaa      >  2:1217827/5‑101 (MQ=255)
                                                                                                gcccccGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAg    >  2:1435079/3‑101 (MQ=255)
                                                                                                  cccgGCTCCGCCGGCGTATTTTCGCCAAGAAAAAGCCGGGGGGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAt  <  2:3890950/101‑1 (MQ=255)
                                                                                                  ccccGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAt  >  1:621185/1‑101 (MQ=255)
                                                                                                  |                                                                                                    
CCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.G.......CGGCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAT  >  NC_014034/2330089‑2330279

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG