New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 1640006 =58 (0.270)57 (0.260) 9/184 30.1 0.57 intergenic (+86/‑59) mccB/mccA methylcrotonoyl‑CoA carboxylase subunit beta/methylcrotonoyl‑CoA carboxylase subunit alpha
?NC_014034 2070642 = 28 (0.120)intergenic (+57/+115) dcd/guaB deoxycytidine triphosphate deaminase/inosine‑5'‑monophosphate dehydrogenase
Rejected: Position hash score below cutoff.

CCCGACAGGCGATCTCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGC....................................................................................................  <  NC_014034/1640113‑1640006
....................................................................................................GGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCGGGGTTTCCCCCGGCCCTTGTGAT  >  NC_014034/2070642‑2070741
                                                                                                                                                                                                                
         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:3827010/101‑1 (MQ=1)
               CGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCC                                                                                              <  1:2907105/101‑1 (MQ=1)
               CGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCC                                                                                              <  1:733607/101‑1 (MQ=1)
               CGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCC                                                                                              <  1:991799/101‑1 (MQ=1)
               CGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCC                                                                                              <  2:112525/101‑1 (MQ=1)
               CGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCC                                                                                              <  1:2394281/101‑1 (MQ=1)
               CGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCC                                                                                              <  1:1753373/101‑1 (MQ=1)
               CGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCC                                                                                              <  2:4586190/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‑2 (MQ=1)
                    CTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCC                                                                                         <  2:4437627/101‑1 (MQ=1)
                    CTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCC                                                                                         <  1:483894/101‑1 (MQ=1)
                     TGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCC                                                                                        <  2:1400708/101‑1 (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:1210783/101‑1 (MQ=1)
                        TGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCC                                                                                     <  1:2626144/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:645818/101‑1 (MQ=1)
                               CAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATAT                                                                              <  2:449998/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‑1 (MQ=1)
                                     CTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCTGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTgggg                                                                        <  2:3654005/101‑5 (MQ=1)
                                     CTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTgggg                                                                        <  1:2541920/101‑5 (MQ=1)
                                       TTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTGGGGCA                                                                      <  1:1549628/101‑1 (MQ=1)
                                           AAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTGGGGCAAGAT                                                                  <  1:1349631/101‑1 (MQ=1)
                                            AACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTGGGGCAAGATG                                                                 <  1:3058043/101‑1 (MQ=1)
                                            AACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTGGGGCAAGATG                                                                 <  2:3897908/101‑1 (MQ=1)
                                            AACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTGGGGCAAGATG                                                                 <  2:4469830/101‑1 (MQ=1)
                                             ACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTGGGGCAAGATGA                                                                <  2:3346597/101‑1 (MQ=1)
                                             ACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTGGGGCAAGATGA                                                                <  2:3558206/101‑1 (MQ=1)
                                             ACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTGGGGCAAGATGA                                                                <  2:3796993/101‑1 (MQ=1)
                                             ACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTGGGGCAAGATGA                                                                <  1:3649654/101‑1 (MQ=1)
                                               ATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTGGGGCAAGATGAAa                                                              <  1:444767/101‑2 (MQ=1)
                                                                                             cccgcCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCGGGGTTTCCC                <  1:3742480/96‑1 (MQ=14)
                                                                                               cgcCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCGGGGTTTCCCCC              <  1:568403/98‑1 (MQ=0)
                                                                                                gcCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCGGGGTTTCCCCCG             <  1:628897/99‑1 (MQ=11)
                                                                                                                                                                                                                
CCCGACAGGCGATCTCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGC....................................................................................................  <  NC_014034/1640113‑1640006
....................................................................................................GGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCGGGGTTTCCCCCGGCCCTTGTGAT  >  NC_014034/2070642‑2070741

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