New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 1640006 =51 (0.250)49 (0.240) 9/184 27.9 0.62 intergenic (+86/‑59) mccB/mccA methylcrotonoyl‑CoA carboxylase subunit beta/methylcrotonoyl‑CoA carboxylase subunit alpha
?NC_014034 2070642 = 8 (0.040)intergenic (+57/+115) dcd/guaB deoxycytidine triphosphate deaminase/inosine‑5'‑monophosphate dehydrogenase
Rejected: Position hash score below cutoff.

CCCGACAGGCGATCTCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGC..................................................................................................  <  NC_014034/1640113‑1640006
....................................................................................................GGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCGGGGTTTCCCCCGGCCCTTGTG  >  NC_014034/2070642‑2070739
                                                                                                                                                                                                              
        GCGATCTCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCT                                                                                                   <  2:1402525/101‑1 (MQ=1)
        GCGATCTCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCT                                                                                                   <  1:428419/101‑1 (MQ=1)
        GCGATCTCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCT                                                                                                   <  1:4175358/101‑1 (MQ=1)
        GCGATCTCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCT                                                                                                   <  1:2505123/101‑1 (MQ=1)
          GATCTCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGC                                                                                                 <  1:1204484/101‑1 (MQ=1)
           ATCTCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCC                                                                                                <  2:1286528/101‑1 (MQ=1)
           ATCTCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCC                                                                                                <  1:3182644/101‑1 (MQ=1)
             CTCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCC                                                                                              <  1:1106617/101‑1 (MQ=1)
             CTCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCC                                                                                              <  2:2964591/101‑1 (MQ=1)
              TCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCC                                                                                             <  2:3170201/101‑1 (MQ=1)
              TCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCC                                                                                             <  2:3942957/101‑1 (MQ=1)
               CGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCC                                                                                            <  2:3458200/101‑1 (MQ=1)
               CGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCC                                                                                            <  2:2764940/101‑1 (MQ=1)
               CGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCC                                                                                            <  1:2067289/101‑1 (MQ=1)
               CGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCC                                                                                            <  2:1149489/101‑1 (MQ=1)
                GGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCG                                                                                           <  1:2777056/101‑1 (MQ=1)
                 GCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGG                                                                                          <  2:2692596/101‑1 (MQ=1)
                 GCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGG                                                                                          <  2:2557792/101‑1 (MQ=1)
                   CCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGgc                                                                                        <  1:729367/101‑3 (MQ=1)
                        TGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCC                                                                                   <  1:1215935/101‑1 (MQ=1)
                        TGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCC                                                                                   <  1:1648436/101‑1 (MQ=1)
                         GGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCG                                                                                  <  1:479916/101‑1 (MQ=1)
                         GGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCG                                                                                  <  2:496413/101‑1 (MQ=1)
                          GCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGG                                                                                 <  1:2873755/101‑1 (MQ=1)
                          GCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGG                                                                                 <  2:1621041/101‑1 (MQ=1)
                            GATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGA                                                                               <  2:1498276/101‑1 (MQ=1)
                             ATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGAT                                                                              <  1:4387768/101‑1 (MQ=1)
                               CAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATAT                                                                            <  1:3411153/101‑1 (MQ=1)
                                  GATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTg                                                                         <  2:4001707/101‑2 (MQ=1)
                                  GATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTCCCCCCCGGGCCCCCCGGGATATTTg                                                                         <  2:3054581/101‑2 (MQ=1)
                                    TCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTggg                                                                       <  1:4190192/101‑4 (MQ=1)
                                     tTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATctgggg                                                                      <  1:865870/100‑7 (MQ=1)
                                     CTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTgggg                                                                      <  2:2437084/101‑5 (MQ=1)
                                      TTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTggggc                                                                     <  1:1805490/101‑6 (MQ=1)
                                       TTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTGGGGCA                                                                    <  1:4352643/101‑1 (MQ=1)
                                        TTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTGGGGCAA                                                                   <  2:309034/101‑1 (MQ=1)
                                         TGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTGGGGCAAG                                                                  <  2:1631071/101‑1 (MQ=1)
                                         TGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTGGGGCAAG                                                                  <  2:473330/101‑1 (MQ=1)
                                          GAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTGGGGCAAGA                                                                 <  2:1831399/101‑1 (MQ=1)
                                          GAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTGGGGCAAGA                                                                 <  1:3410791/101‑1 (MQ=1)
                                          GAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTGGGGCAAGA                                                                 <  2:3063283/101‑1 (MQ=1)
                                           AAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTGGGGCAAGAT                                                                <  2:650917/101‑1 (MQ=1)
                                           AAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTGGGGCAAGAT                                                                <  2:3199038/101‑1 (MQ=1)
                                           AAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTGGGGCAAGAT                                                                <  1:2541996/101‑1 (MQ=1)
                                              CATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTGGGGCAAGATGAA                                                             <  2:1132183/101‑1 (MQ=1)
                                               ATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTGGGGCAAGATGAAa                                                            <  2:2949176/101‑2 (MQ=1)
                                               ATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGCTGCCCCCCGGGCCCCCCGGGATATTTGGGGCAAGATGAAa                                                            <  1:1348068/101‑2 (MQ=1)
                                                                                          ccccccgcCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCGGGGTTT                 <  1:1954065/93‑1 (MQ=255)
                                                                                                   GGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCGGGGTTTCCCCCGGCC        <  1:1698604/101‑1 (MQ=255)
                                                                                                                                                                                                              
CCCGACAGGCGATCTCGGCCCTGTTGGCGATCAGGATCTTTTGAAACATCGCGGTCTCCTTCGTCAACGGCCGGGGCGCTGCCCCGTGGCGGAAGTGGCGGGGGGCGC..................................................................................................  <  NC_014034/1640113‑1640006
....................................................................................................GGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCGGGGTTTCCCCCGGCCCTTGTG  >  NC_014034/2070642‑2070739

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG