New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 1640006 =69 (0.320)47 (0.220) 10/184 29.3 0.42 intergenic (+86/‑59) mccB/mccA methylcrotonoyl‑CoA carboxylase subunit beta/methylcrotonoyl‑CoA carboxylase subunit alpha
?NC_014034 = 2590196 62 (0.260)intergenic (‑95/‑161) glnE/RCAP_rcc02415 glutamate‑ammonia‑ligase adenylyltransferase/hypothetical protein
Rejected: Position hash score below cutoff.

CCCGACAGGCGATCTCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGC....................................................................................................  <  NC_014034/1640113‑1640006
....................................................................................................GGGGGCGCTGCCCCCCGCACCCCCCGGGATATTCGCGCCAAGATGACAGGGGCGGGGCATTTGCGCGCCTTGACCCGGGGCTTGCGGCGGCTAGTCTGGCGCCATGAG  <  NC_014034/2590196‑2590097
                                                                                                                                                                                                                
         CGATCTCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTG                                                                                                    <  2:3318101/101‑1 (MQ=1)
         CGATCTCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTG                                                                                                    <  2:1192324/101‑1 (MQ=1)
          GATCTCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGC                                                                                                   <  1:325783/101‑1 (MQ=1)
           ATCTCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGGGCGGCC                                                                                                  >  1:1973800/1‑101 (MQ=0)
            TCTCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCC                                                                                                 <  1:2811728/101‑1 (MQ=1)
             CTCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGGGGGGCCCC                                                                                                >  1:1775819/1‑101 (MQ=0)
               CGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCC                                                                                              <  1:733607/101‑1 (MQ=1)
               CGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCC                                                                                              <  1:991799/101‑1 (MQ=1)
               CGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCC                                                                                              <  1:1753373/101‑1 (MQ=1)
               CGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCC                                                                                              <  2:4586190/101‑1 (MQ=1)
               CGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCC                                                                                              <  2:112525/101‑1 (MQ=1)
               CGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCC                                                                                              <  1:2394281/101‑1 (MQ=1)
               CGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCC                                                                                              <  1:3827010/101‑1 (MQ=1)
               CGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCC                                                                                              <  1:2907105/101‑1 (MQ=1)
                GGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGGGGGGCCCCCCc                                                                                             >  1:2944500/1‑100 (MQ=0)
                GGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCG                                                                                             <  1:2995333/101‑1 (MQ=1)
                GGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCG                                                                                             <  2:938051/101‑1 (MQ=1)
                  CCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGgg                                                                                           <  2:2302728/101‑3 (MQ=1)
                    CTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGggcc                                                                                         <  2:4437627/101‑5 (MQ=1)
                    CTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGggcc                                                                                         <  1:483894/101‑5 (MQ=1)
                     TGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGggccc                                                                                        <  2:1400708/101‑6 (MQ=1)
                      GTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCC                                                                                       <  1:804821/101‑1 (MQ=1)
                       TTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGGGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCC                                                                                      <  1:4478553/101‑1 (MQ=1)
                       TTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCC                                                                                      <  1:2475538/101‑1 (MQ=1)
                        TGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCC                                                                                     <  1:2626144/101‑1 (MQ=1)
                        TGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCC                                                                                     <  1:1210783/101‑1 (MQ=1)
                         GGCGATCAGGATCTTTTNAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCG                                                                                    <  2:892503/101‑1 (MQ=1)
                         GGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCG                                                                                    <  2:4353945/101‑1 (MQ=1)
                         GGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCG                                                                                    <  1:4382835/101‑1 (MQ=1)
                          GCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGG                                                                                   <  1:1907525/101‑1 (MQ=1)
                           CGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGG                                                                                  <  1:3028055/101‑1 (MQ=1)
                            GATCAGGATCTTTTGAAACATCGCGGTCTCCTTCCTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGA                                                                                 <  1:4002554/101‑1 (MQ=1)
                             gTCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCACCGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGAT                                                                                <  2:3661330/100‑1 (MQ=1)
                             ATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGAT                                                                                <  2:2788189/101‑1 (MQ=1)
                             ATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGAT                                                                                <  1:1884789/101‑1 (MQ=1)
                              TCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATA                                                                               <  1:3437014/101‑1 (MQ=1)
                               CAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATAT                                                                              <  2:449998/101‑1 (MQ=1)
                               CAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATAT                                                                              <  2:645818/101‑1 (MQ=1)
                               CAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATAT                                                                              <  2:3511336/101‑1 (MQ=1)
                                AGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATT                                                                             <  1:695372/101‑1 (MQ=1)
                                AGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATT                                                                             <  2:2849596/101‑1 (MQ=1)
                                 GGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTt                                                                            <  1:1182454/101‑2 (MQ=1)
                                     CTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTtgggg                                                                        <  1:2541920/101‑6 (MQ=1)
                                                                                            acccGGCGGGGGGCGCTGCCCCCCGCACCCCCCGGGATATTCGCGCCAAGATGACAGGGGCGGGGCATTTGCGCGCCTTGACCCGGGGCTTGCGGCGGCTA                 <  2:219530/97‑1 (MQ=14)
                                                                                             cccGGCGGGGGGCGCTGCCCCCCGCACCCCCCGGGATATTCGCGCCAAGATGACAGGGGCGGGGCATTTGCGCGCCTTGACCCGGGGCTTGCGGCGGCTAG                <  1:4527328/98‑1 (MQ=17)
                                                                                                   GGGGGGCGCTGCCCCCCGCACCCCCCGGGATATTCGCGCCAAGATGACAGGGGCGGGGCATTTGCGCGCCTTGACCCGGGGCTTGCGGGGGGTAGTCTTGC          >  2:1380444/1‑101 (MQ=2)
                                                                                                   GGGGGGCGCTGCCCCCCGCACCCCCCGGGATATTCGCGCCAAGATGACAGGGGCGGGGCATTTGCGCGCCTTGACCCGGGGCTTGCGGCGGCTAGTCTGGg          >  1:2198707/1‑100 (MQ=11)
                                                                                                                                                                                                                
CCCGACAGGCGATCTCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGC....................................................................................................  <  NC_014034/1640113‑1640006
....................................................................................................GGGGGCGCTGCCCCCCGCACCCCCCGGGATATTCGCGCCAAGATGACAGGGGCGGGGCATTTGCGCGCCTTGACCCGGGGCTTGCGGCGGCTAGTCTGGCGCCATGAG  <  NC_014034/2590196‑2590097

Base quality scores:  ATCG  < 3 ≤  ATCG  < 21 ≤  ATCG  < 29 ≤  ATCG  < 33 ≤  ATCG  < 37 ≤  ATCG